私は、推定値とUNIXの標準誤差からp値を計算する方法を探しています。 はRでは、たとえばそれは、次のUnixの推定値と標準誤差からp値を計算する
test
BETA STD
-0.0103312 0.0812054
T=(-0.0103312/0.0812054)^2
p <- pchisq(T,df=1,lower=F)
p
0.8987638
を使用して行うことができます任意の提案をいただき、ありがとうございます。
awkレスキュー！アルゴリズム26.2.16を使用して
Abramowitz and Stegun: Handbook of Mathematical Functions
awk 'function z(x) {return exp(-x*x/2)/sqrt(2*atan2(0,-1))}
function p(x) {t=1/(1+0.33267*(x>0?x:-x))
return 1-z(x)*t*(0.4361836+t*(-0.1201676+t*0.9372980))}
{print 2*(1-p($1/$2))}'
からあなたに両側検定のp値が得られます。あなたの入力-0.0103312 0.0812054
については
は0.898743
を実装するために時間がかかりますより多くの有効桁数とのより良いアルゴリズムがあります提供しますが、あなたは、いくつかの桁の精度が必要な場合、これは十分です。
ありがとうございました:) – ThePooh
これはあなたの役に立つかもしれません：https://stackoverflow.com/questions/577982/how-to-compute-the-p-value-in-hypothesis-testing-linear-regression – jrcastle90
[ask] then簡潔でテスト可能なサンプル入力と期待される出力と "p値"の計算方法の説明を加えた[mcve]を含む質問を[編集]してください。 「R」を読んで要件を把握する方法を知っていると仮定している場合、あるいは「p値」が何であるかを知っているとすれば、awk専門家とR/p値セット専門家はあなたが思うほど大きくないかもしれません。 –
"in unix"とはどういう意味ですか？なぜあなたはRをUNIX上で動かすことができるのではないのですか？ – Emer