affairs | R Documentation |
Data from Fair (1978). Although Fair used a tobit model with the data, the outcome measure can be modeled as a count. In fact, Greene (2003) modeled it as Poisson, but given the amount of overdispersion in the data, employing a negative binomial model is an appropriate strategy. The data is stored in the affairs data set. Naffairs is the response variable, indicating the number of affairs reported by the participant in the past year.
data(affairs)
A data frame with 601 observations on the following 18 variables.
naffairs
number of affairs within last year
kids
1=have children;0= no children
vryunhap
(1/0) very unhappily married
unhap
(1/0) unhappily married
avgmarr
(1/0) average married
hapavg
(1/0) happily married
vryhap
(1/0) very happily married
antirel
(1/0) anti religious
notrel
(1/0) not religious
slghtrel
(1/0) slightly religious
smerel
(1/0) somewhat religious
vryrel
(1/0) very religious
yrsmarr1
(1/0) >0.75 yrs
yrsmarr2
(1/0) >1.5 yrs
yrsmarr3
(1/0) >4.0 yrs
yrsmarr4
(1/0) >7.0 yrs
yrsmarr5
(1/0) >10.0 yrs
yrsmarr6
(1/0) >15.0 yrs
rwm5yr is saved as a data frame. Count models use naffairs as response variable. 0 counts are included.
Fair, R. (1978). A Theory of Extramarital Affairs, Journal of Political Economy, 86: 45-61. Greene, W.H. (2003). Econometric Analysis, Fifth Edition, New York: Macmillan.
Hilbe, Joseph M (2011), Negative Binomial Regression, Cambridge University Press Hilbe, Joseph M (2009), Logistic regression Models, Chapman & Hall/CRC
data(affairs) glmaffp <- glm(naffairs ~ kids + yrsmarr2 + yrsmarr3 + yrsmarr4 + yrsmarr5, family = poisson, data = affairs) summary(glmaffp) exp(coef(glmaffp)) require(MASS) glmaffnb <- glm.nb(naffairs ~ kids + yrsmarr2 + yrsmarr3 + yrsmarr4 + yrsmarr5, data=affairs) summary(glmaffnb) exp(coef(glmaffnb))