Rob Forsyth
2013-Jun-10 21:41 UTC
[R] padding specific missing values with NA to allow cbind
Dear list Getting very frustrated with this simple-looking problem> m1 <- lm(x~y, data=mydata) > outliers <- abs(stdres(m1))>2 > plot(x~y, data=mydata)I would like to plot a simple x,y scatter plot with labels giving custom information displayed for the outliers only, i.e. I would like to define a column mydata$labels for the mydata dataframe so that the command> text(mydata$y, mydata$x, labels=mydata$labels)will label those rows where outliers[i] = TRUE with text but is otherwise blank The first problem I have is that due to some NAs in mydata, nrows(outliers) < nrows(mydata) and I'm getting in a tangle trying to pad the appropriate rows of outliers Thanks Rob
William Dunlap
2013-Jun-11 03:59 UTC
[R] padding specific missing values with NA to allow cbind
Try adding the argument
na.action = na.exclude
to your call to lm(). See help("na.exclude") for details.
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at
r-project.org] On Behalf
> Of Rob Forsyth
> Sent: Monday, June 10, 2013 2:42 PM
> To: r-help at r-project.org
> Subject: [R] padding specific missing values with NA to allow cbind
>
> Dear list
>
> Getting very frustrated with this simple-looking problem
>
> > m1 <- lm(x~y, data=mydata)
> > outliers <- abs(stdres(m1))>2
> > plot(x~y, data=mydata)
>
> I would like to plot a simple x,y scatter plot with labels giving custom
information
> displayed for the outliers only, i.e. I would like to define a column
mydata$labels for the
> mydata dataframe so that the command
>
> > text(mydata$y, mydata$x, labels=mydata$labels)
>
> will label those rows where outliers[i] = TRUE with text but is otherwise
blank
>
> The first problem I have is that due to some NAs in mydata, nrows(outliers)
<
> nrows(mydata) and I'm getting in a tangle trying to pad the appropriate
rows of outliers
>
> Thanks
>
> Rob
>
> ______________________________________________
> 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.