如何从 R 输出绘制直方图？

Question

我试图从不是数据框的 R 输出绘制直方图。 下面是我的代码和输出。

x <- replicate(1000, 
               {y <- rpois(200, 1)
               {lambda0 <- 1
for(i in 1:1) 
  {
  if( i == 1 ) cat( sprintf("%15s %15s %15s %15s\n", "LogL", "Score", "Information", "New Estimate"))
  logL        <- sum((-lambda0) + y*(log(lambda0)))
  score       <- sum((y/lambda0)-1)
  information <- sum(y/(lambda0)^2)
  lambda1      <- lambda0 + score/information
  cat( sprintf("%15.4f %15.4f %15.4f %15.5f\n", logL, score, information, lambda1))
  lambda0 <- lambda1
}  
} 
               })

下面是我的输出

我正在尝试从输出中获取新的估计值并创建直方图。 你能帮忙吗？

谢谢你。

Answer 1

您需要在循环期间存储 New Estimate 的值。 这样您就可以在循环完成后检索结果。 通常在使用循环时，您预先指定一个变量，您可以在其中保存每次迭代的结果。 例如：

numbers <- 1:3

result <- list(length = length(numbers)

for (i in seq_along(numbers){

result[[i]] <- numbers[[i]] + 1

}

在这个例子中，有一个由三个数字组成的向量，你想给每个数字加一个并保存结果。 您可以通过创建长度为 3 的列表（添加长度更好，但不是必需的）来完成此操作，并且对于每个第i次迭代，您将结果保存在列表的第i个元素中。

完成循环后，您可以从result变量中检索结果。 您可以使用方括号检索第i个结果： result[[i]] 。