Jika apa yang Anda perlihatkan berfungsi untuk lmer
rumus untuk istilah efek acak maka Anda harus dapat menggunakan fungsi dari paket splines yang datang dengan R untuk mengatur fungsi basis yang relevan.
require("lme4")
require("splines")
lmer(counts ~ dependent_variable + (bs(t) | ID), family="poisson")
Bergantung pada apa yang ingin Anda lakukan, Anda juga harus melihat paket gamm4 dan paket mgcv . Yang pertama pada dasarnya memformalkan bs()
bit dalam lmer()
panggilan di atas dan memungkinkan pemilihan kelancaran untuk dilakukan sebagai bagian dari analisis. Yang terakhir dengan fungsi gam()
memungkinkan untuk beberapa tingkat fleksibilitas dalam model pas seperti ini (jika saya mengerti apa yang Anda coba lakukan). Sepertinya Anda ingin tren terpisah di dalam ID
? Pendekatan efek yang lebih pasti akan seperti:
gam(counts ~ dependent_variable + ID + s(t, by = ID) , family="poisson")
Efek acak dapat dimasukkan dalam gam()
model menggunakan s(foo, bs = "re")
istilah jenis di mana foo
akan berada ID
dalam contoh Anda. Apakah masuk akal untuk menggabungkan by
istilah ide dengan efek acak adalah sesuatu untuk dipikirkan dan bukan sesuatu yang saya memenuhi syarat untuk berkomentar.
Pasang kembali Monica - G. Simpson
sumber