我有一個連續的回答變量yld和一個分類預測器check（有3個級別）。我做了單因子方差分析和事後檢驗，看看哪些級別彼此不同。 mdl<-aov(sqrt(var$yld) ~ var$check); summary(mdl);TukeyHSD(mdl)
Df Sum Sq Mean Sq F value Pr(>F)
var$check 2 5162 2581.2 13.5
當我在重複測量和varIdent中使用lme時，遇到了奇怪的結果。任何與此有關的幫助將非常感謝！ 我正在測試沿時間序列的葉子的13C信號是否在2種（A和B）之間不同。我基本上對物種之間的整體差異感興趣，而不是特定的時間點。 這裏是我的數據集： Block Species time X13C
1 B 2 0.775040865
2 B 2 0.343913792
3 B 2 0.3810536
我試圖對time_of_day進行所有數據框列的方差分析，這是一個因素。我的其餘專欄都是雙打和同等長度的。 x = 0
pdf("Time_of_Day.pdf")
for (i in names(data_in)){
if(x > 9){
test <- aov(paste(i, "~ time_of_day"), data = data_in)
}