2017-08-18 7 views
0

DRCパッケージ内で、私は自分のDRMプロットに配置したテキストボックスにED関数を実行して出力を追加しようとしています。どのように私はテキスト関数に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)) 
+1

(https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible [再現可能な例]のお手伝いをするために容易になるだろう-example)にサンプル入力データ、これまでに試した基本コード、および望ましい出力の明確な説明が含まれています。 – MrFlick

+0

##これまでに試したコードは次のとおりです。## 'dayone < - subset(acute、hours == 24)drmdayone < - drm(cummort/total〜treat、data = dayone、fct = LL。 (DRM - 24時間)、xlab = "[Cl(2)]、" type = "binomial" (1、c(1500,3000,4000,5000,6000,10000))軸(2、c(0,0.1,0.5、 - )(mg/L) "、ylab ="死亡率 "、軸= F) 1)) ED(drmdayone、c(10,50,90)) ' –

+0

ここにEDの機能が表示されます。私は[リンク]を使用しているデータ(https://docs.google.com/spreadsheets/d/1S6TDgwP8wW9Fr7ppyMloCec664xJZ4EXLoM9KGI3saU/edit#gid=0) –

答えて

1

わからない:ここで

は、私がこれまで試したコードです。

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

enter image description here

+0

これはまさに私が探していたものです。ありがとうございました! –

関連する問題