upperT: An utility function to calculate lower limit of T statistic

In marp: Model-Averaged Renewal Process

An utility function to calculate lower limit of T statistic

Description

An utility function to calculate lower limit of T statistic

Usage

upperT(up, hat, sigmasq, Tstar, weights, B, alpha)

Arguments

up	upper limit
hat	estimates
sigmasq	variance
Tstar	T statistics estimated from bootstrap samples
weights	model weights
B	number of bootstraps
alpha	confidence level

Value

returns lower limit of T statistic

Examples

# set some parameters
up <- 100 # upper bound
hat <- rep(150, 6) # estimates obtained from each model
sigmasq <- 10 # variance
Tstar <- matrix(rep(100,600),6,100) # T statistics estimated from bootstrap samples
weights <- rep(1/6, 6) # model weights
B <- 100 # number of bootstrapped samples
alpha <- 0.05 # confidence level

# calculate the upper limit of T statistics
res <-  marp::upperT(up, hat, sigmasq, Tstar, weights, B, alpha)

# print result
cat("res = ", res, "\n")

marp documentation built on Aug. 11, 2022, 5:10 p.m.