「Linear Quantile mixed models」、Geraci M. and Bottai M. (2013) で公開されたいくつかの結果を複製しようとしています。データは以下からダウンロードできます: http://www.hsph.harvard.edu/fitzmaur/ala/tlc.txt
自然な 3 次スプライン (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 と分位点混合モデルにかなり慣れていないので、これが私のせいなのか、コーディング中なのか、知識不足なのかわかりません!