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.

EMPIRICAL BAYES TEST PROBLEMS OF VARIANCE COMPONENTS IN RANDOM EFFECTS MODEL

Bayes decision rule of variance components for one-way random effects model is derived and empirical Bayes (EB) decision rules are constructed by kernel estimation method. Under suitable conditions, it is shown that the proposed EB decision rules are asymptotically optimal with convergence rates near O(n-1/2). Finally, an example concerning the main result is given.

TESTING FOR VARYING DISPERSION OF LONGITUDINAL BINOMIAL DATA IN NONLINEAR LOGISTIC MODELS WITH RANDOM EFFECTS

In this paper, it is discussed that two tests for varying dispersion of binomial data in the framework of nonlinear logistic models with random effects, which are widely used in analyzing longitudinal binomial data. One is the individual test and power calculation for varying dispersion through testing the randomness of cluster effects, which is extensions of Dean(1992) and Commenges et al (1994). The second test is the composite test for varying dispersion through simultaneously testing the randomness of cluster effects and the equality of random-effect means. The score test statistics are constructed and expressed in simple, easy to use, matrix formulas. The authors illustrate their test methods using the insecticide data (Giltinan, Capizzi & Malani (1988)).

INFLUENCE ANALYSIS ON EXPONENTIAL NONLINEAR MODELS WITH RANDOM EFFECTS

This paper presents a unified diagnostic method for exponential nonlinear models with random effects based upon the joint likelihood given by Robinson in 1991. The authors show that the case deletion model is equivalent to mean shift outlier model. From this point of view, several diagnostic measures, such as Cook distance, score statistics are derived. The local influence measure of Cook is also presented. Numerical example illustrates that our method is available.

INFLUENCE ANALYSIS IN NONLINEAR MODELS WITH RANDOM EFFECTS

In this paper,a unified diagnostic method for the nonlinear models with random effects based upon the joint likelihood given by Robinson in 1991 is presented.It is shown that the case deletion model is equivalent to the mean shift outlier model.From this point of view,several diagnostic measures,such as Cook distance,score statistics are derived.The local influence measure of Cook is also presented. A numerical example illustrates that the method is available.

Bootstrap inference of the skew-normal two-way classification random effects model with interaction

In this paper,we consider the statistical inference problems for the fixed effect and variance component functions in the two-way classification random effects model with skewnormal errors.Firstly,the exact test statistic for the fixed effect is constructed.Secondly,using the Bootstrap approach and generalized approach,the one-sided hypothesis testing and interval estimation problems for the single variance component,the sum and ratio of variance components are discussed respectively.Further,the Monte Carlo simulation results indicate that the exact test statistic performs well in the one-sided hypothesis testing problem for the fixed effect.And the Bootstrap approach is better than the generalized approach in the one-sided hypothesis testing problems for variance component functions in most cases.Finally,the above approaches are applied to the real data examples of the consumer price index and value-added index of three industries to verify their rationality and effectiveness.

Bias and Mean Square Error of Reliability Estimators under the One and Two Random Effects Models: The Effect of Non-Normality

The coefficient of reliability is often estimated from a sample that includes few subjects. It is therefore expected that the precision of this estimate would be low. Measures of precision such as bias and variance depend heavily on the assumption of normality, which may not be tenable in practice. Expressions for the bias and variance of the reliability coefficient in the one and two way random effects models using the multivariate Taylor’s expansion have been obtained under the assumption of normality of the score (Atenafu et al. [1]). In the present paper we derive analytic expressions for the bias and variance, hence the mean square error when the measured responses are not normal under the one-way data layout. Similar expressions are derived in the case of the two-way data layout. We assess the effect of departure from normality on the sample size requirements and on the power of Wald’s test on specified hypotheses. We analyze two data sets, and draw comparisons with results obtained via the Bootstrap methods. It was found that the estimated bias and variance based on the bootstrap method are quite close to those obtained by the first order approximation using the Taylor’s expansion. This is an indication that for the given data sets the approximations are quite adequate.

Comparison of REML and MINQUE for Estimated Variance Components and Predicted Random Effects

Linear mixed model (LMM) approaches have been widely applied in many areas of research data analysis because they offer great flexibility for different data structures and linear model systems. In this study, emphasis is placed on comparing the properties of two LMM approaches: restricted maximum likelihood (REML) and minimum norm quadratic unbiased estimation (MINQUE) with and without resampling techniques being included. Bias, testing power, Type I error, and computing time were compared between REML and MINQUE approaches with and without Jackknife technique based on 500 simulated data sets. Results showed that MINQUE and REML methods performed equally regarding bias, Type I error, and power. Jackknife-based MINQUE and REML greatly improved power compared to non-Jackknife based linear mixed model approaches. Results also showed that MINQUE is more time-saving compared to REML, especially with the use of resampling techniques and large data set analysis. Results from the actual cotton data analysis were in agreement with our simulated results. Therefore, Jackknife-based MINQUE approaches could be recommended to achieve desirable power with reduced time for a large data analysis and model simulations.

Additive Hazards Regression with Random Effects for Clustered Failure Times

Additive hazards model with random effects is proposed for modelling the correlated failure time data when focus is on comparing the failure times within clusters and on estimating the correlation between failure times from the same cluster, as well as the marginal regression parameters. Our model features that, when marginalized over the random effect variable, it still enjoys the structure of the additive hazards model. We develop the estimating equations for inferring the regression parameters. The proposed estimators are shown to be consistent and asymptotically normal under appropriate regularity conditions. Furthermore, the estimator of the baseline hazards function is proposed and its asymptotic properties are also established. We propose a class of diagnostic methods to assess the overall fitting adequacy of the additive hazards model with random effects. We conduct simulation studies to evaluate the finite sample behaviors of the proposed estimators in various scenarios. Analysis of the Diabetic Retinopathy Study is provided as an illustration for the proposed method.

Testing for Random Effects in Linear Mixed Models for Longitudinal Data under Moment Conditions

In this paper, we consider whether the random effect exists in linear mixed models （LMMs） when only moment conditions are assumed. Based on the estimators of parameters and their asymptotic properties, a Wald-type test is constructed. It is consistent against global alternatives and is sensitive to the local alternatives converging to the null hypothesis at parametric rates, a fastest possibly rate for goodness-of-fit testing. Moreover, a simulation study shows the performance of the test is good. The procedure also applies to a real data.

Bootstrap Inference on the Variance Component Functions in the Two-Way Random Effects Model with Interaction

In this paper,using the Bootstrap approach and generalized approach,the authors consider the one-sided hypothesis testing problems for variance component functions in the two-way random effects model.Firstly,the test statistics and confidence intervals for the sum of variance components are constructed.Next,the one-sided hypothesis testing problems for the ratio of variance components are also discussed.The Monte Carlo simulation results indicate that the Bootstrap approach is better than the generalized approach in most cases.Finally,the above approaches are applied to the real data examples of mice blood p H and molded plastic part’s dimensions.

Mixed-effects modeling for tree height prediction models of Oriental beech in the Hyrcanian forests

Height–diameter relationships are essential elements of forest assessment and modeling efforts.In this work,two linear and eighteen nonlinear height–diameter equations were evaluated to find a local model for Oriental beech（Fagus orientalis Lipsky） in the Hyrcanian Forest in Iran.The predictive performance of these models was first assessed by different evaluation criteria： adjusted R^2（R^2_（adj））,root mean square error（RMSE）,relative RMSE（%RMSE）,bias,and relative bias（%bias） criteria.The best model was selected for use as the base mixed-effects model.Random parameters for test plots were estimated with different tree selection options.Results show that the Chapman–Richards model had better predictive ability in terms of adj R^2（0.81）,RMSE（3.7 m）,%RMSE（12.9）,bias（0.8）,%Bias（2.79） than the other models.Furthermore,the calibration response,based on a selection of four trees from the sample plots,resulted in a reduction percentage for bias and RMSE of about 1.6–2.7%.Our results indicate that the calibrated model produced the most accurate results.

Optimal Credibility Estimation of Random Parameters in Hierarchical Random Effect Linear Model

In the hierarchical random effect linear model, the Bayes estimator of random parameter are not only dependent on specific prior distribution but also it is difficult to calculate in most cases. This paper derives the distributed-free optimal linear estimator of random parameters in the model by means of the credibility theory method. The estimators the authors derive can be applied in more extensive practical scenarios since they are only dependent on the first two moments of prior parameter rather than on specific prior distribution. Finally, the results are compared with some classical models and a numerical example is given to show the effectiveness of the estimators.

Hydrodynamic Analysis of Elastic Floating Collars in Random Waves

As the main load-bearing component of fish cages, the floating collar supports the whole cage and undergoes large deformations. In this paper, a mathematical method is developed to study the motions and elastic deformations of elastic floating collars in random waves. The irregular wave is simulated by the random phase method and the statistical approach and Fourier transfer are applied to analyze the elastic response in both time and frequency domains. The governing equations of motions are established by Newton＇s second law, and the governing equations of deformations are obtained based on curved beam theory and modal superposition method. In order to validate the numerical model of the floating collar attacked by random waves, a series of physical model tests are conducted. Good relationship between numerical simulation and experimental observations is obtained. The numerical results indicate that the transfer function of out-of-plane and in-plane deformations increase with the increasing of wave frequency. In the frequency range between 0.6 Hz and 1.1 Hz, a linear relationship exists between the wave elevations and the deformations. The average phase difference between the wave elevation and out-of-plane deformation is 60° with waves leading and the phase between the wave elevation and in-plane deformation is 10° with waves lagging. In addition, the effect of fish net on the elastic response is analyzed. The results suggest that the deformation of the floating collar with fish net is a little larger than that without net.

Ionizing radiation effect on single event upset sensitivity of ferroelectric random access memory

The impact of ionizing radiation effect on single event upset（SEU） sensitivity of ferroelectric random access memory（FRAM） is studied in this work. The test specimens were firstly subjected to ^60Co γ-ray and then the SEU evaluation was conducted using ^209Bi ions. As a result of TID-induced fatigue-like and imprint-like phenomena of the ferroelectric material, the SEU cross sections of the post-irradiated devices shift substantially. Different trends of SEU cross section with elevated dose were also found, depending on whether the same or complementary test pattern was employed during the TID exposure and the SEU measurement.

Seismic loss assessment of RC high-rise buildings designed according to Eurocode 8

A probabilistic seismic loss assessment of RC high-rise(RCHR)buildings designed according to Eurocode 8 and located in the Southern Euro-Mediterranean zone is presented herein.The loss assessment methodology is based on a comprehensive simulation approach which takes into account ground motion(GM)uncertainty,and the random effects in seismic demand,as well as in predicting the damage states(DSs).The methodology is implemented on three RCHR buildings of 20-story,30-story and 40-story with a core wall structural system.The loss functions described by a cumulative lognormal probability distribution are obtained for two intensity levels for a large set of simulations(NLTHAs)based on 60 GM records with a wide range of magnitude(M),distance to source(R)and different site soil conditions(SS).The losses expressed in percent of building replacement cost for RCHR buildings are obtained.In the estimation of losses,both structural(S)and nonstructural(NS)damage for four DSs are considered.The effect of different GM characteristics(M,R and SS)on the obtained losses are investigated.Finally,the estimated performance of the RCHR buildings are checked to ensure that they fulfill limit state requirements according to Eurocode 8.

二元随机效应模型的相关性问题研究

文章讨论二元随机效应模型的相关性问题。对于已给数据,拟合二元随机效应模型中两个随机效应相关系数的估计值ρ,能产生对第i次观测随机效应的估计值ui1、ui2,并计算样本相关系数r。结果显示,r和ρ有时差距较大,甚至ρ不显著,而r却是显著的,出现"不一致"现象,特别是"独立随机效应模型"(ρ恒为0)的r可能显著。有时优化过程中ρ非常接近-1,从而使优化过程非正常终止,此时可以采用RE(τ)模型。

Parameter Estimates in Random Intercept Mixed Effects Model for Repeated Measures

In this article the following random intercept mixed effects model will be considered： yij = vi =v^τijβ＋ εij,i=1,…,m;j=1,2,…,ni, where {vi} are i.i.d, random effects with mean α 2. 2 and finite variance σ^2 v, {εij} are i.i.d, random errors with finite variance ε^2 ε. Here we will estimate α,σ^2 v,σ^2 ε,β and study their large sample properties, such as strong consistency, strong convergence rates and asymptotic normality.

TESTING FOR VARYING DISPERSION IN DISCRETE EXPONENTIAL FAMILY NONLINEAR MODELS

It is necessary to test for varying dispersion in generalized nonlinear models.Wei,et al(1998) developed a likelihood ratio test,a score test and their adjustments to test for varying dispersion in continuous exponential family nonlinear models.This type of problem in the framework of general discrete exponential family nonlinear models is discussed.Two types of varying dispersion,which are random coefficients model and random effects model,are proposed,and corresponding score test statistics are constructed and expressed in simple,easy to use,matrix formulas.

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.