小编eve*_*lyn的帖子

无法使用 emmeans 进行 arcsin 反向转换

我正在使用我的响应变量作为百分比(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)

r gam lsmeans emmeans

5
推荐指数
1
解决办法
163
查看次数

标签 统计

emmeans ×1

gam ×1

lsmeans ×1

r ×1