核fisher matlab,使用Fisher法matlab结合P值?

我不认为这是建立在MATLAB费舍尔的组合概率测试,但它不应该是很难实现:

P = [0.43 0.12 0.0021 0.05 0.017 0.001 0.025 0.038 0.35 0.29];

k = length(P);

首先我们要造一个配偶帮助矩阵总结中的元素数p作为我们想:

% the following matrix is used to sun each n elements in a row:

n = 3;

summer = diag(ones(k,1));

for d = 1:n-1

summer = summer + diag(ones(k-d,1),-d);

end

如果我们运行P*summer,我们得到:

ans =

Columns 1 through 6

0.5521 0.1721 0.0691 0.068 0.043 0.064

Columns 7 through 10

0.413 0.678 0.64 0.29

接下来,我们首先采取一切P的LN计算统计,比它们求和3的(和-2乘):

% compute the combine fisher statistic:

X = -2.*log(P(:).')*summer;

结果:

X =

Columns 1 through 6

18.26 22.564 26.472 27.956 29.342 27.734

Columns 7 through 10

16.018 11.116 4.5754 2.4757

最后我们计算P-

% get the p-values for all combinations:

p_vals = chi2cdf(X(1:end-n+1),6,'upper');

而我们得到:与

2*3 = 6 DF从卡方分布值

p_vals =

Columns 1 through 6

0.005614 0.00095661 0.00018177 9.577e-05 5.2399e-05 0.00010546

Columns 7 through 8

0.013659 0.084865

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值