我正在尝试按如下方式执行 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