首页 > 解决方案 > 难以可视化混合效应模型

问题描述

我的模型拟合似乎与我的数据不匹配。这与对数转换有关吗?随着缩放?还是有其他问题?

我试图可视化距离对 PAI 的影响,这对于我的数据集中的不同类型(T1、T2 和 T3)。我想保持不变的其他预测变量的效果(但如果这不可能,他们可以参与可视化)。

m1PAI <- lmer(PAI~ (scale(HMSL)  +Type + scale(log(Distance)))^2 + scale(ST_Trees)+ (1|Region:Transect),na.action = na.omit, data=DF, REML=FALSE)

m1line <- predict(m1PAI, se.fit=TRUE, re.form=NA,  type= "response")$fit

p1PAI <- ggplot()+geom_point(data=DF,aes(x=Distance, y=PAI,color=Type))+geom_smooth(aes(x=DF$PAI,y=m1line))
p1PAI  

> dput(head(DF,20))
structure(list(Region = c("BE", "BE", "BE", "BE", "BE", "BE", 
"BE", "BE", "BE", "BE", "BE", "BE", "BE", "BE", "BE", "BE", "BE", 
"BE", "BE", "BE"), Elevation = c("HI", "HI", "HI", "HI", "HI", 
"HI", "HI", "HI", "HI", "HI", "HI", "HI", "HI", "HI", "HI", "LO", 
"LO", "LO", "LO", "LO"), Type = c("T1", "T1", "T1", "T1", "T1", 
"T2", "T2", "T2", "T2", "T2", "T3", "T3", "T3", "T3", "T3", "T1", 
"T1", "T1", "T1", "T1"), plot = c("P1", "P2", "P3", "P4", "P5", 
"P1", "P2", "P3", "P4", "P5", "P1", "P2", "P3", "P4", "P5", "P1", 
"P2", "P3", "P4", "P5"), Distance = c(1.5, 4.5, 12.5, 35.5, 99.5, 
1.5, 4.5, 12.5, 35.5, 99.5, 1.5, 4.5, 12.5, 35.5, 99.5, 1.5, 
4.5, 12.5, 35.5, 75), PAI = c(3.97, 5.43, 6.08, 6.09, 6.28, 5.47, 
5.04, 6.1, 6.12, 6.05, 4.38, 4.37, 3.98, 4.03, 4.97, 5.25, 4.46, 
5.33, 6.26, 6.74), ST_Trees = c(3.97, 4.36333333333333, 2.819375, 
2.45, 2.54333333333333, 3.49241379310345, 3.16764705882353, 3.0252380952381, 
3.57545454545455, 2.73, 2.73, 2.73, 2.86071428571429, 2.73, 2.73, 
2.66266666666667, 2.58523178807947, 2.53111111111111, 3.19662790697674, 
3.3348275862069), Transect = c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 
2L, 2L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L)), row.names = c(NA, 
20L), class = "data.frame")

p1PAI:可视化模型和数据点

我希望我的模型预测与我的数据点一致。秤好像有问题。我该如何解决这个问题?

标签: rggplot2visualizationmixed-models

解决方案


推荐阅读