首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 0 毫秒
1.
    
Likelihood analysis for regression models with measurement errors in explanatory variables typically involves integrals that do not have a closed-form solution. In this case, numerical methods such as Gaussian quadrature are generally employed. However, when the dimension of the integral is large, these methods become computationally demanding or even unfeasible. This paper proposes the use of the Laplace approximation to deal with measurement error problems when the likelihood function involves high-dimensional integrals. The cases considered are generalized linear models with multiple covariates measured with error and generalized linear mixed models with measurement error in the covariates. The asymptotic order of the approximation and the asymptotic properties of the Laplace-based estimator for these models are derived. The method is illustrated using simulations and real-data analysis.  相似文献   

2.
This paper is motivated from the analysis of neuroscience data in a study of neural and muscular mechanisms of muscle fatigue. Multidimensional outcomes of different natures were obtained simultaneously from multiple modalities, including handgrip force, electromyography (EMG), and functional magnetic resonance imaging (fMRI). We first study individual modeling of the univariate response depending on its nature. A mixed‐effects beta model and a mixed‐effects simplex model are compared for modeling the force/EMG percentages. A mixed‐effects negative‐binomial model is proposed for modeling the fMRI counts. Then, I present a joint modeling approach to model the multidimensional outcomes together, which allows us to not only estimate the covariate effects but also to evaluate the strength of association among the multiple responses from different modalities. A simulation study is conducted to quantify the possible benefits by the new approaches in finite sample situations. Finally, the analysis of the fatigue data is illustrated with the use of the proposed methods.  相似文献   

3.
    
Neuhaus JM  McCulloch CE  Boylan R 《Biometrics》2011,67(2):654-6; disucssion 656-60
Litière, Alonso, and Molenberghs (2007, Biometrics, 63, 1038-1044) presented the results of simulation studies that they claimed showed that misspecification of the shape of the random effects distribution can produce marked increases in Type II error (decreases in power) of tests based on fits of generalized linear mixed models. However, the article contains a logical fallacy that invalidates this claim. We present logically correct simulation studies that demonstrate little increase in Type II error, consistent with the earlier work that shows little effect due to misspecification.  相似文献   

4.
Bayesian inference for variance components using only error contrasts   总被引:6,自引:0,他引:6  
HARVILLE  DAVID A. 《Biometrika》1974,61(2):383-385
  相似文献   

5.
6.
A note on the logistic link function   总被引:1,自引:0,他引:1  
Kagan  Abram 《Biometrika》2001,88(2):599-601
  相似文献   

7.
A note on estimation for gamma and stable processes   总被引:1,自引:0,他引:1  
BASAWA  I. V.; BROCKWELL  P. J. 《Biometrika》1980,67(1):234-236
  相似文献   

8.
  总被引:1,自引:0,他引:1  
A fundamental challenge to understanding patterns in ecological systems lies in employing methods that can analyse, test and draw inference from measured associations between variables across scales. Hierarchical linear models (HLM) use advanced estimation algorithms to measure regression relationships and variance–covariance parameters in hierarchically structured data. Although hierarchical models have occasionally been used in the analysis of ecological data, their full potential to describe scales of association, diagnose variance explained, and to partition uncertainty has not been employed. In this paper we argue that the use of the HLM framework can enable significantly improved inference about ecological processes across levels of organization. After briefly describing the principals behind HLM, we give two examples that demonstrate a protocol for building hierarchical models and answering questions about the relationships between variables at multiple scales. The first example employs maximum likelihood methods to construct a two-level linear model predicting herbivore damage to a perennial plant at the individual- and patch-scale; the second example uses Bayesian estimation techniques to develop a three-level logistic model of plant flowering probability across individual plants, microsites and populations. HLM model development and diagnostics illustrate the importance of incorporating scale when modelling associations in ecological systems and offer a sophisticated yet accessible method for studies of populations, communities and ecosystems. We suggest that a greater coupling of hierarchical study designs and hierarchical analysis will yield significant insights on how ecological processes operate across scales.  相似文献   

9.
    
In addition to the processes structuring free‐living communities, host‐associated microbiota are directly or indirectly shaped by the host. Therefore, microbiota data have a hierarchical structure where samples are nested under one or several variables representing host‐specific factors, often spanning multiple levels of biological organization. Current statistical methods do not accommodate this hierarchical data structure and therefore cannot explicitly account for the effect of the host in structuring the microbiota. We introduce a novel extension of joint species distribution models (JSDMs) which can straightforwardly accommodate and discern between effects such as host phylogeny and traits, recorded covariates such as diet and collection site, among other ecological processes. Our proposed methodology includes powerful yet familiar outputs seen in community ecology overall, including (a) model‐based ordination to visualize and quantify the main patterns in the data; (b) variance partitioning to assess how influential the included host‐specific factors are in structuring the microbiota; and (c) co‐occurrence networks to visualize microbe‐to‐microbe associations.  相似文献   

10.
    
Analysis of longitudinal data with excessive zeros has gained increasing attention in recent years; however, current approaches to the analysis of longitudinal data with excessive zeros have primarily focused on balanced data. Dropouts are common in longitudinal studies; therefore, the analysis of the resulting unbalanced data is complicated by the missing mechanism. Our study is motivated by the analysis of longitudinal skin cancer count data presented by Greenberg, Baron, Stukel, Stevens, Mandel, Spencer, Elias, Lowe, Nierenberg, Bayrd, Vance, Freeman, Clendenning, Kwan, and the Skin Cancer Prevention Study Group[New England Journal of Medicine 323 , 789–795]. The data consist of a large number of zero responses (83% of the observations) as well as a substantial amount of dropout (about 52% of the observations). To account for both excessive zeros and dropout patterns, we propose a pattern‐mixture zero‐inflated model with compound Poisson random effects for the unbalanced longitudinal skin cancer data. We also incorporate an autoregressive of order 1 correlation structure in the model to capture longitudinal correlation of the count responses. A quasi‐likelihood approach has been developed in the estimation of our model. We illustrated the method with analysis of the longitudinal skin cancer data.  相似文献   

11.
A note on the difference between profile and modified profile likelihood   总被引:1,自引:0,他引:1  
COX  D. R.; REID  N. 《Biometrika》1992,79(2):408-411
  相似文献   

12.
13.
14.
    
Alternative parameterizations and problems of identification and estimation of multivariate random effects models for categorical responses are investigated. The issues are illustrated in the context of the multivariate binomial logit-normal (BLN) model introduced by Coull and Agresti (2000, Biometrics 56, 73-80). We demonstrate that the BLN model is poorly identified unless proper restrictions are imposed on the parameters. Moreover, estimation of BLN models is unduly computationally complex. In the first application considered by Coull and Agresti, an identification problem results in highly unstable, highly correlated parameter estimates and large standard errors. A probit-normal version of the specified BLN model is demonstrated to be underidentified, whereas the BLN model is empirically underidentified. Identification can be achieved by constraining one of the parameters. We show that a one-factor probit model is equivalent to the probit version of the specified BLN model and that a one-factor logit model is empirically equivalent to the BLN model. Estimation is greatly simplified by using a factor model.  相似文献   

15.
    
Auxiliary covariate data are often collected in biomedical studies when the primary exposure variable is only assessed on a subset of the study subjects. In this study, we investigate a semiparametric‐estimated likelihood estimation for the generalized linear mixed models (GLMM) in the presence of a continuous auxiliary variable. We use a kernel smoother to handle continuous auxiliary data. The method can be used to deal with missing or mismeasured covariate data problems in a variety of applications when an auxiliary variable is available and cluster sizes are not too small. Simulation study results show that the proposed method performs better than that which ignores the random effects in GLMM and that which only uses data in the validation data set. We illustrate the proposed method with a real data set from a recent environmental epidemiology study on the maternal serum 1,1‐dichloro‐2,2‐bis(p‐chlorophenyl) ethylene level in relationship to preterm births.  相似文献   

16.
    
Summary .   Motivated by the spatial modeling of aberrant crypt foci (ACF) in colon carcinogenesis, we consider binary data with probabilities modeled as the sum of a nonparametric mean plus a latent Gaussian spatial process that accounts for short-range dependencies. The mean is modeled in a general way using regression splines. The mean function can be viewed as a fixed effect and is estimated with a penalty for regularization. With the latent process viewed as another random effect, the model becomes a generalized linear mixed model. In our motivating data set and other applications, the sample size is too large to easily accommodate maximum likelihood or restricted maximum likelihood estimation (REML), so pairwise likelihood, a special case of composite likelihood, is used instead. We develop an asymptotic theory for models that are sufficiently general to be used in a wide variety of applications, including, but not limited to, the problem that motivated this work. The splines have penalty parameters that must converge to zero asymptotically: we derive theory for this along with a data-driven method for selecting the penalty parameter, a method that is shown in simulations to improve greatly upon standard devices, such as likelihood crossvalidation. Finally, we apply the methods to the data from our experiment ACF. We discover an unexpected location for peak formation of ACF.  相似文献   

17.
    
Summary In recent years, nonlinear mixed‐effects (NLME) models have been proposed for modeling complex longitudinal data. Covariates are usually introduced in the models to partially explain intersubject variations. However, one often assumes that both model random error and random effects are normally distributed, which may not always give reliable results if the data exhibit skewness. Moreover, some covariates such as CD4 cell count may be often measured with substantial errors. In this article, we address these issues simultaneously by jointly modeling the response and covariate processes using a Bayesian approach to NLME models with covariate measurement errors and a skew‐normal distribution. A real data example is offered to illustrate the methodologies by comparing various potential models with different distribution specifications. It is showed that the models with skew‐normality assumption may provide more reasonable results if the data exhibit skewness and the results may be important for HIV/AIDS studies in providing quantitative guidance to better understand the virologic responses to antiretroviral treatment.  相似文献   

18.
基于rDNA ITS序列对绒泡菌目黏菌系统发育的探讨   总被引:1,自引:0,他引:1  
李倩  闫淑珍  陈双林 《菌物学报》2015,34(3):424-433
绒泡菌目Physarida是黏菌纲Myxogastria最大的一个目,对其系统发育关系的研究一直是根据形态特征。为了从分子水平探讨绒泡菌目乃至黏菌纲的系统发育关系,以黏菌r DNA ITS通用引物对绒泡菌目5属8种黏菌的r DNA ITS进行扩增和测序,结合Gen Bank中已有的黏菌r DNA ITS序列,利用贝叶斯推断法(Bayesian inference,BI)和最大似然法(Maximum likelihood,ML)构建系统发育树。结果表明:绒泡菌目不同物种的r DNA ITS区在碱基组成和长度上差异明显,长度为777–1 445bp,G+C mol%在53.4%–61.9%之间。绒泡菌目与发网菌目Stemonitida聚类为两个明显的分支,在绒泡菌目分支上,绒泡菌科Physaraceae和钙皮菌科Didymiaceae各聚为一支,支持了形态学上以孢丝是否具有石灰质为依据区分这两个科的观点。由多份不同地理来源的鳞钙皮菌Didymium squamulosum材料组成的钙皮菌科又形成3个分支,证实了这个形态种是由地域来源广泛、繁殖亲和性各异和遗传变异较大的不同生物种组成的复合体。  相似文献   

19.
20.
A note on the generalized information criterion for choice of a model   总被引:3,自引:0,他引:3  
ATKINSON  A. C. 《Biometrika》1980,67(2):413-418
  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号