GPRMortalitySummary {GPRMortality}

summary table for Gaussian Process Regression results

Description

This function summarize the results by percentiles.

This function works for both child and age-sex specific mortality rates.

Usage

GPRMortalitySummary(model,percentile=c(0.025,0.5,0.975))

Arguments

model

a model object.

percentile

a vector of uncertainty interval percentile. The default is 0.025, 0.5 and 0.975.

Value

a matrix of GPR result including percentile.

Author(s)

Parinaz Mehdipour, Ali Ghanbari, Iman Navidi

Examples

 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