I congratulate both teams for these welcome contributions on modelling the Covid-19 pandemic. To produce results of such quality within exacting timescales is a genuine achievement. Both studies infer a time-varying reproduction number R t from summary data by construct-ing hierarchical Bayesian frameworks embodying R t as an intrinsic parameter. Observations arise as noisy, time-shifted representations of an autoregressive infection process with weights specified by generation-time probabilities and moderated by R t . With a common root in Flaxman et al. (2020), the papers differ in their treatment of temporal effects and spatial cou-pling (with Teh et al. (2022) adopting an explicitly spatio-temporal Gaussian process for log R t while Mishra et al. (2022) use a random walk prior), in their use of data, and in underlying assumptions. Neither study, in the prior for R t , incorporates foreseeable effects such as step changes follow-ing interventions, the impact of improved testing on track-and-trace measures, or the expected decline in R t due to susceptible depletion. Incidentally, the presentation of the infection model in Mishra et al. (2022) seems confusing, with R t between equations (1) and (2) changing from an instantaneous reproduction number to a ‘raw’ reproduction number, subsequently re-scaled by the susceptible proportion before reporting. The papers’ general approach is arguably the ‘image analyst’s take’ on epidemic modelling, where the objective is to recover a ‘true’ R t from a noisy image, with prior distributions providing regularisation rather than capturing mechanistic thinking. This approach differs
我祝贺这两个团队在模拟Covid-19大流行方面做出的可喜贡献。在严格的时间尺度内产生如此高质量的结果是一项真正的成就。两项研究都通过构建层次贝叶斯框架,从汇总数据中推断出时变的再现数R t $$ {R}_t $$R t $$ {R}_t $$作为内在参数。观察结果是自回归感染过程的噪声时移表示,其权重由代时间概率指定,并由R t $$ {R}_t $$调节。与Flaxman等人(2020)的共同根源,这两篇论文在处理时间效应和空间耦合方面有所不同(Teh等人(2022)对log R t采用了明确的时空高斯过程$$ log {R}_t $$,而Mishra等人则采用了不同的方法。(2022)使用随机漫步先验),在数据的使用和潜在的假设中。在R t $$ {R}_t $$之前的研究中,这两项研究都没有纳入可预见的影响,例如干预后的阶跃变化,改进测试对跟踪和跟踪措施的影响,或R t $$ {R}_t $$由于易感耗竭的预期下降。顺便提一下,Mishra等人(2022)对感染模型的描述似乎令人困惑,方程(1)和(2)之间的R t $$ {R}_t $$从瞬时繁殖数变为“原始”繁殖数,随后由报告前的敏感比例重新缩放。论文的一般方法可以说是“图像分析师对流行病建模的看法”,其目标是从噪声图像中恢复“真实”的R t $$ {R}_t $$,先验分布提供正则化,而不是捕获机械思维。这种方法不同于植物或动物病原体建模者经常采用的方法,后者旨在估计控制传播过程不同方面的参数,例如接触率和空间核函数,然后将“机制”理解外推到其他环境。R t $$ {R}_t $$这个可以定义的量是传播过程和假定的监测和控制策略的副产品,而不是一个内在参数。 , 2019)?例如,当模拟具有繁殖矩阵R t $$ {mathbf{R}}_t $$的结构化种群的类似数据时,输入的R t $$ {R}_t $$能成功地跟踪真实R的最大特征值吗T $$ {mathbf{R}}_t $$,或者它可能低估了这个数量,因为组间的感染分布可能与相应的特征向量不匹配?与更简单的平滑方法进行比较也是受欢迎的。这两篇论文强调了流行病统计建模的一个重要挑战——对更复杂的机制模型的统计推断,这些模型可能为有针对性的控制策略的设计提供信息。这就要求可用数据的丰富程度与模型的复杂性更好地匹配;实现这样的匹配本身就是一个重大挑战。这些论文的作者有效地利用了现有数据,他们的建模是理解空间相互作用影响的重要一步。探索他们的框架是否延伸到其他异质性将是有趣的,例如年龄结构引起的异质性,其重要性在其他研究中已经得到强调(例如Lau等人,2020)。
{"title":"Gavin J. Gibson's invited discussion contribution to the papers in Session 2 of the Royal Statistical Society's Special Topic Meeting on Covid-19 Transmission: 11 June 2021","authors":"Gavin J. Gibson","doi":"10.1111/rssa.12972","DOIUrl":"10.1111/rssa.12972","url":null,"abstract":"I congratulate both teams for these welcome contributions on modelling the Covid-19 pandemic. To produce results of such quality within exacting timescales is a genuine achievement. Both studies infer a time-varying reproduction number R t from summary data by construct-ing hierarchical Bayesian frameworks embodying R t as an intrinsic parameter. Observations arise as noisy, time-shifted representations of an autoregressive infection process with weights specified by generation-time probabilities and moderated by R t . With a common root in Flaxman et al. (2020), the papers differ in their treatment of temporal effects and spatial cou-pling (with Teh et al. (2022) adopting an explicitly spatio-temporal Gaussian process for log R t while Mishra et al. (2022) use a random walk prior), in their use of data, and in underlying assumptions. Neither study, in the prior for R t , incorporates foreseeable effects such as step changes follow-ing interventions, the impact of improved testing on track-and-trace measures, or the expected decline in R t due to susceptible depletion. Incidentally, the presentation of the infection model in Mishra et al. (2022) seems confusing, with R t between equations (1) and (2) changing from an instantaneous reproduction number to a ‘raw’ reproduction number, subsequently re-scaled by the susceptible proportion before reporting. The papers’ general approach is arguably the ‘image analyst’s take’ on epidemic modelling, where the objective is to recover a ‘true’ R t from a noisy image, with prior distributions providing regularisation rather than capturing mechanistic thinking. This approach differs","PeriodicalId":49983,"journal":{"name":"Journal of the Royal Statistical Society Series A-Statistics in Society","volume":"185 S1","pages":"S96-S98"},"PeriodicalIF":2.0,"publicationDate":"2022-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://rss.onlinelibrary.wiley.com/doi/epdf/10.1111/rssa.12972","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41366942","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"John Kingman’s invited discussion contribution to the papers in Session 1 of the Royal Statistical Society’s Special Topic Meeting on COVID-19 Transmission: 9 June 2021","authors":"John Kingman","doi":"10.1111/rssa.12886","DOIUrl":"10.1111/rssa.12886","url":null,"abstract":"","PeriodicalId":49983,"journal":{"name":"Journal of the Royal Statistical Society Series A-Statistics in Society","volume":"185 S1","pages":"S41-S43"},"PeriodicalIF":2.0,"publicationDate":"2022-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41911124","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
When disaggregation of national estimates in several domains or areas is required, direct survey estimators, which use only the domain-specific survey data, are usually design-unbiased even under complex survey designs (at least approximately) and require no model assumptions. Nevertheless, they are appropriate only for domains or areas with sufficiently large sample size. For example, when estimating poverty in a domain with a small sample size (small area), the volatility of a direct estimator might make that area seems like very poor in one period and very rich in the next one. Small area (or indirect) estimators have been developed in order to avoid such undesired instability. Small area estimators borrow strength from the other areas so as to improve the precision and therefore obtain much more stable estimators. However, the usual model-based assumptions, which include some kind of area homogeneity, may not hold in real applications. A more flexible model based on multivariate mixtures of normal distributions that generalises the usual nested error linear regression model is proposed for estimation of general parameters in small areas. This flexibility makes the model adaptable to more general situations, where there may be areas with a different behaviour from the other ones, making the model less restrictive (hence, more close to nonparametric) and more robust to outlying areas. An expectation-maximisation (E-M) method is designed for fitting the proposed mixture model. Under the proposed mixture model, two different new predictors of general small area indicators are proposed. A parametric bootstrap method is used to estimate the mean squared errors of the proposed predictors. Small sample properties of the new predictors and of the bootstrap procedure are analysed by simulation studies and the new methodology is illustrated with an application to poverty mapping in Palestine.
{"title":"Multivariate mixture model for small area estimation of poverty indicators","authors":"Agne Bikauskaite, Isabel Molina, Domingo Morales","doi":"10.1111/rssa.12965","DOIUrl":"10.1111/rssa.12965","url":null,"abstract":"<p>When disaggregation of national estimates in several domains or areas is required, direct survey estimators, which use only the domain-specific survey data, are usually design-unbiased even under complex survey designs (at least approximately) and require no model assumptions. Nevertheless, they are appropriate only for domains or areas with sufficiently large sample size. For example, when estimating poverty in a domain with a small sample size (small area), the volatility of a direct estimator might make that area seems like very poor in one period and very rich in the next one. Small area (or indirect) estimators have been developed in order to avoid such undesired instability. Small area estimators borrow strength from the other areas so as to improve the precision and therefore obtain much more stable estimators. However, the usual model-based assumptions, which include some kind of area homogeneity, may not hold in real applications. A more flexible model based on multivariate mixtures of normal distributions that generalises the usual nested error linear regression model is proposed for estimation of general parameters in small areas. This flexibility makes the model adaptable to more general situations, where there may be areas with a different behaviour from the other ones, making the model less restrictive (hence, more close to nonparametric) and more robust to outlying areas. An expectation-maximisation (E-M) method is designed for fitting the proposed mixture model. Under the proposed mixture model, two different new predictors of general small area indicators are proposed. A parametric bootstrap method is used to estimate the mean squared errors of the proposed predictors. Small sample properties of the new predictors and of the bootstrap procedure are analysed by simulation studies and the new methodology is illustrated with an application to poverty mapping in Palestine.</p>","PeriodicalId":49983,"journal":{"name":"Journal of the Royal Statistical Society Series A-Statistics in Society","volume":"185 S2","pages":"S724-S755"},"PeriodicalIF":2.0,"publicationDate":"2022-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://rss.onlinelibrary.wiley.com/doi/epdf/10.1111/rssa.12965","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45493631","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
I congratulate: Parag, Thompson, and Donnelly; Jewell and Lewnard; and Coffeng and de Vlas on their papers which highlight both the benefits and potential pitfalls associated with statistics such as the doubling time