fm <- lm(Gas ~ Insul/zs(Temp), whiteside)
gm <- lm(Gas ~ Insul/zu(Temp), whiteside)
hm <- lm(Gas ~ Insul/Temp, whiteside)
c(fm = unname(predict(fm, data.frame(Insul = "Before", Temp = 0.0))),
gm = unname(predict(gm, data.frame(Insul = "Before", Temp = 0.0))),
hm = unname(predict(hm, data.frame(Insul = "Before", Temp = 0.0))))
rm(fm, gm, hm)
Run the code above in your browser using DataLab