从Unix中的估计和标准错误计算p值

-1

我正在寻找一种方法来根据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

感谢您的任何建议来完成。

来源

2017-05-31 ThePooh

这可能对你有用：https：//stackoverflow.com/questions/577982/how-to-compute-the-p-value-in-hypothesis-testing-linear-regression – jrcastle90

阅读[问]然后你的问题要包含一个简洁，可测试的样本输入和期望输出的[mcve]加上如何计算“p值”的描述。如果你假设我们都知道如何读出“R”来找出你的需求和/或我们知道什么是“p值”，那么YMMV--这组awk专家和R/p值集之间的交集专家可能不如你想象的那么好。 –

“unix”中的含义是什么？你为什么不使用R，因为它可以在UNIX上运行？ – Emer

awk来救援！

使用算法16年2月26日从 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

有与将需要更长的时间来实现更显著的数字更好的算法，但如果你需要一些十进制数字精度，这是不够好。

来源

2017-06-01 17:24:27 karakfa

谢谢你的帮助:) – ThePooh

从Unix中的估计和标准错误计算p值

回答

相关问题