#############################################
#######Commands used to produce STA242/ENV255 lecture plots, March 1 & 6, 2000
#######This script file produces diagnostic plots of the pollen data.
#############################################
attach(pollen)
code_c(rep(0,35),rep(1,12))
#############################################
###Tentative model is parallel lines model###
#############################################
########coded scatterplot
par(mfrow=c(1,1))
plot(log(duration),log(removed/(1-removed)),ylab="logit(p)",xlab="log(duration)",type="n",main="Parallel Lines Model")
points(log(duration)[code==1],log(removed/(1-removed))[code==1],pch=1)
points(log(duration)[code==0],log(removed/(1-removed))[code==0],pch=2)
########Fit parallel lines model
parallel.mod_lm(log(removed/(1-removed))~log(duration)+code,qr=T)
coef_summary(parallel.mod)$coefficients[,1]
abline(coef[1],coef[2])
abline(coef[1]+coef[3],coef[2])
text(log(duration),log(removed/(1-removed))+.15,labels)
########Now get diagnostics
parallel.mod.diag_ls.diag(parallel.mod)
#############################################
###Note: elements of parallel.mod.diag are:
###[1] "std.dev"      "hat"          "std.res"      "stud.res"     "cooks"       
###[6] "dfits"        "correlation"  "std.err"      "cov.unscaled"
#############################################
###########Residual vs. Fitted Plots
par(mfrow=c(1,1))
plot(parallel.mod$fitted,parallel.mod$residuals,main="Resids vs. Fitted: Parallel Lines Model")
text(parallel.mod$fitted,parallel.mod$residuals+.07,labels)
abline(0,0)
##########Diagnostic Plots
par(mfrow=c(4,1))
plot(seq(1:length(parallel.mod.diag$cooks)),parallel.mod.diag$cooks,main="Cooks Distance, Parallel Lines Model")
plot(seq(1:length(parallel.mod.diag$hat)),parallel.mod.diag$hat,main="Leverage, Parallel Lines Model")
plot(seq(1:length(parallel.mod.diag$stud.res)),abs(parallel.mod.diag$stud.res),main="Abs. Student. Res, Parallel Lines Model")
plot(seq(1:length(parallel.mod.diag$dfits)),parallel.mod.diag$dfits,main="DFFITS, Parallel Lines Model")
##################