Hi, all. I've written a function that returns the survival function for a Gompertz mortality model. I've specified the two model parameters. Using a simple integration, I can calculate the life expectancy at any age. Is there a way I can use a loop with the integration that will quickly return life expectancy over a range of ages, say 0 to 80, so that I don't have to manually type in the age in which I'm interested? Please see the code below. Thanks so much.
--Trey hk.bothsex_Gompsurv <- function (t) { x=c(0.02342671, 0.05837508) a3<-x[1] b3<-x[2] shift<-15 S.t<-exp(a3/b3*(1-exp(b3*(t-shift)))) return<-S.t } integrate(hk.bothsex_Gompsurv,0,Inf)$value/hk.bothsex_Gompsurv(0) # life expectancy at birth (change lower limit of integral and corresponding "t" in denominator to calculate life expectancy at any age ----- Trey Batey---Anthropology Instructor Division of Social Sciences Mt. Hood Community College Gresham, OR 97030 Alt. Email: trey.batey[at]mhcc[dot]edu -- View this message in context: http://r.789695.n4.nabble.com/using-a-loop-with-an-integration-tp4578752p4578752.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.