Re: [R] p-value of the pooled Z score

2012-03-13 Thread Thomas Lumley
On Wed, Mar 14, 2012 at 3:41 AM, cheba meier wrote: > Hello, > > I have to compute the pooled z-value and I would like to know which way is > more appropriate > > > b <- c( -0.205,1.040,0.087) > s <- c(0.449,0.167,0.241) > n <- c(310, 342, 348) > z <- b/s > > Z <- sum(z)/sqrt(length(n)) > P <- 2*(

[R] p-value of the pooled Z score

2012-03-13 Thread cheba meier
Hello, I have to compute the pooled z-value and I would like to know which way is more appropriate b <- c( -0.205,1.040,0.087) s <- c(0.449,0.167,0.241) n <- c(310, 342, 348) z <- b/s Z <- sum(z)/sqrt(length(n)) P <- 2*(1-pnorm(abs(Z))) P w <- sqrt(n) Zw <- sum(w * z)/sqrt(sum(w^2)) Pw <- 1 -