r - proper form for sm.ancova -
i'm trying run non-parametric ancova tukey-kramer post hoc. suggested me sm.ancova info sheet on proc confusing. "group" function? have 4 sites measured on multiple years. year co-variate.
statistical questions i'm addressing: main effects: 1) slopes different? post-hoc: 2) sites different each other?
any appreciated! i'm lost thing.
current code not work:
data1 library(sm) attach(data1) x<-year x y<-c(data1[,2:5]) y as.vector(y) sm.ancova(x,y,group=x)
sm.ancova(year,data,site)
i believe group
independent variable in case sites (or treatment group).
you need specify model equal or parallel.
Comments
Post a Comment