# NOT RUN {
#Two groups
n <- 3
A <- c(1,10,30)
beta <- diag(1,n)
beta[1,2] <- -0.32
beta[2,1] <- 1/beta[1,2]
B <- compute.B.from.beta(beta)
L_Phi <- class_group(beta)
correl <- "RegNeg"
#apply function of all groups
A_app <- unlist(lapply(L_Phi,apparent.activities.Aq,A,B,correl))
# }
Run the code above in your browser using DataLab