R/ANOVAposthoc.R

In DeducerANOVA: Factorial ANOVA for Deducer

ANOVAposthoc <- function(data,dv,between) {
	btw <- data[,match(between,names(data))]
	y <- data[,match(dv,names(data))]
	data$tall <- with(data, interaction(btw))
	postHoc <- glht(lm(y ~ tall, data = data),linfct = mcp(tall = "Tukey"))
	return(postHoc)
	}