data(tire)
C <- c(0,1,-1,0,10,-10)
for ( x in seq(15,70,5) ) { C <- rbind( C,c(0,1,-1,0,x,-x) ) }
# numerate the contrasts
rownames(C) <- paste("C", 1:nrow(C), sep="")
# simultaneous confidence intervals for difference of regression functions
summary(simint(cost ~ make*mph, data=tire,
cmatrix=C, eps=0.001))
Run the code above in your browser using DataLab