R/pPval.r

In selectMeta: Estimation of Weight Functions in Meta Analysis

Defines functions pPval

Documented in pPval

pPval <- function(q, u, theta, sigma2){
    res <- rep(NA, length(q))
    for (j in 1:length(q)){
	   res[j] <- try(integrate(dPval, lower = 0, upper = q[j], u, theta, sigma2)$value, silent = TRUE)
	   if(is.numeric(res[j]) == FALSE){res[j] <- 0.5}
	  }
    return(res)
}

Try the selectMeta package in your browser