我正在使用我的响应变量作为百分比(0-100)进行 GAM。我使用反正弦变换来改进模型拟合(asin(sqrt(myvariable/100)))。我现在想评估原始规模上解释因子变量水平之间的对比。我一直在尝试使用 emmeans 并按照转换和链接函数小插图中的步骤来设置我的模型,并以 emmeans 可以读取的格式进行转换。但是,当我运行 emmeans 函数时,出现以下错误: link$mu.eta(object@bhat[estble]) 中的错误:尝试应用非函数。
我这样设置转换对象:
tran <- make.tran("asin.sqrt", 100)
Run Code Online (Sandbox Code Playgroud)
我相信 bit 可以工作,因为当我用 emmeans 在线性模型上尝试它时,它工作了:
warp.t <- with(tran, lm(linkfun(breaks)~wool*tension, warpbreaks))
emmeans(warp.t, ~wool|tension, type="response")
tension = L:
wool response SE df lower.CL upper.CL
A 44.2 4.00 48 36.3 52.3
B 27.7 3.61 48 20.8 35.3
tension = M:
wool response SE df lower.CL upper.CL
A 23.5 3.41 48 17.0 30.7
B 28.4 3.63 48 21.4 35.9
tension = H:
wool response SE df lower.CL …Run Code Online (Sandbox Code Playgroud)