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")) 
Comments
Post a Comment