在for循环R中从`curve（）`保存x和y？答案

【问题标题】：Saving x, and y from `curve()` in a for loop R?在for循环R中从`curve（）`保存x和y？
【发布时间】：2017-12-07 16:59:28
【问题描述】：

我有一个名为 all.priors 的函数（请参阅下面的 R 代码）。我的目标是从for 循环内的curve() 调用中获取x 和y，并将这些xs 和ys 保存为对象h。

（我想在h 中有101 行和2*length(d) 列。这样，每2 列包含来自curve() 的x 和y 在for 循环中运行.)

问题：

如何正确保存来自curve() 调用的xs 和ys？ [我收到错误：incorrect number of subscripts on matrix]

all.priors = function(a, b, lo, hi, d, Bi = 55, n = 1e2){

h = matrix(NA, 101, 2*length(d)) 

for(i in 1:length(d)){
         p = function(x) get(d[i])(x, a, b)
     prior = function(x) p(x)/integrate(p, lo, hi)[[1]]
likelihood = function(x) dbinom(Bi, n, x)
 posterior = function(x) prior(x)*likelihood(x)
     h[i,] = curve(posterior, ty = "n", ann = FALSE, yaxt = "n", xaxt = "n", add = i!= 1, bty = "n")
     }
}
#Example of use:
all.priors(lo = 0, hi = 1, a = 2, b = 3, d = c("dgamma", "dnorm", "dcauchy", "dlogis"))

【问题讨论】：

我不明白你想要的输出是什么。在h 的每一行和每一列中你想要什么？ curve() 返回一个包含组件x 和y 的列表。您不能直接将列表存储在矩阵中。
@MrFlick，我要保存x和y的每一组，然后分别绘制每组？
这并不能更好地解释事情。您希望从该函数返回的数据的结构究竟如何？
@MrFlick，我明白了，所以我想在h 中有101 行和2*length(d) 列？这样，每 2 列包含来自 curve() 的 x 和 y 在 for 循环中运行。
@MrFlick，行数和列数有意义吗？您认为以这种方式存储x 和y 是否适合单独绘制它们？

标签： r function for-loop curve

【解决方案1】：

您只需要小心地将值放入矩阵中，然后从您的函数中返回矩阵。试试这个

all.priors = function(a, b, lo, hi, d, Bi = 55, n = 1e2){

  h = matrix(NA, 101, 2*length(d)) 

  for(i in 1:length(d)){
    p = function(x) get(d[i])(x, a, b)
    prior = function(x) p(x)/integrate(p, lo, hi)[[1]]
    likelihood = function(x) dbinom(Bi, n, x)
    posterior = function(x) prior(x)*likelihood(x)
    cv <- curve(posterior, ty = "n", ann = FALSE, yaxt = "n", xaxt = "n", add = i!= 1, bty = "n")
    h[,i*2-1] <- cv$x
    h[,i*2] <- cv$y
  }
  h
}
all.priors(lo = 0, hi = 1, a = 2, b = 3, d = c("dgamma", "dnorm", "dcauchy", "dlogis"))

【讨论】：

【解决方案2】：

解决此问题的另一种方法可能是将答案保存在列表而不是矩阵中。我认为您的函数使正在发生的事情变得复杂，因此我将使用一个更简单的示例。

h = list()
for(i in 1:5) {
    h[i] = list(curve(sin(i*x), xlim=c(0,6.3))) }

生成的数据结构应该易于使用。

【讨论】：

@morouzian 是的。你想要一个包含所有曲线的图还是多个图？