2

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

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

例如

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()在某些问题上阅读仍然是“建设阶段”让我想知道我是否代表雄心勃勃gamm4()?

或者我可能只是完全错了，并且在一些死胡同中远足......

4

0 回答 0