inst/snippet/lm-star-dfbeta.R

In fastR: Foundations and Applications of Statistics Using R

star.plot3 <- xyplot(dfbeta(star.model2)[,'temp']~index,
    data=hotstar,
    ylab="DFBETA",
    panel=function(x,y,...) {
        ids <- which(abs(y) > 0.5)
        panel.xyplot(x,y,...)
        grid.text(
            x=x[ids]+1.5, y = y[ids],
            as.character(ids), default.units="native")
    })
coef(lm(light~temp,hotstar))
coef(lm(light~temp,hotstar[-7,]))