R/fRatioStat.R

#
# F-ratio test statistic for effect of groups.
#
fRatioStat<-function(resp,groups) {
  
  grMeans<-tapply(resp,groups,mean)
  grNums<-tapply(resp,groups,length)
  grandMean<-mean(resp)
  n<-length(resp)
	
	grpFactor<-as.factor(groups)
	MSTr <- sum( grNums * (grMeans - grandMean)^2)
	diffs <- resp - grMeans[grpFactor]
	MSE = sum(diffs^2)
	
  (MSTr/(nlevels(grpFactor)-1))/(MSE/(n-nlevels(grpFactor)))
}

Try the EffectsRelBaseline package in your browser

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

EffectsRelBaseline documentation built on May 2, 2019, 9:44 a.m.