Displaying 2 results from an estimated 2 matches for "p_pj".
Did you mean:
p_j
2005 Apr 07
2
about mantelhaen.test (PR#7779)
...S_GMH <- matrix(0, nc = df_GMH, nr = df_GMH)
S_SMH <- matrix(0, nc = df_SMH, nr = df_SMH)
S_CSMH <- matrix(0, nc = df_CSMH, nr = df_CSMH)
for(k in 1:K) {
V <- NULL
f <- x[, , k]
ntot <- sum(f)
p_ip <- apply(f, 1, sum) / ntot
p_pj <- apply(f, 2, sum) / ntot
m <- p_ip %x% p_pj * ntot
V <- ntot^2 * ((diag(p_ip) - p_ip %*% t(p_ip)) %x% (diag(p_pj) - p_pj
%*% t(p_pj))) / (ntot-1)
Y_GMH <- Y_GMH + A_GMH %*% (c(t(f)) - m)
Y_SMH <- Y_SMH + A_SMH %*% (c(t(f)) - m)
Y_CSMH <- Y...
2006 Oct 31
0
about mantelhaen.test (PR#7779)
...- matrix(0, nc =3D df_GMH, nr =3D df_GMH)
S_SMH <- matrix(0, nc =3D df_SMH, nr =3D df_SMH)
S_CSMH <- matrix(0, nc =3D df_CSMH, nr =3D df_CSMH)
for(k in 1:K) {
V <- NULL
f <- x[, , k]
ntot <- sum(f)
p_ip <- apply(f, 1, sum) / ntot
p_pj <- apply(f, 2, sum) / ntot
m <- p_ip %x% p_pj * ntot
V <- ntot^2 * ((diag(p_ip) - p_ip %*% t(p_ip)) %x% (diag(p_pj) -
p_pj %*% t(p_pj))) / (ntot-1)
Y_GMH <- Y_GMH + A_GMH %*% (c(t(f)) - m)
Y_SMH <- Y_SMH + A_SMH %*% (c(t(f)) - m)
Y_CSMH <- Y...