Inline. -- Bert
On Tue, Apr 26, 2016 at 2:25 PM, André Luis Neves <andrl...@ualberta.ca> wrote: > Ok. I`m trying to run a Poisson glmm with an observation-level random > intercept. But I`m getting the following error for the 'Baci' variable: > > 'Error: (maxstephalfit) PIRLS step-halvings failed to reduce deviance in > pwrssUpdate'. I guess this message is because the baci variable is not a > an integer, Why would you "guess" that? Given your magnitudes, I suspect it's a scaling issue. Try centering and scaling your Baci variable before fitting. Cheers, Bert and cannot be transformed into an integer as R has a threshold > of > 2x10^9 even in 64 bit R. > > It runs fine for the fungii variable. > > If you guys want to run the data (attached), the full command is below. > > Thanks. > > --------------------------------------------- > > ##Import data: > > qPCR <- read.delim(file.choose(), > header = TRUE, > dec = ".") > > ##Load package > > library(lme4) > > ##Other steps: > > qPCR$obs <- 1:nrow(qPCR) > qPCR$fID<-as.factor(qPCR$ID) > qPCR$fDiet<-as.factor(qPCR$Diet) > > ##Run the model: > > M1 <- glmer (Baci ~ fDiet + Crossover + (1|fID:Crossover) + (1|obs), > family = poisson, data=qPCR) > > > > Andre > > > > On Tue, Apr 26, 2016 at 2:36 PM, jim holtman <jholt...@gmail.com> wrote: > >> Can you explain why you need them as 'integer', A floating point >> representation can hold a value upto ~4.5e15 as an "integer" keeping the >> precision that you might need. >> >> >> Jim Holtman >> Data Munger Guru >> >> What is the problem that you are trying to solve? >> Tell me what you want to do, not how you want to do it. >> >> On Tue, Apr 26, 2016 at 1:11 PM, André Luis Neves <andrl...@ualberta.ca> >> wrote: >> >>> Dear all: >>> >>> I converted the columns (Baci, Meti, Fungii, Protozoai) into integers >>> (using excel) and then imported the data (.txt) into R. Interestingly, the >>> other three variables were loaded as INT, but the 'Baci' one continued as >>> Num. >>> >>> I imported the data using the following command line: >>> >>> X <- read.delim(file.choose(), >>> header = TRUE, >>> dec = ".") >>> >>> Here is the structure of X: >>> >>> > str(X) >>> 'data.frame': 115 obs. of 5 variables: >>> $ ID : Factor w/ 61 levels "107ZRED","112BLKW",..: 8 12 15 18 26 27 >>> 29 31 32 36 ... >>> $ Baci : num 2.90e+12 5.55e+11 9.46e+11 8.13e+11 4.06e+11 ... >>> $ Meti : int 352645997 334146268 767208656 171567266 462747405 >>> 414905627 237010514 387480048 214671355 328813226 ... >>> $ Fungii : int 43645 19009 15998 2189 8972 8240 3133 17922 6156 13746 >>> ... >>> $ Protozoai: int 3220523 1851891 3252462 1665675 34123768 23175015 >>> 203685 >>> 4261780 43110492 69802572 ... >>> >>> >>> I need Baci as an integer, and tried to convert it using as.integer >>> function, but was not successful. >>> >>> >>> Could anyone please help me to solve this problem. >>> >>> Thanks, >>> >>> >>> >>> -- >>> Andre >>> >>> ______________________________________________ >>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >>> https://stat.ethz.ch/mailman/listinfo/r-help >>> PLEASE do read the posting guide >>> http://www.R-project.org/posting-guide.html >>> and provide commented, minimal, self-contained, reproducible code. >>> >> >> > > > -- > Andre > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
