problem with zero-weighted observations in predict.lm?
Prof Brian Ripley wrote:
I think you will find that 'n' is used in several ways in predict.lm, and since NA-handling was introduced in R 1.8.0 they may differ in value. So the safest route seems to be to change just 'n' in df <- n - p
Yes, that seems to fix things. Will commit to R-devel shortly. -p
Peter Dalgaard Center for Statistics, Copenhagen Business School Phone: (+45)38153501 Email: pd.mes at cbs.dk Priv: PDalgd at gmail.com