R/AveS20w.R

"AveS20w" <-
function(v,args){
Beta <- args$Beta; n <- args$n
nr   <- length(args$rr); nv <- length(v); A <- rep(0,nv)
for (i in 1:nv) {
rs   <- args$rr/v[i] 
tmp  <- rs*(exp(rs)-1)
A[i] <- sum(tmp)/n-Beta}; A}

Try the RobustAFT package in your browser

Any scripts or data that you put into this service are public.

RobustAFT documentation built on Aug. 21, 2023, 5:13 p.m.