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正在分析的干预措施数量

...

ķ   网格近似值支付意愿的值向量

...

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 )

将绘图数据从R导出到Excel

1 个答案: