R/repeated_measures_function.R

Defines functions RManovaNov

Documented in RManovaNov

RManovaNov<-function(r,k,f,n){
  tm = gl(k, 1, n*k, unique(factor(f)))   # matching treatment
  tm

  blk = gl(n, k, k*n)           # blocking factor
  blk

  av = aov(r ~ blk+tm)
  summary(av)

  qqnorm(av$residuals)
  plot(av$fitted.values,av$residuals)

  TukeyHSD(av)$tm
}

Try the NMVANOVA package in your browser

Any scripts or data that you put into this service are public.

NMVANOVA documentation built on Oct. 6, 2019, 5:05 p.m.