Skip to content

outlier identify in qqplot

4 messages · agent dunham, R. Michael Weylandt

agent dunham
# 
Sorry, but I need the same and i don't understand your help. 

So, after fitting my lm model, how can i identify my data? I was trying the
following, but it doesn't work.

/identify(qnorm(c(0.25,
0.75)),quantile(rstandard(mymodel)[!is.na(rstandard(mymodel))], c(0.25,
0.75)),row.names(mydata))
warning: no point within 0.25 inches/

/In identify.default(qnorm(c(0.25, 0.75)),
quantile(rstandard(lmt2t)[!is.na(rstandard(lmt2t))],  :
  more 'labels'  than points/

Thanks in advance, show user at host.com

--
View this message in context: http://r.789695.n4.nabble.com/outlier-identify-in-qqplot-tp4076587p4623088.html
Sent from the R help mailing list archive at Nabble.com.
R. Michael Weylandt
# 
Reproducible example? https://github.com/hadley/devtools/wiki/Reproducibility

Michael
On Thu, May 10, 2012 at 6:31 AM, agent dunham <crosspide at hotmail.com> wrote:
agent dunham
# 
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. Michael Weylandt
# 
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: