################################################################## ### Example 9.7: Linearization ################################################################## # Generate Fig 9.4 n<-1000 uvals<-(1:n)/(n+1) par(mfrow=c(2,2)) qqplot(1+0.02*qnorm(uvals),exp(0.02*qnorm(uvals)),xlab="",ylab="",main="",type="l") lines(1+0.02*qnorm(uvals),1+0.02*qnorm(uvals),lty=3) qqplot(1+0.3*qnorm(uvals),exp(0.3*qnorm(uvals)),xlab="",ylab="",main="",type="l") lines(1+0.3*qnorm(uvals),1+0.3*qnorm(uvals),lty=3) qqplot(1+0.02/sqrt(3)*qt(uvals,3),exp(0.02/sqrt(3)*qt(uvals,3)),xlab="",ylab="",main="",type="l") lines(1+0.02/sqrt(3)*qt(uvals,3),1+0.02/sqrt(3)*qt(uvals,3),lty=3) qqplot(1+0.3/sqrt(3)*qt(uvals,3),exp(0.3/sqrt(3)*qt(uvals,3)),xlab="",ylab="",main="",type="l") lines(1+0.3/sqrt(3)*qt(uvals,3),1+0.3/sqrt(3)*qt(uvals,3),lty=3)