similar to: p values in lme4 package

Hi, I have not been able to calculate residuals in the lme4 package. I've been trying the resid() function after I ran a GLMM with the lmer() function, but I get an error message that says "residuals are not inserted yet". I looked it up in the "help" history and I realized that several people have had this problem in the past, related to some bug in this function and

Hi all, I was wondering if somebody could give me advice regarding the dispersion parameter in GLMM's. I'm a beginner in R and basically in GLMM's. I've ran a GLMM with a poisson family and got really nice results that conform with theory, as well with results that I've obtained previously with other analysis and that others have obtained in similar studies. But the

library(glmmADMB) #Example for glmm.admb data(epil2) glmm.admb(y~Base*trt+Age +Visit,random=~Visit,group="subject",data=epil2,family="nbinom") Gives: Error in glmm.admb(y ~ Base * trt + Age + Visit, random = ~Visit, group = "subject", : The function maximizer failed ****************** R version 2.4.1 RC (2006-12-14 r40181) powerpc-apple-darwin8.8.0 locale: C

Dear R users I have built several glmm models using glmmPQL in the following structure: m1<-glmmPQL(dev~env*har*treat+dens, random = ~1|pop/rep, family = Gamma) (full script below, data attached) I have tried all the methods I can find to obtain some sort of model fit score or to compare between models using following the deletion of terms (i.e. AIC, logLik, anova.lme(m1,m2)), but I

Hi, Does anyone know if: with R can you take a set of numbers and aggregate them like you can in SPSS? For example, could you calculate the percentage of people who smoke based on a dataset like the following: smoke = 1 non-smoke = 2 variable 1 1 1 2 2 1 1 1 2 2 2 2 2 2 When aggregated, SPSS can tell you what percentage of persons are smokers based on the frequency of 1's and 2's. Can

Hello, I am trying to use the lmer function from the lme4 package I have installed today (lme4_0.999375-32.zip; R-2.10.1). According to the information, I should be able to use a generalized linear mixed model. However when I try to fit a model with Gamma distribution of the errors, it gives me the following error model1<-lmer(Cmic~1+(1|FOREST/DEPTH),data=DATOS,family=Gamma) Error en

Dear community, I'm trying to track down a quote, but can't recall the source or the exact structure - not very helpful, I know - something along the lines that: 80% of [applied] statistics is linear regression ... ? Does this ring a bell for anyone? Thanks, Andrew -- Andrew Robinson Department of Mathematics and Statistics Tel: +61-3-8344-9763 University of

Hello, I posted a question yesterday but I got no replies, so I'll try to reformulate it in a more concise way. I have the following data, summarizing approval ratings on two different surveys for a random sample of 1600 individuals: > ## Example: Ratings of prime minister (Agresti, Table 12.1, p.494) > rating <- matrix(c(794, 86, 150, 570), 2, 2) > dimnames(rating) <-

Hi there, I've been wanting to fit a GLMM and I'm not completely sure I'm doing things right. As I said in a previous message my response variable is continuous with many zeros, so I was having a hard time finding an appropriate error distribution. I read some previous help mails given to other people advising them to use the Tweedie distribution. I'm still not sure if this

Using the weight argument with a variance function in lme (nlme), you can allow for heteroscedasticity of the within-group error. Is there a way to do this for the other variance components? For example, suppose you had subjects, days nested within subjects, and visits nested within days within subjects (a fully nested two-way design) and you had, say men and women subjects. Could you allow for

Hi all, Does anybody have the experience of using optim to estimate variables with integral forms? here the code: trun.mean<- function(x) # t is the threshold { mu=x[1]; sigma=x[2]; t=x[3]; f <- function(x) (1/(sigma*sqrt(2*pi)))*exp(-(x-mu)^2/(2*sigma^2)); pdf.fun <- function(x) x*f(x); integrate(f,thre,upper=Inf)$value/integrate(pdf.fun,thre,upper=Inf)$value ; } when I

Dear All, I tried to use the boot function, provided in the boot package, in such a simple task as to create a bootstrap distribution of the mean of a vector x. I wrote: b <- boot(x, mean, R=200) Well, it doesn't work. I suspect it has something to do with what is called "second argument" of the "statistic" in the help page of boot. What is the second argument of

Hi, I have to build a logistic regression model on a data set that I have. I have three input variables (x1, x2, x3) and one output variable (y). The syntax of lrm function looks like this lrm(formula, data, subset, na.action=na.delete, method="lrm.fit", model=FALSE, x=FALSE, y=FALSE, linear.predictors=TRUE, se.fit=FALSE, penalty=0, penalty.matrix, tol=1e-7,

Hi everyone, I am fitting a bivariate smoothing model by using gam. But I got an error message like this: "Error in eigen(hess1, symmetric = TRUE) : 0 x 0 matrix" If anyone know how to figure it out, pleaselet me know. Thanks very much. [[alternative HTML version deleted]]

Hi everybody, I am using predict.gam() now. I but it seems there is no such option to get standard errors of the predicted values. I tried to set se=T or se.fit=T but no use. If you know anything about that please let me know. Thanks very much. Kevin. [[alternative HTML version deleted]]

[Originally sent this to r-help at lists.R-projects.org, but in case that's the wrong list I'm re-posting. Apologies if this becomes a re-post] -------------- next part -------------- An embedded message was scrubbed... From: Jonathon Kopecky <jkopecky at umich.edu> Subject: Need to fit a regression line using orthogonal residuals Date: Mon, 29 Jan 2007 14:52:24 -0500 Size: 1138

I'm investigating a number of dependent variables using mixed models, e.g. data.lmer45 = lmer(ampStopB ~ (type + stress + MorD)^3 + (1|speaker) + (1|word), data=data) The p-values for some of the 2-way and 3-way interactions are significant at a 0.05 level and I have been trying to find out how to understand the exact nature of the interactions. Does anyone know if it is possible to run

Hi, I try to make a model using lmer, but the weigths is not accept. m1<-lmer(ocup/total~tempo+(tempo|estacao),family=binomial,weights=total) Erro em lmer(ocup/total ~ tempo + (tempo | estacao), family = binomial, : object `weights' of incorrect type I dont understand why this error, with glm this work. the total object is a vector. Any idea? Thanks Ronaldo -- God is subtle, but

I`m doing some functions on C that gives me the x and y coordinates. I`d like to now how I can get these coordinates (both are a vector of number) on R to that I can make a graphic. I`ve already made a package with my functions, so I just wanna how about how to get the coordinates. Thanks, Heloise.

hello i have a mixed effect model which gives slope and intercept terms for 6 groups (diagnosis (3 levels) by risk group(2 levels)). the fixed part of the model is -- brain volume ~ Diagnosis + Risk Group + (Risk Group * age : Diagnosis) - 1 thus allowing risk group age/slope terms to vary within diagnosis and omitting a nonsignificant diagnosis by risk group intercept (age was centered)