calgspn(
xm = qnorm(matrix(rep(c(0.03,0.04,0.01),times=2),ncol=3,nrow=2)),
alpham = matrix(rep(c(0.02,0.03,0.05),each=2),ncol=3,nrow=2),
critm = matrix(rep(qnorm(c(0.02,0.03,0.05)),each=2),ncol=3,nrow=2),
matrix.list = list(diag(3),diag(3)),
sided = rep(-1,2)
)
Run the code above in your browser using DataLab