my.2x2.dev <- function(p,q,r){ #tmp <- r*sqrt(p*q*(1-p)*(1-q)) x <- c(p*q,p*(1-q),(1-p)*q,(1-p)*(1-q)) tmp <- r*min(x) matrix(x+c(tmp,-tmp,-tmp,tmp),byrow=TRUE,2,2) } my.2x2x2.dev <- function(p,q,r){ tmp <- r*sqrt(prod(p)*prod(q)*r) tmp1 <…