[R] From NUM to INT

[Bert Gunter]

Inline.

-- Bert


On Tue, Apr 26, 2016 at 2:25 PM, André Luis Neves <andrluis at 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 <jholtman at 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 <andrluis at 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 at 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 at 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.