lqmm - package: Error in f(arg, …) : NA/NaN/Inf in foreign function call (arg 1)
Hi Ben, thank you for your time and effort, I will contact the project maintainer as you suggested. I provided only a simplified version of my model and dataset to reproduce the error. My original model has twenty variables and the dataset has 42 subjects with around 400k observations. The plots you created look beautiful but I am not sure whether I understand your code completely. Especially those two lines: with(so2,tapply(gsDeviationMio,list(vpName,aoi),quantile,0.15)) with(so2,table(vpName,aoi)) However I think they will be hard to absorb when I do the same procedure for all of my 42 test subjects. So here is some background information: Being a psychology student and still a novice to R and statistics in general, I am not even sure whether the "big gun" is the right method here. So here is some context. The data was obtained in a flight simulator. The gsDeviationMio variable is a measurement for how far a pilot deviates from his original landing course. 0 is no deviation, so perfect landing course - higher or lower numbers indicate deviation. I multiplied all values with 1 million as the values usually oscillate between 0 and 1 and I had some issues with rounding going on when calculating quantile regressions. Now I want to investigate the relationship between a pilot's eye movements and the flying performance. Variable Aoi are areas of interest so in other words where the pilot is looking at. First I calculated linear models and quantile regressions for RMSD course deviations, now I want to investigate single observations. Therefore I learned I need to delve into mixed models. I soon calculated my first linear mixed model. (I guess I should apply a transformation here, so how about square rooting everything to get rid of the negative values of gsDeviationMio?). The problem is theoretically average flight performance is not that interesting, I am much more interested in predictors for either very good or very bad flight performance. Therefore I wanted to use lqmm. I hope this makes sense? Hope you have a nice Sunday, Florian
On Sun, Aug 17, 2014 at 12:45 AM, Ben Bolker <bbolker at gmail.com> wrote:
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