适用于sm.ancova的形式

Question

我正在尝试使用Tukey-Kramer post hoc运行非参数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是您的自变量，因此在这种情况下您的网站（或治疗组）。

您还需要将模型指定为相等或平行。