【发布时间】：2023-01-25 08:00:36
【问题描述】：
我正在用 brm 函数（brms 包）做 HGAM。现在，我保留了默认的 k = 10，但我想知道是否应该将 k 设置得更高。我知道 mgcv 包中的函数 k.check 可以进行诊断，但我找不到 brms 的任何等效项。
brms 是否存在 mgcv::k.check 的等价物？ 如果不是，我怎么知道我的 k 是否足够？我可以更改它并检查是否合适，但我的时间和计算机能力有限。
谢谢！
【问题讨论】：
-
您可以获取残差样本并将它们传递到
k.check()的本地黑客版本中，以对许多样本运行测试，以获得测试统计的后验，然后总结该后验。 “本地黑客”是指将k.check的内容剥离到您自己的函数中，该函数可以获取所需的输入以及残差样本。 -
完美谢谢！我会试试的。顺便说一下@GavinSimpson，我的所有 GAM 都基于您撰写或合着的论文以及您举办的会议。感谢您的超级有用的工作（我是粉丝！）。