similar to: Cause of error message in cov function?

Hi, I'm still having problems putting the variance components of my model in to a data frame, it is a continuation of this discussion, http://r.789695.n4.nabble.com/ANOVA-problem-td4609062.html, but now focussed on the problem of extracting variance components. I have got my mixed effects model now /narrow$line<-as.factor(narrow$line) rg.lmer <- lapply(split(narrow,

Hello, I am trying to use varcomp to decompose the variance across multiple nested levels on a lme object. I am able to successfully do this and when I view the varcomp object I can see the individual values / estimates for the variance at different levels. However, I want to be able to extract each of them separately, as I need to build a confidence interval using bootstrapping on the sample

Hi try str(varcompobject) to see structure of this object. You can extract parts by standard R means. Cheers Petr > -----Original Message----- > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Sharada > Ramadass > Sent: Tuesday, August 8, 2017 3:33 PM > To: r-help at r-project.org > Subject: [R] how to extract individual values from varcomp? > >

Hello, I have some data covering contaminant concentrations in fish over a time period of ~35 years. Each year, multiple samples of fish were taken (with varying sample sizes each year). Ultimately, I want an estimation of the variance between years, and the variance within years + random effects. I used a linear mixed model to estimate these variances, but after reading a number of different

I am trying to use varcomp to obtain the variance partitioning across different nested levels of random effects (say x,y and z). I get the three variance components (for each of my along with an additional one called 'within' from varcomp output. I am using the 'scale total variance to 1' option and though the within component is small, it does form a part of what explains the

Hi Everyone, I am using ASReml to fit a spatial model. I do not have all the components of ASReml when I call; names(summary()) e.g. names(summary(fcov.asr2)) [1] "call" "loglik" "nedf" "sigma" "varcomp" I am trying to get the coefs but I get "NULL". Does anybody know the reason? Any help would be much appreciated. Regards

Hi, I'm trying to fit a model using variance components in R, but if very new on it, so I'm asking for your help. I have imported the SPSS database onto R, but I don't know how to convert the commands... the SPSS commands I'm trying to convert are: VARCOMP RATING BY CHAIN SECTOR RESP ASPECT ITEM /RANDOM = CHAIN SECTOR RESP ASPECT ITEM /METHOD = MINQUE (1) /DESIGN

Dear all, I'm using zeroinfl() from the pscl-package for zero inflated Poisson regression. I would like to calculate (aproximate) prediction intervals for the fitted values. The package itself does not provide them. Can this be calculated analyticaly? Or do I have to use bootstrap? What I tried until now is to use bootstrap to estimate these intervals. Any comments on the code are welcome.

I'm running simulation using lme and sometimes the estimated variance-covariance matrix is not positive definite so that the intervals function won't work for the random effect coefficients. I've tried varcomp from the ape package but this does not return all the coefficients. How can I extract the random effect coefficients without using intervals? Rick B.

Dear all, I'd like to remove the space between the x-axis and the 0 of the y-axis in the hist function. I saw the previous post https://stat.ethz.ch/pipermail/r-help/2012-September/324177.html suggesting adding the x-axis at position 0 (or the lowest value of y) after making the histogram. Unfortunately, it doesn't work in my case as I want to use it with a histogram with a ylim starting

Dear All, I need to calculate the variance components in a mixed effect model (one fixed and one random effect) with REML (maximizing the proportion of the likelihood that does not depend on the fixed effects). In S+ there is the varcomp function, but I would like to do it in R. Is there a way to do that? Thanks! Katalin ___ Katalin Csillery Division of Biological Sciences University of

Hello, I just installed R-2.15.0 on windows XP and cannot load package tcltk2 (which I just downloaded from CRAN as tcltk2_1.2-1.zip; package install reported no problems): > library(tcltk2) Carico il pacchetto richiesto: tcltk Loading Tcl/Tk interface ... done Error : .onLoad failed in loadNamespace() for 'tcltk2', details: call: system("cat /etc/issue", intern = TRUE,

I am trying to understand better an analysis mean RT in various conditions in a within subjects design with the overall mean RT / subject as one of the factors. LME seems to be the right way to do this. using something like m<- lme(rt~ a *b *subjectRT, random= ~1|subject) and then anova(m,type = "marginal"). My understanding is that lme is an easy interface for dummy coding

Hi all I am trying to generate a normal unbalanced data to estimate the coefficients of LM, LMM, GLM, and GLMM and their standard errors. Also, I am trying to estimate the variance components and their standard errors. Further, I am trying to use the likelihood ratio test to test H0: sigma^2_b = 0 (random effects variance component), and the t-test to test H0:mu=0 (intercept of the model Yij = mu

Hello, I am having trouble loading Matrix (0.96-3) in self-compiled R-2.1.1- beta, and self compiled R.app (got it a couple of days ago via: svn co https://svn.r-project.org/R-packages/trunk/Mac-GUI Mac-GUI ) on Mac OS 10.4.1. The problem I get when I try to load Matrix is the following: (I know I do not need to worry about the warnings, even though it might be nice to know why they

Hi, apparently I sent my question about using R and C++ to the wrong list, ironically seeing as that list was called Rcpp. Anyway, I was directed to post my question here. To summarize my current question, I have found two commands that I want to be able to put into a package. The commands are 'R CMD SHLIB X.cc X_main.cc' and

Hello all R-users! In Jun-2005, I find the follow discussion about using of LME function ( in NLME library ) for fitting non-replicate data The thread: ANOVA vs REML approach to variance component estimation http://tolstoy.newcastle.edu.au/R/help/05/06/6498.html Someone expose the follow problem: # non-replicate data y <- c(2.2, -1.4, -0.5, -0.3, -2.1, 1.5, 1.3, -0.3, 0.5, -1.4,

I have some time data and which is in seconds time <-c( 126230400 126252000 126273600 126295200 126316800 126338400) now I wanted to convert this time to Y M D H M S format I have tried following codes but it does not give me the out put in Y M D H M S time_t1 <- as.POSIXlt(time, origin="2005-01-01", tz="GMT") & time_f <- as.POSIXct(time,

Dear List, Is there a way to extract variance components from lmeObjects or summary.lme objects without using intervals()? For my purposes I don't need the confidence intervals which I'm obtaining using parametric bootstrap. Thanks, Mike [[alternative HTML version deleted]]

Hello all, This may be aimed for r-devel, but I encountered this as an R-user and not an R-developer so I start here (having said that, please direct me to R-devel if you think this is appropriate. I am not cross- posting, as I believe this is bad netiquette). I am a recent, but extremely happy R-user (especially after getting my own copy of MASS 2002). My adventures started when I wanted