1

我正在尝试按如下方式执行 lme 函数

mydata <- read.table(
    "H:/edu/Multivariat/HCMpart2.TXT", header=TRUE, sep="\t",
    na.strings="*", dec=",", strip.white=TRUE
)
mydata = data.frame(mydata)

summary(lme(mydata$x~1+mydata$grp+mydata$var, random~1|mydata$id))

其中 x 包含我的值,grp 和 var 表示导致 x 值的组和变量,id 是患者的 id。

其中 HCMpart2.txt 包含一个带有“id grp var x”的标题,中间有标签以及所有这些的相应值。我尝试使用“as.numeric”函数将因子转换为数字因子,但它没有完成我的问题。

当我尝试执行 lme 功能时,我得到以下信息

Error in as.data.frame.default(data) : 
cannot coerce class '"formula"' into a data.frame

任何人都可以帮忙吗?我的印象是我做的一切都是正确的……问候 Cenderze

4

1 回答 1

4

三件事:

就像影子说的那样,您在随机参数中缺少等号。

lme有一个data论点可以阻止您编写带有很多$符号的丑陋代码。

拦截是隐含的。

summary(model <- lme(x ~ grp + var, mydata, random = ~ 1 | id))
于 2013-10-16T11:52:42.457 回答