建立研究队列
从ANZARD数据库筛选出符合目标试验条件的患者,即首次刺激ART周期且受精后仅有单个合子的未产患者。
回顾性分析2009-2022年澳大利亚和新西兰的ANZARD数据。筛选条件包括:使用自己的精子和卵母细胞、首次刺激ART周期、单合子、未产。排除标准包括:完全受精失败、使用PGT、使用捐赠配子或代孕、进行胚胎捐赠、冻胚后90天内未解冻、在无治疗变异的诊所就诊。最终纳入11,163名患者。
Cleavage stage versus blastocyst stage transfers in patients with a single zygote: an emulated target trial
基于大型回顾性注册数据库(ANZARD),采用目标试验模拟框架和g计算,包含统计建模和敏感性分析,但缺乏胚胎质量等关键混杂因素数据,且未提供直接实验操作详情。
基于ANZARD注册数据库的回顾性队列(2009-2022年,n=11,163)
患者人群限定为首次刺激ART周期且仅有单个合子的未产患者 治疗分配(卵裂期或囊胚期移植)通过多项式混合模型推断,利用治疗年份作为外部变异 主要结局为活产,按WHO定义(≥20周妊娠或≥400g出生体重) 使用g计算模拟理想化目标试验,估计意向性治疗效果 使用bootstrap(B=1000)计算平均治疗效果和条件平均治疗效果的置信区间
提炼研究问题、关键发现与证据,快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
从ANZARD数据库筛选出符合目标试验条件的患者,即首次刺激ART周期且受精后仅有单个合子的未产患者。
回顾性分析2009-2022年澳大利亚和新西兰的ANZARD数据。筛选条件包括:使用自己的精子和卵母细胞、首次刺激ART周期、单合子、未产。排除标准包括:完全受精失败、使用PGT、使用捐赠配子或代孕、进行胚胎捐赠、冻胚后90天内未解冻、在无治疗变异的诊所就诊。最终纳入11,163名患者。
估计基于患者特征的胚胎存活到卵裂期或囊胚期的概率,以及患者被分配接受卵裂期或囊胚期移植的概率,以解决选择性/生存偏倚。
使用多项式混合模型,以治疗年份、诊所随机效应和患者及治疗协变量预测合子存活到适合移植或冷冻的卵裂期或囊胚期胚胎的概率。模型以可用卵裂期胚胎、可用囊胚期胚胎和无可用胚胎的患者计数为结局,每个结局建模为治疗分配和每个阶段可用胚胎概率的子模型的混合。协变量包括女性年龄、男性年龄、不孕诊断、取卵数、受精方法、受精率和使用睾丸精子。模型形式通过AIC进行变量选择。
估计在给定的卵裂期或囊胚期移植下,实现活产或临床妊娠的概率。
使用广义加性模型(GAM)对二元结局(活产、临床妊娠)的概率进行建模,给定协变量。模型包括对连续变量(女性年龄和男性年龄)的平滑项和分类变量(不孕诊断、取卵数、受精方法、使用睾丸精子)的收缩项。ART诊所作为随机截距。通过广义交叉验证(GCV)确定适当的收缩程度,使用R包mgcv。
模拟理想化的随机对照试验,估计意向性治疗效果(卵裂期 vs 囊胚期移植)对活产和临床妊娠的影响,并评估效应的异质性。
使用g计算,基于开发的胚胎发育和移植结局模型,在满足目标试验标准的患者队列上模拟理想化RCT。对于每位患者,预测在两种意向分配(卵裂期或囊胚期移植)下的潜在结局,计算平均治疗效果(ATE),包括风险差(RD)和相对风险(RR),以囊胚期为参照组。通过拟合以女性年龄为平滑项的GAM,计算条件平均治疗效果(CATE)。使用bootstrap(B=1000)计算置信区间。
评估估计的治疗效果对不同胚胎发育率和移植结局参数的敏感性,以检验结果的稳健性。
绘制结果图,改变平均囊胚发育率(30%-70%)和平均囊胚移植活产率(5%-25%),同时将卵裂期发育和活产率固定为研究中的平均值,以说明治疗效果的敏感性。
按研究目的归类文中使用的方法,便于定位所需技术。
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 数据分析 | R包mgcv | -- | -- |
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 队列定义 | 患者必须是首次刺激ART周期、仅形成单个合子、未产;排除使用PGT、捐赠配子或代孕等情况 阅读提示:Materials and methods - Data and cohort |
| 因果模型 | 需要验证多项式混合模型和GAM模型的具体设定,包括协变量、随机效应、变量选择方法(AIC、GCV) 阅读提示:Materials and methods - Embryo development and treatment assignment models, Embryo transfer outcome models, and Supplementary Data File S1 |
| 治疗效果估计 | g计算的具体实现,包括模拟次数、bootstrap样本量(B=1000) 阅读提示:Materials and methods - Causal estimation |
| 敏感性分析 | 需要核对敏感性分析中设定的囊胚发育率和活产率的范围(30%-70%,5%-25%) 阅读提示:Materials and methods - Sensitivity analyses and Figure 3 legend |