我已經建立了 followinf 模型
model2 <- lmerTest::lmer(LPP2POz ~ 1 COND (1|ID), data = dataLPP2POz)
如果我嘗試運行以下函式，它會回傳此錯誤：
printCoefmat(summary(model2)$tTable,
has.Pvalue = T, P.values = T)
這是我正在處理的資料集中的一個簡短資料框
dput(head(dataLPP2POz))
structure(list(ID = structure(c(1L, 1L, 1L, 2L, 2L, 2L), .Label = c("01",
"04", "06", "07", "08", "09", "10", "11", "12", "13", "15", "16",
"17", "18", "19", "21", "22", "23", "25", "27", "28", "30", "44",
"46", "49"), class = "factor"), GR = c("RP", "RP", "RP", "RP",
"RP", "RP"), SES = c("V", "V", "V", "V", "V", "V"), COND = structure(c(1L,
2L, 3L, 1L, 2L, 3L), .Label = c("NEG-CTR", "NEG-NOC", "NEU-NOC"
), class = "factor"), LPP2POz = c(7.91468942320841, 9.94838815736199,
10.2186482048953, 1.07455889922813, 1.65917850515029, 3.22422743232682
)), row.names = c(NA, 6L), class = "data.frame")
Error in printCoefmat(summary(model2)$tTable, has.Pvalue = T, P.values = T) :
'x' must be coefficient matrix/data frame
任何人都能夠理解錯誤是什么？
uj5u.com熱心網友回復：
跟進并澄清@aosmith 的評論。
對于不同的混合模型包，提取系數表的作業方式不同。較新的默認值（適用于lme4, lmerTest, glmmTMB）是可以將系數表提取為summary(model)$coefficients（在幕后，這意味著該summary()方法回傳一個串列，其中將系數表存盤為$coefficients）元素。對于這些包，coef(summary(model))是更好的做法。
對于nlme包中的物件（即lme），匯總表存盤為summary(model)$tTable（這很不幸，但nlme比 R 本身更舊......
它不會給出與 完全相同的結果printCoefmat，但您也可以查看一些用于漂亮列印 輸出的選項broom.mixed::tidy()，其目的是構建一個兼容層，這樣您就不必記住所有這些東西......
轉載請註明出處，本文鏈接：https://www.uj5u.com/gongcheng/338773.html
上一篇：R-每列整個資料框的百分比
下一篇：從函式回傳基本圖和ggplot