Research article Special Issues

Modeling of daily confirmed Saudi COVID-19 cases using inverted exponential regression

  • Received: 14 December 2020 Accepted: 24 February 2021 Published: 08 March 2021
  • The coronavirus disease 2019 (COVID-19) pandemic caused by the coronavirus strain has had massive global impact, and has interrupted economic and social activity. The daily confirmed COVID-19 cases in Saudi Arabia are shown to be affected by some explanatory variables that are recorded daily: recovered COVID-19 cases, critical cases, daily active cases, tests per million, curfew hours, maximal temperatures, maximal relative humidity, maximal wind speed, and maximal pressure. Restrictions applied by the Saudi Arabia government due to the COVID-19 outbreak, from the suspension of Umrah and flights, and the lockdown of some cities with a curfew are based on information about COVID-15. The aim of the paper is to propose some predictive regression models similar to generalized linear models (GLMs) for fitting COVID-19 data in Saudi Arabia to analyze, forecast, and extract meaningful information that helps decision makers. In this direction, we propose some regression models on the basis of inverted exponential distribution (IE-Reg), Bayesian (BReg) and empirical Bayesian regression (EBReg) models for use in conjunction with inverted exponential distribution (IE-BReg and IE-EBReg). In all approaches, we use the logarithm (log) link function, gamma prior and two loss functions in the Bayesian approach, namely, the zero-one and LINEX loss functions. To deal with the outliers in the proposed models, we apply Huber and Tukey's bisquare (biweight) functions. In addition, we use the iteratively reweighted least squares (IRLS) algorithm to estimate Bayesian regression coefficients. Further, we compare IE-Reg, IE-BReg, and IE-EBReg using some criteria, such as Akaike's information criterion (AIC), Bayesian information criterion (BIC), deviance (D), and mean squared error (MSE). Finally, we apply the collected data of the daily confirmed from March 23 - June 21, 2020 with the corresponding explanatory variables to the theoretical findings. IE-EBReg shows good model for the COVID-19 cases in Saudi Arabia compared with the other models

    Citation: Sarah R. Al-Dawsari, Khalaf S. Sultan. Modeling of daily confirmed Saudi COVID-19 cases using inverted exponential regression[J]. Mathematical Biosciences and Engineering, 2021, 18(3): 2303-2330. doi: 10.3934/mbe.2021117

    Related Papers:

    [1] EYK Ng, Leonard Jun Cong Looi . Numerical analysis of biothermal-fluids and cardiac thermal pulse of abdominal aortic aneurysm. Mathematical Biosciences and Engineering, 2022, 19(10): 10213-10251. doi: 10.3934/mbe.2022479
    [2] Guoyi Ke, Chetan Hans, Gunjan Agarwal, Kristine Orion, Michael Go, Wenrui Hao . Mathematical model of atherosclerotic aneurysm. Mathematical Biosciences and Engineering, 2021, 18(2): 1465-1484. doi: 10.3934/mbe.2021076
    [3] Rongyanqi Wang, Lan Li, Man Chen, Xiaojuan Li, Yueyun Liu, Zhe Xue, Qingyu Ma, Jiaxu Chen . Gene expression insights: Chronic stress and bipolar disorder: A bioinformatics investigation. Mathematical Biosciences and Engineering, 2024, 21(1): 392-414. doi: 10.3934/mbe.2024018
    [4] Tingting Chen, Wei Hua, Bing Xu, Hui Chen, Minhao Xie, Xinchen Sun, Xiaolin Ge . Robust rank aggregation and cibersort algorithm applied to the identification of key genes in head and neck squamous cell cancer. Mathematical Biosciences and Engineering, 2021, 18(4): 4491-4507. doi: 10.3934/mbe.2021228
    [5] Miao Zhu, Tao Yan, Shijie Zhu, Fan Weng, Kai Zhu, Chunsheng Wang, Changfa Guo . Identification and verification of FN1, P4HA1 and CREBBP as potential biomarkers in human atrial fibrillation. Mathematical Biosciences and Engineering, 2023, 20(4): 6947-6965. doi: 10.3934/mbe.2023300
    [6] Jing Cao, Hui Gan, Han Xiao, Hui Chen, Dan Jian, Ning Jiang, Xuan Zhai . Key protein-coding genes related to microglia in immune regulation and inflammatory response induced by epilepsy. Mathematical Biosciences and Engineering, 2021, 18(6): 9563-9578. doi: 10.3934/mbe.2021469
    [7] Ahmed Hammad, Mohamed Elshaer, Xiuwen Tang . Identification of potential biomarkers with colorectal cancer based on bioinformatics analysis and machine learning. Mathematical Biosciences and Engineering, 2021, 18(6): 8997-9015. doi: 10.3934/mbe.2021443
    [8] Chenyang Jiang, Weidong Jiang . AGTR1, PLTP, and SCG2 associated with immune genes and immune cell infiltration in calcific aortic valve stenosis: analysis from integrated bioinformatics and machine learning. Mathematical Biosciences and Engineering, 2022, 19(4): 3787-3802. doi: 10.3934/mbe.2022174
    [9] Yan Li, Yuzhang Zhu, Guiping Dai, Dongjuan Wu, Zhenzhen Gao, Lei Zhang, Yaohua Fan . Screening and validating the core biomarkers in patients with pancreatic ductal adenocarcinoma. Mathematical Biosciences and Engineering, 2020, 17(1): 910-927. doi: 10.3934/mbe.2020048
    [10] Kewei Ni, Gaozhong Sun . The identification of key biomarkers in patients with lung adenocarcinoma based on bioinformatics. Mathematical Biosciences and Engineering, 2019, 16(6): 7671-7687. doi: 10.3934/mbe.2019384
  • The coronavirus disease 2019 (COVID-19) pandemic caused by the coronavirus strain has had massive global impact, and has interrupted economic and social activity. The daily confirmed COVID-19 cases in Saudi Arabia are shown to be affected by some explanatory variables that are recorded daily: recovered COVID-19 cases, critical cases, daily active cases, tests per million, curfew hours, maximal temperatures, maximal relative humidity, maximal wind speed, and maximal pressure. Restrictions applied by the Saudi Arabia government due to the COVID-19 outbreak, from the suspension of Umrah and flights, and the lockdown of some cities with a curfew are based on information about COVID-15. The aim of the paper is to propose some predictive regression models similar to generalized linear models (GLMs) for fitting COVID-19 data in Saudi Arabia to analyze, forecast, and extract meaningful information that helps decision makers. In this direction, we propose some regression models on the basis of inverted exponential distribution (IE-Reg), Bayesian (BReg) and empirical Bayesian regression (EBReg) models for use in conjunction with inverted exponential distribution (IE-BReg and IE-EBReg). In all approaches, we use the logarithm (log) link function, gamma prior and two loss functions in the Bayesian approach, namely, the zero-one and LINEX loss functions. To deal with the outliers in the proposed models, we apply Huber and Tukey's bisquare (biweight) functions. In addition, we use the iteratively reweighted least squares (IRLS) algorithm to estimate Bayesian regression coefficients. Further, we compare IE-Reg, IE-BReg, and IE-EBReg using some criteria, such as Akaike's information criterion (AIC), Bayesian information criterion (BIC), deviance (D), and mean squared error (MSE). Finally, we apply the collected data of the daily confirmed from March 23 - June 21, 2020 with the corresponding explanatory variables to the theoretical findings. IE-EBReg shows good model for the COVID-19 cases in Saudi Arabia compared with the other models



    Abdominal aortic aneurysm (AAA) has a high prevalence in western countries and leads to a large number of cardiovascular deaths [1]. Aneurysm is defined as an increase in artery diameter by no less than 1.5 times. Aneurysm can appear in any part of arteries in the human body, with AAA being the most common subtype [2]. AAA will cause artery rupture and continuous bleeding, eventually leading to deaths. The mortality rate of AAA after artery rupture is 85~90% [3]. If the diameter of artery is more than 5.5 cm in AAA, surgery is usually performed to prevent artery rupture. AAA is a disease with complex pathological mechanisms, therefore, it is challenging to develop efficient therapeutic strategies to prevent or treat AAA [4]. Although several lines of evidence have suggested that AAA is a multi-gene and multi-factor disease, its pathogenesis is not completely clear [5].

    The very association of arteriosclerosis (AS) and AAA may be up for debate. Previous studies have also indicated that arteriosclerosis (AS) is likely negatively associated with the growth of AAA but positively correlated with its presence [68]. It is widely believed that AAA is merely a consequence/manifestation of AS, but this notion is yet to be fully substantiated. Although revascularization and risk factor control methods are extensively utilized for the treatment of AS-related AAA clinically, the disease is still the most common cause of death. Therefore, in the era of precision medicine, novel strategies are urgently needed to prevent and treat AS-related AAA. However, the biomarkers and therapeutic targets for AS and AAA are yet to be sufficiently explored. Besides, a comprehensive analysis of the associations between AS and AAA at the transcriptomic level is still lacking.

    Herein, we utilized gene expression profiles retrieved from the Gene Expression Omnibus (GEO) database to identify differentially-expressed genes (DEGs) in both AS and AAA patients. The physiopathological functions of these DEGs were then characterized bioinformatically to discover novel diagnostic biomarkers and treatment targets for AS-related AAA.

    This research was supervised and supported by the research ethics committee of the Tianjin Medical University General Hospital. Moreover, this study was conducted in accordance with the declaration of Helsinki. The participants were granted informed consent prior to their enrollments in our study. Artery samples were collected from popliteal arteries of 4 AS and abdominal aortas of 4 AAA patients treated with lower limb amputation and artificial vessel replacement, respectively. Normal artery specimens from visceral arteries of 4 healthy donors were also collected as controls. After arterial resection, total RNA was extracted from tissues rapidly frozen in liquid nitrogen.

    Two GEO (https://www.ncbi.nlm.nih.gov/geo/) datasets, namely GSE100927 (deposited by Marja) and GSE7084 (deposited by Gerard) were employed in this research. The GSE100927 dataset contains 104 samples, including 69 atherosclerotic peripheral artery samples and 35 control samples; while the GSE7084 dataset includes gene expression profiles of 15 abdominal aortas samples obtained from 7 AAA patients and 8 healthy controls.

    We utilized the "affy", "affyPLM", and "limma" packages in R (http://www.bioconductor.org/packages/release/bioc/html/affy.html) [9] to analyze the GSE100927 and GSE7084 datasets. The data were background-adjusted, quantile-normalized, subjected to probe summarization, and log2 transformed, to generate a robust multi-array average (RMA) and log-transformed perfect match (PM) and mismatch (MM) probes. The initial p-values were adjusted by employing the Benjamini-Hochberg method, and fold changes (FCs) were computed on the basis of the false discovery rate (FDR). A gene was considered as a DEG if it had a |log2 FC| > 1 and an adjusted p value < 0.05. Subsequently, DEGs shared by the GSE100927 and GSE7084 datasets were isolated, based on which a Venn diagram was established using Venny 2.1 (http://bioinfogp.cnb.csic.es/tools/venny/index.html).

    To annotate the functions of DEGs shared by the GSE100927 and GSE7084 datasets, we performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses on them using R. We also used the "clusterProfiler" software in R for gene set enrichment analysis (GSEA) [10]. All visualizations were processed by the "ggplot2" package in R.

    The comprehensive information of the hub proteins retrieved from STRING (https://string-db.org/), an online database for the retrieval of inter-gene and inter-protein interactions, was utilized to establish a PPI network, which was visualized using the Cytoscape software (Version 3.7.2) [11]. Subsequently, most important modules within the PPI network were identified by the "Molecular Complex Detection (MCODE)" function of Cytoscape with the following identification criteria: an MCODE score of more than 5, degree and node score cutoff values of 2 and 0.2, respectively, and a maximum depth of 100.

    Based on the cutoff criteria of degree computed by the "cytoHubba" function of Cytoscape, genes that played important roles in the PPI network were the hub genes. The DEGs were scored by the 12 scoring methods of the "cytoHubba" function, and the top ten genes with the highest scores in each scoring method were registered. Based on the 12 ranking lists obtained, upset charts were constructed and the number of times that a gene was registered was counted.

    We extracted total RNA from the above-mentioned artery samples using the TRIzol reagent (TransGen Biotech, China). The total RNA was then reversely transcribed into cDNA using a real-time PCR mRNA detection kit (TIANGEN Biotech, China) before PCR amplification on a Roche Photorecycler 480 Real-time PCR system. The reaction consisted of the following steps: one step of 95˚C for 30 sec and 40 rounds of 95˚C for 5 sec and 60˚C for 20 sec. GAPDH was included as the internal reference. Relative transcript levels of target genes were determined by the 2-ΔΔCt method. The PCR primers were designed and synthesized by AoKe Dingsheng (Beijing, China) and their sequences are shown in Table 1.

    Table 1.  Primer sequences for qRT-PCR.
    Gene Primer sequence (5′→3′)
    SPI1 F: GTGCCCTATGACACGGATCTA R: AGTCCCAGTAATGGTCGCTAT
    TYROBP F: ACTGAGACCGAGTCGCCTTAT R: ATACGGCCTCTGTGTGTTGAG
    TLR2 F: ATCCTCCAATCAGGCTTCTCT R: GGACAGGTCAAGGCTTTTTACA
    FCER1G F: AGCAGTGGTCTTGCTCTTACT R: TGCCTTTCGCACTTGGATCTT
    MMP9 F: TGTACCGCTATGGTTACACTCG R: GGCAGGGACAGTTGCTTCT
    GAPDH F: ATGGCTACAGCAACAGGGT R: TTATGGGGTCTGGGATGG

     | Show Table
    DownLoad: CSV

    GraphPad Prism 8 was used to compare data from test and control samples. To analyze the clinical characteristics of AS patients, AAA patients and healthy controls, the quantitative data are expressed as means ± standard deviations (SDs). The Student's t-test and Wilcoxon rank sum test were implemented for analysis of data as appropriate. P values less than 0.05 were indicative of statistical significance. All statistical analyses in this study were conducted with the SPSS 26.0 software (IBM, USA).

    Microarray data from both GEO datasets were normalized and subjected to principal component analysis (PCA) for evaluating differences in biological characteristics between samples (Figure 1A, B), thereby ensuring data accuracy. Genes with a |log2 FC| > 1 and an adjusted p value < 0.05 were identified using the "Limma" package in R. Compared with controls, 324 and 477 overexpressed and 113 and 370 underexpressed genes were identified in the peripheral artery samples of AS and AAA patients in the GSE100927 and GSE7084 datasets, respectively (Figure 2). The DEGs are also presented by volcano plots and heat maps shown in Figure 1CF. The DEGs shared by the GSE100927 and GSE7084 datasets were exhibited by a Venn diagram (Figure 3A), including 169 upregulated and 37 downregulated DEGs.

    Figure 1.  Identification of differentially expressed mRNAs analysis. Principal component analysis (PCA) graph in (A) GSE7084 and (B) GSE100927. Volcano plot of differentially expressed mRNAs in (C) GSE7084 and (D) GSE100927. The red dot represents upregulated mRNAs and the blue dot represents downregulated mRNAs. A heatmap of 85 differentially expressed genes between (E) AAA patients and normal individuals, (F) AS patients and normal individuals. Red represents upregulated genes, and green represents downregulated genes.
    Figure 2.  Flow diagram of the study approach. mRNA microarray analyses were performed on aortic specimens obtained from GSE100927 (69 AS patients and 35 control group) and GSE7084 (7 AAA patients and 8 control group). From among differentially-expressed genes (DEGs), network analysis of gene expression both in AAA and AS identifies 206 co-DEGs. Afterwards, functional and pathway enrichment analysis were performed for the co-DEGs. The top 5 genes with the highest scores in cytoHubba module were taken as hub genes. At last, validation quantitative polymerase chain reaction (qRT-PCR) was performed for the 5 hub genes in samples obtained from our institute, including 4 organ donors, 4 AAA patients and 4 AS patients.
    Figure 3.  Enrichment analysis of the DEGs. (A) A Venn-diagram of DEGs between GSE100927 and GSE7084. (B) GO functional analysis results illustrating the significantly enriched terms for the co-DEGs. (C) Results of the KEGG enrichment analysis of the co-DEGs. DEGs, differentially expressed genes.

    The functions of the identified DEGs were explored via GO and KEGG enrichment analyses in R. The results suggested that the DEGs shared by the GSE100927 and GSE7084 datasets were significantly enriched with biological processes (BP) terms of "neutrophil activation", "immune response", and "neutrophil degranulation", a molecular function (MF) term of "actin binding", as well as a cellular components (CC) term of "secretory granule membrane" in the GO analysis (Figure 3B). Our KEGG analysis suggested that these shared DEGs were associated with 10 pathways, such as "lipid and atherosclerosis", "phagosome", and "osteoclast differentiation" (Figure 3C).

    To characterize the interactions between proteins encoded by these shared DEGs, a PPI network was built by STRING (Version 11.0) and visualized by Cytoscape (Figure 4A). The network consisted of 190 nodes (representing 190 proteins) and 1715 edges (representing 1715 interactions between the proteins). In addition, 3 essential modules were identified from the network by MCODE in Cytoscape (Table 2). According to the upset chart, five genes (namely SPI1, TYROBP, TLR2, FCER1G and MMP9) were registered more than 5 times and were regarded as hub genes (Figure 4B). Interestingly, these genes were also included in the modules with the highest MCODE scores (Figure 4C).

    Figure 4.  PPI network of co-DEGs constructed by the STRING database. (A) The PPI network of co-DEGs was constructed using Cytoscape. (B) Among the 12 scoring methods of the cytoHubba module, the intersection of the top ten genes in each group. (C) Identification of a sub-network using MCODE in Cytoscape software. Diamond nodes, hub genes; Red nodes, up-regulated genes; Blue nodes, down-regulated genes.
    Table 2.  Module analysis of DEGs using Cytoscape.
    Cluster Score Nodes Edges Node IDs
    1 19.6 21 196 LAPTM5, FCER1G, LY86, CCR1, CD53, FCGR2A, FCGR3A, CSF1R, TYROBP, C1QA, SPI1, IL10RA, LILRB2, CYBB, CTSS, LCP2, MNDA, CD74, MPEG1, AIF1, ITGAM
    2 8.5 21 85 CCL5, ITGAL, CCL3, PARVG, CCR5, CCR7, IGSF6, BTK, CD83, TLR2, IL1B, CSF3R, HLA-DRA, VAV1, CD68, HK3, FGR, CD14, MYO1F, CCL4, NCF2
    3 6.207 30 90 EVI2B, MS4A6A, CCRL2, IL7R, VAMP8, SYK, APOE, LILRB3, STK10, ADAM8, HCST, LYN, HMOX1, INPP5D, SLA, ALOX5, CXCR4, CORO1A, MMP9, CXCL16, NCF4, PTPN6, TREM1, IL1RN, HCLS1, SLC2A5, AQP9, PLAU, RAC2, VAV3

     | Show Table
    DownLoad: CSV

    Differential expression of the DEGs, as revealed by bioinformatic methods, was further confirmed by qRT-PCR using cDNA obtained from artery samples from AS patients (n = 4), AAA patients (n = 4), and healthy controls (n = 4). As shown in Figure 5, SPI1, TYROBP, TLR2, FCER1G and MMP9 were significantly upregulated both in AS and AAA samples compared with normal controls.

    Figure 5.  Validation of the gene expression levels of hub genes between patients and healthy control group. (A–E) The expression of SPI1, TYROBP, TLR2, FCER1G and MMP9 in ASO and AAA arteries compared with normal arteries. The data are presented as the mean ± SEM, **P < .01 vs control group; ***P < .001 vs control group.

    A large number of clinical investigations h.ave comprehensively exhibited an enhanced risk of AAA in patients with AS. However, the mechanism underlying the association between AS and AAA remains elusive, although it has been speculated to be associated with chronic inflammation, degradation of extracellular matrix, apoptosis of vascular smooth muscle cells, and thrombosis [12], physiological characteristics that also contribute to the formation of atherosclerotic plaques, suggesting that AS, which frequently affects the aneurysm wall, is a crucial risk factor for AAA. We herein discovered 5 key genes implicated in AS-related AAA, namely SPI1, TYROBP, TLR2, FCER1G and MMP9. These genes have the potential to become novel biomarkers and treatment targets for AS-related AAA. Our findings not only provide further insights into the pathogenic mechanisms of AS-related AAA, but also have clinical significance—they suggest several potential early diagnostic markers to prevent the progression of AAA.

    Among the identified hub genes, SPI1 encodes PU.1, an ETS-domain transcription factor that is required for the development of myeloid cells and can regulate intercellular communications to facilitate immune responses [13]. The role of the SPI1 in AAA progression is unclear. According to previous findings, the expression of SPI1 gene is upregulated during the differentiation of myeloid cells and remains persistently high in B cells, monocytes, mast cells, circulating neutrophils, and other granulocytes in human [14]. Pan et al. found that SPI1 expression was significantly increased in Tibetan minipigs subjected to a high-fat/cholesterol diet [15], further implying that the activation of myeloid cells is critical for AS development. TYROBP, which has a synonym of DNAX-activating protein of 12 kDa (DAP12), encodes a transmembrane receptor ubiquitously expressed in neutrophils, natural killer cells, and monocytes/macrophages that transduces immune signals via a tyrosine-based activation motif [16]. DAP12 has been extensively implicated in the survival, reproduction, polarization, and differentiation of multiple types of immune cells, particularly the monocytes/macrophages [16]. Wang et al. observed that DAP12 was expressed at a high level in APOE mice, which facilitated the development of AS plaques via the TREM-1/DAP12 pathway [17]. Hinterseher et al. observed that TYROBP was highly expressed in human AAA and was involved in the pathogenesis of AAA through the NK cell-mediated cytotoxicity pathway [18]. TLR2 encodes a cell membrane protein that identifies microorganism invasions and regulates innate immune reactions [19]. Jabłońska et al. found significantly higher levels of TNF-α, IL-4 and TLR2 in AAA patients compared with healthy controls, and suggested that TLR2 may be involved in the aortic and systemic inflammatory responses in AAA patients [20]. Lee et al. showed that TLR2 could activate p38 and ERK1/2 signaling pathways, selectively regulate IL-6–mediated RANKL upregulation and OPG downregulation, and induce the differentiation of VSMCs into chondrocytes, thereby resulting in vascular calcification during AS development [21]. FCER1G (CD23) encodes a high-affinity cell membrane receptor for IgE, which is upregulated during aging [22]. FCER1G was demonstrated to be able to regulate hypertension and platelet hyperaggregability during thrombus formation [2325]. Given that AAA is tightly associated with these diseases, it is reasonable that FCER1G may be implicated in the development of AAA. MMP9 encodes a matrix metalloproteinase that reorganizes and degrades extracellular matrix during collagen catabolic processes, leading to the instability and rupture of atherosclerotic plaques [26,27]. This molecule activates the immune responses that lead to atherosclerotic plaque destabilization by facilitating the degranulation, transendothelial migration, and differentiation of leukocytes, as well as activating interleukin signaling pathways. However, its role in AAA progression remains controversial. Hurks et al. found that the MMP9 levels remained stable in a study of inflammation in AAA artery wall [28]. Another study found that MMP9 level in AAA was correlated with intraluminal thrombus amount, implying that thrombus activated MMP9 [29].

    In this study, we identified DEGs shared by the AS and AAA cohorts and discovered 5 hub genes (namely SPI1, TYROBP, TLR2, FCER1G and MMP9) by PPI network construction and identification of key modules. The relationship between these genes, AS and AAA needs to be further investigated in future studies.

    By analyzing GEO datasets using bioinformatic methods and verifying the results by in vitro experiments, we not only identified 206 DEGs shared by both AS and AAA cohorts, but also detected 5 hub genes (SPI1, TYROBP, TLR2, FCER1G and MMP9) linking AS and AAA, which may be used as potential biomarkers for the diseases. More importantly, these hub genes may become potential treatment targets for AS-related AAA, which will facilitate the establishment of new therapeutic strategies, including early prevention and precision treatment.

    Our study has some limitations. Firstly, our research was based on microarray data retrieved from a public database. Secondly, since transcript abundances of genes may not be able to accurately reflect the corresponding protein levels, for example, we do not know which type of cells are mainly involved in the expression of these hub genes, and whether the complex in vivo environment will mask some of more meaningful genes, additional in vitro and in vivo studies are required to characterize the functions of the key genes. Therefore, large-scale prospective clinical investigations can, to some extent, better support our findings.

    We thank Yudi Liu for technical assistance.

    This study was supported by the National Natural Science Foundation of China [Grant Numbers 82070489].

    This study was approved on June, 2021 by the Institutional Review Board of the Tianjin Medical University General Hospital with an approval number of IRB2021-WZ-080.

    The authors declare that they have no competing interests.



    [1] E. Tello-Leal, B. A. Macias-Hernandez, Association of environmental and meteorological factors on the spread of COVID-19 in Victoria, Mexico, and air quality during the lockdown, Environ. Res., (2020), 110442.
    [2] S. Kodera, E. A. Rashed, A. Hirata, Correlation between COVID-19 morbidity and mortality rates in Japan and local population density, temperature, and absolute humidity, Int. J. Env. Res. Pub. He., 17 (2020), 5477. doi: 10.3390/ijerph17155477
    [3] S. A. Meo, A. A. Abukhalaf, A. A. Alomar, N. M. Alsalame, T. Al-Khlaiwi, A. M. Usmani, Effect of temperature and humidity on the dynamics of daily new cases and deaths due to COVID-19 outbreak in Gulf countries in Middle East Region, Eur. Rev. Med. Pharmacol. Sci., 24 (2020), 7524-7533.
    [4] L. A. Casado-Aranda, J. Sanchez-Fernandez, M. I. Viedma-del-Jesus, Analysis of the scientific production of the effect of COVID-19 on the environment: A bibliometric study, Environ. Res., (2020), 110416.
    [5] B. Dogan, M. B. Jebli, K. Shahzad, T. H. Farooq, U. Shahzad, Investigating the effects of meteorological parameters on COVID-19: Case study of New Jersey, United States, Environ. Res., 191 (2020), 110148. doi: 10.1016/j.envres.2020.110148
    [6] S. A. Meo, A. A. Abukhalaf, A. A. Alomar, O. M. Alessa, W. Sami, D. C. Klonoff, Effect of environmental pollutants PM-2.5, carbon monoxide, and ozone on the incidence and mortality of SARS-COV-2 infection in ten wildfire affected counties in California, Sci. Total Environ., 757 (2021), 143948. doi: 10.1016/j.scitotenv.2020.143948
    [7] J. Yuan, Y. Wu, W. Jing, J. Liu, M. Du, Y. Wang, et al., Non-linear correlation between daily new cases of COVID-19 and meteorological factors in 127 countries, Environ. Res., 193 (2021), 110521. doi: 10.1016/j.envres.2020.110521
    [8] P. McCullagh, J. A. Nelder, Generalized Linear Models, 1st edition, Chapman and Hall, London, 1983.
    [9] J. A. Nelder, D. Pregibon, An extended quasi-likelihood function, Biometrika, 74 (1987), 221-232. doi: 10.1093/biomet/74.2.221
    [10] J. A. Nelder, R. W. M. Wedderburn, Generalized linear models, J. R. Stat. Soc. Ser. A, 135 (1972), 370-384.
    [11] K. H. Yuan, P. M. Bentler, Improving the convergence rate and speed of Fisher-scoring algorithm: ridge and anti-ridge methods in structural equation modeling, Ann. Inst. Stat. Math., 69 (2017), 571-597. doi: 10.1007/s10463-016-0552-2
    [12] P. De Jong, G. Z. Heller, Generalized linear models for insurance data, 1st edition, Cambridge Books, 2008.
    [13] T. F. Liao, Interpreting Probability Models: Logit, Probit, and Other Generalized Linear Models, No 07-101, SAGE Publications, Thousand Oaks, 1994.
    [14] R. Richardson, B. Hartman, Bayesian nonparametric regression models for modeling and predicting healthcare claims, Insur. Math. Econ., 83 (2018), 1-8. doi: 10.1016/j.insmatheco.2018.06.002
    [15] C. Song, Y. Wang, X. Yang, Y. Yang, Z. Tang, X. Wang, et al., Spatial and Temporal Impacts of Socioeconomic and Environmental Factors on Healthcare Resources: A County-Level Bayesian Local Spatiotemporal Regression Modeling Study of Hospital Beds in Southwest China, Int. J. Env. Res. Pub. He., 17 (2020), 5890. doi: 10.3390/ijerph17165890
    [16] Y. Mohamadou, A. Halidou, P. T. Kapen, A review of mathematical modeling, artificial intelligence and datasets used in the study, prediction and management of COVID-19, Appl. Intell., 50 (2020), 3913-3925. doi: 10.1007/s10489-020-01770-9
    [17] T. A. Trunfio, A. Scala, A. D. Vecchia, A. Marra, A. Borrelli, Multiple Regression Model to Predict Length of Hospital Stay for Patients Undergoing Femur Fracture Surgery at "San Giovanni di Dio e Ruggi d'Aragona" University Hospital, In European Medical and Biological Engineering Conference, Springer, Cham, (2020), 840-847.
    [18] A. Z. Keller, A. R. R. Kamath, U. D. Perera, Reliability analysis of CNC machine tools, Reliab. Eng., 3 (1982), 449-473. doi: 10.1016/0143-8174(82)90036-1
    [19] Y. Abdel-Aty, A. Shafay, M. M. M. El-Din, M. Nagy, Bayesian inference for the inverse exponential distribution based on pooled type-II censored samples, J. Stat. Appl. Pro., 4 (2015), 235.
    [20] S. Dey, Inverted exponential distribution as a life distribution model from a Bayesian viewpoint, Data Sci. J., 6 (2007), 107-113. doi: 10.2481/dsj.6.107
    [21] S. K. Singh, U. Singh, A. S. Yadav, P. K. Vishwkarma, On the estimation of stress strength reliability parameter of inverted exponential distribution, IJSW, 3 (2015), 98-112. doi: 10.14419/ijsw.v3i1.4329
    [22] L. Fahrmeir, G. Tutz, Multivariate Statistical Modelling Based on Generalized Linear Models, 2nd edition, Springer Science and Business Media, Berlin/Heidelberg, 2013.
    [23] E. Cepeda, D. Gamerman, Bayesian methodology for modeling parameters in the two parameter exponential family, Rev. Estad., 57 (2015), 93-105.
    [24] D. K. Dey, S. K. Ghosh, B. K. Mallick, Generalized Linear Models: A Bayesian Perspective, 1st edition, CRC Press, New York, 2000.
    [25] U. Olsson, Generalized Linear Models, An Applied Approach, 1st edition, Student Litteratur Lund., Sweden, 2002.
    [26] N. Sano, H. Suzuki, M. Koda, A robust ensemble learning using zero-one loss function, J. Oper. Res. Soc. Japan, 51 (2008), 95-110.
    [27] H. Robbins, An empirical Bayes approach to statistics, In Breakthroughs in statistics, Springer, (1955), 388-394.
    [28] L. Wei, Empirical Bayes test of regression coefficient in a multiple linear regression model, Acta Math. Appl. Sin-E, 6 (1990), 251-262. doi: 10.1007/BF02019151
    [29] R. S. Singh, Empirical Bayes estimation in a multiple linear regression model, Ann. Inst. Stat. Math., 37 (1985), 71-86. doi: 10.1007/BF02481081
    [30] W. M. Houston, D. J. Woodruff, Empirical Bayes Estimates of Parameters from the Logistic Regression Model, ACT Res. Report Ser., (1997), 97-96.
    [31] S. L. Wind, An empirical Bayes approach to multiple linear regression, Ann. Stat., 1 (1973), 93-103. doi: 10.1214/aos/1193342385
    [32] S. Y. Huang, Empirical Bayes testing procedures in some nonexponential families using asymmetric Linex loss function, J. Stat. Plan. Infer., 46 (1995), 293-305. doi: 10.1016/0378-3758(94)00112-9
    [33] R. J. Karunamuni, Optimal rates of convergence of empirical Bayes tests for the continuous one-parameter exponential family, Ann. Stat., (1996), 212-231.
    [34] M. Yuan, Y. Lin, Efficient empirical Bayes variable selection and estimation in linear models, J. Am. Stat. Assoc., 100 (2005), 1215-1225. doi: 10.1198/016214505000000367
    [35] L. S. Chen, Empirical Bayes testing for a nonexponential family distribution, Commun. Stat., Theor. M., 36 (2007), 2061-2074. doi: 10.1080/03610920601143675
    [36] B. Efron, Large-scale inference: Empirical Bayes methods for estimation, testing, and prediction, 1st, Cambridge University Press, 2012.
    [37] M. Shao, An empirical Bayes test of parameters for a nonexponential distribution family with Negative Quadrant Dependent random samples, In 2013 10th International Conference on Fuzzy Systems and Knowledge Discovery (FSKD), IEEE, (2013), 648-652.
    [38] J. E. Kim, D. A. Nembhard, Parametric empirical Bayes estimation of individual time-pressure reactivity, Int. J. Prod. Res., 56 (2018), 2452-2463. doi: 10.1080/00207543.2017.1380321
    [39] K. Jampachaisri, K. Tinochai, S. Sukparungsee, Y. Areepong, Empirical Bayes Based on Squared Error Loss and Precautionary Loss Functions in Sequential Sampling Plan, IEEE Access, 8 (2020), 51460-51465. doi: 10.1109/ACCESS.2020.2979872
    [40] Y. Li, L. Hou, Y. Yang, J. Tong, Huber's M-Estimation-Based Cubature Kalman Filter for an INS/DVL Integrated System, Math. Probl. Eng., (2020), 2020.
    [41] B. Sinova, S. Van Aelst, Advantages of M-estimators of location for fuzzy numbers based on Tukey's biweight loss function, Int. J. Approx. Reason., 93 (2018), 219-237. doi: 10.1016/j.ijar.2017.10.032
    [42] P. McCullagh, J. A. Nelder, Generalized Linear Models, 2nd edition, Chapman and Hall/CRC, 1985.
    [43] S. Das, D. K. Dey, On Bayesian analysis of generalized linear models using the Jacobian technique, Am. Stat., 60 (2006), 264-268. doi: 10.1198/000313006X128150
    [44] S. Ferrari, F. Cribari-Neto, Beta regression for modelling rates and proportions, J. Appl. Stat., 31 (2004), 799-815. doi: 10.1080/0266476042000214501
    [45] S. Das, D. K. Dey, On Bayesian analysis of generalized linear models: A new perspective, Technical Report, Statistical and Applied Mathematical Sciences Institute, Research Triangle Park, (2007), 33.
    [46] P. J. Huber, Robust estimation of a location parameter, Ann. Math. Stat., 35 (1964), 73-101. doi: 10.1214/aoms/1177703732
    [47] P. J. Rousseeuw, A. M. Leroy, Robust Regression and Outlier Detection, 1st edition, John Wiley and Sons, NY, 1987.
    [48] L. Chang, B. Hu, G. Chang, A. Li, Robust derivative-free Kalman filter based on Huber's M-estimation methodology, J. Process Control, 23 (2013), 1555-1561. doi: 10.1016/j.jprocont.2013.05.004
    [49] P. J. Huber, Robust Statistics, 1st edition, John Wiley and Sons, NY, 1981.
    [50] R. A. Maronna, R. D. Martin, V. J. Yohai, Robust Statistics: Theory and Methods, 1st edition, John Wiley and Sons, West Sussex, 2006.
    [51] F. Wen, W. Liu, Iteratively reweighted optimum linear regression in the presence of generalized Gaussian noise, In 2016 IEEE International Conference on Digital Signal Processing (DSP), IEEE, (2016), 657-661.
    [52] H. Kikuchi, H. Yasunaga, H. Matsui, C. I. Fan, Efficient privacy-preserving logistic regression with iteratively Re-weighted least squares, In 2016 11th Asia Joint Conference on Information Security (AsiaJCIS), IEEE, (2016), 48-54.
    [53] J. Tellinghuisen, Least squares with non-normal data: Estimating experimental variance functions, Analyst, 133(2) (2008), 161-166.
    [54] R. M. Leuthold, On the use of Theil's inequality coefficients, Am. J. Agr. Econ., 57 (1975), 344-346. doi: 10.2307/1238512
    [55] T. Niu, L. Zhang, B. Zhang, B. Yang, S. Wei, An Improved Prediction Model Combining Inverse Exponential Smoothing and Markov Chain, Math. Probl. Eng., 2020 (2020), 11.
    [56] J. J. Faraway, Extending the linear model with R: Generalized linear, mixed effects and nonparametric regression models, 2nd edition, CRC press, 2016.
    [57] J. Fox, S. Weisberg, An R companion to applied regression, 3rd edition, Sage publications, Inc., 2018.
    [58] E. Dikici, F. Orderud, B. H. Lindqvist, Empirical Bayes estimator for endocardial edge detection in 3D+ T echocardiography, In 2012 9th IEEE International Symposium on Biomedical Imaging (ISBI), IEEE, (2012), 1331-1334.
    [59] A. Coluccia, F. Ricciato, Improved estimation of instantaneous arrival rates via empirical Bayes, In 2014 13th Annual Mediterranean Ad Hoc Networking Workshop, IEEE, (2014), 211-216.
  • This article has been cited by:

    1. Xuezhi Du, Han Qi, Wenbin Ji, Peiyuan Li, Run Hua, Wenliang Hu, Feng Qi, Kamal Shah, Construction of a Colorectal Cancer Prognostic Risk Model and Screening of Prognostic Risk Genes Using Machine-Learning Algorithms, 2022, 2022, 1748-6718, 1, 10.1155/2022/9408839
    2. Jan Vincent B. Beltran, Fang-Ping Lin, Chaw-Liang Chang, Tai-Ming Ko, Single-Cell Meta-Analysis of Neutrophil Activation in Kawasaki Disease and Multisystem Inflammatory Syndrome in Children Reveals Potential Shared Immunological Drivers, 2023, 148, 0009-7322, 1778, 10.1161/CIRCULATIONAHA.123.064734
    3. Xiaohong Song, Zhaoran Chen, Yuehong Zheng, Jianqiang Wu, Urine proteomics in cardiovascular diseases: advances in biomarker discovery and clinical applications, 2024, 1478-9450, 10.1080/14789450.2024.2436401
    4. Jianqing Deng, Lei Tian, Haitao Chi, Lei Chen, Junhui Wang, Yan Xue, Qiang Zhao, Nan Zheng, Jie Dong, Jiaying Li, Wei Guo, Cangsong Xiao, Ming Yang, Establishment of a murine chronic proximal thoracic aortic aneurysm model by combining periaortic elastase application with oral BAPN administration, 2025, 999, 00142999, 177678, 10.1016/j.ejphar.2025.177678
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3504) PDF downloads(186) Cited by(1)

Figures and Tables

Figures(6)  /  Tables(8)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog