cricketers {DAAGxtras} | R Documentation |
Data were derived from the Who's who of cricketers, with year of birth up to 1960, lifespan or lifetime up to 1992.
data(cricketers)
A data frame with 6172 observations on the following 8 variables.
left
right
left
year
life
dead
acd
kia
inbed
cause
alive
acd
(accidental death) inbed
(died in bed)Note that those 'killed in action' (mostly during World Wars I and II) form a subset of those who died by accident.
John Aggleton, Martin Bland.
Bland M and Altman D. 2005. Do the left-handed die young? Significance 2, 166-170.
Aggleton JP, Bland JM, Kentridge RW, Neave NJ 1994. Handedness and longevity: an archival study of cricketers. British Medical Journal 309, 1681-1684.
data(cricketers) cricket <- subset(cricketers, year>=1840 & year <= 1960) numLH <- xtabs(~ left+year, data=cricket) propLH <- prop.table(numLH, margin=2)[2,] yr <- as.numeric(colnames(numLH)) plot(propLH ~ yr) cricket$lh <- unclass(cricket$left)-1 left2.hat <- fitted(lm(lh ~ poly(year,2), data=cricket)) ord <- order(cricket$year) lines(left2.hat[ord] ~ cricket$year[ord]) library(splines) ns3.hat <- fitted(lm(lh ~ ns(year,3), data=cricket)) lines(ns3.hat[ord] ~ cricket$year[ord], col="red") library(survival) summary(coxph(Surv(life, kia) ~ bs(year,3) +left, data=cricket)) cricket$notacdDead <- with(cricket, {dead[acd==1]<-0; dead}) summary(coxph(Surv(life, notacdDead) ~ ns(year,2) +left, data=cricket))