Description Usage Arguments Value Author(s) Examples
View source: R/GPRMortalitySummary.R
This function summarize the results by percentiles.
This function works for both child and age-sex specific mortality rates.
1 | GPRMortalitySummary(model,percentile=c(0.025,0.5,0.975))
|
model |
a model object. |
percentile |
a vector of uncertainty interval percentile. The default is 0.025, 0.5 and 0.975. |
a matrix of GPR result including percentile.
Parinaz Mehdipour, Ali Ghanbari, Iman Navidi
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | library("rstan")
library("GPRMortality")
head(data.mortality)
head(data.mean)
mortality <- data.mortality[data.mortality$location%in%c(0,5) &
data.mortality$age_cat%in%c(1,10) & data.mortality$sex%in%c(0,1),]
mean <- data.mean[data.mean$location%in%c(0,5) &
data.mean$age_cat%in%c(1,10) & data.mean$sex%in%c(0,1),]
# WARNING: The following code will take a long time to run
fit = GPRMortality(mortality,mean,minYear = 1990,maxYear = 2015,
nu = 2,rho_ =0.4 ,product = 0.1 ,verbose=TRUE)
####### summary
fit_sum = GPRMortalitySummary(fit)
fit_sum
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.