spT.pCOVER: Nominal Coverage

In jmarca/spTimer: Spatio-Temporal Bayesian Modelling

Description

This function is used to obtain nominal coverage.

Usage

spT.pCOVER(z=NULL,zup=NULL,zlow=NULL,zsample=NULL,level=95)

Arguments

z	The original values (matrix or vector).
zup	The predicted values for upper interval (matrix or vector).
zlow	The predicted values for lower interval (matrix or vector).
zsample	Predicted MCMC samples.
level	Level of coverages.

See Also

spT.validation.

Examples

##

# Create `x': the true values.
# Create `yup': the upper interval.
# Create `ylow': the lower interval.

x <- rnorm(1000,5,0.1)
yup <- rnorm(1000,7,2)
ylow <- rnorm(1000,3,2)
	
# The pCOVER is:

spT.pCOVER(z=x, zup=yup, zlow=ylow)

# create predicted MCMC samples

y <- matrix(rnorm(1000*5000,5,1),1000,5000)

# The pCOVER is:

spT.pCOVER(z=x, zsample=y)
spT.pCOVER(z=x, zsample=y, level=50)

##  

jmarca/spTimer documentation built on May 19, 2019, 1:51 p.m.