1

我有一个在R中使用5个变量开发的模型。选择线性混合建模方法来开发具有嵌套效果的模型。

我用于模型开发的 R 代码如下:

model1 <- lmer(Reduction.factor ~ (1|PAI:Open.wind) + 
     (1|PAI:Temp) + (1|PAI:Height)+ (1|PAI:Density)+ PAI  ,
    data = model)

PAI 嵌套了 4 个参数。

为了了解模型如何在不同条件下进行预测,我为 5 个参数创建了不同的数字作为数据框(称为“案例研究”)。

它看起来像这样:

Temp    Height  Density PAI Open wind
20.000  0.041   0.033   1.960   30.000
20.000  0.082   0.061   1.960   30.000
20.000  0.122   0.059   1.960   30.000
20.000  0.163   0.061   1.960   30.000
20.000  0.204   0.043   1.960   30.000
20.000  0.245   0.048   1.960   30.000
20.000  0.286   0.052   1.960   30.000
40.000  0.082   0.061   1.960   40.000
40.000  0.122   0.059   1.960   40.000
40.000  0.163   0.061   1.960   40.000
40.000  0.204   0.043   1.960   40.000
40.000  0.245   0.048   1.960   40.000
40.000  0.286   0.052   1.960   40.000

我想通过使用 5 个参数的值来预测“Reduction.factor”。

因此,我p1 <- predict(model1,case study)查看了预测..但只有(函数(x,n)中的错误:newdata中检测到的新级别<-出现了此消息..

如何通过将这些值用于 5 个参数来预测“Reduction.factor”?而且..如果我想有置信区间(95%),我应该输入'interval ='confidence''吗?

先感谢您。

4

0 回答 0