如何用 R 求联合概率

新生刚接触 R，想问下如何用 R 求联合概率。
已知X 服从正态分布，
Y 服从 gamma分布，

求解联合概率使得 ：X 不小于 a且Y 不小于 b的概率
以及 X 不大于 a，Y 不大于 b 的联合概率。

还请赐教代码。

参数随意，只想知道代码怎么写。谢谢大家

由边际分布不能决定联合分布，得要假设条件分布
假设X~N(0,1),Y服从条件分布Y|X ~ Gamma(|x|, b)。

1. set.seed(123)
   
2. 
   
3. a <- 1.5
   
4. b <- 1
   
5. 
   
6. # generate 1 million obs for x and y|x
   
7. len <- 1e6
   
8. x <- rnorm(len, mean = 0, sd = 1)
   
9. y <- numeric(len)
   
10. for (i in 1:len) {
    
11.     y[i] <- rgamma(1, shape = abs(x[i]), rate = 1)
    
12. }
    
13. 
    
14. prob1 <- sum(x>a & y>b)/len
    
15. prob2 <- sum(x<a & y<b)/len

复制代码

您好，请问您的问题解决了吗？我现在也困在这一步了