cricketer | R Documentation |
Year and birth, lifespan, etc, of British first class cricketers, born 1840-1960, whose handedness could be determined from information in the Who's who of cricketers. The status (alive=0, dead =1), and lifetime or lifespan, is for 1992.
data(cricketer)
A data frame with 5960 observations on the following 8 variables.
left
a factor with levels right
left
year
numeric, year of birth
life
numeric, lifetime or lifespan to 1992
dead
numeric (0 = alive (censored), 1 = dead, in 1992)
acd
numeric (0 = not accidental or not dead, 1 = accidental death)
kia
numeric (0 = not killed in action, 1 = killed in action)
inbed
numeric (0 = did not die in bed, 1 = died in bed)
cause
a factor with levels 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. Data were collated as described in Aggleton et al.
Aggleton JP, Bland JM, Kentridge RW, Neave NJ 1994. Handedness and longevity: an archival study of cricketers. British Medical Journal 309, 1681-1684.
Bailey P, Thorne P, Wynne-Thomas P. 1993. Who's Who of Cricketers. 2nd ed, London, Hamlyn.
Bland M and Altman D. 2005. Do the left-handed die young? Significance 2, 166-170.
earlycrcktr
.
data(cricketer)
numLH <- xtabs(~ left+year, data=cricketer)
propLH <- prop.table(numLH, margin=2)[2,]
yr <- as.numeric(colnames(numLH))
plot(propLH ~ yr)
cricketer$lh <- unclass(cricketer$left)-1
left2.hat <- fitted(lm(lh ~ poly(year,2), data=cricketer))
ord <- order(cricketer$year)
lines(left2.hat[ord] ~ cricketer$year[ord])
library(splines)
ns3.hat <- fitted(lm(lh ~ ns(year,3), data=cricketer))
lines(ns3.hat[ord] ~ cricketer$year[ord], col="red")
require(survival)
summary(coxph(Surv(life, kia) ~ bs(year,3) +left, data=cricketer))
cricketer$notacdDead <- with(cricketer, {dead[acd==1]<-0; dead})
summary(coxph(Surv(life, notacdDead) ~ ns(year,2) +left, data=cricketer))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.