用R估计两个样本均值之间的差异

时间:2015-11-15 19:06:19

标签: r static static-methods

我正处于语言R编程研究的开始阶段,

我尝试编写自己的函数来构造两个成对样本均值之间差异的置信区间估计。

这是我的代码

::msg.en.hello=Hello
::msg.fi.hello=Päivää
powershell -c "Get-Content '%~f0' -encoding UTF8 | Select-String -Pattern '^::msg\.' | ForEach-Object {$_ -replace '::msg.',''} | Set-Content '%temp%\messages.txt' -encoding Unicode"

我需要显示结果并用它来书写说明,但我不能。

感谢您的帮助

1 个答案:

答案 0 :(得分:1)

你离得太远了。

R按照您编写的顺序运行代码。因此,您必须先创建n,然后才能在计算中使用它。

My.Function=function(X1,X2,con.int){
  n = length(X1) # this comes first!!
  X1bar=sum(X1)/n
  X2bar=sum(X2)/n
  XD = X2-X1
  XDbar=sum(X2-X1)/n
  Sd =sqrt((sum((XD-XDbar)^2))/(n-1))
  Alpha = 1- con.int
  T = abs(qt(Alpha/2, n-1))
  Lower.B = (X2bar - X1bar)- T*Sd/sqrt(n)
  Uper.B = (X2bar - X1bar)+ T*Sd/sqrt(n)
  print(c(Lower.B,"X2bar - X1bar",Uper.B))
}

我还取得了T分数的绝对值,这样当qt()返回负数时,你的上限和下限就不会混淆。

您可能还希望使用cat功能打印结果。 https://stat.ethz.ch/R-manual/R-devel/library/base/html/cat.html