data(UPBdata)
weightData <- neWeight(negaff ~ att + gender + educ + age,
data = UPBdata, nRep = 2)
## extract regression weights for natural effect model
head(weights(weightData))
## assess normality
qqnorm(residuals(weightData))
## assess model adequacy
library(car)
residualPlots(weightData)
Run the code above in your browser using DataLab