sm.ancova的正確格式

Question

我正在嘗試使用Tukey-Kramer事件運行非參數ANCOVA。有人向我建議sm.ancova可以做到這一點，但他們關於這個過程的信息表很混亂。什麼是「組」功能？ 我有四個網站測量多年。年是我的合作伙伴。

我正在處理的統計問題： 主要影響：1）斜率不同？ 事後：2）哪些網站彼此不同？

任何幫助，將不勝感激！這件事我很迷茫。

當前代碼不起作用：

data1 
    library(sm) 
    attach(data1) 
    x<-year 
    x 
    y<-c(data1[,2:5]) 
    y 
    as.vector(y) 
    sm.ancova(x,y,group=x) 

Answer 1

sm.ancova(year,data,site) 

相信group是你的獨立變量，因此在這種情況下，您的網站（或治療組）。

您還需要將模型指定爲相等或平行。