2

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

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)

非常感谢你的帮助！

最好的，弗洛里安

4

1 回答 1

3

这是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))
于 2014-08-18T09:25:19.047 回答