我正在尝试执行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
发布于 2013-10-16 19:52:42
三件事:
就像shadow说的,你在随机参数中遗漏了一个等号。
lme有一个data参数,它可以让你不必编写带有很多$符号的丑陋代码。
截取被隐式包括在内。
summary(model <- lme(x ~ grp + var, mydata, random = ~ 1 | id))https://stackoverflow.com/questions/19402429
复制相似问题