R语言使用runif、rnorm结合初始值计算需求与市场份额增长率

原代码错误原因 #

你之前编写的MarketSizeGrowth<-cumprod(c(runif(1000,1.d,1.e),rep(1.f,2))存在维度匹配错误：代码将1000个第一年市场份额样本和2个年度增长率拼接为长度1002的向量，cumprod会对整个序列做累积乘积，无法实现每个模拟样本单独计算3年增长的需求。

实现思路 #

• 每次模拟为独立样本，分别计算对应3年的需求、市场份额数值
• 固定增长率下，后一年数值 = 前一年数值 * 增长乘数，无需重复抽样初始分布
• 最终输出矩阵格式结果：每一行对应1次模拟，每一列对应1个年份，方便后续统计分析

可运行代码示例 #

基础循环写法（可读性更高） #

# 可根据实际业务值替换以下参数
n_sim <- 1000       # 模拟次数
a <- 800            # 第一年需求均值
b <- 150            # 第一年需求标准差
c_pct <- 6          # 需求年增长率c%，示例为6%
d_pct <- 15         # 第一年市场份额下限d%，示例为15%
e_pct <- 25         # 第一年市场份额上限e%，示例为25%
f_pct <- 7          # 市场份额年增长率f%，示例为7%
n_year <- 3         # 计算周期为3年

# 转换增长乘数
demand_grow <- 1 + c_pct / 100
share_grow <- 1 + f_pct / 100

# 计算3年需求
demand <- matrix(NA, nrow = n_sim, ncol = n_year)
demand[, 1] <- rnorm(n_sim, mean = a, sd = b)
for (y in 2:n_year) {
  demand[, y] <- demand[, y-1] * demand_grow
}

# 计算3年市场份额乘数（你提到的加1格式，如需实际份额减1即可）
share_multi <- matrix(NA, nrow = n_sim, ncol = n_year)
share_multi[, 1] <- runif(n_sim, min = 1 + d_pct/100, max = 1 + e_pct/100)
for (y in 2:n_year) {
  share_multi[, y] <- share_multi[, y-1] * share_grow
}

向量化写法（运行效率更高） #

如果不需要保留中间过程，可直接用replicate + cumprod实现同等效果：

# 3年需求增长乘数序列
demand_seq <- cumprod(c(1, rep(demand_grow, n_year - 1)))
demand <- t(replicate(n_sim, rnorm(1, a, b) * demand_seq))

# 3年市场份额增长乘数序列
share_seq <- cumprod(c(1, rep(share_grow, n_year - 1)))
share_multi <- t(replicate(n_sim, runif(1, 1 + d_pct/100, 1 + e_pct/100) * share_seq))

结果说明 #

• 输出的demand和share_multi均为1000行3列的矩阵，第1/2/3列分别对应第1/2/3年的数值
• 如需计算每年实际销售额，直接对应位置相乘即可：annual_sales <- demand * (share_multi - 1)，减1是将加1后的份额乘数还原为实际占比

内容的提问来源于stack exchange，提问作者PythonSOS