r - How to add ED function output from the console to a DRM plot -


within drc package, trying add output running ed function text box have placed on drm plot. ideas on how write script can include ed function in text function , automatically place results onto plot?

here code i've tried far:

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 want ed function (below) results displayed", font=2)  ed(drmdayone, c(10,50,90)) 

not sure how want look, can capture text output function , draw directly.

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

enter image description here


Comments

Popular posts from this blog

What is happening when Matlab is starting a "parallel pool"? -

angular - DownloadURL return null in below code -

php - Cannot override Laravel Spark authentication with own implementation -