从诊断测试数据计算95％置信区间

Question

我正在尝试根据Matlab中的一系列矩阵计算95％置信区间：

我知道如何计算所需的敏感度，特异性，阴性预测值和阳性预测值，但是我不确定，如果这些数据，如何计算95％置信区间。

非常感谢任何帮助。

Answer 1

您需要整个分布来计算置信区间。 获得后，将其命名为“data”并运行以下代码。 该代码使用数据中的nans，使用向量（Nx1）或矩阵（NxM）。只要每个分布由列

表示

  %%% calculate CI
    alpha_lvl = [2.5 97.5]; % lower and upper boundaries for the CI 
    SEM = nanstd(data)./sqrt(sum(~isnan(data )));   % Standard Error
    ts  = tinv(alpha_lvl/100,length(data)-1);     % T-Score
    CI  = bsxfun(@plus,nanmean(data ),bsxfun(@times,SEM,ts')); % confidence intervals