In this article, the problem of estimating the covariance matrix in general linear mixed models is considered. Two new classes of estimators obtained by shrinking the eigenvalues towards the origin and the arithmetic ...In this article, the problem of estimating the covariance matrix in general linear mixed models is considered. Two new classes of estimators obtained by shrinking the eigenvalues towards the origin and the arithmetic mean, respectively, are proposed. It is shown that these new estimators dominate the unbiased estimator under the squared error loss function. Finally, some simulation results to compare the performance of the proposed estimators with that of the unbiased estimator are reported. The simulation results indicate that these new shrinkage estimators provide a substantial improvement in risk under most situations.展开更多
Today, Linear Mixed Models (LMMs) are fitted, mostly, by assuming that random effects and errors have Gaussian distributions, therefore using Maximum Likelihood (ML) or REML estimation. However, for many data sets, th...Today, Linear Mixed Models (LMMs) are fitted, mostly, by assuming that random effects and errors have Gaussian distributions, therefore using Maximum Likelihood (ML) or REML estimation. However, for many data sets, that double assumption is unlikely to hold, particularly for the random effects, a crucial component </span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">in </span></span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">which assessment of magnitude is key in such modeling. Alternative fitting methods not relying on that assumption (as ANOVA ones and Rao</span></span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">’</span></span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">s MINQUE) apply, quite often, only to the very constrained class of variance components models. In this paper, a new computationally feasible estimation methodology is designed, first for the widely used class of 2-level (or longitudinal) LMMs with only assumption (beyond the usual basic ones) that residual errors are uncorrelated and homoscedastic, with no distributional assumption imposed on the random effects. A major asset of this new approach is that it yields nonnegative variance estimates and covariance matrices estimates which are symmetric and, at least, positive semi-definite. Furthermore, it is shown that when the LMM is, indeed, Gaussian, this new methodology differs from ML just through a slight variation in the denominator of the residual variance estimate. The new methodology actually generalizes to LMMs a well known nonparametric fitting procedure for standard Linear Models. Finally, the methodology is also extended to ANOVA LMMs, generalizing an old method by Henderson for ML estimation in such models under normality.展开更多
Stable water isotopes are natural tracers quantifying the contribution of moisture recycling to local precipitation,i.e.,the moisture recycling ratio,but various isotope-based models usually lead to different results,...Stable water isotopes are natural tracers quantifying the contribution of moisture recycling to local precipitation,i.e.,the moisture recycling ratio,but various isotope-based models usually lead to different results,which affects the accuracy of local moisture recycling.In this study,a total of 18 stations from four typical areas in China were selected to compare the performance of isotope-based linear and Bayesian mixing models and to determine local moisture recycling ratio.Among the three vapor sources including advection,transpiration,and surface evaporation,the advection vapor usually played a dominant role,and the contribution of surface evaporation was less than that of transpiration.When the abnormal values were ignored,the arithmetic averages of differences between isotope-based linear and the Bayesian mixing models were 0.9%for transpiration,0.2%for surface evaporation,and–1.1%for advection,respectively,and the medians were 0.5%,0.2%,and–0.8%,respectively.The importance of transpiration was slightly less for most cases when the Bayesian mixing model was applied,and the contribution of advection was relatively larger.The Bayesian mixing model was found to perform better in determining an efficient solution since linear model sometimes resulted in negative contribution ratios.Sensitivity test with two isotope scenarios indicated that the Bayesian model had a relatively low sensitivity to the changes in isotope input,and it was important to accurately estimate the isotopes in precipitation vapor.Generally,the Bayesian mixing model should be recommended instead of a linear model.The findings are useful for understanding the performance of isotope-based linear and Bayesian mixing models under various climate backgrounds.展开更多
Taking the nonlinear nature of runoff system into account,and combining auto-regression method and multi-regression method,a Nonlinear Mixed Regression Model (NMR) was established to analyze the impact of temperature ...Taking the nonlinear nature of runoff system into account,and combining auto-regression method and multi-regression method,a Nonlinear Mixed Regression Model (NMR) was established to analyze the impact of temperature and precipitation changes on annual river runoff process. The model was calibrated and verified by using BP neural network with observed meteorological and runoff data from Daiying Hydrological Station in the Chaohe River of Hebei Province in 1956–2000. Compared with auto-regression model,linear multi-regression model and linear mixed regression model,NMR can improve forecasting precision remarkably. Therefore,the simulation of climate change scenarios was carried out by NMR. The results show that the nonlinear mixed regression model can simulate annual river runoff well.展开更多
Impacts of the minimum purchase price policy for grain on the planting area of rice in Hubei Province were analyzed based on a mixed linear model.After the indicator system containing the minimum purchase price policy...Impacts of the minimum purchase price policy for grain on the planting area of rice in Hubei Province were analyzed based on a mixed linear model.After the indicator system containing the minimum purchase price policy and other factors influencing the planting area of rice was constructed,principal component analysis of the system was conducted,and then a mixed linear model where the planting area of rice was as the dependent variable was established.The results show that after the exclusion of the interference from other factors,the minimum purchase price policy for grain had a positive impact on the planting area of rice in Hubei Province.That is,the minimum purchase price policy significantly stimulated the growth of rice planting area in Hubei Province.展开更多
A linear mixed model is used to determine the explaining infant mortality rate data of United Nations countries. The HDI （human development index） has a significant negative linear relationship with infant mortality...A linear mixed model is used to determine the explaining infant mortality rate data of United Nations countries. The HDI （human development index） has a significant negative linear relationship with infant mortality rate. United Nations data shows that the infant mortality rate has a descending trend over the period 1990-2010. This study aims to assess the value of the HDI as a predictor of infant mortality rate. Findings in the paper suggest that significant percentage reductions in infant mortality might be possible for countries for controlling the HDI.展开更多
We focus on the development of model selection criteria in linear mixed models. In particular, we propose the model selection criteria following the Mallows’ Conceptual Predictive Statistic (Cp) [1] [2] in linear mix...We focus on the development of model selection criteria in linear mixed models. In particular, we propose the model selection criteria following the Mallows’ Conceptual Predictive Statistic (Cp) [1] [2] in linear mixed models. When correlation exists between the observations in data, the normal Gauss discrepancy in univariate case is not appropriate to measure the distance between the true model and a candidate model. Instead, we define a marginal Gauss discrepancy which takes the correlation into account in the mixed models. The model selection criterion, marginal Cp, called MCp, serves as an asymptotically unbiased estimator of the expected marginal Gauss discrepancy. An improvement of MCp, called IMCp, is then derived and proved to be a more accurate estimator of the expected marginal Gauss discrepancy than MCp. The performance of the proposed criteria is investigated in a simulation study. The simulation results show that in small samples, the proposed criteria outperform the Akaike Information Criteria (AIC) [3] [4] and Bayesian Information Criterion (BIC) [5] in selecting the correct model;in large samples, their performance is competitive. Further, the proposed criteria perform significantly better for highly correlated response data than for weakly correlated data.展开更多
For the linear mixed model with skew-normal random effects, this paper gives the density function, moment generating function and independence conditions. The noncentral skew chi-square distribution is defined and its...For the linear mixed model with skew-normal random effects, this paper gives the density function, moment generating function and independence conditions. The noncentral skew chi-square distribution is defined and its density function is shown. The necessary and sufficient conditions under which a quadratic form is distributed as noncentral skew chi-square distribution are obtained. Also, a version of Cochran＇s theorem is given~ which modifies the result of Wang et al. （2009） and is used to set up exact tests for fixed effects and variance components of the proposed model. For illustration, our main results are applied to a real data problem.展开更多
Over the past decades,the expansion of natu-ral secondary forests has played a crucial role in offsetting the loss of primary forests and combating climate change.Despite this,there is a gap in our understanding of ho...Over the past decades,the expansion of natu-ral secondary forests has played a crucial role in offsetting the loss of primary forests and combating climate change.Despite this,there is a gap in our understanding of how tree species’growth and mortality patterns vary with eleva-tion in these secondary forests.In this study,we analyzed data from two censuses(spanning a five-year interval)conducted in both evergreen broadleaved forests(EBF)and temperate coniferous forests(TCF),which have been recovering for half a century,across elevation gradients in a subtropical mountain region,Mount Wuyi,China.The results indicated that the relative growth rate(RGR)of EBF(0.028±0.001 cm·cm^(-1)·a^(-1))and the mortality rate(MR)(20.03%±1.70%)were 27.3%and 16.4%higher,respec-tively,than those of TCF.Interestingly,the trade-off between RGR and MR in EBF weakened as elevation increased,a trend not observed in TCF.Conversely,TCF consistently showed a stronger trade-off between RGR and MR compared to EBF.Generalized linear mixed models revealed that ele-vation influences RGR both directly and indirectly through its interactions with slope,crown competition index(CCI),and tree canopy height(CH).However,tree mortality did not show a significant correlation with elevation.Additionally,DBH significantly influenced both tree growth and mortal-ity,whereas and CH and CCI had opposite effects on tree growth between EBF and TCF.Our study underscores the importance of elevation in shaping the population dynamics and the biomass carbon sink balance of mountain forests.These insights enhance our understanding of tree species’life strategies,enabling more accurate predictions of forest dynamics and their response to environmental changes.展开更多
The rapid evolution of unmanned aerial vehicle(UAV)technology and autonomous capabilities has positioned UAV as promising last-mile delivery means.Vehicle and onboard UAV collaborative delivery is introduced as a nove...The rapid evolution of unmanned aerial vehicle(UAV)technology and autonomous capabilities has positioned UAV as promising last-mile delivery means.Vehicle and onboard UAV collaborative delivery is introduced as a novel delivery mode.Spatiotemporal collaboration,along with energy consumption with payload and wind conditions play important roles in delivery route planning.This paper introduces the traveling salesman problem with time window and onboard UAV(TSPTWOUAV)and emphasizes the consideration of real-world scenarios,focusing on time collaboration and energy consumption with wind and payload.To address this,a mixed integer linear programming(MILP)model is formulated to minimize the energy consumption costs of vehicle and UAV.Furthermore,an adaptive large neighborhood search(ALNS)algorithm is applied to identify high-quality solutions efficiently.The effectiveness of the proposed model and algorithm is validated through numerical tests on real geographic instances and sensitivity analysis of key parameters is conducted.展开更多
In this paper, the problem of estimating the covariance matrix in general linear mixed models is considered. A new class of estimators is proposed. It is shown that this new estimator dominates the analysis of varianc...In this paper, the problem of estimating the covariance matrix in general linear mixed models is considered. A new class of estimators is proposed. It is shown that this new estimator dominates the analysis of variance estimate under two squared loss functions. Finally, some simulation results to compare the performance of the proposed estimator with that of the analysis of variance estimate are reported. The simulation results indicate that this new estimator provides a substantial improvement in risk under most situations.展开更多
In this paper, we propose a bias-corrected empirical likelihood （BCEL） ratio to construct a goodness- of-fit test for generalized linear mixed models. BCEL test maintains the advantage of empirical likelihood that i...In this paper, we propose a bias-corrected empirical likelihood （BCEL） ratio to construct a goodness- of-fit test for generalized linear mixed models. BCEL test maintains the advantage of empirical likelihood that is self scale invariant and then does not involve estimating limiting variance of the test statistic to avoid deteri- orating power of test. Furthermore, the bias correction makes the limit to be a process in which every variable is standard chi-squared. This simple structure of the process enables us to construct a Monte Carlo test proce- dure to approximate the null distribution. Thus, it overcomes a problem we encounter when classical empirical likelihood test is used, as it is asymptotically a functional of Gaussian process plus a normal shift function. The complicated covariance function makes it difficult to employ any approximation for the null distribution. The test is omnibus and power study shows that the test can detect local alternatives approaching the null at parametric rate. Simulations are carried out for illustration and for a comparison with existing method.展开更多
The genus Beta encompasses economically important root crops such as sugar and table beet.A Beta diversity set including the wild relative B.vulgaris ssp.maritima was grown in the field,and a large phenotypic diversit...The genus Beta encompasses economically important root crops such as sugar and table beet.A Beta diversity set including the wild relative B.vulgaris ssp.maritima was grown in the field,and a large phenotypic diversity was observed.The genomes of 290 accessions were sequenced,and more than 10 million high-quality SNPs were employed to study genetic diversity.A genome-wide association study was performed,and marker-trait associations were found for nine phenotypic traits.The candidate gene within the M locus controlling monogermity on chromosome 4 was previously unknown.The most significant association for monogermity was identified at the end of chromosome 4.Within this region,a non-synonymous mutation within the zinc-finger domain of the WIP2 gene co-segregated with monogermity.This gene plays a regulatory role in AGL8/FUL in Arabidopsis.Intriguingly,commercial hybrids are in a heterozygous state at this position.Thus,the long-sought gene for monogermity was identified in this study.Red and yellow pigmentation due to betalain accumulation in shoots and roots is an important characteristic of table and leaf beets.The strongest associations were found upstream or downstream of two genes encoding Cytochrome P450 and anthocyanin MYB-like transcription factor proteins involved in betalain biosynthesis.Significant associations for Cercospora leaf spot resistance were identified on chromosomes 1,2,7,and 9.The associated regions harbor genes encoding proteins with leucinerich repeats and nucleotide binding sites whose homologs are major constituents of plant-pathogen defense.展开更多
This study investigated the factors contributing to intravenous admixture preparation errors(IAPEs)within Pharmacy Intravenous Admixture Services(PIVAS).A retrospective analysis was conducted on IAPEs documented in th...This study investigated the factors contributing to intravenous admixture preparation errors(IAPEs)within Pharmacy Intravenous Admixture Services(PIVAS).A retrospective analysis was conducted on IAPEs documented in the PIVAS unit of a large multi-specialty hospital in China,which houses over 2000 beds,covering the period from January 1,2015 to December 31,2022.Drug preparation records were examined using a generalized linear mixed model(GLMM)to identify both univariate and multivariate factors associated with IAPE occurrences.A total of 824 IAPE cases were recorded during the study period,yielding an overall error rate of 0.018%.Univariate analysis identified drug categories(general drugs,anti-infective drugs,and antineoplastic drugs),preparation time(workdays),and years of work experience as significant determinants(P<0.05).Multivariate analysis further confirmed that drug categories(general and antineoplastic drugs),preparation time(workdays),and work experience remained statistically significant predictors of IAPE incidence(P<0.05).IAPEs in PIVAS were influenced by multiple factors,predominantly those related to personnel and drug characteristics.Targeted interventions,informed by multivariate analysis,are essential to mitigating these errors and enhancing medication safety.展开更多
The main purpose of this paper is to investigate D-optimal population designs in multi-response linear mixed models for longitudinal data.Observations of each response variable within subjects are assumed to have a fi...The main purpose of this paper is to investigate D-optimal population designs in multi-response linear mixed models for longitudinal data.Observations of each response variable within subjects are assumed to have a first-order autoregressive structure,possibly with observation error.The equivalence theorems are provided to characterise theD-optimal population designs for the estimation of fixed effects in the model.The semi-Bayesian D-optimal design which is robust against the serial correlation coefficient is also considered.Simulation studies show that the correlation between multi-response variables has tiny effects on the optimal design,while the experimental costs are important factors in the optimal designs.展开更多
Territory risk analysis has played an important role in the decision-making of auto insurance rate regulation.Due to the optimality of insurance loss data groupings,clustering methods become the natural choice for suc...Territory risk analysis has played an important role in the decision-making of auto insurance rate regulation.Due to the optimality of insurance loss data groupings,clustering methods become the natural choice for such territory risk classification.In this work,spatially constrained clustering is first applied to insurance loss data to form rating territories.The generalized linear model(GLM)and generalized linear mixed model(GLMM)are then proposed to derive the risk relativities of obtained clusters.Each basic rating unit within the same cluster,namely Forward Sortation Area(FSA),takes the same risk relativity value as its cluster.The obtained risk relativities from GLM or GLMM are used to calculate the performance metrics,including RMSE,MAD,and Gini coefficients.The spatially constrained clustering and the risk relativity estimate help obtain a set of territory risk benchmarks used in rate filings to guide the rate regulation process.展开更多
The health hazard of mercury （Hg） compounds is internationally recognized, and the main pathways for methylmercury （MeHg） intake in humans are through consumption of food, especially fish. Given the large releases...The health hazard of mercury （Hg） compounds is internationally recognized, and the main pathways for methylmercury （MeHg） intake in humans are through consumption of food, especially fish. Given the large releases of Hg to the environment in China, combined with the fast development of hydropower, this issue deserves attention. Provided similar mobilization pathways of Hg in China as seen in reservoirs in North America and Europe one should expect increased Hg contamination in relation to future hydropower reservoir construction in this country. This study presents total Hg （THg） concentrations in wild fish from six Guizhou reservoirs, China. The THg concentrations in fish were generally low despite high background levels in the bedrock and depositions from local point sources. The over all mean ± SD concentration of THg was （0.066 ± 0.078） μg/g （n = 235）. After adjusting for among-reservoir variation in THg, there were significant differences in THg among functional groups of the fish, assumed to re？ect trophic levels. Predicted THg- concentration ratios, retrieved from a mixed linear model, between the functional groups were 9：4：4：1 for carnivorous, omnivorous, planktivorous and herbivorous fish. This result indicated that MeHg accumulation may prevail even under circumstances with short food chains as in this Chinese water system. No fish exceeded recommended maximum THg limit for human consumption set by World Health Organization and the Standardization Administration of China （0.5 μg/g fish wet weight （ww））. Only six fish （2.5%） exceeded the maximum THg limit set by US Environmental Protection Agency （0.3 μg/g fish ww）.展开更多
Eleven evaluating parameters for rice core collection were assessed based on genotypic values and molecular marke＇ information. Monte Carlo simulation combined with mixed linear model was used to eliminate the interf...Eleven evaluating parameters for rice core collection were assessed based on genotypic values and molecular marke＇ information. Monte Carlo simulation combined with mixed linear model was used to eliminate the interference from environment in order to draw more reliable results. The coincidence rate of range （CR） was the optimal parameter. Mean Simpson index （MD）, mean Shannon-Weaver index of genetic diversity （M1） and mean polymorphism information content （MPIC） were important evaluating parameters. The variable rate of coefficient of variation （VR） could act as an important reference parameter for evaluating the variation degree of core collection. Percentage of polymorphic loci （p） could be used as a determination parameter for the size of core collection. Mean difference percentage （MD） was a determination parameter for the reliability judgment of core collection. The effective evaluating parameters for core collection selected in the research could be used as criteria for sampling percentage in different plant germplasm populations.展开更多
Dissecting the genetic architecture of complex traits is an ongoing challenge for geneticists.Two complementary approaches for genetic mapping,linkage mapping and association mapping have led to successful dissection ...Dissecting the genetic architecture of complex traits is an ongoing challenge for geneticists.Two complementary approaches for genetic mapping,linkage mapping and association mapping have led to successful dissection of complex traits in many crop species.Both of these methods detect quantitative trait loci(QTL) by identifying marker–trait associations,and the only fundamental difference between them is that between mapping populations,which directly determine mapping resolution and power.Based on this difference,we first summarize in this review the advances and limitations of family-based mapping and natural population-based mapping instead of linkage mapping and association mapping.We then describe statistical methods used for improving detection power and computational speed and outline emerging areas such as large-scale meta-analysis for genetic mapping in crops.In the era of next-generation sequencing,there has arisen an urgent need for proper population design,advanced statistical strategies,and precision phenotyping to fully exploit high-throughput genotyping.展开更多
One hundred and sixty-eight genotypes of cotton from the same growing region were used as a germplasm group to study the validity of different genetic distances in constructing cotton core subset. Mixed linear model a...One hundred and sixty-eight genotypes of cotton from the same growing region were used as a germplasm group to study the validity of different genetic distances in constructing cotton core subset. Mixed linear model approach was employed to unbiasedly predict genotypic values of 20 traits for eliminating the environmental effect. Six commonly used genetic distances(Euclidean,standardized Euclidean,Mahalanobis,city block,cosine and correlation distances) combining four commonly used hierarchical cluster methods(single distance,complete distance,unweighted pair-group average and Ward's methods) were used in the least distance stepwise sampling(LDSS) method for constructing different core subsets. The analyses of variance(ANOVA) of different evaluating parameters showed that the validities of cosine and correlation distances were inferior to those of Euclidean,standardized Euclidean,Mahalanobis and city block distances. Standardized Euclidean distance was slightly more effective than Euclidean,Mahalanobis and city block distances. The principal analysis validated standardized Euclidean distance in the course of constructing practical core subsets. The covariance matrix of accessions might be ill-conditioned when Mahalanobis distance was used to calculate genetic distance at low sampling percentages,which led to bias in small-sized core subset construction. The standardized Euclidean distance is recommended in core subset construction with LDSS method.展开更多
基金supported by the Funding Project for Academic Human Resources Development in Institutions of Higher Learning Under the Jurisdiction of Beijing Municipality (0506011200702)National Natural Science Foundation of China+2 种基金Tian Yuan Special Foundation (10926059)Foundation of Zhejiang Educational Committee (Y200803920)Scientific Research Foundation of Hangzhou Dianzi University(KYS025608094)
文摘In this article, the problem of estimating the covariance matrix in general linear mixed models is considered. Two new classes of estimators obtained by shrinking the eigenvalues towards the origin and the arithmetic mean, respectively, are proposed. It is shown that these new estimators dominate the unbiased estimator under the squared error loss function. Finally, some simulation results to compare the performance of the proposed estimators with that of the unbiased estimator are reported. The simulation results indicate that these new shrinkage estimators provide a substantial improvement in risk under most situations.
文摘Today, Linear Mixed Models (LMMs) are fitted, mostly, by assuming that random effects and errors have Gaussian distributions, therefore using Maximum Likelihood (ML) or REML estimation. However, for many data sets, that double assumption is unlikely to hold, particularly for the random effects, a crucial component </span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">in </span></span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">which assessment of magnitude is key in such modeling. Alternative fitting methods not relying on that assumption (as ANOVA ones and Rao</span></span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">’</span></span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">s MINQUE) apply, quite often, only to the very constrained class of variance components models. In this paper, a new computationally feasible estimation methodology is designed, first for the widely used class of 2-level (or longitudinal) LMMs with only assumption (beyond the usual basic ones) that residual errors are uncorrelated and homoscedastic, with no distributional assumption imposed on the random effects. A major asset of this new approach is that it yields nonnegative variance estimates and covariance matrices estimates which are symmetric and, at least, positive semi-definite. Furthermore, it is shown that when the LMM is, indeed, Gaussian, this new methodology differs from ML just through a slight variation in the denominator of the residual variance estimate. The new methodology actually generalizes to LMMs a well known nonparametric fitting procedure for standard Linear Models. Finally, the methodology is also extended to ANOVA LMMs, generalizing an old method by Henderson for ML estimation in such models under normality.
基金This study was supported by the National Natural Science Foundation of China(42261008,41971034)the Natural Science Foundation of Gansu Province,China(22JR5RA074).
文摘Stable water isotopes are natural tracers quantifying the contribution of moisture recycling to local precipitation,i.e.,the moisture recycling ratio,but various isotope-based models usually lead to different results,which affects the accuracy of local moisture recycling.In this study,a total of 18 stations from four typical areas in China were selected to compare the performance of isotope-based linear and Bayesian mixing models and to determine local moisture recycling ratio.Among the three vapor sources including advection,transpiration,and surface evaporation,the advection vapor usually played a dominant role,and the contribution of surface evaporation was less than that of transpiration.When the abnormal values were ignored,the arithmetic averages of differences between isotope-based linear and the Bayesian mixing models were 0.9%for transpiration,0.2%for surface evaporation,and–1.1%for advection,respectively,and the medians were 0.5%,0.2%,and–0.8%,respectively.The importance of transpiration was slightly less for most cases when the Bayesian mixing model was applied,and the contribution of advection was relatively larger.The Bayesian mixing model was found to perform better in determining an efficient solution since linear model sometimes resulted in negative contribution ratios.Sensitivity test with two isotope scenarios indicated that the Bayesian model had a relatively low sensitivity to the changes in isotope input,and it was important to accurately estimate the isotopes in precipitation vapor.Generally,the Bayesian mixing model should be recommended instead of a linear model.The findings are useful for understanding the performance of isotope-based linear and Bayesian mixing models under various climate backgrounds.
基金Under the auspices of National Natural Science Foundation of China (No. 50809004)
文摘Taking the nonlinear nature of runoff system into account,and combining auto-regression method and multi-regression method,a Nonlinear Mixed Regression Model (NMR) was established to analyze the impact of temperature and precipitation changes on annual river runoff process. The model was calibrated and verified by using BP neural network with observed meteorological and runoff data from Daiying Hydrological Station in the Chaohe River of Hebei Province in 1956–2000. Compared with auto-regression model,linear multi-regression model and linear mixed regression model,NMR can improve forecasting precision remarkably. Therefore,the simulation of climate change scenarios was carried out by NMR. The results show that the nonlinear mixed regression model can simulate annual river runoff well.
基金Supported by the Humanities and Social Sciences Foundation for Young Scholars of Ministry of Education of China(11y3jc630197)
文摘Impacts of the minimum purchase price policy for grain on the planting area of rice in Hubei Province were analyzed based on a mixed linear model.After the indicator system containing the minimum purchase price policy and other factors influencing the planting area of rice was constructed,principal component analysis of the system was conducted,and then a mixed linear model where the planting area of rice was as the dependent variable was established.The results show that after the exclusion of the interference from other factors,the minimum purchase price policy for grain had a positive impact on the planting area of rice in Hubei Province.That is,the minimum purchase price policy significantly stimulated the growth of rice planting area in Hubei Province.
文摘A linear mixed model is used to determine the explaining infant mortality rate data of United Nations countries. The HDI （human development index） has a significant negative linear relationship with infant mortality rate. United Nations data shows that the infant mortality rate has a descending trend over the period 1990-2010. This study aims to assess the value of the HDI as a predictor of infant mortality rate. Findings in the paper suggest that significant percentage reductions in infant mortality might be possible for countries for controlling the HDI.
文摘We focus on the development of model selection criteria in linear mixed models. In particular, we propose the model selection criteria following the Mallows’ Conceptual Predictive Statistic (Cp) [1] [2] in linear mixed models. When correlation exists between the observations in data, the normal Gauss discrepancy in univariate case is not appropriate to measure the distance between the true model and a candidate model. Instead, we define a marginal Gauss discrepancy which takes the correlation into account in the mixed models. The model selection criterion, marginal Cp, called MCp, serves as an asymptotically unbiased estimator of the expected marginal Gauss discrepancy. An improvement of MCp, called IMCp, is then derived and proved to be a more accurate estimator of the expected marginal Gauss discrepancy than MCp. The performance of the proposed criteria is investigated in a simulation study. The simulation results show that in small samples, the proposed criteria outperform the Akaike Information Criteria (AIC) [3] [4] and Bayesian Information Criterion (BIC) [5] in selecting the correct model;in large samples, their performance is competitive. Further, the proposed criteria perform significantly better for highly correlated response data than for weakly correlated data.
基金supported by National Natural Science Foundation of China(Grant No.11401148)Ministry of Education of China,Humanities and Social Science Projects(Grant Nos.14YJC910005,10YJC790184)+2 种基金Zhejiang Provincial Natural Science Foundation of China(Grant No.LY14A010030)Zhejiang Provincial Philosophy and Social Science Planning Project of China(Grant No.13NDJC089YB)Houji Scholar Fund of Northwest A and F University,China
文摘For the linear mixed model with skew-normal random effects, this paper gives the density function, moment generating function and independence conditions. The noncentral skew chi-square distribution is defined and its density function is shown. The necessary and sufficient conditions under which a quadratic form is distributed as noncentral skew chi-square distribution are obtained. Also, a version of Cochran＇s theorem is given~ which modifies the result of Wang et al. （2009） and is used to set up exact tests for fixed effects and variance components of the proposed model. For illustration, our main results are applied to a real data problem.
基金funded by the National Natural Science Foundation of China(Grant No.32271872).
文摘Over the past decades,the expansion of natu-ral secondary forests has played a crucial role in offsetting the loss of primary forests and combating climate change.Despite this,there is a gap in our understanding of how tree species’growth and mortality patterns vary with eleva-tion in these secondary forests.In this study,we analyzed data from two censuses(spanning a five-year interval)conducted in both evergreen broadleaved forests(EBF)and temperate coniferous forests(TCF),which have been recovering for half a century,across elevation gradients in a subtropical mountain region,Mount Wuyi,China.The results indicated that the relative growth rate(RGR)of EBF(0.028±0.001 cm·cm^(-1)·a^(-1))and the mortality rate(MR)(20.03%±1.70%)were 27.3%and 16.4%higher,respec-tively,than those of TCF.Interestingly,the trade-off between RGR and MR in EBF weakened as elevation increased,a trend not observed in TCF.Conversely,TCF consistently showed a stronger trade-off between RGR and MR compared to EBF.Generalized linear mixed models revealed that ele-vation influences RGR both directly and indirectly through its interactions with slope,crown competition index(CCI),and tree canopy height(CH).However,tree mortality did not show a significant correlation with elevation.Additionally,DBH significantly influenced both tree growth and mortal-ity,whereas and CH and CCI had opposite effects on tree growth between EBF and TCF.Our study underscores the importance of elevation in shaping the population dynamics and the biomass carbon sink balance of mountain forests.These insights enhance our understanding of tree species’life strategies,enabling more accurate predictions of forest dynamics and their response to environmental changes.
基金Fundamental Research Funds for the Central Universities(2024JBZX038)National Natural Science F oundation of China(62076023)。
文摘The rapid evolution of unmanned aerial vehicle(UAV)technology and autonomous capabilities has positioned UAV as promising last-mile delivery means.Vehicle and onboard UAV collaborative delivery is introduced as a novel delivery mode.Spatiotemporal collaboration,along with energy consumption with payload and wind conditions play important roles in delivery route planning.This paper introduces the traveling salesman problem with time window and onboard UAV(TSPTWOUAV)and emphasizes the consideration of real-world scenarios,focusing on time collaboration and energy consumption with wind and payload.To address this,a mixed integer linear programming(MILP)model is formulated to minimize the energy consumption costs of vehicle and UAV.Furthermore,an adaptive large neighborhood search(ALNS)algorithm is applied to identify high-quality solutions efficiently.The effectiveness of the proposed model and algorithm is validated through numerical tests on real geographic instances and sensitivity analysis of key parameters is conducted.
基金This research is supported by National Natural Science Foundation of China, Tian Yuan Special Foundation under Grant No. 10926059 and Zhejiang Provincial Natural Science Foundation of China under Grant No. Y6100053.
文摘In this paper, the problem of estimating the covariance matrix in general linear mixed models is considered. A new class of estimators is proposed. It is shown that this new estimator dominates the analysis of variance estimate under two squared loss functions. Finally, some simulation results to compare the performance of the proposed estimator with that of the analysis of variance estimate are reported. The simulation results indicate that this new estimator provides a substantial improvement in risk under most situations.
基金Supported by the National Natural Science Foundation of China(No.10901109)a grant(HKBU2030/07P)from the Research Grants Council of Hong Kong,Hong Kong,China
文摘In this paper, we propose a bias-corrected empirical likelihood （BCEL） ratio to construct a goodness- of-fit test for generalized linear mixed models. BCEL test maintains the advantage of empirical likelihood that is self scale invariant and then does not involve estimating limiting variance of the test statistic to avoid deteri- orating power of test. Furthermore, the bias correction makes the limit to be a process in which every variable is standard chi-squared. This simple structure of the process enables us to construct a Monte Carlo test proce- dure to approximate the null distribution. Thus, it overcomes a problem we encounter when classical empirical likelihood test is used, as it is asymptotically a functional of Gaussian process plus a normal shift function. The complicated covariance function makes it difficult to employ any approximation for the null distribution. The test is omnibus and power study shows that the test can detect local alternatives approaching the null at parametric rate. Simulations are carried out for illustration and for a comparison with existing method.
基金the financial support from the German Research Foundation(Deutsche Forschungsgemeinschaft,DFG)—Project Number 400993799(Project 2 within the Research Training Group 2501 Translational Evolutionary Research,https://gepris.dfg.de/gepris/projekt/400993799)supported by the BMBF-funded de.NBI Cloud,part of the German Network for Bioinformatics(de.NBI)。
文摘The genus Beta encompasses economically important root crops such as sugar and table beet.A Beta diversity set including the wild relative B.vulgaris ssp.maritima was grown in the field,and a large phenotypic diversity was observed.The genomes of 290 accessions were sequenced,and more than 10 million high-quality SNPs were employed to study genetic diversity.A genome-wide association study was performed,and marker-trait associations were found for nine phenotypic traits.The candidate gene within the M locus controlling monogermity on chromosome 4 was previously unknown.The most significant association for monogermity was identified at the end of chromosome 4.Within this region,a non-synonymous mutation within the zinc-finger domain of the WIP2 gene co-segregated with monogermity.This gene plays a regulatory role in AGL8/FUL in Arabidopsis.Intriguingly,commercial hybrids are in a heterozygous state at this position.Thus,the long-sought gene for monogermity was identified in this study.Red and yellow pigmentation due to betalain accumulation in shoots and roots is an important characteristic of table and leaf beets.The strongest associations were found upstream or downstream of two genes encoding Cytochrome P450 and anthocyanin MYB-like transcription factor proteins involved in betalain biosynthesis.Significant associations for Cercospora leaf spot resistance were identified on chromosomes 1,2,7,and 9.The associated regions harbor genes encoding proteins with leucinerich repeats and nucleotide binding sites whose homologs are major constituents of plant-pathogen defense.
基金The National Natural Science Foundation of China(Grant No.72474013)the Beijing Health Technology Promotion Project(Grant No.BHTPP2024007)。
文摘This study investigated the factors contributing to intravenous admixture preparation errors(IAPEs)within Pharmacy Intravenous Admixture Services(PIVAS).A retrospective analysis was conducted on IAPEs documented in the PIVAS unit of a large multi-specialty hospital in China,which houses over 2000 beds,covering the period from January 1,2015 to December 31,2022.Drug preparation records were examined using a generalized linear mixed model(GLMM)to identify both univariate and multivariate factors associated with IAPE occurrences.A total of 824 IAPE cases were recorded during the study period,yielding an overall error rate of 0.018%.Univariate analysis identified drug categories(general drugs,anti-infective drugs,and antineoplastic drugs),preparation time(workdays),and years of work experience as significant determinants(P<0.05).Multivariate analysis further confirmed that drug categories(general and antineoplastic drugs),preparation time(workdays),and work experience remained statistically significant predictors of IAPE incidence(P<0.05).IAPEs in PIVAS were influenced by multiple factors,predominantly those related to personnel and drug characteristics.Targeted interventions,informed by multivariate analysis,are essential to mitigating these errors and enhancing medication safety.
基金partly supported by the National Natural Science Foundation of China(Nos.11971318,11871143)Shanghai Rising-Star Program(No.20QA1407500).
文摘The main purpose of this paper is to investigate D-optimal population designs in multi-response linear mixed models for longitudinal data.Observations of each response variable within subjects are assumed to have a first-order autoregressive structure,possibly with observation error.The equivalence theorems are provided to characterise theD-optimal population designs for the estimation of fixed effects in the model.The semi-Bayesian D-optimal design which is robust against the serial correlation coefficient is also considered.Simulation studies show that the correlation between multi-response variables has tiny effects on the optimal design,while the experimental costs are important factors in the optimal designs.
文摘Territory risk analysis has played an important role in the decision-making of auto insurance rate regulation.Due to the optimality of insurance loss data groupings,clustering methods become the natural choice for such territory risk classification.In this work,spatially constrained clustering is first applied to insurance loss data to form rating territories.The generalized linear model(GLM)and generalized linear mixed model(GLMM)are then proposed to derive the risk relativities of obtained clusters.Each basic rating unit within the same cluster,namely Forward Sortation Area(FSA),takes the same risk relativity value as its cluster.The obtained risk relativities from GLM or GLMM are used to calculate the performance metrics,including RMSE,MAD,and Gini coefficients.The spatially constrained clustering and the risk relativity estimate help obtain a set of territory risk benchmarks used in rate filings to guide the rate regulation process.
文摘The health hazard of mercury （Hg） compounds is internationally recognized, and the main pathways for methylmercury （MeHg） intake in humans are through consumption of food, especially fish. Given the large releases of Hg to the environment in China, combined with the fast development of hydropower, this issue deserves attention. Provided similar mobilization pathways of Hg in China as seen in reservoirs in North America and Europe one should expect increased Hg contamination in relation to future hydropower reservoir construction in this country. This study presents total Hg （THg） concentrations in wild fish from six Guizhou reservoirs, China. The THg concentrations in fish were generally low despite high background levels in the bedrock and depositions from local point sources. The over all mean ± SD concentration of THg was （0.066 ± 0.078） μg/g （n = 235）. After adjusting for among-reservoir variation in THg, there were significant differences in THg among functional groups of the fish, assumed to re？ect trophic levels. Predicted THg- concentration ratios, retrieved from a mixed linear model, between the functional groups were 9：4：4：1 for carnivorous, omnivorous, planktivorous and herbivorous fish. This result indicated that MeHg accumulation may prevail even under circumstances with short food chains as in this Chinese water system. No fish exceeded recommended maximum THg limit for human consumption set by World Health Organization and the Standardization Administration of China （0.5 μg/g fish wet weight （ww））. Only six fish （2.5%） exceeded the maximum THg limit set by US Environmental Protection Agency （0.3 μg/g fish ww）.
基金the National Natural Science Foundation of China （Grant No. 30270759） the Science and Technology Department of Zhejiang Province （Grant No. 2005C32001）.
文摘Eleven evaluating parameters for rice core collection were assessed based on genotypic values and molecular marke＇ information. Monte Carlo simulation combined with mixed linear model was used to eliminate the interference from environment in order to draw more reliable results. The coincidence rate of range （CR） was the optimal parameter. Mean Simpson index （MD）, mean Shannon-Weaver index of genetic diversity （M1） and mean polymorphism information content （MPIC） were important evaluating parameters. The variable rate of coefficient of variation （VR） could act as an important reference parameter for evaluating the variation degree of core collection. Percentage of polymorphic loci （p） could be used as a determination parameter for the size of core collection. Mean difference percentage （MD） was a determination parameter for the reliability judgment of core collection. The effective evaluating parameters for core collection selected in the research could be used as criteria for sampling percentage in different plant germplasm populations.
基金supported by the Priority Academic Program Development of Jiangsu Higher Education Institutionthe National Natural Science Foundation of China(Nos.91535103,31391632,and 31200943)+4 种基金the National High Technology Research and Development Program of China(No.2014AA10A601-5)the Natural Science Foundation of Jiangsu Province(No.BK2012261)the Natural Science Foundation of Jiangsu Higher Education Institution(No.14KJA210005)the Postgraduate Research and Innovation Project in Jiangsu Province(No.KYLX151368)the Innovative Research Team of University in Jiangsu Province
文摘Dissecting the genetic architecture of complex traits is an ongoing challenge for geneticists.Two complementary approaches for genetic mapping,linkage mapping and association mapping have led to successful dissection of complex traits in many crop species.Both of these methods detect quantitative trait loci(QTL) by identifying marker–trait associations,and the only fundamental difference between them is that between mapping populations,which directly determine mapping resolution and power.Based on this difference,we first summarize in this review the advances and limitations of family-based mapping and natural population-based mapping instead of linkage mapping and association mapping.We then describe statistical methods used for improving detection power and computational speed and outline emerging areas such as large-scale meta-analysis for genetic mapping in crops.In the era of next-generation sequencing,there has arisen an urgent need for proper population design,advanced statistical strategies,and precision phenotyping to fully exploit high-throughput genotyping.
基金Project supported by the National Natural Science Foundation of China (No. 30270759)the Cooperation Project in Science and Technology between China and Poland Governments (No. 32-38)the Scientific Research Foundation for Doctors in Shandong Academy of Agricultural Sciences (No. [2007]20), China
文摘One hundred and sixty-eight genotypes of cotton from the same growing region were used as a germplasm group to study the validity of different genetic distances in constructing cotton core subset. Mixed linear model approach was employed to unbiasedly predict genotypic values of 20 traits for eliminating the environmental effect. Six commonly used genetic distances(Euclidean,standardized Euclidean,Mahalanobis,city block,cosine and correlation distances) combining four commonly used hierarchical cluster methods(single distance,complete distance,unweighted pair-group average and Ward's methods) were used in the least distance stepwise sampling(LDSS) method for constructing different core subsets. The analyses of variance(ANOVA) of different evaluating parameters showed that the validities of cosine and correlation distances were inferior to those of Euclidean,standardized Euclidean,Mahalanobis and city block distances. Standardized Euclidean distance was slightly more effective than Euclidean,Mahalanobis and city block distances. The principal analysis validated standardized Euclidean distance in the course of constructing practical core subsets. The covariance matrix of accessions might be ill-conditioned when Mahalanobis distance was used to calculate genetic distance at low sampling percentages,which led to bias in small-sized core subset construction. The standardized Euclidean distance is recommended in core subset construction with LDSS method.
