marp (version 0.1.0)

upperT: 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)

Value

returns lower limit of T statistic

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

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")