Displaying 1 result from an estimated 1 matches for "narepwt".
2009 Mar 05
1
hatvalues?
...e help in understanding the theory. First, I take the simplest case of a single variant. For this I turn o John Fox's book, "Applied Regression Analysis and Generalized Linear Models, p 245 and generate this 'R' code:
> library(car)
> attach(Davis)
# remove the NA's
> narepwt <- repwt[!is.na(repwt)]
> meanrw <- mean(narepwt)
> drw <- narepwt - meanrw
> ssrw <- sum(drw * drw)
> h <- 1/length(narepwt) + (drw * drw)/ssrw
> h
This gives me a array of values the largest of which is
> order(h, decreasing=TRUE)
[1] 21 52 17 93 30 62 1...