荟萃分析中风险比率的概貌似然检验改进方法

0 下载量 73 浏览量 更新于2024-07-15 收藏 764KB PDF 举报
"这篇论文研究了在荟萃分析中如何更有效地检验风险比率的异质性,特别是针对Cochran Q检验在处理少量研究时的局限性。文章提出了两种改进的检验方法,即PL1和PL2,它们基于轮廓最大似然估计量,并在标准卡方检验统计量中应用风险比对数的方差公式。此外,还介绍了一个简单的天真测试(SIM)作为比较。通过模拟实验评估了这些测试在不同研究数量、样本量和风险比率情况下的性能。结果显示,对于中等到大规模的研究(k≥16)和中等到大样本量(≥50),PL1和PL2测试在控制I型错误率和提高检测能力方面表现最优。然而,当样本量较大(≥500)时,标准的Cochran Q测试表现最佳。对于小样本量(≤10)的情况,没有哪种测试表现出色。论文强调了这些发现对于实际应用中的选择具有指导意义。" 这篇论文探讨了在荟萃分析研究中,尤其是在涉及临床疗效评估时,如何改进风险比率的异质性检验方法。传统的Cochran Q检验由于合并研究数量有限,其检测异质性的能力相对较弱。为了解决这个问题,研究者提出PL1和PL2两种新的测试方法,这两种方法的核心在于将风险比的最大似然估计量引入到卡方检验的统计量中,以提高检验的精确度。同时,SIM测试被引入作为对比,以评估这三种方法的性能。 通过大量的模拟实验,研究者系统地考察了这三种方法在不同参数设置下的表现,包括显著性水平、研究数量、治疗组和对照组的样本量以及真实的风险比率。实验结果表明,PL1和PL2测试在研究数量较大(至少16个研究)且样本量适中(50或100)的情况下,具有最高的功效,能够更好地控制I型错误率。而当样本量进一步增加(如500以上)时,标准的Cochran Q检验则表现出更好的性能。然而,无论研究规模如何,对于小样本量研究(样本量≤10),所有测试方法的表现都不理想。 这些发现对于实际的荟萃分析工作具有重要的实践指导价值。研究人员可以根据研究的具体条件,如预期的样本规模和研究数量,来选择最合适的统计检验方法。对于大型研究项目,PL1和PL2可能是更优的选择,而Cochran Q测试则适用于样本量较大的情况。这篇论文为优化荟萃分析中的统计检验提供了理论依据和实用工具。

clear %清除内存 load('1797 b007_0.mat') %根据实际需要更改地址、路径 sig=X118_DE_time(1:12000); fs=12000; N=12000; Ts=1/fs; sig=sig(1:N);%设置取样频率fs,取样数N t=0:Ts:(N-1)*Ts;%时间轴 t sig=(sig-mean(sig))/std(sig,1);%对 sig 进行归一化 subplot(211);plot(t,sig);%绘制 sig 波形 xlabel('时间 t/s'); ylabel('振动加速度/V'); nfft=fs/2; % 16384 S=pspectrum(sig,nfft);%对 sig 做功率谱 subplot(212); plot((0:nfft/2 -1)/nfft*fs,S(1:nfft/2));% 绘制功率谱 xlabel('频率 f/Hz'); ylabel('功率谱 P/W') [c,l]=wavedec(sig,3,'db2');%利用 db2 对 sig 进行 3 级小波分解 c3=wrcoef ('a',c ,l,'db2',3); d3=wrcoef('d',c,l,'db2',3); d2 =wrcoef ('d',c,l,'db2',2); d1 =wrcoef('d',c,l,'db2',1);%重构第 1-3 层细节 d1~d3 和第 3 层概貌 c3 figure; subplot(414); plot(t,c3); ylabel('c3');%绘制 c3 subplot(413); plot(t,d3); ylabel('d3');%绘制 d3 subplot(412); plot(t,d2); ylabel('d2');%绘制 d2 subplot(411); plot(t,d1); ylabel('d1');%绘制 d1 y=hilbert(d1); %对 d1 进行 Hilbert 变换,得y ydata=abs(y); %ydata=|y| ydata=ydata-mean(ydata);%对 ydata 去均值(目的是去除幅度较大的直流分量) P=pspectrum(ydata,nfft);%ydata 的功率谱为 P figure; plot((0:nfft/2-1)/nfft*fs,P(1:nfft/2)); xlabel('频率 f/Hz');%绘出 d1 的 Hilbert 包络谱 P=P(1:nfft/2); [M,f1]=max(P); f1=f1*fs/nfft-1 %故障频率 f1为包络谱中幅度最大处的频率 将代码由利用db2进行3级小波分解改为利用db10进行5级小波分解

2023-05-09 上传