r - 在 R 中绘制顺序统计

Question

我需要绘制由以下定义的顺序统计的二元随机向量的联合密度函数：

f(x,y) = (n!)*f(x)*f(y) for x < y and 0 otherwise

其中f是单变量随机变量的密度函数。例如，f(x)=e^(-x)。

我不知道我是否必须使用persp以及如何修复有条件的支撑x<y。

非常感谢。胡里奥

Answer 1

这种绘图的典型工作流程是定义一个 x 和 y 网格，创建您的函数（在您的情况下，您将使用它ifelse来检查您的x < y条件），将函数应用到网格上outer，最后使用 绘制persp。

f <- function(x) exp(-x)
g <- function(x, y) ifelse(x < y, f(x)*f(y), 0)
x <- seq(0, 1, length.out=50)
y <- seq(0, 1, length.out=50)
z <- outer(x, y, g)
persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "lightblue")

在你的例子中，n！似乎只是一个标量，所以它不会影响情节的形状。