# residuals and predicted values for survreg
# for more information, try help(residuals.survreg)
 # help(predict.survreg)

insa.wei <- survreg(Surv(t,cens)~kV,dist="weib") # -300.5  -339.7
insa.exp <- survreg(Surv(t,cens)~kV,dist="exp") # -304.9  -424.9
insa.ll <- survreg(Surv(t,cens)~kV,dist="logl") # -303.4 -332.7
insa.ln <- survreg(Surv(t,cens)~kV,dist="logn")  # -303.7 -331.1

 dresid.ins.exp <- residuals(insa.exp, type="deviance")
# get the deviance residuals 

  postscript("Qins.exp.ps")  # qq plot, with added line. 
qqnorm( dresid.ins.exp)  ; qqline(dresid.ins.exp)
 dev.off()

 dresid.ins.wei <- residuals(insa.wei, type="deviance")
 postscript("Qins.wei.ps"); qqnorm( dresid.ins.wei, main="QQ normal plot for Weibull fit, insulation")  
 qqline(dresid.ins.wei)
dev.off()

pred.ins.wei <-  predict(insa.wei,type="lp")
# lienar predictor or risk score.

 postscript("Rvfit.wei.ps")
plot(pred.ins.wei,dresid.ins.wei, main="residuals vs linear predictor, insulation ")
dev.off()

predr.ins.wei <- predict(insa.wei,type="response")

plot(predr.ins.wei,dresid.ins.wei, main="residuals vs time, insulation")

 postscript("Rv-kV.wei.ps")
plot(kV,dresid.ins.wei, main="residuals vs kV , insulation ")