镇静药种类与呛咳反应的卡方检验 - Matlab代码实现

假设药物B组的观察值为O1和O2，药物R组的观察值为O3和O4，无不良反应的期望频数为E1和E2，呛咳反应的期望频数为E3和E4。\r\n根据卡方检验的原理，计算出观察值和期望值的差异：\r\nO1-E1, O2-E2, O3-E3, O4-E4\r\n然后计算卡方统计量：\r\nX^2 = (O1-E1)^2/E1 + (O2-E2)^2/E2 + (O3-E3)^2/E3 + (O4-E4)^2/E4\r\n将上述数据代入matlab代码进行计算：\r\nobserved = [7, 468; 46, 724];\r\nexpected = (sum(observed,2)sum(observed))/sum(sum(observed));\r\nchi2 = sum(sum((observed-expected).^2./expected));\r\n根据卡方分布的自由度计算p值：\r\ndf = (size(observed,1)-1)(size(observed,2)-1);\r\np = 1 - chi2cdf(chi2, df);\r\n最后，根据p值判断镇静药种类对是否发生呛咳反应是否存在显著性差异。如果p值小于设定的显著性水平（通常为0.05），则可以拒绝原假设，即存在显著性差异。