Question

我收集了一组数据，其中包含一个时间序列，其中每个y值是通过取30个葡萄簇重量样本的平均值得到的。

我希望用相同数量的x和y值来模拟更多数据，这样我就可以进行一些贝叶斯分析来找到数据的后验分布。我有数据，我知道增长遵循Gompertz曲线，公式为：

[y = a*exp(-exp(-(x-x0)/b))]，a = 88.8，b = 11.7和x0 = 15.1。

我的数据是

x = c(0, 28, 36, 42, 50, 58, 63, 71, 79, 85, 92, 99, 106, 112)

y = c(0, 15, 35, 55, 62, 74, 80, 96, 127, 120, 146, 160, 177, 165).

任何帮助将不胜感激

Answer 1

*将在给出更多信息时进行编辑**

我对你的问题感到有些困惑。我已经编写了你写入R的内容。请详细说明我，以便我可以帮助你：

require 'csv'

total_sales = 0

CSV.foreach('sales-data.csv', headers: true, converters: :all) do |row|

  if (row[2]) == "Music"

    # TODO: check if category is "Music" (row[2])

    # TODO: if it is music, add total_sales + the row's sales (row[6])

    puts total_sales + (row[6])

  end

puts total_sales.round(2)

模拟R中Gompertz曲线的数据

1 个答案: