请求验证R语言客户动态模拟代码的正确性

需求场景 #

• 第1天初始有100个客户，编号为customer1至customer100
• 任意一天，有0.5的概率新增2个客户，0.5的概率流失2个现有客户
• 客户流失后，其编号永不复用
• 新增客户的编号为当前已使用的最大客户编号+1

提交的模拟代码 #

初始化参数 #

set.seed(123)
n <- 100
p_exit <- 0.5
p_join <- 0.5
customers <- sort(paste0("customer", 1:n))
max_customer_number <- n

创建存储对象 #

customers_per_day <- list()

模拟循环代码 #

for (day in 1:100) {

    customers_per_day[[day]] <- customers[order(as.numeric(gsub("customer", "", customers)))]

    if (runif(1) < p_exit) {
        customers[sample(n, 2)] <- NA
    }

    if (runif(1) < p_join) {
        new_customers <- sort(paste0("customer", max_customer_number + 1:max_customer_number + 2))
        customers[is.na(customers)][1:2] <- new_customers
        max_customer_number <- max_customer_number + 2
    }
}
customers_per_day

代码逻辑验证与问题指出 #

你的代码整体思路方向正确，但存在几处不符合需求的逻辑错误：

1. 流失客户的逻辑错误
   你用customers[sample(n, 2)] <- NA处理流失，但n是固定的初始值100，不会随客户变化更新。这会导致你从固定的100个位置中随机选2个，可能选中已经是NA的无效位置，实际没有流失有效客户，不符合“流失2个现有客户”的要求。
   正确做法是从当前非NA的有效客户位置中随机选2个，同时加判断避免有效客户不足2个时报错：

   if (runif(1) < p_exit) {
       valid_indices <- which(!is.na(customers))
       if (length(valid_indices) >= 2) {
           customers[sample(valid_indices, 2)] <- NA
       }
   }
   

2. 新增客户的编号生成错误
   你写的new_customers <- sort(paste0("customer", max_customer_number + 1:max_customer_number + 2))完全不符合需求。R中运算符优先级是先计算1:max_customer_number，比如当max_customer_number=100时，会生成1到100的序列，再加100和2后变成103到202，这显然不是“当前最大编号+1、+2”的逻辑。
   正确写法应为：

   new_customers <- paste0("customer", max_customer_number + 1:2)
   

   这里sort可以省略，因为1:2是递增序列，生成的编号自然有序。

3. 新增客户的赋值逻辑隐患
   你用customers[is.na(customers)][1:2] <- new_customers填充新增客户，但如果当前customers中NA位置不足2个（比如所有客户都是有效的），这行代码会报错，因为目标向量长度小于2，无法完成赋值。
   更合理的方式是直接维护有效客户向量，不用NA占位：

   # 初始化改为仅存有效客户
   customers <- paste0("customer", 1:n)
   # 流失时直接移除对应客户
   if (runif(1) < p_exit && length(customers)>=2) {
       customers <- customers[-sample(length(customers), 2)]
   }
   # 新增时直接追加
   if (runif(1) < p_join) {
       new_customers <- paste0("customer", max_customer_number + 1:2)
       customers <- c(customers, new_customers)
       max_customer_number <- max_customer_number + 2
   }
   

   这种方式更直观，也避免了NA占位带来的各种问题。

4. 每日客户存储的顺序逻辑
   你用customers[order(as.numeric(gsub("customer", "", customers)))]排序存储的逻辑是正确的，若换用有效客户向量的方式，每次存储前执行排序即可。

修正后核心逻辑符合需求 #

修正上述错误后，代码可满足所有需求：初始100个客户、流失客户编号永不复用、每日按概率执行新增/流失操作、新增客户编号基于已使用的最大编号递增。

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