錯誤使用LME功能

Question

當我試圖如下

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執行LME功能「不能強迫類的‘公式’'成data.frame」代表導致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

Answer 1

三兩件事：

就像影子說，你錯過了一個等於在隨機參數的跡象。

lme有一個data的說法，停止你不得不用很多$標誌寫難看的代碼。

隱含地包含截距。

summary(model <- lme(x ~ grp + var, mydata, random = ~ 1 | id))