similar to: glmmPQL: Na/NaN/Inf in foreign function call

Hi, I'm having trouble with glmmPQL from the MASS package. I'm trying to fit a model with a binary response variable, two fixed and two random variables (nested), with a sample of about 200,000 data points. Unfortunately, I'm getting an error message that is difficult to understand without knowing the internals of the glmmPQL function. > model <- glmmPQL(primed ~

After upgrading to R 2.6.0, I'm having trouble running lmer: model <- lmer(primed ~ log(dist.time)*role + 1|target.utt, data=data.utts) Error in UseMethod("as.logical") : no applicable method for "as.logical" So I thought I'd upgrade lme4 to the latest version, but unfortunately the compilation fails - perhaps there's a missing #include: R CMD INSTALL

Dear all, I try to fit a glmmPQL on a huge data with 384189 individuals id=1:384189: working in 1520 establishments est:1:1516. The minimum number of individuals in every establishment is 30. This works for a subsample excluding establishemnet cells smaller than 100, but fail when we include smaller cells: R> summary(glmmPQL(count ~ + I( age-ave(age,est) )* ave(age,est) + + I(

I'm trying to use "predict" with a linear mixed-effects logistic regression model fitted with nlmmPQL from the MASS library. Unfortunately, I'm getting an error "non-conformable arguments" in predict.lme, and I would like to understand why. I have used the same call to "predict" with "glm" models without problems. I assume I'm doing

--Apple-Mail-44--797532055 Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes Content-Transfer-Encoding: 7bit I believe read.table may report misleading errors. In this example, where a header line in a file has an incorrect number of row names (28 instead of 29), I get the error message "duplicate row.names are not allowed". However, I cannot not find any

Dear R users, I'm attempting to fit a GLM with random effects using the tweedie family for the error structure. I'm getting the error: iteration 1 Error in logLik.reStruct(object, conLin) : NA/NaN/Inf in foreign function call (arg 3) I'm running V2.1.0 I notice from searching the lists that the same error was reported in May 2004 by Spencer Graves, but no-one was able to

Dear listers, First let me appologize if the same mail arrives multiple times. Recently I had some probelms sending my e-mails to the list. I encountered a problem when running glmmPQL procuedure doing multilevel modeling with a dichotomous outcome. Those are the two error messages I usually get: Error in logLik.reStruct(object, conLin) : NA/NaN/Inf in foreign function call (arg 3)

Can anyone shed any light on this? > doubt.demographic.pql<-glmmPQL(random = ~ 1 | groupid/participantid, + fixed = r.info.doubt ~ + realage + minority + female + education + income + scenario, + data = fgdata.df[coded.resource,], + na.action=na.omit, +

I've got a problem with "mcmcsamp" used with glmer objects produced with "lmer" from the lme4 package. When calling mcmcsamp, I get the error Error in if (var(y) == 0) { : missing value where TRUE/FALSE needed This does not occur with all models, but I can't find anything wrong with the dataset. If the error is in my data, can someone tell me what I am looking

Hi all, I'm writing a script to do some basic text analysis in R. Let's assume I have a data frame named data which contains a column named 'utt' which contains strings. Is there a straightforward way to achieve something like this: data$ContainsThe <- ifelse(startsWith(data$Utt,"the"),"y","n") or data$ContainsThe <-

Asterisk v1.0 is running on RH 9. I installed festival RPM (festival-1.4.2-16.i386.rpm) and edited the festival.scm file to add: (define (tts_textasterisk string mode) "(tts_textasterisk STRING MODE) Apply tts to STRING. This function is specifically designed for use in server mode so a single function call may synthesize the string. This function name may be added to the server safe

Dear R-ers, I am having trouble understanding why I am getting an error using glmmPQL (library MASS). I am getting the following error: iteration 1 Error in MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 The long story: I have data from an experiment on pairwise comparisons between 3 treatments (a, b, c). So a typical run of an experiment

I'm getting a strange error from glmmPQL. Consider the following sample code: set.seed(8) N. <- 1000 z <- rnorm(N.) pr.good <- exp(-1e-4*exp(2+2*z)) quantile(pr.good) DF. <- data.frame(yield=rbinom(N., N., pr.good)/N., Offset=rep(-10, N.), nest=1:N.) fit <- glmmPQL(fixed=1-yield~offset(Offset), random=~1|nest, family=binomial(link="cloglog"),

I know it's conventional to report bugs to the maintainer, but I'm not sure which package actually contains this bug(s), so I apologize for sending this to the list at large. I see the bug under both R 2.1.1, and R 2.2.1. (I sent a related message a while ago, but this one has more detail.) library(MASS) library(nlme) fit.model <- function(il, model.family) { cs <-

Hello, I am having a Mac 10.6.4 (Snow Leopard). I have compiled and built Asterisk 1.6.2.9 and Festival 2.0.95:beta on my machine. Asterisk is working fine with SIP channels without Festival. I have written following context in extension.conf: [connect-to-me] exten => s,1,Answer Exten => s,n,SayDigits(?1?) exten => s,n,Festival(hello john) exten => s,n,Hangup I use call files to

Dear listers, glmmPQL (package MASS) is given to work by repeated call to lme. In the classical outputs glmmPQL the Variance Structure is given as " fixed weights, Formula: ~invwt". The script shows that the function varFixed() is used, though the place where 'invwt' is defined remains unclear to me. I wonder if there is an easy way to specify another variance

I'm a bit puzzled by how to write out additive random effects in glmmPQL. In my situation, I have a factorial design on two (categorical) random factors, A and B. At each combination, I have a binary response, y, and two binary fixed covariates, C and D. If everything were fixed, I would use glm(y ~ A + B + C + D, family = binomial) My first thought was to use glmmPQL(y ~ A + B, random

Dear R-help, I recently tried implementing glmmPQL in 2.3.1, and I discovered a few differences as compared to 2.2.1. I am fitting a regression with fixed and random effects with Gamma error structure. First, 2.3.1 gives different estimates than 2.2.1, and 2.3.1, takes more iterations to converge. Second, when I try using the anova function it says, "'anova' is not available

Hello, I am using the R function glmmPQL to fit a logistic GLMM, with weights. I am finding that I get fairly different parameter estimates in glmmPQL from fitting the full dataset (with no "weight" statement) and an equivalent, shorter dataset with the weights statement. I am using the weights statement in the 'glmmPQL' function exactly as in the 'glm' function. I

Hello R users, I am trying to run a model with a binary response variable (nesting success: 0 failure, 1 success) and 8 fixed terms. Nesting success was examined in 72 cases in 34 territories (TER) during a 6 study years. Territories are nested within 14 patches (PATCH). I want to run a model taking into account these nested factors and repeated observation. To do this, I assume that the best