详细的说明请参考下面的文章:
http://www.sciencenet.cn/blog/user_content.aspx?id=305762
在知道了分布函数和相应的统计量后怎样得到p-value呢?用R可以很容易的实现,如下(来自<统计建模和R软件>)
P_value <- function(cdf,x,paramet,side=0){
n <-length(paramet)
p <-switch(n+1,
cdf(x),
cdf(x,paramet),
cdf(x,paramet[1],paramet[2]),
cdf(x,paramet[1],paramet[2],paramet[3])
)
if(side <0) p
else if(side >0) 1-p
else
if(p<1/2) 2*p
else 2*(1-p)
}
其中cdf是分布函数(下面会列出一些分布),
x是计算出来的统计量的值,如t分布中的t值,f分布中的f值,
paramet是分布的参数,如正态分布需要指定(mu,sigma),f分布需要指定(df1,df2),
side是指计算p值的类型,是单侧,还是双侧.