r - 如何绘制模型 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()-设备窗口放在前面查看图层。需要使用下拉菜单将其拉到最前面。

我（现在）将在代码中演示如何 a）计算 R^2 和 b）在 png() 设备上注释基本图形图形。

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()