##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function (residuals, MSE)
{
w <- 1:length(residuals)
r <- cor(sort(residuals), sqrt(MSE) * (qnorm((w - 0.375)/(length(residuals) +
0.25))))
return(r)
}
Run the code above in your browser using DataLab