R/medciMeeker.R

In RMediation: Mediation Analysis Confidence Intervals

medciMeeker <-
function(mu.x,mu.y,se.x,se.y,rho=0,alpha=.05,...){
  p <- alpha/2
  q.l <- qprodnormalMeeker(p, mu.x=mu.x, mu.y=mu.y, se.x=se.x, se.y=se.y, rho=rho, lower.tail=TRUE)
  q.u <- qprodnormalMeeker(p, mu.x=mu.x, mu.y=mu.y, se.x=se.x, se.y=se.y, rho=rho, lower.tail=FALSE)
  CI <- c(q.l[[1]],q.u[[1]]) #confidence interval
  quantMean <- mu.x*mu.y + se.x*se.y*rho
  quantSE <- sqrt(se.y^2*mu.x^2+se.x^2*mu.y^2+2*mu.x*mu.y*rho*se.x*se.y+se.x^2*se.y^2+se.x^2*se.y^2*rho^2);
  res <- list(CI,quantMean,quantSE)
  names(res) <- c(paste((1-alpha)*100,"% ","CI",sep=""), "Estimate", "SE")
  return(res)
}