在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！似乎只是一个标量，所以它不会影响情节的形状。