Question

我试图扩展一个类似于Gavin Simpson描述的模型，这里也包括随机效应。 https://stats.stackexchange.com/questions/32730/how-to-include-an-interaction-term-in-gam

loc被4级因子取代。固定效应是两个分类变量（A和B）和两个连续的变量（C1和C2），具有所有成对相互作用。响应是1s和0s。

e.g。

gamm4(y~A*C1 + A*C2 + B*C1 + B*C2    
                s(yDay, bs = "cc", k = 12) + 
                s(yDay, bs = "cc", by = A, k = 12, m = 1) + 
                s(Hour, bs = "cc", k = 12) + 
                s(Hour, bs = "cc", by = Loc, k = 12, m = 1) +
                t2(Hour, yDay, by = A, bs = rep("cc",2)),
                random = ~ (1 | ID),
                family = "binomial"(link="logit"),
                data = dat)

试图运行它我得到了错误
Error in X %*% diag(diagU[indi]) : non-conformable arguments

从潜伏在gamm4()和gamm()的代码中，它似乎与惩罚随机效应的每个级别内的平滑点有关？

我目前正在运行该模型，其中包含4个人的子集（最终可能大约40个人），并且最终还应该在几年内嵌套个人。

可能是我的问题出在特定情况上（读：我代表我的数据过于野心勃勃），但是阅读gamm4()仍然是＃34;建设阶段＆＃34;在某些问题上让我想知道我是否代表gamm4()雄心勃勃？

或者我可能只是错误而且在一些死胡同里徒步旅行......

gamm4中的张力平滑

0 个答案: