| 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.
leftright leftyearlifedeadacdkiainbedcausealive 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))