將繪圖數據從R導出到Excel

Question

我有這個圖，我從這段代碼生成它：

m <- bcea(e=effects,c=costs, ref=2, interventions=treatments, Kmax=50000)

情節是：

evi.plot(m)

現在，我需要在excel文件中導出這個evpi.plot（m），而不是創建的jpeg，但是數據和它一起，我的意思是創建了X和Y軸。

我一直在使用這樣的東西，但不是這種情況

write.table( thresholds, 'clipboard', sep='\t', row.names=FALSE, col.names=FALSE )

Answer 1

在函數的文檔bcea從包BCEA你可以看到你對象的結構：

值

包含以下元素的類“bcea”的對象

n.sim貝葉斯模型產生的模擬數

n.comparators正在分析的干預措施的數量

...

k網格近似值支付意願的向量

...

evi信息的預期價值的值向量，作為支付意願的函數

如果你看一下evi.plot的函數定義，你會發現你的x和y值是名為k和evi的元素：

> evi.plot
function (he) 
{
    options(scipen = 10)
    plot(he$k, he$evi, t = "l", xlab = "Willingness to pay", 
        ylab = "EVPI", main = "Expected Value of Information")
    if (length(he$kstar) > 0) {
        points(rep(he$kstar, 3), c(-10000, he$evi[he$k == he$kstar]/2, 
            he$evi[he$k == he$kstar]), t = "l", lty = 2, col = "dark grey")
        points(c(-10000, he$kstar/2, he$kstar), rep(he$evi[he$k == 
            he$kstar], 3), t = "l", lty = 2, col = "dark grey")
    }
}
<environment: namespace:BCEA>

所以：

res <- cbind(m$k, m$evi)
write.table(res, file="bcea.csv", sep=',', row.names=FALSE, col.names=FALSE )