Stata绘制带置信区间条形图报错求助及显著性值添加咨询

一、解决"too few variables specified"报错 #

报错核心原因是graph twoway的bar和rcap命令缺少x轴分组变量，且by(T_Meat)会生成子图，不符合你要的并列条形图需求。修改后的完整代码如下：

collapse (mean) meaninfo_avoiding_ba= info_avoiding_ba (sd) sdinfo_avoiding_ba=info_avoiding_ba (count) n=info_avoiding_ba, by(T_Meat)
generate hiwrite = meaninfo_avoiding_ba + invttail(n-1,0.025)*(sdinfo_avoiding_ba / sqrt(n))
generate lowrite = meaninfo_avoiding_ba - invttail(n-1,0.025)*(sdinfo_avoiding_ba / sqrt(n))
// 修正绘图命令：指定x轴为T_Meat，取消by()改为同一图展示
graph twoway (bar meaninfo_avoiding_ba T_Meat, barwidth(0.6)) ///
            (rcap hiwrite lowrite T_Meat), ///
            ytitle("info_avoiding_ba均值") ///
            xtitle("T_Meat分组") ///
            legend(off) ///
            title("分组均值与95%置信区间")

关键调整说明：

• bar meaninfo_avoiding_ba T_Meat明确指定y轴为均值、x轴为分组变量T_Meat
• 移除by(T_Meat)，避免生成两个独立子图，改为同一图内并列展示两组条形
• 添加barwidth、ytitle等参数优化图表美观度

二、在条形上方添加显著性值 #

需先通过t检验获取两组差异的显著性标记，再将标记添加到图中，完整流程如下：

步骤1：对原始数据做t检验并生成显著性标记 #

不要直接覆盖原始数据，先执行检验并生成标记：

// 独立样本t检验
ttest info_avoiding_ba, by(T_Meat)
// 提取p值并生成显著性标记
local p_val = r(p)
gen sig_mark = ""
if `p_val' < 0.001 replace sig_mark = "***"
else if `p_val' < 0.01 replace sig_mark = "**"
else if `p_val' < 0.05 replace sig_mark = "*"
else replace sig_mark = "ns"

步骤2：将标记合并到分组统计数据中 #

// 给每组匹配对应的显著性标记
bysort T_Meat: gen sig = sig_mark[1]
// 重新执行collapse，保留显著性标记
collapse (mean) meaninfo_avoiding_ba= info_avoiding_ba ///
         (sd) sdinfo_avoiding_ba=info_avoiding_ba ///
         (count) n=info_avoiding_ba ///
         (first) sig=sig, by(T_Meat)
// 重新生成置信区间（若之前已生成可跳过）
generate hiwrite = meaninfo_avoiding_ba + invttail(n-1,0.025)*(sdinfo_avoiding_ba / sqrt(n))
generate lowrite = meaninfo_avoiding_ba - invttail(n-1,0.025)*(sdinfo_avoiding_ba / sqrt(n))

步骤3：绘图时添加显著性标记 #

修改绘图命令，通过text将标记放在条形顶端：

graph twoway (bar meaninfo_avoiding_ba T_Meat, barwidth(0.6)) ///
            (rcap hiwrite lowrite T_Meat) ///
            (text meaninfo_avoiding_ba T_Meat, text(sig) offset(0.2) size(medium)), ///
            ytitle("info_avoiding_ba均值") ///
            xtitle("T_Meat分组") ///
            legend(off) ///
            title("分组均值与95%置信区间") ///
            note("*p<0.05, **p<0.01, ***p<0.001, ns=不显著")

参数说明：

• offset(0.2)调整文本与条形顶端的距离，避免重叠
• note()添加显著性规则注释，提升图表可读性

三、额外美观优化建议 #

• 给条形设置不同颜色：bar(..., color(navy maroon))区分两组
• 自定义y轴刻度：ylabel(0(0.5)3)（根据你的数据范围调整）
• 调整置信区间线条样式：rcap(..., lwidth(medium) lcolor(black))

内容的提问来源于stack exchange，提问作者Bénédicte