spTimer (version 0.02)

spT.pCOVER: Nominal Coverage

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, zsamples=y)
spT.pCOVER(z=x, zsamples=y, level=50)

##