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 general...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.展开更多
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 t...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.展开更多
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...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.展开更多
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 stati...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.展开更多
In this article, we develop a statistical inference technique for the unknown coefficient functions in the varying coeffi- cient model with random effect. A residual-adjusted block empirical likelihood (RABEL) method ...In this article, we develop a statistical inference technique for the unknown coefficient functions in the varying coeffi- cient model with random effect. A residual-adjusted block empirical likelihood (RABEL) method is suggested to inves- tigate the model by taking the within-subject correlation into account. Due to the residual adjustment, the proposed RABEL is asymptotically chi-squared distribution. We illustrate the large sample performance of the proposed method via Monte Carlo simulations and a real data application.展开更多
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 de...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.展开更多
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 t...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.展开更多
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. O...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)).展开更多
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 equivale...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.展开更多
We consider the problem of variable selection for the single-index random effects models with longitudinal data. An automatic variable selection procedure is developed using smooth-threshold. The proposed method share...We consider the problem of variable selection for the single-index random effects models with longitudinal data. An automatic variable selection procedure is developed using smooth-threshold. The proposed method shares some of the desired features of existing variable selection methods: the resulting estimator enjoys the oracle property;the proposed procedure avoids the convex optimization problem and is flexible and easy to implement. Moreover, we use the penalized weighted deviance criterion for a data-driven choice of the tuning parameters. Simulation studies are carried out to assess the performance of our method, and a real dataset is analyzed for further illustration.展开更多
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 Orient...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.展开更多
Based on empirical likelihood method and QR decomposition technique, an orthogonality empirical likelihood based estimation method for the fixed effects in linear mixed effects models is proposed. Under some regularit...Based on empirical likelihood method and QR decomposition technique, an orthogonality empirical likelihood based estimation method for the fixed effects in linear mixed effects models is proposed. Under some regularity conditions, the proposed empirical log-likelihood ratio is proved to be asymptotically chi-squared, and then the confidence intervals for the fixed effects are constructed. The proposed estimation procedure is not affected by the random effects,and then the resulting estimator is more effective. Some simulations and a real data application are conducted for further illustrating the performances of the proposed method.展开更多
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 t...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.展开更多
文摘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.
基金The project is partly supported by NSFC (19971085)the Doctoral Program Foundation of the Institute of High Education and the Special Foundation of Chinese Academy of Sciences.
文摘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.
基金supported by the National Science Foundation of China under Grant Nos.71361015,71340010,71371074the Jiangxi Provincial Natural Science Foundation under Grant No.20142BAB201013+2 种基金China Postdoctoral Science Foundation under Grant No.2013M540534China Postdoctoral Fund special Project under Grant No.2014T70615Jiangxi Postdoctoral Science Foundation under Grant No.2013KY53
文摘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.
基金supported by National Social Science Foundation of China(21BTJ068)。
文摘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.
文摘In this article, we develop a statistical inference technique for the unknown coefficient functions in the varying coeffi- cient model with random effect. A residual-adjusted block empirical likelihood (RABEL) method is suggested to inves- tigate the model by taking the within-subject correlation into account. Due to the residual adjustment, the proposed RABEL is asymptotically chi-squared distribution. We illustrate the large sample performance of the proposed method via Monte Carlo simulations and a real data application.
文摘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.
基金The research project supported by NSFC(1 9631 0 4 0 ) and NSFJ
文摘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.
基金The project supported by NNSFC (19631040), NSSFC (04BTJ002) and the grant for post-doctor fellows in SELF.
文摘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)).
文摘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.
文摘We consider the problem of variable selection for the single-index random effects models with longitudinal data. An automatic variable selection procedure is developed using smooth-threshold. The proposed method shares some of the desired features of existing variable selection methods: the resulting estimator enjoys the oracle property;the proposed procedure avoids the convex optimization problem and is flexible and easy to implement. Moreover, we use the penalized weighted deviance criterion for a data-driven choice of the tuning parameters. Simulation studies are carried out to assess the performance of our method, and a real dataset is analyzed for further illustration.
基金This research received no specific grant from any funding agency in the public,commercial,or not-for-profit sectors
文摘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.
基金Supported by the National Social Science Foundation of China(Grant No.18BTJ035).
文摘Based on empirical likelihood method and QR decomposition technique, an orthogonality empirical likelihood based estimation method for the fixed effects in linear mixed effects models is proposed. Under some regularity conditions, the proposed empirical log-likelihood ratio is proved to be asymptotically chi-squared, and then the confidence intervals for the fixed effects are constructed. The proposed estimation procedure is not affected by the random effects,and then the resulting estimator is more effective. Some simulations and a real data application are conducted for further illustrating the performances of the proposed method.
文摘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.
