0

我正在尝试运行线性混合效果模型(包 nlme),但我反复收到错误:'closure' 类型的对象不是子集。

> apoeht <- read.csv("apoeht.csv")  
> library(nlme)  
> model.a <- lme(Timmrec ~ age, data = apoeht, random = ~ age | pathid, 
+                na.exclude)  
Error: object of type 'closure' is not subsettable

谢谢。

4

1 回答 1

3

问题是您正在将一个函数传递na.exclude()给. 实际上,您的电话是:correlationlme()

model.a <- lme(Timmrec ~ age, data = apoeht, random = ~ age | pathid, 
               correlation = na.exclude)

处理correlation参数的代码做出了某些假设,但它肯定不会期望传递一个不相关的函数。

您可能想要使用该na.action参数,但如果您不提供其他参数,则必须为其命名。你要

model.a <- lme(Timmrec ~ age, data = apoeht, random = ~ age | pathid, 
               na.action = na.exclude)
于 2014-01-17T02:40:53.477 回答