我需要绘制由以下各项定义的阶数统计量的二元随机向量的联合密度函数:
f(x,y) = (n!)*f(x)*f(y) for x < y and 0 otherwise
其中f是单变量随机变量的密度函数。例如,f(x)=e^(-x)
。
我不知道我是否必须使用persp
该条件以及如何修复该支撑x<y
。
非常感谢你。朱利奥
此类绘图的典型工作流程是定义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!似乎只是一个标量,因此不会影响图的形状。
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句