如何使用 R 的 Luminescence 包中的 plot_KDE 存储模拟结果并将所有结果绘制在一个图中

Question

我正在使用随机数模拟创建模拟。这给出了 100 组 45 个值的错误。

首先，我想存储这些模拟的结果。 然后，我需要将这些模拟的结果绘制在一张图上。我需要制作的图使用 Luminescence 包，属于 KDE 类型。 我已经设法生成了单独的实体，但在存储结果和生成带有所有模拟的图时都在努力。 到目前为止，我已经创建了模拟：

Simulation <- function() {
    RNC <- rescale (SFMT(45, dim=1, mexp=216091, 
                            usepset=T, withtorus= F, usetime=T), 
                    c(0.01,130))
    RNC_error <- RNC*0.15
    df <-data.frame(RNC,RNC_error)
}

我要创建的情节使用以下内容：

library("Luminescence")
plot_KDE(data=df, na.rm = TRUE,
         values.cumulative = TRUE, order = TRUE, 
         boxplot = F, rug = F,
         summary.method = "MCM", bw = "nrd0", 
         output = TRUE)

对于我的最终结果，我需要存储所有模拟的数值结果以及包含所有模拟结果的单个 KDE 图。

Answer 1

把你的问题分成两部分。

1. 存储结果。你有一个数据框，df，所以只需使用write.csv() 将结果存储到 CSV 文件，即

   write.csv(df, file="some_file.csv")
   

2. 存储你的情节。显然你不能使用 csv 文件，所以我们将使用 pdf 或 png，例如

   # Open the file
   pdf("figure_file.pdf")
   plot_KDE(data=df, na.rm = TRUE,
       values.cumulative = TRUE, order = TRUE, 
       boxplot = F, rug = F,
       summary.method = "MCM", bw = "nrd0", 
       output = TRUE)
   # Close the file
   dev.off()
   

   要保存为png，请使用png() 而不是pdf()