| Class | Analysis and data manipulation command |
| Name | kaplan-meier |
| Arguments | <age-at-onset> <censor> [residuals] |
Prints the product-limit estimator for the survivor
function for the quantitative trait age-at-onset, where censor
is the binary outcome trait, which is affected when
age-at-onset represents the age at which the individual first expressed
the trait. The age-at-onset is replaced by a nonparametric residual
when requested. If affected, this is:
sgn(1-H(t)).(-2(1-H(t)+log(H(t))))1/2
If unaffected:
-(-2H(t))1/2
where H(t) is the Nelson-Aalen estimate of the integrated
hazard function at that age t.
Example:
>> include liuex.in
>> kaplan-meier onset ad
------------------------------------------------
Kaplan-Meier survivor function for "onset"
------------------------------------------------
"ad" is outcome (censoring) trait.
Age-at-onset Failed Riskset H(t) S(t) ase
---------------------------------------------------
52.0000 1 108 0.0093 0.9907 0.0092
57.0000 1 95 0.0198 0.9803 0.0138
58.0000 1 92 0.0307 0.9697 0.0173
59.0000 1 91 0.0416 0.9590 0.0201
60.0000 4 88 0.0871 0.9154 0.0287
61.0000 2 81 0.1118 0.8928 0.0321
62.0000 1 76 0.1249 0.8811 0.0338
65.0000 5 71 0.1954 0.8190 0.0412
67.0000 4 64 0.2579 0.7678 0.0459
68.0000 1 56 0.2757 0.7541 0.0471
69.0000 3 55 0.3303 0.7130 0.0502
70.0000 4 49 0.4119 0.6548 0.0539
72.0000 5 36 0.5508 0.5638 0.0598
73.0000 1 29 0.5853 0.5444 0.0608
76.0000 5 23 0.8027 0.4260 0.0668
79.0000 2 15 0.9360 0.3692 0.0689
80.0000 1 12 1.0193 0.3385 0.0697
83.0000 1 7 1.1622 0.2901 0.0746
H(t) = Nelson-Aalen estimator of integrated hazard
S(t) = Kaplan-Meier estimator of survivor function
43 affecteds and 119 unaffecteds used
| << (impute) | Up to index | >> (rank) |