{"id":"61e25c9f-3e9e-483d-b86e-04c81035faf9","arxiv_id":"2606.23146","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"Introduces a bias-corrected two-stage approach for joint modeling of multivariate ordinal longitudinal HRQoL and survival data that approximates full joint Bayesian estimation with lower computation time.","lead":"The paper proposes a slope-corrected two-stage (SC2S) estimator for jointly modeling multidimensional ordinal longitudinal HRQoL data and survival outcomes. It aims to reduce bias relative to standard two-stage methods while remaining computationally lighter than full joint Bayesian estimation.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"No significant objection identified","rationale":"The reader's weakest assumption is the very point the simulations are designed to test. Because the full manuscript supplies those simulations and they support the claim under the conditions examined, the assumption does not appear to be load-bearing.","tokens_in":1704,"tokens_out":229,"duration_ms":13844,"concrete_test":"Re-run the main simulation design (Section 4) with the association parameter set to zero; confirm that SC2S and the naive two-stage estimator coincide within Monte Carlo error while the joint estimator remains unbiased.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on the SC2S procedure (informative priors on random effects from stage 1 plus re-estimation of slope parameters) being sufficient to recover most of the bias reduction of full joint estimation. The manuscript supplies simulation evidence under the data-generating model and a real-data application; both show the method tracks the joint Bayesian benchmark while remaining faster. No internal inconsistency, untested regime, or hidden assumption that would invalidate the headline result is apparent.","agreement_with_reader":"disagree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript proposes a slope-corrected two-stage (SC2S) approach for jointly modeling multidimensional longitudinal ordinal HRQoL data and survival outcomes in a latent trait framework. The SC2S method obtains random-effect estimates in stage 1, feeds them as informative priors into the survival submodel, and additionally re-estimates longitudinal slope parameters to mitigate bias induced by ignoring the dependence between the processes. Simulation studies under the data-generating model and a real-data application to progressive glioblastoma patients are presented to show that SC2S substantially reduces bias relative to standard two-stage procedures while closely approximating fully joint Bayesian estimation at lower computational cost.","tokens_in":1787,"tokens_out":414,"duration_ms":18378,"significance":"If the performance claims hold, the SC2S procedure supplies a practical, scalable alternative for joint modeling when the number of latent dimensions and random effects renders full joint likelihood maximization or MCMC prohibitive. By preserving most of the efficiency of two-stage methods while recovering much of the bias reduction of joint estimation, the approach could enable routine incorporation of multidimensional HRQoL trajectories into survival analyses in oncology settings where computational resources are limited.","major_comments":[],"minor_comments":[{"comment":"Abstract: the statement that SC2S 'substantially reduces bias' would be strengthened by a parenthetical reference to the specific bias or MSE reductions reported in the simulation tables.","section":"Abstract"},{"comment":"Section 4 (Simulation studies): clarify whether the re-estimation of slope parameters is performed inside the second-stage survival model or as a separate post-processing step, and state the convergence criteria used for the MCMC chains.","section":"Section 4"},{"comment":"Table 2 or equivalent results table: include the effective sample size or computation time per replication so readers can directly compare the efficiency gain against the joint Bayesian benchmark.","section":"Simulation results"}],"recommendation":"minor_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the positive assessment of our manuscript, the accurate summary of the SC2S approach, and the recommendation for minor revision. No major comments were raised in the report.","responses":[],"tokens_in":1255,"tokens_out":57,"duration_ms":15746,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main thing here is a two-stage method that cuts bias in joint models of multi-dimensional longitudinal HRQoL scores and survival by feeding stage-one results through informative priors on the random effects and then re-estimating the slopes.\n\nWhat is new is the slope-correction step combined with those priors inside a multidimensional latent-trait model. Earlier two-stage work is cited, but this specific pairing for ordinal multi-dimensional data does not appear in the references they give.\n\nThe paper shows through simulations that the approach reduces bias in both the longitudinal and survival parts relative to plain two-stage fitting, and the glioblastoma application demonstrates it runs faster than full joint Bayesian estimation while staying close to that benchmark.\n\nThe evidence is simulation-based under the data-generating model plus one real-data case. That is standard for this kind of work, though it leaves open how the method behaves under misspecification. The central assumption—that the priors plus slope re-estimation recover most of the dependence without the full joint likelihood—holds up in the reported results, so it is not a load-bearing flaw.\n\nThis is for biostatisticians who routinely fit joint models to oncology data with multiple HRQoL dimensions and need something faster than full joint estimation. Readers working on similar computational shortcuts would get direct value from the method and the timing comparisons.\n\nI would send it for peer review. The computational claim is concrete and the evidence is proportionate to what is claimed.","headline":"The SC2S procedure gives a workable bias fix for multidimensional HRQoL-survival joint models that tracks full Bayesian estimation at lower cost.","tokens_in":2242,"tokens_out":367,"would_cite":false,"duration_ms":17359,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"A slope-corrected two-stage method reduces bias in joint models of multidimensional quality-of-life trajectories and survival times.","keywords":["joint models","longitudinal data","survival analysis","HRQoL","bias correction","two-stage estimation","ordinal data","latent trait model"],"falsifier":"A simulation study in which the SC2S estimates show no reduction in bias relative to a standard two-stage fit when the true dependence between longitudinal trajectories and survival is known and nonzero.","tokens_in":2611,"feed_emoji":"📊","tokens_out":601,"duration_ms":17008,"temperature":0.7,"pith_summary":"Joint models link repeated health-related quality-of-life scores to survival risk so that dropout from death is properly accounted for, yet full joint fitting grows too slow once multiple questionnaire dimensions and random effects enter the picture. The paper introduces a slope-corrected two-stage procedure that first fits the longitudinal model, then feeds its information into the survival model through informative priors on the random effects and re-estimates the slope parameters. This keeps most of the speed of separate two-stage fitting while cutting the bias that arises when the dependence between the processes is ignored. Simulation studies show the resulting estimates stay close to those obtained from fully joint Bayesian estimation. An application to glioblastoma patients confirms that the method delivers usable accuracy gains without the full computational load.","feed_headline":"Slope correction cuts bias in joint HRQoL-survival models","feed_subtitle":"The two-stage method matches full joint estimation accuracy at far lower computational cost for multidimensional quality-of-life data.","key_machinery":"The slope-corrected two-stage (SC2S) procedure, which transfers longitudinal information via informative priors on random effects and re-estimates slope parameters to correct bias from ignored dependence.","core_discovery":"The slope-corrected two-stage (SC2S) approach propagates longitudinal information to the survival model through informative priors on the random effects while additionally re-estimating longitudinal slope parameters; this substantially reduces bias in both the longitudinal and survival submodels and closely approximates fully joint Bayesian estimation while retaining much of the computational efficiency of two-stage procedures.","pith_inferences":["The same information-passing step could be tested in other joint modeling settings where full likelihoods become intractable.","Researchers facing large questionnaires might adopt the procedure instead of ignoring dependence or accepting long run times.","The approach suggests that targeted re-estimation of a few key parameters can substitute for complete joint likelihood evaluation in bias correction.","If the method scales to higher-dimensional data, it could open routine joint analysis for trials that currently use only separate models."],"forward_implications":["Bias drops in both the longitudinal HRQoL submodel and the survival submodel.","Computational cost stays far below that of full joint estimation even with multiple latent dimensions.","The method handles multivariate ordinal longitudinal data within a latent trait framework.","Results in the glioblastoma application match full joint estimates closely enough for practical use."],"fun_headline_variants":["SC2S reduces bias in multidimensional HRQoL-survival models","Slope-corrected two-stage approximates full joint estimation","Two-stage approach corrects bias for HRQoL and survival data","Bias-corrected two-stage for longitudinal HRQoL-survival analysis"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"That passing longitudinal information through informative priors on random effects and then re-estimating slopes is sufficient to remove bias caused by ignoring dependence between the processes.","fun_headline_variants_meta":{"raw":{"variants":["SC2S reduces bias in multidimensional HRQoL-survival models","Slope-corrected two-stage approximates full joint estimation","Two-stage approach corrects bias for HRQoL and survival data","Bias-corrected two-stage for longitudinal HRQoL-survival analysis"]},"model":"grok-4.3","cost_usd":0.004339,"raw_usage":{"total_tokens":2167,"prompt_tokens":648,"num_sources_used":0,"completion_tokens":66,"cost_in_usd_ticks":43387000,"prompt_tokens_details":{"text_tokens":648,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":1453,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":648,"tokens_out":66,"duration_ms":11504,"temperature":1.0,"reasoning_tokens":1453,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-26T07:35:49.358846+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A simulation study in which the SC2S estimates show no reduction in bias relative to a standard two-stage fit when the true dependence between longitudinal trajectories and survival is known and nonzero.","supporting_citations":[],"review_version":1}