Marginal Conceptual Predictive Statistic for Mixed Model Selection

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.

WOMBAT—A tool for mixed model analyses in quantitative genetics by restricted maximum likelihood (REML)

WOMBAT is a software package for quantitative genetic analyses of continuous traits, fitting a linear, mixed model; estimates of covariance components and the resulting genetic parameters are obtained by restricted maximum likelihood. A wide range of models, comprising numerous traits, multiple fixed and random effects, selected genetic covariance structures, random regression models and reduced rank estimation are accommodated. WOMBAT employs up-to-date numerical and computational methods. Together with the use of efficient compilers, this generates fast executable programs, suitable for large scale analyses. Use of WOMBAT is illustrated for a bivariate analysis. The package consists of the executable program, available for LINUX and WINDOWS environments, manual and a set of worked example, and can be downloaded free of charge from http：//agbu. une.edu.au/-kmeyer/wombat.html

Analysis of Mixed Model in Gear Transmission Based on ADAMS

The traditional method of mechanical gear driving simulation includes gear pair method and solid to solid contact method. The former has higher solving efficiency but lower results accuracy; the latter usually obtains higher precision of results while the calculation process is complex, also it is not easy to converge. Currently, most of the researches are focused on the description of geometric models and the definition of boundary conditions. However, none of them can solve the problems fundamentally. To improve the simulation efficiency while ensure the results with high accuracy, a mixed model method which uses gear tooth profiles to take the place of the solid gear to simulate gear movement is presented under these circumstances. In the process of modeling, build the solid models of the mechanism in the SolidWorks firstly; Then collect the point coordinates of outline curves of the gear using SolidWorks API and create fit curves in Adams based on the point coordinates; Next, adjust the position of those fitting curves according to the position of the contact area; Finally, define the loading conditions, boundary conditions and simulation parameters. The method provides gear shape information by tooth profile curves; simulates the mesh process through tooth profile curve to curve contact and offer mass as well as inertia data via solid gear models. This simulation process combines the two models to complete the gear driving analysis. In order to verify the validity of the method presented, both theoretical derivation and numerical simulation on a runaway escapement are conducted. The results show that the computational efficiency of the mixed model method is 1.4 times over the traditional method which contains solid to solid contact. Meanwhile, the simulation results are more closely to theoretical calculations. Consequently, mixed model method has a high application value regarding to the study of the dynamics of gear mechanism.

IMPROVED ESTIMATES OF THE COVARIANCE MATRIX IN GENERAL LINEAR MIXED MODELS

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.

Software applications for providing comprehensive computing capabilities to problems related to mixed models in animal breeding

Several computer packages have been developed to accomplish improved programs for animal breeding and genetic selection. This paper described most of the currant software and provided suggestions for improved software. Khon Kaen University, Thailand, will provide free of charge the new software developed at Khon Kaen University by the author of this paper. The contact for requesting the software is listed: monchai@kku.ac.th.

A Mixed Model Analysis of a Fertilizer Experiment on Oil Palm in Nigeria

Mixed model analysis procedure was used to analyze the effect of fertilizer application on the Fresh Fruit Bunch (FFB) yield of oil palm. This was with a view to achieve the most appropriate and a robust model for analyzing yield response for fertilizer application in oil palm. In this study, a mixed model analysis procedure was used to analyze yield data obtained from a fertilizer trial conducted between 1997 and 2005. In mixed effect model, replicates and years were used as block. In contrast the fixed effect ANOVA model usually lumped up replicates and years as a random error. In the model replicates were used as block with no block interaction, replicates as block with allowance for block-fertilizer interaction, years as block with allowance for block-fertilizer interaction, and years and replicates as block with allowance for year fertilizer and replicate-fertilizer interaction. Mixed model theory was also used to provide the explicit description of the design matrices in the models. Also, hypotheses relevant to each model were formulated and used to test for specific effects in the models such as, fixed part, random part and interacting parts using appropriate error terms as determined by the derived Expected Mean Squares (EMS). The results revealed that at 5% significant level (p

Nonparametric Estimation in Linear Mixed Models with Uncorrelated Homoscedastic Errors

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 in which assessment of magnitude is key in such modeling. Alternative fitting methods not relying on that assumption (as ANOVA ones and Rao’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.

Linear Mixed Model Analysis of Worldwide Longitudinal Infant Mortality Rate Data and Association with Human Development Index

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.

Combined Selection in Backcross Population of Papaya(Carica papaya L.)by the Mixed Model Methodology

The selection of superior genotypes based on the simultaneous response to several characteristics of agronomic importance is a key strategy to overcome the scarcity of available varieties of papaya. This study aimed to apply the combined selection by using distinct selection indexes based on both the genetic values obtained by the REML/BLUP methodology and the real measured values to select agronomically superior genotypes of papaya within backcross progenies. The combined selection was carried out based on genetic and phenotypic values, original and standardized, multiplied by the agronomic weights. The results of the analysis of genetic parameters indicate that the evaluated progenies have expressive genetic variability for the considered traits, and that there are real possibilities of genetic progress with the selection. Among the analyzed indexes, the one based on standardized genetic value presented greater consistency in the ranking of genetic material, demonstrating the advantage of data standardization. Five progenies belonging to the BC1 generation, and five to the BC3 generation were selected using this index. A total of 27 plants ag-ronomically superior were selected within the top five progenies and recommended for generation advance, 23 being selected by combined selection and 4 using the direct selection for the four mainly characters in papaya breeding program: production, pulp and fruit firmness and soluble solids. Beyond the selection of superior genotypes for the development of future inbred lines, this study also allowed defining the best strategy to apply the combined selection in papaya using pre-dicted breeding values obtained by BLUP. This strategy may allow higher accuracy in the selection process, thus increasing the chances of success of the breeding programs.

Solving Level Scheduling in Mixed Model Assembly Line by Simulated Annealing Method

This paper presents an application of the simulated annealing algorithm to solve level schedules in mixed model assembly line. Solving production sequences with both number of setups and material usage rates to the minimum rate will optimize the level schedule. Miltenburg algorithm (1989) is first used to get seed sequence to optimize further. For this the utility time of the line and setup time requirement on each station is considered. This seed sequence is optimized by simulated annealing. This investigation helps to understand the importance of utility in the assembly line. Up to 15 product sequences are taken and constructed by using randomizing method and find the objective function value for this. For a sequence optimization, a meta-heuristic seems much more promising to guide the search into feasible regions of the solution space. Simulated annealing is a stochastic local search meta-heuristic, which bases the acceptance of a modified neighboring solution on a probabilistic scheme inspired by thermal processes for obtaining low-energy states in heat baths. Experimental results show that the simulated annealing approach is favorable and competitive compared to Miltenburg’s constructive algorithm for the problems set considered. It is proposed to found 16,985 solutions, the time taken for computation is 23.47 to 130.35, and the simulated annealing improves 49.33% than Miltenberg.

Classification of territory risk by generalized linear and generalized linear mixed models

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 relationship between compartment models and their stochastic counterparts:A comparative study with examples of the COVID-19 epidemic modeling

Deterministic compartment models(CMs)and stochastic models,including stochastic CMs and agent-based models,are widely utilized in epidemic modeling.However,the relationship between CMs and their corresponding stochastic models is not well understood.The present study aimed to address this gap by conducting a comparative study using the susceptible,exposed,infectious,and recovered(SEIR)model and its extended CMs from the coronavirus disease 2019 modeling literature.We demonstrated the equivalence of the numerical solution of CMs using the Euler scheme and their stochastic counterparts through theoretical analysis and simulations.Based on this equivalence,we proposed an efficient model calibration method that could replicate the exact solution of CMs in the corresponding stochastic models through parameter adjustment.The advancement in calibration techniques enhanced the accuracy of stochastic modeling in capturing the dynamics of epidemics.However,it should be noted that discrete-time stochastic models cannot perfectly reproduce the exact solution of continuous-time CMs.Additionally,we proposed a new stochastic compartment and agent mixed model as an alternative to agent-based models for large-scale population simulations with a limited number of agents.This model offered a balance between computational efficiency and accuracy.The results of this research contributed to the comparison and unification of deterministic CMs and stochastic models in epidemic modeling.Furthermore,the results had implications for the development of hybrid models that integrated the strengths of both frameworks.Overall,the present study has provided valuable epidemic modeling techniques and their practical applications for understanding and controlling the spread of infectious diseases.

Adaptive Random Effects/Coefficients Modeling

Adaptive fractional polynomial modeling of general correlated outcomes is formulated to address nonlinearity in means, variances/dispersions, and correlations. Means and variances/dispersions are modeled using generalized linear models in fixed effects/coefficients. Correlations are modeled using random effects/coefficients. Nonlinearity is addressed using power transforms of primary (untransformed) predictors. Parameter estimation is based on extended linear mixed modeling generalizing both generalized estimating equations and linear mixed modeling. Models are evaluated using likelihood cross-validation (LCV) scores and are generated adaptively using a heuristic search controlled by LCV scores. Cases covered include linear, Poisson, logistic, exponential, and discrete regression of correlated continuous, count/rate, dichotomous, positive continuous, and discrete numeric outcomes treated as normally, Poisson, Bernoulli, exponentially, and discrete numerically distributed, respectively. Example analyses are also generated for these five cases to compare adaptive random effects/coefficients modeling of correlated outcomes to previously developed adaptive modeling based on directly specified covariance structures. Adaptive random effects/coefficients modeling substantially outperforms direct covariance modeling in the linear, exponential, and discrete regression example analyses. It generates equivalent results in the logistic regression example analyses and it is substantially outperformed in the Poisson regression case. Random effects/coefficients modeling of correlated outcomes can provide substantial improvements in model selection compared to directly specified covariance modeling. However, directly specified covariance modeling can generate competitive or substantially better results in some cases while usually requiring less computation time.

Generalized Height-Diameter Models for Pinus montezumae Lamb. and Pinus pseudostrobus Lindl. Plantations in Michoacan, Mexico

Tree height (H) in a natural stand or forest plantation is a fundamental variable in management, and the use of mathematical expressions that estimate H as a function of diameter at breast height (d) or variables at the stand level is a valuable support tool in forest inventories. The objective was to fit and propose a generalized H-d model for Pinus montezumae and Pinus pseudostrobus established in forest plantations of Nuevo San Juan Parangaricutiro, Michoacan, Mexico. Using nonlinear least squares (NLS), 10 generalized H-d models were fitted to 883 and 1226 pairs of H-d data from Pinus montezumae and Pinus pseudostrobus, respectively. The best model was refitted with the maximum likelihood mixed effects model (MEM) approach by including the site as a classification variable and a known variance structure. The Wang and Tang equation was selected as the best model with NLS;the MEM with an additive effect on two of its parameters and an exponential variance function improved the fit statistics for Pinus montezumae and Pinus pseudostrobus, respectively. The model validation showed equality of means among the estimates for both species and an independent subsample. The calibration of the MEM at the plot level was efficient and might increase the applicability of these results. The inclusion of dominant height in the MEM approach helped to reduce bias in the estimates and also to better explain the variability among plots.

Simulation analysis of wear characteristics of connecting rod bearing bush based on improved mixed lubrication model

The failure rate of crankpin bearing bush of diesel engine under complex working conditions such as high temperature,dynamic load and variable speed is high.After serious wear,it is easy to deteriorate the stress state of connecting rod body and connecting rod bolt,resulting in serious accidents such as connecting rod fracture and body damage.Based on the mixed lubrication characteristics of connecting rod big endbearing shell of diesel engine under high explosion pressure impact load,an improved mixed lubrication mechanism model is established,which considers the influence of viscoelastic micro deformation of bearing bush material,integrates the full film lubrication model and dry friction model,couples dynamic equation of connecting rod.Then the actual lubrication state of big end bearing shell is simulated numerically.Further,the correctness of the theoretical research results is verified by fault simulation experiments.The results show that the high-frequency impact signal with fixed angle domain characteristics will be generated after the serious wear of bearing bush and the deterioration of lubrication state.The fault feature capture and alarm can be realized through the condition monitoring system,which can be applied to the fault monitoring of connecting rod bearing bush of diesel engine in the future.

Research on Mixed Logic Dynamic Modeling and Finite Control Set Model Predictive Control of Multi-Inverter Parallel System

Parallel connection of multiple inverters is an important means to solve the expansion,reserve and protection of distributed power generation,such as photovoltaics.In view of the shortcomings of traditional droop control methods such as weak anti-interference ability,low tracking accuracy of inverter output voltage and serious circulation phenomenon,a finite control set model predictive control(FCS-MPC)strategy of microgrid multiinverter parallel system based on Mixed Logical Dynamical(MLD)modeling is proposed.Firstly,the MLD modeling method is introduced logical variables,combining discrete events and continuous events to form an overall differential equation,which makes the modeling more accurate.Then a predictive controller is designed based on the model,and constraints are added to the objective function,which can not only solve the real-time changes of the control system by online optimization,but also effectively obtain a higher tracking accuracy of the inverter output voltage and lower total harmonic distortion rate(Total Harmonics Distortion,THD);and suppress the circulating current between the inverters,to obtain a good dynamic response.Finally,the simulation is carried out onMATLAB/Simulink to verify the correctness of the model and the rationality of the proposed strategy.This paper aims to provide guidance for the design and optimal control of multi-inverter parallel systems.

Comparison of isotope-based linear and Bayesian mixing models in determining moisture recycling ratio

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.

A unified Minorization-Maximization approach for estimation of general mixture models

The mixed distribution model is often used to extract information from heteroge-neous data and perform modeling analysis.When the density function of mixed distribution is complicated or the variable dimension is high,it usually brings challenges to the parameter es-timation of the mixed distribution model.The application of MM algorithm can avoid complex expectation calculations,and can also solve the problem of high-dimensional optimization by decomposing the objective function.In this paper,MM algorithm is applied to the parameter estimation problem of mixed distribution model.The method of assembly and decomposition is used to construct the substitute function with separable parameters,which avoids the problems of complex expectation calculations and the inversion of high-dimensional matrices.

STUDY ON MIXED MODEL OF NEURAL NETWORK FOR FARMLAND FLOOD/DROUGHT PREDICTION

The paper concerns a flood/drought prediction model involving the continuation of time series of a predictand and the physical factors influencing the change of predictand.Attempt is made to construct the model by the neural network scheme for the nonlinear mapping relation based on multi-input and single output.The model is found of steadily higher predictive accuracy by testing the output from one and multiple stepwise predictions against observations and comparing the results to those from a traditional statistical model.

A compressed variance component mixed model for detecting QTNs and QTN-by-environment and QTN-by-QTN interactions in genome-wide association studies

Although genome-wide association studies are widely used to mine genes for quantitative traits,the effects to be estimated are confounded,and the methodologies for detecting interactions are imperfect.To address these issues,the mixed model proposed here first estimates the genotypic effects for AA,Aa,and aa,and the genotypic polygenic background replaces additive and dominance polygenic backgrounds.Then,the estimated genotypic effects are partitioned into additive and dominance effects using a one-way analysis of variance model.This strategy was further expanded to cover QTN-by-environment interactions(QEIs)and QTN-by-QTN interactions(QQIs)using the same mixed-model framework.Thus,a three-variance-component mixed model was integrated with our multi-locus random-SNP-effect mixed linear model(mrMLM)method to establish a new methodological framework,3VmrMLM,that detects all types of loci and estimates their effects.In Monte Carlo studies,3VmrMLM correctly detected all types of loci and almost unbiasedly estimated their effects,with high powers and accuracies and a low false positive rate.In re-analyses of 10 traits in 1439 rice hybrids,detection of 269 known genes,45 known gene-by-environment interactions,and 20 known gene-by-gene interactions strongly validated 3VmrMLM.Further analyses of known genes showed more small(67.49%),minor-allele-frequency(35.52%),and pleiotropic(30.54%)genes,with higher repeatability across datasets(54.36%)and more dominance loci.In addition,a heteroscedasticity mixed model in multiple environments and dimension reduction methods in quite a number of environments were developed to detect QEIs,and variable selection under a polygenic background was proposed for QQI detection.This study provides a new approach for revealing the genetic architecture of quantitative traits.