Subject: Re: [R] (-1 as index) OR (envelope for QQ)
From: Peter Dalgaard BSA (p.dalgaard@biostat.ku.dk)
Date: Fri 25 Feb 2000 - 06:17:39 EST
Message-ID: <x2k8junyn0.fsf@blueberry.kubism.ku.dk>
"Dan E. Kelley" <kelley@Phys.Ocean.Dal.CA> writes:
> Content-Type: TEXT/PLAIN; charset=US-ASCII
>
> I'm new to R (and to S) and am wondering about code from pages 72 and
> 83 of MASS (Venables+Ripley, 3rd edition), to draw an envelope on a QQ
> plot. Copying from the book, I've got:
>
> #... code whose gist is "a.fit <- nls(..."
> num.points <- length(resid(a.fit))
> qqnorm(residuals(a.fit)) # illustrate data-model residuals
> qqline(residuals(a.fit))
> samp <- cbind(residuals(a.fit), matrix(rnorm(num.points*19),num.points,19))
> samp <- apply(scale(samp), 2, sort)
> rs <- samp[,1]
> xs <- qqnorm(rs, plot=FALSE)$x
> env <- t(apply(samp[,-1], 1, range)) ###########################3
> xyul <- par("usr")
> smidge <- min(diff(c(xyul[1], xs, xyul[2])))/2
> segments(xs-smidge,env[,1], xs+smidge, env[,1])
> segments(xs-smidge,env[,2], xs+smidge, env[,2])
>
....
> 1) You idiot! QQ plots are supposed to have standardized residuals of
> the data-model misfit. Do that and the QQ envelop is ok.
>
> 2) Don't copy S stuff into R code. The S code "samp[,-1]" from VR
> should be written "..." for use in R.
>
> 3) Don't be so rude as to ask such questions (with length attachments,
> no less) on a group like this. Give us a break!
>
> ... but I'm hoping that the mode won't be in the final category!
Well, 2) is *not* the answer... Try removing the plot=F from the 2nd
qqnorm and I think you'll see the light. Notice the scale on the y
axis.
-- O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk) FAX: (+45) 35327907 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
This archive was generated by hypermail 2b25 : Mon 17 Jul 2000 - 12:33:09 EST