#library(DPpackage) No longer supported in R! library(gam) #data(orings) #?orings Flight=1:23 ThermalDistress=c(0,1,0,0,0,0,0,0,1,1,1,0,0,2,0,0,0,0,0,0,2,0,1) NoThermalDistress=6-ThermalDistress Temperature=c(66,70,69,68,67,72,73,70,57,63,70,78,67,53,67,75,70,81,76,79,75,76,58) Pressure=c(rep(50,6),100,100,rep(200,15)) Challenger=data.frame(cbind(Flight,ThermalDistress,NoThermalDistress,Temperature,Pressure)) plot(Challenger) Challenger$Pressure=as.factor(Challenger$Pressure) fit=gam(cbind(ThermalDistress,NoThermalDistress)~s(Temperature)+Pressure+s(Flight), family=binomial(link=logit),data=Challenger) par(mfrow=c(2,2)) plot(fit,se=TRUE) summary(fit)