Displaying 2 results from an estimated 2 matches for "hadza".
2011 Jun 28
1
plotting survival curves with model parameters
...) # reproduces Figure 7.4 from Wood et al. (2002)
}
siler()
#########################################################################
How can I modify the function so that I can plot curves using
published Siler parameters I have already compiled into a dataframe
(below)?
names<-c("Hadza","Ache")
a1<-c(0.351,0.157)
b1<-c(0.895,0.721)
a2<-c(0.011,0.013)
a3<-c(0.0000067,0.000048)
b3<-c(0.125,0.103)
sil.anthro<-data.frame(cbind(names,a1,b1,a2,a3,b3))
For example, how can I modify the function above to produce a curve
for the a specific group (e.g., Ha...
2011 Jul 05
3
plotting survival curves (multiple curves on single graph)
...return(S.t)
#return(h.t)
}
t<-seq(0,90,1)
plot(t,sil(t),ylim=c(0,1),type='l',cex.lab=0.8,cex.axis=0.75,ylab='S(t)',xlab='Age
(years)')
}
with(hazanth[1,3:7],silsurv(a1=a1,b1=b1,a2=a2,a3=a3,b3=b3));title(main=hazanth[1,1],cex.main=0.9)
# plot for Hadza
with(hazanth[2,3:7],silsurv(a1=a1,b1=b1,a2=a2,a3=a3,b3=b3));title(main=hazanth[2,1],cex.main=0.9)
# plot for Ache
with(hazanth[3,3:7],silsurv(a1=a1,b1=b1,a2=a2,a3=a3,b3=b3));title(main=hazanth[3,1],cex.main=0.9)
# plot for Hiwi
with(hazanth[4,3:7],silsurv(a1=a1,b1=b1,a2=a2,a3=a3,b3=b3));title(mai...