r - nlme 包：augPred 用 log='x' 改变 x 轴

Question

您好我正在尝试将 x 轴上的比例更改为由生成的图上的对数比例

model2 <- nlme(consumption~Q*10^(5*(exp(-alpha*Q*price)-1)), fixed=Q+alpha~1,`random=Q+alpha~1|ratid, groups = ~ treatment,`data = oxydata, start = c(Q=3, alpha=0.001))

plot(augPred(model2), log='x', xlim = c(1,3000), xlab = "Price (resp/mg", ylab = "Consumption")

但是，log='x'orlog="x"不会改变 x 轴刻度。我可以更改 xlim 和 ylim，添加标题等，但我似乎无法将 x 轴更改为对数刻度。非常感谢任何可以或无法完成的帮助或验证。

丑陋的情节，数据全部塞满

标签： rplotnlme

Answer 1

感谢任何试图解决这个问题的人，但我明白了。此函数使用 xyplot()，它需要 scale = list...

`plot(augPred(model2), scale = list(x = list(log = 10)), xlim = c(1,3000), xlab = "Price (resp/mg", ylab = "Consumption")