你这是一个两样本自由度为4的t统计量的检验，查一下他们计算公式，也就是把两个样本的均值相减以后除以他们联合样本标表准差，算出以后针对t分布的表就可以查出P值了。随便找本概率统计的书都有这个公式的。
- 好像是散的，不能整体下载，麻烦！
- 这是一个多重检验问题，比较每一种情况下两种系统的好坏，最后做一综合评价啊。
- 谢谢两位指点！
- a=1125;b=501;
> c=b/a
> for(i in 1:13){
+ c=c*(b+i)/(a+i)}
> d=1125*c
> d
[1] 0.01499478
>
这个应该是正确答案！ - 本人又模拟了一下，代码如下：
m=100000
c<-rep(0,m)
for(j in 1:m)
{
x<-1:1138
y<-sample(x,514)
z<-sort(y)
for(i in 1:501){ U<-z[(i+1):(i+13)] ;
V<-z[i:(i+12)]+1 ;
temp<-(U==V)
if(sum(temp)==13) {c[j]=1;break}
}
}
sum(c)/m
这个计算要花点时间，运行结果是 0.00793，好像模拟次数越大，结果越小！（m=10000时为0.0088），大家可以试一下。 - a=1125;b=501;
> c=b/a
> for(i in 1:13){
+ c=c*(b+i)/(a+i)}
> d=1125*c
> d
[1] 0.01499478
>
这个应该是正确答案！ - 补充说明一下，分子应该是从1125中任选一个，乘以从1124中任选500个，再乘以501的阶乘（即捆绑后的全排列）。
分母应该是从1138个中任选514进行全排列。 - 我用捆绑法给出一种思路，大家看行不行？
首先1138个数里14个连号的共有1125个，而这14个号在514个号里占了14各位，就有501个位置任其选，剩下的500个号来自剩下的1124个数字。按这种思路，即先从1125个数任选一个放在514个位置中的任意一个位置，然后再从1124个数中任选500个数全排剩下的500个位置，则运用排列组合计算式应该是：
这里没有编辑器，又复制不过来，不好意思！有了这个思想写起来很简单，各位自己写吧！ - 各位同仁：
本人最近需要学习空间统计，不知道有没有比较好的教材可以提供的？非常感谢！ - 最近需要应用R解决保序回归的问题，但是带有一些条件限制，比如限制估计范围等，不知道有没有现成的命令可以调用的？希望各位高手指点一二，非常感谢！
- 看一下
- 是不是只有校园网才能打开啊？我是校外联通的，进不去啊
- FDR是控制同时进行多个检验时犯错误的概率，至于生物上的解释就不清楚了！
- 你问的是多重检验的问题，FDR及local FDR是控制检验错误的度量，前者是由Benjiamini（1995）年提出，后者是指局部度量，由Efron提出，相关的细节要查它的文章，q-value是由Storey（2002）提出的与之有关的概念，记不清了，好像什么强度吧。
- 看看
- 想看一下！
- 想看一看！
- 可以安装包bootstrap，直接调用就行了
- 希望有一本针对R中运行出错问题的解答，最好是中文版的啊！