我想计算线性分位数混合模型,但总是出现以下错误
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)
非常感谢您的帮助!
最好, 弗洛里安
最佳答案
这是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))
关于r - 线性分位数混合模型 [R] lqmm - 包 : Error in f(arg, ...) : 外部函数调用中的 NA/NaN/Inf (arg 1),我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/25341816/