来自fMultivar的随机变量不'；t匹配密度函数

我试图使用fMultivar包绘制t分布的直方图，并显示相应的密度函数与直方图匹配

问题是直方图似乎与密度函数不匹配

下面是一个简单的例子

library(fMultivar)

n = 50000
qtile = seq(-5,5,by=0.05)

xtt = rmvst(n,dim=1,df=3)
hist(xtt,300,xlim=c(-5,5),probability=TRUE)
lines(qtile,dmvst(qtile,dim=1,df=3),col='red')

xn = rmvsnorm(n,dim=1)
hist(xn,300,xlim=c(-5,5),probability=TRUE)
lines(qtile,dmvsnorm(qtile,dim=1),col='red')

xt = rt(n,df=3)
hist(xt,300,xlim=c(-5,5),probability=TRUE)
lines(qtile,dt(qtile,df=3),col='red')

很明显，xtt与密度函数不匹配。 xn和xt都有

---

我在这里做错了什么吗？

我不相信，至少在随机抽签是否给出了不同的分布方面。柱状图看起来是一样的，我在几次运行中查看了'xt'和'xtt'的

qqplot

，它们似乎排列得很好。将密度线与直方图进行比较确实表明密度估计值可能不正确，当我进行本次检查时，我同意密度估计值之间存在系统性差异：

plot( density(dt(qtile,df=3)- dmvst(qtile,dim=1,df=3) )  )

---

分布是双峰的，dmvst密度系统地大于尾部的dt密度。然而，我不明白为什么t分布必须与斜t分布相匹配。您不应该将其与包中的dst{in-package:sn}进行比较吗？（这是最好向包作者提出的问题。）

它看起来像是

dmvst

中的一个bug：在维度1中，忽略了

df

参数，而是使用高斯分布……Vincent，我没有想到看代码。谢谢你的建议。我想我应该。Dwin，我确实联系了其中一位作者，因为没有列出维护人员。希望我能收到他的回信。