我记得有人告诉我要生成变量 Marshall Olkin 分布，一个过程如下：$p$

对于，生成，然后 其中和Y_2=\max(X_1,X_3)。$p=2$$X_1,X_2,X_3\sim exp(\lambda)$$(Y_1,Y_2)\sim MMO$$Y_1=\max(X_1,X_2)$$Y_2=\max(X_1,X_3)$

但我一定是弄错了定义，因为这看起来不像“MMO”（因为沿$Y_1=Y_2$有很多观察，我不记得是这种情况）：

library(matrixStats)
n<-100
p<-2
a1<-matrix(rexp(n*(p+1),2),n,p+1)
a2<-cbind(rowMaxs(a1[,1:2]),rowMaxs(a1[,c(1,3)]))
plot(a2)

任何人都可以帮助生成 MMO 的算法吗？