Displaying 1 result from an estimated 1 matches for "alargerthan3".
2012 May 03
1
NA's when subset in a dataframe
...aving this silly problem.
I've a linear model. After fixing it, I wanted to know which data had
studentized residuals larger than 3, so i tried this:
d1 <- cooks.distance(lmmodel)
r <- sqrt(abs(rstandard(lmmodel)))
rstu <- abs(rstudent(lmmodel))
a <- cbind( mydata, d1, r,rstu)
alargerthan3 <- a[rstu >3, ]
And suddenly a[rstu >3, ] has 17 rows, 7 of them are "new rows", where all
the entries are NA's, even its rownames.
Because of this I'm not sure of the dimension of a[rstu >3, ] (Do I only
have 8 entries?)
Has this happened to anybody before?...