# NOT RUN {
x = ModelMatrix(uptake ~ Type, CO2)
y = model.frame(uptake ~ Type, CO2)[,1]
rx = lfit(x, y)
est(t(c(0, -1, 1)), rx) # Quevec - Mississippi
t.test(uptake ~ Type, CO2) # compare with the above
# }
Run the code above in your browser using DataLab