使用预先指定的平均值生成t分布密度图

时间:2014-06-27 16:42:35

标签: r distribution

我想创建一个学生的分布密度图,平均值为0.02而不是0.这可能吗?

分配应该有2个自由度。

这可能吗?

我尝试了以下内容:

X<-rnorm(100000,mean=0.02, sd=(1/sqrt(878)))
pop.mean<-mean(X)
t<-sapply(1:10000, function(x) (mean(sample(X,100))-pop.mean)/(1/sqrt(878)))

plot(density(t))
  1. 这种做法是否正确?

  2. 如果它是正确的,我怎样才能获得真正的密度,而不仅仅是近似值?

1 个答案:

答案 0 :(得分:2)

你的陈述和例子在某种程度上相互矛盾。

你想要一个基于平均值为0.02的非中心t分布吗?这就是你的例子所暗示的,但请注意,非中心t不仅仅是一个偏移的t,它现在是偏斜的。

如果您想要非中心t,那么您可以使用如下命令绘制它:

 curve(dt(x,2,0.02), from=-5, to=6)

或者,您想要转移t分布吗?在0.02附近对称的分布,具有t分布的形状?

您可以使用如下命令绘制移动的曲线:

curve(dt(x-0.02,2), from=-5, to=6 )

curve函数有一个add参数,如果你想比较它们,可以用它们在同一个图上绘制(在这种情况下差别不大),改变其中一个他们会被建议。