Matlab FDR校正的使用

Posted 2023-04-18

参考技术A 分析脑电在涉及到多次比较的时候，往往要进行多重比较校正，常用的是有FWE校正和FDR，还有NBS校正，以下的是关于FDR校正的实现，资料来自网上

%最简单的实现方式，基于Storey procedure ( introduced by Storey, 2002)，适用于P值数量>1000的情况，否则原则上会崩溃。我用MATLAB测试过，会出现warning，但不会报错。严格程度较低，如果你的 ttest P值不是特别显著（0.01-0.05），可以用这个试试，或许可以过FDR校正.

%基于linear step-up (LSU) procedure (introduced by Benjamini and Hochberg, 1995)。最常见的FDR校正方式，严格程度较高，但比
Bonferroni校正低，适用于 ttest P值显著（<0.01）。

注：

BHFDR的计算过程

公式：FDR = p * (n/i)， p是pvalue，n是p值个数，最大的P值的i值为n，第二大则是n-1，依次至最小为1。

%指定调整参数λ，用于估计零假设为真的先验概率，Lambda Value可以是：（1）0-1内任意一值 （2）4个以上序列，或以示例的矩阵形式：[first:incr:last]，mafdr函数自动选择最优参数。严格程度比1略低，我用第一个不能过校正但是用这个居然可以过，为FDR苦恼得同志们可以被解救了。

%对选择最优参数lambda的方法进行选择，'bootstrap' (default)，'polynomial'。第3种方法的延伸。

MATLAB 图像校正，希望路过的大神抽点宝贵的时间出来为我解决一下

把图像里的A4纸的梯形模型校正成长方形的，多点注释，谢谢（没有财富，sorry）
就是本来是长方形的，由于拍的角度造成看起来像梯形的了，现在怎么把它校正过来

参考技术A

仿射变换

tform=maketform('projective',[0 0;n 0;n m;0 m],[a b;c d;e f;g h]);
img2=imtransform(img,tform);

追问

看不懂坐标，A B C D E F F G H ，该怎么给值，还有n m 怎么知道它的坐标

追答

[m,n,d]=size(img); %图像是彩色的，所以是三维

abcdefg，可以通过imshow(img)窗口中 data cursor按钮得到坐标。

参考技术B 1，二值化，抠出纸张区域.
2，imresize, 参数采用双线性内插。 help imresize追问

额，不懂，怎么抠区域，