[R] outlier identify in qqplot

R. Michael Weylandt michael.weylandt at gmail.com
Thu May 10 19:54:36 CEST 2012


Do you want a qqnorm() instead of a qqplot() ? [Reproducibility also
involves posting the code you used that lead to the error / warning]

This seems to work for me:

# mydata <- source("http://r.789695.n4.nabble.com/file/n4623493/mydata.txt")[[1]]
# Have to drop visible attribute

lmmodel <- lm(log(vdep) ~ v1 + sqrt(v2) + v3 +v5 + v6 + v7 + v8 + v9 +
v10, data = mydata)

qqnorm(residuals(lmmodel))

# Or if interactive:

qqnormInt <- function(..., IDENTIFY = TRUE){
     qqnorm(...) -> X
     if(IDENTIFY) return(identify(X))
     invisisble(X)
}


qqnormInt(residuals(lmmodel))

Michael


On Thu, May 10, 2012 at 9:20 AM, agent dunham <crosspide at hotmail.com> wrote:
> Find the data attached,
>
> http://r.789695.n4.nabble.com/file/n4623493/mydata.txt mydata.txt
>
> The model would be /lmmodel <- lm(log(vdep) ~ v1 + sqrt(v2) + v3 +v5 + v6 +
> v7 + v8 + v9 + v10, data = mydata)/
>
> Thanks again,
>
>
> user at host.com
>
> --
> View this message in context: http://r.789695.n4.nabble.com/outlier-identify-in-qqplot-tp4076587p4623493.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.



More information about the R-help mailing list