The data set lrb.csv has data on 61 married couples.

The dependent variable is HAPPYMAR, respondent's marital happiness (1 = Very Happy, 0 = Otherwise).

CHURCH, Church attendance (1 = Often attends, 0 = other), FEMALE (1 = female, 0 = male), and EDUC, Years of education, are the explanatory variables.

mydata=read.csv("lrb.csv")

names(mydata) ##

[1] "happymar" "church" "female" "educ"

a) Run the logistic regression of happymar on church, female and educ

b) Which of the three explanatory variables are needed in the model? Peform the appropriate hypothesis tests.

c) Based on the signs of the logistic regression coefficients (positive or negative), interpret their impact on the probability of a happy marriage. Any surprises?

d) Determine the probability of marital happiness for: (a) a male with 8 years of education who is not a regular churchgoer (b) a male with 8 years of education who is a regular churchgoer (c) a female with 16 years of education who is not a regular churchgoer (d) a female with 16 years of education who is a regular churchgoer.