Displaying 1 result from an estimated 1 matches for "paleodemographers".
2012 Apr 16
0
Gompertz-Makeham hazard models---test for significant difference
...),.Dimnames=list(NULL,c("col1","col2","col3")))
> naga
col1 col2 col3
[1,] 15 20 46
[2,] 20 35 43
[3,] 35 50 26
[4,] 50 Inf 14
So above, the first two columns are the lower and upper limits of the
age-range categories (typically used by paleodemographers), and the third
column is the number of skeletons assigned to each group. I've co-opted
some R script for fitting a Gompertz-Makeham hazard model to this data...
##############################
GM.naga <- function(x,deaths=naga)
{
a2=x[1]
a3=x[2]
b3=x[3]
shift<-15
nrow<-NRO...