Displaying 1 result from an estimated 1 matches for "xyul".
Did you mean:
xul
2000 Feb 24
2
(-1 as index) OR (envelope for QQ)
...ata-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])
where I've marked a confusing line with ########. From what I gather
(from section 2.1, or page 5, of "R comp...