如何将控制台的ED函数输出添加到DRM绘图

如何将控制台的ED函数输出添加到DRM绘图,r,drc,R,Drc,在DRC包中,我试图将运行ED函数得到的输出添加到一个文本框中,该文本框已放置在DRM绘图中。关于如何编写脚本以便在文本函数中包含ED函数并自动将结果放置到绘图上,有什么想法吗 以下是我迄今为止尝试过的代码: dayone <- structure(list(hours = c(24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24

在DRC包中,我试图将运行ED函数得到的输出添加到一个文本框中,该文本框已放置在DRM绘图中。关于如何编写脚本以便在文本函数中包含ED函数并自动将结果放置到绘图上,有什么想法吗

以下是我迄今为止尝试过的代码:

dayone <- structure(list(hours = c(24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 
    24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 
    24L, 24L, 24L, 24L, 24L, 24L, 24L), jarnum = c(NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_, NA_character_, NA_character_, NA_character_, 
    NA_character_, NA_character_), treat = c(0L, 0L, 0L, 0L, 1500L, 
    1500L, 1500L, 1500L, 3000L, 3000L, 3000L, 3000L, 4000L, 4000L, 
    4000L, 4000L, 5000L, 5000L, 5000L, 5000L, 6000L, 6000L, 6000L, 
    6000L, 10000L, 10000L, 10000L, 10000L), rep = c(1L, 2L, 3L, 4L, 
    1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 
    1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L), total = c(5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L, 5L), mort = c(0L, 0L, 0L, 0L, 0L, 0L, 
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 2L, 1L, 1L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L), cummort = c(0L, 0L, 0L, 0L, 0L, 0L, 
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 2L, 1L, 1L, 5L, 5L, 5L, 5L, 5L, 5L, 
    5L, 5L, 5L, 5L, 5L, 5L)), .Names = c("hours", "jarnum", "treat", 
    "rep", "total", "mort", "cummort"), row.names = c(NA, -28L), class = c("tbl_df", 
    "tbl", "data.frame"))

drmdayone <- drm(cummort/total~treat, data=dayone, fct=LL.2(), type="binomial")
plot(drmdayone, type="all", pch=16, lwd=2,log="x", main="Acute DRC - 24 Hours",xlab="[Cl-] (mg/L)", ylab="Mortality", axes=F) 
axis(1, c(1500,3000,4000,5000,6000,10000))
axis(2, c(0,0.1,0.5,1))
text(500, 0.6, "Here is where I want the ED function (below) results to be displayed", font=2)

ED(drmdayone, c(10,50,90))

dayone不确定您希望它看起来如何,但是您可以捕获函数的文本输出并直接绘制它

text(500, .9, paste(capture.output(ED(drmdayone, c(10,50,90))), collapse="\n"))

这将更容易帮助您获得输入数据示例、您迄今为止尝试过的基本代码以及所需输出的清晰描述。#######
第一天我尝试过的代码以及我正在使用的数据。很抱歉,注释中的代码格式不正确,我是stackoverflow的新手,还不知道如何使用它。您应该编辑问题以包含此信息,以便正确格式化。还可以看到我提供的关于如何在问题本身中包含数据的链接,而不是感谢外部链接。这正是我想要的。非常感谢。