Displaying 2 results from an estimated 2 matches for "g1g2".
Did you mean:
g12
2010 Mar 04
1
ifthen() question
...ue when length of unique id == 1 but I cannot get it to
work. e.g.:
#function to aggregate effect sizes:
aggs <- function(g, n.1, n.2, cor = .50) {
n.1 <- mean(n.1)
n.2 <- mean(n.2)
N_ES <- length(g)
corr.mat <- matrix (rep(cor, N_ES^2), nrow=N_ES)
diag(corr.mat) <- 1
g1g2 <- cbind(g) %*% g
PSI <- (8*corr.mat + g1g2*corr.mat^2)/(2*(n.1+n.2))
PSI.inv <- solve(PSI)
a <- rowSums(PSI.inv)/sum(PSI.inv)
var.g <- 1/sum(PSI.inv)
g <- sum(g*a)
out<-cbind(g,var.g, n.1, n.2)
return(out)
}
# automating this procedure for all rows of df. Thi...
2010 Mar 04
0
ifthen() question -- whoops--ifelse()
...it to
> work. e.g.:
>
> #function to aggregate effect sizes:
> aggs <- function(g, n.1, n.2, cor = .50) {
> n.1 <- mean(n.1)
> n.2 <- mean(n.2)
> N_ES <- length(g)
> corr.mat <- matrix (rep(cor, N_ES^2), nrow=N_ES)
> diag(corr.mat) <- 1
> g1g2 <- cbind(g) %*% g
> PSI <- (8*corr.mat + g1g2*corr.mat^2)/(2*(n.1+n.2))
> PSI.inv <- solve(PSI)
> a <- rowSums(PSI.inv)/sum(PSI.inv)
> var.g <- 1/sum(PSI.inv)
> g <- sum(g*a)
> out<-cbind(g,var.g, n.1, n.2)
> return(out)
> }
>
>
>...