如何在STATA中计算两组间的几何最小二乘均值比（含95%CI）

你已经通过对数转换+混合效应模型完成了蛋白尿相对基线的几何最小二乘均值百分比变化计算，接下来计算两组间的几何均值比值及95%CI可以通过以下步骤实现：

原理说明 #

由于你对结局变量percentupcr做了对数转换（生成lnpercent），两组的几何均值比值等于exp(HTN组对数尺度预测均值 - 非HTN组对数尺度预测均值)——这是对数转换的核心性质：对数的差值对应原尺度的比值。

Stata代码实现 #

1. 先确保模型已加载（如果之前的模型不在内存中，重新运行）：

gen lnpercent = ln(percentupcr)
xtmixed lnpercent htn##year || no:, var

2. 按年份计算HTN组 vs 非HTN组的几何均值比值及95%CI：

margins, over(year) contrast(htn) expression(exp(xb))

代码解释 #

• over(year)：指定按随访年份分层计算，因为你的模型包含htn与year的交互项，不同年份的组间比值可能存在差异。
• contrast(htn)：计算HTN组（假设htn=1）与非HTN组（htn=0）在对数尺度上的均值差。
• expression(exp(xb))：将对数尺度的均值差转换为原尺度的几何均值比值，同时Stata会自动对置信区间做指数转换，得到对应原尺度的95%CI。

如果需要不按年份的整体比值（仅当交互项无统计学意义时推荐），可以去掉over(year)参数：

margins, contrast(htn) expression(exp(xb))

内容的提问来源于stack exchange，提问作者donghoon kang