{"id":"d77d619f-c704-4b8e-895b-12d481f589d9","arxiv_id":"2606.00795","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"Extends MDBNs to non-Markovian queues via phase-type approximation, enabling accurate PCQs with orders-of-magnitude inference speedup on G/M/1 queue.","lead":"This paper extends modular dynamic Bayesian networks to non-Markovian queues by approximating non-exponential distributions with phase-type distributions. A smart generalist might read it to see how metamodels can enable faster causal analysis of realistic queueing systems.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.3","headline":"Phase-type approximation's accuracy-tractability tradeoff lacks quantitative validation for PCQ accuracy","rationale":"The reader's weakest_assumption directly identifies the same load-bearing point extracted from the abstract. Because the full text is referenced but the assessment remains abstract-limited, no stronger internal inconsistency or alternative concern can be isolated; the phase-count tradeoff remains the least-secured step in the argument.","tokens_in":1679,"tokens_out":271,"duration_ms":19150,"concrete_test":"For the reported G/M/1 experiments, recompute all PCQ results while varying the phase-type order from 2 to 8 phases and plot absolute error versus direct simulation; if error fails to drop below 5% or plateaus before tractability is lost, the balancing assumption does not support the accuracy claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on approximating non-exponential distributions via phase-type distributions to extend MDBNs to non-Markovian queues. The abstract explicitly flags the challenge of balancing metamodeling accuracy and tractability when choosing the number of phases, yet reports no analysis or sensitivity results on how phase count affects PCQ error in the G/M/1 experiments. Without this, the assertion of 'accurate answers' cannot be tied to the approximation method.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper extends modular dynamic Bayesian networks (MDBNs) for causal metamodeling from Markovian to non-Markovian queueing systems. It approximates general (non-exponential) distributions via phase-type distributions to enable estimation of probabilistic and causal queries (PCQs) from a single trained model. The work identifies three challenges—selecting the number of phases to balance accuracy and tractability, efficient parameter learning, and choosing a sampling interval for discrete-time approximation of continuous-time dynamics—and offers preliminary solutions. Experiments on a G/M/1 queue are reported to show that the resulting MDBN yields accurate PCQ answers while providing orders-of-magnitude speedup over direct simulation.","tokens_in":1750,"tokens_out":363,"duration_ms":19150,"significance":"If the central claims hold, this would constitute the first causal metamodeling approach applicable to non-Markovian systems, addressing a clear limitation of prior MDBN work. The ability to answer a range of PCQs from one model without repeated expensive simulations could be valuable for queueing analysis in operations research and performance modeling, where general distributions are the norm.","major_comments":[{"comment":"Experiments section (G/M/1 queue results): the manuscript reports that the MDBN produces 'accurate answers' to PCQs but contains no sensitivity analysis, error tables, or quantitative results showing how PCQ error varies with the number of phases used in the phase-type approximation. Because the abstract itself identifies the accuracy-tractability tradeoff as a novel challenge and the central claim rests on this approximation, the absence of such validation leaves the accuracy assertion unsupported.","section":"Experiments section (G/M/1 queue results)"}],"minor_comments":[],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their thoughtful review and for highlighting the need for stronger validation of the phase-type approximation. We address the major comment below and will revise the manuscript accordingly.","responses":[{"response":"We agree that the current experiments section provides only a qualitative demonstration of accuracy on the G/M/1 queue without quantitative sensitivity analysis or error tables versus the number of phases. The manuscript is framed as initiating the extension with preliminary solutions to the identified challenges, which is why the validation was kept at a high level. To strengthen the central claim, we will add in the revised version a new subsection (or expanded experiments) that includes (i) error tables for representative PCQs as a function of phase count (e.g., 2, 5, 10 phases) and (ii) a plot or table showing the accuracy-tractability tradeoff. This will directly address the referee's concern while preserving the preliminary nature of the work.","revision_made":"yes","referee_comment":"[Experiments section (G/M/1 queue results)] Experiments section (G/M/1 queue results): the manuscript reports that the MDBN produces 'accurate answers' to PCQs but contains no sensitivity analysis, error tables, or quantitative results showing how PCQ error varies with the number of phases used in the phase-type approximation. Because the abstract itself identifies the accuracy-tractability tradeoff as a novel challenge and the central claim rests on this approximation, the absence of such validation leaves the accuracy assertion unsupported."}],"tokens_in":1322,"tokens_out":328,"duration_ms":9000,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main advance is adapting modular dynamic Bayesian networks to non-Markovian systems. Prior MDBN work stayed inside Markovian assumptions, so using phase-type distributions to handle general interarrival or service times is the concrete step forward. The authors flag the practical issues this creates—selecting the number of phases, learning the parameters, and picking a sampling interval—and sketch preliminary fixes for each. The G/M/1 experiments are presented as evidence that the resulting metamodel can answer probabilistic causal queries with large speedups over direct simulation.\n\nThe soft spot is exactly the one the stress-test note flags. The abstract itself notes the accuracy-tractability tradeoff when choosing phases, yet the reported results give no sensitivity numbers, error curves, or ablation on phase count versus PCQ error. Without that, the claim of “accurate answers” floats on an untested assumption about how good the approximation needs to be. If the full paper contains those checks, they need to be front and center; if not, the empirical support stays thin.\n\nThis is aimed at people already working on simulation metamodels or queueing applications who need causal queries rather than just means and variances. A reader in that niche can extract the problem formulation and the list of open challenges even if the current experiments are preliminary.\n\nI would send it to peer review. The direction is sensible and the gap it targets is real; referees can usefully press for the missing validation on the phase-type step.","headline":"The paper takes a first cut at causal metamodeling for non-Markovian queues by swapping in phase-type approximations, but the validation on how well that approximation holds up is missing.","tokens_in":2226,"tokens_out":379,"would_cite":false,"duration_ms":16309,"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":"Approximating general distributions with phase-type distributions extends modular dynamic Bayesian networks to non-Markovian queues.","keywords":["causal metamodeling","modular dynamic Bayesian networks","non-Markovian queues","phase-type distributions","G/M/1 queue","probabilistic causal queries","discrete-event simulation"],"falsifier":"Running the MDBN and direct simulation in parallel on the same G/M/1 queue and finding that the MDBN answers to probabilistic causal queries deviate substantially from the simulation results.","tokens_in":2577,"feed_emoji":"⏱️","tokens_out":607,"duration_ms":20067,"temperature":0.7,"pith_summary":"This paper extends modular dynamic Bayesian networks, which previously worked only for Markovian systems, to non-Markovian queues. It achieves the extension by replacing non-exponential distributions with phase-type approximations. The work identifies and gives initial solutions for the resulting issues of phase count selection, parameter learning, and sampling interval choice. Experiments on a G/M/1 queue confirm that the extended model answers probabilistic and causal queries accurately while delivering large speedups over direct simulation.","feed_headline":"Phase-type approximations extend causal metamodels to non-Markovian queues","feed_subtitle":"MDBNs on G/M/1 queues answer PCQs accurately with orders-of-magnitude faster inference than direct simulation","key_machinery":"Modular dynamic Bayesian networks (MDBNs) extended by phase-type distribution approximations to represent non-Markovian continuous-time dynamics in a discrete-time model.","core_discovery":"By approximating non-exponential distributions using phase-type distributions, modular dynamic Bayesian networks can be extended to non-Markovian queues, yielding the first causal metamodeling technique for such systems. Experiments on a G/M/1 queue demonstrate that the MDBN can produce accurate answers to PCQs with orders-of-magnitude speedup of inference times relative to direct simulation.","pith_inferences":["The speedup may allow real-time causal analysis in applied queueing settings such as service systems or networks.","Similar phase-type approximations could be tested on multi-server or network-of-queues models.","Automatic methods for selecting the number of phases might be developed to reduce manual tuning."],"forward_implications":["A single trained model estimates a range of probabilistic and causal queries on the non-Markovian system.","Inference times are orders of magnitude shorter than those of direct simulation.","The number of phases can be chosen to balance metamodel accuracy against computational cost.","The same framework applies to other non-Markovian queueing models once the phase-type approximation is set."],"fun_headline_variants":["Phase-types extend causal MDBNs to non-Markovian queues","Causal metamodeling for non-Markovian queues using MDBNs","MDBN metamodels for non-Markovian queue causal queries","Phase-type method for MDBN non-Markovian queue metamodels"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"Approximating non-exponential distributions with phase-type distributions preserves enough accuracy for the metamodel to answer probabilistic and causal queries correctly.","fun_headline_variants_meta":{"raw":{"variants":["Phase-types extend causal MDBNs to non-Markovian queues","Causal metamodeling for non-Markovian queues using MDBNs","MDBN metamodels for non-Markovian queue causal queries","Phase-type method for MDBN non-Markovian queue metamodels"]},"model":"grok-4.3","cost_usd":0.010051,"raw_usage":{"total_tokens":4440,"prompt_tokens":624,"num_sources_used":0,"completion_tokens":74,"cost_in_usd_ticks":100512000,"prompt_tokens_details":{"text_tokens":624,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3742,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":624,"tokens_out":74,"duration_ms":42534,"temperature":1.0,"reasoning_tokens":3742,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-28T19:21:32.059755+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Running the MDBN and direct simulation in parallel on the same G/M/1 queue and finding that the MDBN answers to probabilistic causal queries deviate substantially from the simulation results.","supporting_citations":[],"review_version":1}