Example contingency table

conting.R

#--------------------------------------------------
# Contingency table
#--------------------------------------------------
# Productivity Smoking   No smoking
# High         13        44
# Low          25        29

count<-c(13,25,44,29)
smoking<-factor(c("Yes","Yes","No","No"))
productivity<-factor(c("High","High","Low","Low"))
model<- glm(count~ smoking+productivity, poisson)
fitted(model)
countsummary(model)
model2<- update(model,~.-productivity, test="Chi")

1-pchisq(<valor-de-chi>,<df>)
#--------------------------------------------------