Learn R Programming

COUNT (version 1.3.2)

titanicgrp: titanicgrp

Description

The data is an grouped version of the 1912 Titanic passenger survival log,

Usage

data(titanicgrp)

Arguments

Format

A data frame with 12 observations on the following 5 variables.
survive
number of passengers who survived
cases
number of passengers with same pattern of covariates
age
1=adult; 0=child
sex
1=Male; 0=female
class
ticket class 1= 1st class; 2= second class; 3= third class

Source

Found in many other texts

Details

titanicgrp is saved as a data frame. Used to assess risk ratios

References

Hilbe, Joseph M (2014), Modeling Count Data, Cambridge University Press Hilbe, Joseph M (2007, 2011), Negative Binomial Regression, Cambridge University Press Hilbe, Joseph M (2009), Logistic Regression Models, Chapman & Hall/CRC

Examples

Run this code
library(MASS)
library(msme)
data(titanicgrp)
glmlr <- glm(survive ~ age + sex + factor(class) + offset(log(cases)),
             family=poisson, data=titanicgrp)
summary(glmlr)
exp(coef(glmlr))

lcases <- titanicgrp$cases
nb2o <- nbinomial(survive ~ age + sex + factor(class), 
                                        formula2 =~ age + sex,
                                        offset = lcases,
                                        mean.link="log",
                                        scale.link="log_s",
                                        data=titanicgrp)
summary(nb2o)
exp(coef(nb2o))

Run the code above in your browser using DataLab