[s:15]
问题是这样的:
R下可以直接cor.test()对相关系数r进行显著性检验,获得P值.
现在需要直接根据t和df,运用脚本来计算p value,费了大半天,终于找到了t分布的密度函数计算公式 f(t),并在EXCEL下直接计算,公式大致为:
EXP(GAMMALN((df+1)/2))*(1+t^2/df)^(-(df+1)/2)/SQRT(df*PI())/EXP(GAMMALN(df/2))
但是计算结果和EXCEL下的TDIST计算结果不太一致,而TDIST和R的cor.test()计算结果是一致的.
请问f(t)和TDIST的结果(很近似),有么区别,我该选择哪个?
我需要解决的问题,是根据t和df直接运用公式(计算机语言实现)计算p值.
特地注册了帐号,想请高手指点,不胜感谢.