Ошибка NA/NaN/Inf при установке LQMM

Question

Я пытаюсь воспроизвести некоторые результаты, которые были опубликованы в "Смешанных моделях линейного квантиля", 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 и смешанных моделях квантилей, поэтому я не уверен, что это ошибка от моего имени, кодирование или недостаток знаний!