如何绘制模型drc包？

Question

脚本：

lettuce

m0<-drm(weight~conc, data = lettuce, fct = LL.3())

summary(m0)

modelFit(m0)

plot(m0)

我想知道如何绘制conc在轴x中的精确值和它们在y中的相应值，以及如何绘制(r²)。

Answer 1

我正在尽我最大的努力来回答这个问题，但它是基于这样的假设:这个问题可能是由于用户界面的困难造成的，尽管我在评论中做出了努力，但术语仍然没有得到澄清。我假设这是一个用R^2值注释drc图的请求。plot.drc函数是一个基本的图形函数，因此默认情况下会绘制到交互式设备上，在我的机器( Mac)上似乎也是如此，它并不总是将quartz()-device窗口显示在视图层的前面。需要使用下拉菜单将其拉到最前面。

我(现在)将用代码演示如何a)计算R^2和b)在png()-device上注释基本图形图形。

png()      # opens file for graphics output
plot(m0)   # will need dev.off() to complete output
cor( fitted(m0), lettuce$weight) ^2  # calculation of R^2
# [1] 0.8143727 # output to interactive console
text( x=20, y=1.2, labels= bquote(R^2 == .(cor( fitted(m0), lettuce$weight) ^2) ) )
dev.off()

​

​