Paper record
A practical workflow combining Kaplan-Meier and Bayesian accelerated failure time analyses for censored reproductive phenology data in soybean
BMC Plant Biology · 12 Sept 2026 · 10.1186/s12870-026-09851-6
Abstract
Abstract Developmental progression in plants is inherently a time-to-event process, yet plant phenology data are often analyzed as ordinary endpoint traits even when some individuals fail to reach the target stage within the observation window. Such observations are right-censored rather than missing and should be retained in inference. Here, we evaluated 21 soybean cultivars across 10 controlled-environment treatment combinations varying in CO 2 concentration, photoperiod, and LED light quality, and analyzed days from sowing to the R6 stage using a practical workflow that combined Kaplan-Meier survival analysis with a Bayesian Weibull accelerated failure time (AFT) model. Based on 210 cultivar × treatment observations derived from 1,050 plants, treatments resolved into favorable, intermediate, and strongly inhibitory classes. Elevated CO 2 showed the strongest association with reproductive timing: increasing CO 2 from ambient to 1000–1400 ppm reduced median survival time, while the additional gain from 1000 ppm to 1400 ppm was minimal, indicating a saturating response. Blue - Red spectral treatments and 6–8 h photoperiods were associated with rapid and synchronized development within the tested treatment combinations. The Bayesian AFT model included treatment as a random effect to account for the chamber-level experimental structure; the between-treatment standard deviation was estimated at 0.02 (95% credible interval: 0.00–0.06), indicating negligible chamber-to-chamber variation. The highest-ranked predicted combination (1400 ppm, blue: red = 2:1, 6 h) gave a median time to R6 of 66.90 d (mean 72.18 d). This study provides a practical and transferable workflow for censored plant phenology datasets in controlled-environment research, phenotyping, and breeding. The workflow integrates established survival analysis methods and explicitly accounts for chamber-level design, offering a framework for structured plant phenotyping experiments.
Code and data availability
The supplied article blocks contain no data availability, code availability, or supplementary material statements. The paper describes soybean phenotyping (21 cultivars, 10 treatments, Kaplan-Meier and Bayesian Weibull AFT analyses), but no public phenotype dataset, images, analysis code, or trained model is referenced
No evidence-backed public reproduction asset is currently recorded.