R/pb.score.R

In altmeta: Alternative Meta-Analysis Methods

pb.score <- function(y = NA, s2 = NA, n = NA, n00, n01, n10, n11){
  n <- n00 + n01 + n10 + n11
  nSum0 <- n00 + n10
  nSum1 <- n01 + n11
  n0Sum <- n00 + n01
  n1Sum <- n10 + n11
  Z <- n11 - nSum1/n*n1Sum
  V <- n0Sum/n*n1Sum/n*nSum0/(n - 1)*nSum1
  y.reg <- Z/V
  x.reg <- 1/sqrt(V)
  w.reg <- V
  if(all(abs(diff(x.reg[!is.na(x.reg)])) < 1e-5) | all(abs(diff(y.reg[!is.na(y.reg)])) < 1e-5)){
    out <- list(pval = 1, stat = 0, coef = c(0, 0))
  }else{
    out <- lm(y.reg ~ x.reg, weights = w.reg)
    out <- summary(out)$coefficients
    pval <- out["x.reg", "Pr(>|t|)"]
    coef <- as.numeric(out[,1])
    out <- list(pval = pval, stat = as.numeric(out["x.reg", "t value"]), coef = coef)
  }
  return(out)
}