R/MetaFEE2.R

In YandaLang/RankMeta: Rank-based Mata-analysis

#' Meta-analysis for fixed-effect - E2
#' @description Meta-analysis for fixed-effect - E2
#' @usage
#' MetaFEE2(info,alpha=0.05)
#' @param info summary statistics
#' @param alpha significance level
#' @author Yanda Lang, Joseph McKean
#' @note format for info (2SHL,sample size 1,sample size 2)
#' @examples
#' study1 <- c(-80.0, 17, 16)
#' study2 <- c(-89.0, 23, 23)
#' info <- rbind(study1,study2)
#' MetaFEE2(info,alpha=0.05)
#' @keywords robust
#' @keywords Meta-analysis

#' @import stats Rfit
#' @export

MetaFEE2 = function(info,alpha=0.05){
  # Define h studies split into K pair of studies
  K = nrow(info)
  # Define median diff, sample size
  deltahat = info[,1]; n1 = info[,2]; n2 = info[,3]

  # The estimate of delta using median of walsh averages of differences
  deltatilde = median(walsh(deltahat))
  D = sort(walsh(deltahat))

  # proposed combined CI
  c = round(K*(K+1)/4-0.5-qnorm(1-alpha/2)*sqrt(K*(K+1)*(2*K+1)/24))
  lower = D[c+1]
  upper = D[K*(K+1)/2-c]

  # collect values
  combineCI = c(lower, upper)
  return(combineCI)
}