拟合LQMM时出现NA / NaN / Inf误差

Question

我试图复制在＆＃34;线性分位数混合模型＆＃34;，Geraci M.和Bottai M.（2013）中发表的一些结果。数据可以从以下网址下载： http://www.hsph.harvard.edu/fitzmaur/ala/tlc.txt

我使用线性分位数混合模型，具有自然三次样条（df = 3，断点0,4,6）和随机效应。我使用

创建了样条线

t1.t2 <-ns(Lead.r$week, df=3, knots=4)
t1.t2
t1 <- t1.t2[1:400, 2]
t1

分位数混合模型

当我对样条曲线和处理（交互和主效应）进行建模时，我可以设法获得输出，如下所示：

fit.int2 <-lqmm(lead ~ t1.t2*group, random= ~1, group=id, 
            data=Lead.r, tau=c(0.25, 0.5, 0.75))

然而，一旦我尝试输入T1和治疗期限（主效应和相互作用），如下所示，我得到一个错误。

fit.int2 <-lqmm(lead ~ t1.t2*group+t1*group, random= ~1, group=id, 
            data=Lead.r, tau=c(0.25, 0.5, 0.75))
Error: NA/NaN/Inf in foreign function call (arg 1)

我对R和分位数混合模型很陌生，所以我不确定这是代表我的错误，是编码还是缺乏知识！