我想创建一个学生的分布密度图,平均值为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))
这种做法是否正确?
如果它是正确的,我怎样才能获得真正的密度,而不仅仅是近似值?
答案 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
参数,如果你想比较它们,可以用它们在同一个图上绘制(在这种情况下差别不大),改变其中一个他们会被建议。