2

我正在尝试复制发表在“线性分位数混合模型”、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 和分位数混合模型很陌生，所以我不确定这是否代表我的错误、编码或缺乏知识！

4

0 回答 0