[R-sig-ME] [R] lqmm - package: Error in f(arg, …) : NA/NaN/Inf in foreign function call (arg 1)

[Ben Bolker]

I did some poking around and couldn't figure it out -- it fails fairly deep
inside the function (options(error=recover) is useful ...).  I would ask
the maintainer (see maintainer("lqmm")).

  However, I wonder why you're pulling out these big guns for this
problem.  With only 4 random effects groups, you won't get much shrinkage,
and with so many data points all comparisons will probably be statistically
significant.  Why not just compute the quantiles you're interested in
across groups?


require("lqmm")
stackoverflow <- read.table("mixedModelDataSet.txt",
                            sep="\t", header = TRUE ) # import

so2 <- stackoverflow[complete.cases(stackoverflow),]
mixMoGs15 <- lqmm(gsDeviationMio ~ aoi, random =  ~ 1, group = vpName,
                  data = so2, tau = 0.15)

library(ggplot2); theme_set(theme_bw())
ggplot(so2,aes(x=aoi,y=gsDeviationMio,colour=vpName))+geom_point()
ggplot(so2,aes(x=aoi,y=gsDeviationMio,colour=vpName))+geom_boxplot()
ggplot(so2,aes(x=aoi,y=gsDeviationMio,colour=vpName,fill=vpName))+
    geom_violin()
with(so2,tapply(gsDeviationMio,list(vpName,aoi),quantile,0.15))
with(so2,table(vpName,aoi))



On Sat, Aug 16, 2014 at 9:02 AM, florian k. <zournyque at gmail.com> wrote:

> Dear mailing list,
>
> After some requests I finally managed to make my error reproducible -
> please find the dataset here:
>
> https://dl.dropboxusercontent.com/u/79415744/mixedModelDataSet.txt
>
> require("lqmm")
>
> stackoverflow <- read.table("mixedModelDataSet.txt",  sep="\t", header
> = TRUE ) # import
>
> mixMoGs15 <- lqmm(gsDeviationMio ~ aoi, random =  ~ 1, group = vpName,
> data = stackoverflow, tau = 0.15)
>
> Thanks in advance for your help,
> Best,
> Florian
>
> On Sat, Aug 16, 2014 at 5:57 PM, florian k. <zournyque at gmail.com> wrote:
> > Dear mailing list,
> >
> > After some requests I finally managed to make my error reproducible -
> > find the data set attached
> >
> > require("lqmm")
> >
> > stackoverflow <- read.table("mixedModelDataSet.txt",  sep="\t", header
> > = TRUE ) # import
> >
> > mixMoGs15 <- lqmm(gsDeviationMio ~ aoi, random =  ~ 1, group = vpName,
> > data = stackoverflow, tau = 0.15)
> >
> > Thanks in advance for your help,
> > Best,
> > Florian
> >
> > On Fri, Aug 15, 2014 at 11:30 AM, florian k. <zournyque at gmail.com>
> wrote:
> >> Dear mailing list,
> >>
> >> I keep getting the following error when trying to compute linear
> >> quantile mixed models using the lqmm package:
> >>
> >>  "NA/NaN/Inf in foreign function call (arg 1)"
> >>
> >> For better readability I posted my question on stackoverflow:
> >>
> >>
> http://stackoverflow.com/questions/25323880/error-in-farg-na-nan-inf-in-foreign-function-call-arg-1-r-lqmm-pa
> >>
> >> I would be very thankful if any of you could help me! Thank you for
> your time!
> >>
> >> Best Greetings,
> >> Florian
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>

	[[alternative HTML version deleted]]