如何绘制模型 drc 包?
How to plot a model drc package?
脚本:
lettuce
m0<-drm(weight~conc, data = lettuce, fct = LL.3())
summary(m0)
modelFit(m0)
plot(m0)
我想知道如何在轴 x
中绘制 conc
的确切值及其在 y
中的相应值以及如何绘制 (r²
).
我正在尽我最大的努力解决这个问题,但假设这个问题可能是由于用户界面的困难导致的,尽管我在评论中做出了努力,但未说明的用户界面和术语尚未澄清.我假设这是用 R^2 值注释 drc 图的请求。 plot.drc 函数是一个基本的图形函数,因此默认情况下会绘制到交互式设备上,并且在我的机器上似乎是这样做的(Mac),它不会总是带来 quartz()-device window 到查看层的前面。需要使用下拉菜单将其拉到最前面。
我(现在)将在代码中演示如何 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()
脚本:
lettuce
m0<-drm(weight~conc, data = lettuce, fct = LL.3())
summary(m0)
modelFit(m0)
plot(m0)
我想知道如何在轴 x
中绘制 conc
的确切值及其在 y
中的相应值以及如何绘制 (r²
).
我正在尽我最大的努力解决这个问题,但假设这个问题可能是由于用户界面的困难导致的,尽管我在评论中做出了努力,但未说明的用户界面和术语尚未澄清.我假设这是用 R^2 值注释 drc 图的请求。 plot.drc 函数是一个基本的图形函数,因此默认情况下会绘制到交互式设备上,并且在我的机器上似乎是这样做的(Mac),它不会总是带来 quartz()-device window 到查看层的前面。需要使用下拉菜单将其拉到最前面。
我(现在)将在代码中演示如何 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()