sm.ancova的适当形式

Question

我正在尝试使用Tukey-Kramer帖子运行非参数ANCOVA。有人向我建议sm.ancova可以做到这一点，但这个过程的帮助页令人困惑。什么是"group“函数？我有四个多年来测量过的网站。年份是我的协变量。

我正在解决的统计问题:主要影响: 1)坡度是否不同？Post-hoc: 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是你的自变量，所以在这种情况下，你的网站(或治疗组)。

您还需要将model指定为相等或并行。