线性分位数混合模型 [R] lqmm - 包：f(arg, ...) 中的错误：外部函数调用中的 NA/NaN/Inf (arg 1)答案

【问题标题】：linear quantile mixed model [R] lqmm - package: Error in f(arg, ...) : NA/NaN/Inf in foreign function call (arg 1)线性分位数混合模型 [R] lqmm - 包：f(arg, ...) 中的错误：外部函数调用中的 NA/NaN/Inf (arg 1)
【发布时间】：2014-08-16 16:47:54
【问题描述】：

我想计算线性分位数混合模型，但总是出现以下错误

Error in f(arg, ...) : NA/NaN/Inf in foreign function call (arg 1)

要重现，请下载数据集并导入：

https://dl.dropboxusercontent.com/u/79415744/mixedModelDataSet.txt

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

然后尝试计算模型：

require("lqmm")
stack15 <- lqmm(gsDeviationMio ~ aoi, random =  ~ 1, group = vpName, data = stackoverflow, tau = 0.15)

我做错了什么？

计算非分位数混合模型有效：

stackLme <- lme(gsDeviationMio ~ aoi, random =  ~ 1|vpName, data = stackoverflow)

非常感谢您的帮助！

最好，弗洛里安

【问题讨论】：

标签： r mixed-models quantile

【解决方案1】：

这是 Marco Geraci（lqmm 的作者）的答案

响应的规模存在问题。此外，“gs”算法似乎与这个数据集有一些问题。试试下面的

stackoverflow$y <- scale(stackoverflow$gsDeviationMio, center = T, scale = T)

lqmm(y ~ aoi, random =  ~ 1, group = vpName, data = stackoverflow, tau = 0.15, control = lqmmControl(method = "df", UP_max_iter = 200))

【讨论】：