Displaying 1 result from an estimated 1 matches for "veryspeedy".
2008 Mar 07
1
Finding Interaction and main effects contrasts for two-way ANOVA
...contrasts(),
se.contrasts() and fit.contrasts() in package gmodels. Given the url
for a small dataset and the two-way anova model below, I'd like to
reproduce the results from appended SAS code. Thanks. --Dale.
## the dataset (from Montgomery)
twoway <- read.table("http://dsteele.veryspeedy.net/sta501/twoway.txt",
col.names=c('material', 'temp','voltage'),colClasses=c('factor',
'factor', 'numeric'))
## the model
fit <- aov(voltage ~ material*temp, data=twoway)
/* SAS code */
proc glm data=twoway;
class material temp;
model vo...