我正在努力編寫一個腳本,該腳本允許使用lme4或nlme軟件包比較不同的線性混合效果模型。由於我不想調整每個添加或刪除模型的腳本,因此我正在尋找一種動態的方法。這樣做我只需要調整一個包含模型公式的字符串的變量。 這工作得很好,除非anova()進來anova()不接受列表包含適當的類要素: ###### Here is my problem
# comparing models by means
我正試圖運行方差分析,用於名爲MegaData的數據框中的二項分佈數據,其中前4列是具有單位,年份和物種的分類變量。下面是我的R代碼,當我試圖循環訪問數據框的不同列以運行我的模型時,它不起作用。 mod <- list()
for (i in MegaData[,5:36]) {
for(j in length(MegaData[,5:36])){
mod[[j]] <- g