lqmm - package: Error in f(arg, …) : NA/NaN/Inf in foreign function call (arg 1)
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:
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