R语言使用ifelse函数初始化向量时变量异常归零问题咨询

问题核心原因 #

• 错误使用ifelse做流程控制：ifelse是面向向量化运算的函数，返回值长度和判断条件长度完全一致，不能替代if(){...}else{...}做分支逻辑控制。你当前写法中判断条件i==0长度为1，执行w<-x0时只会给w的第一个元素赋值为3，第二个元素保留初始化的0值，此时w为(3,0)；后续循环i不等于0时，要把w赋值为初始值为(0,0)的x，同样因为ifelse的特性，仅给w的第一个元素赋值为0，此时w就变成了(0,0)，后续生成y时均值为0，自然会出现你观察到的变量全为0的情况。
• 依赖对象未明确传入：代码中S1、h(x)、g(x)三个对象没有在函数内定义，也没有作为参数传入，若全局环境无对应定义会直接报错，即使有定义也可能因为作用域问题拿到非预期值。
• 索引逻辑错误：alfa每次循环都被赋值为单个标量，不存在alfa[n]的索引值，运行到return语句会直接报错。

修复后代码示例 #

# 先定义你的实际计算逻辑，此处为示例可替换
h <- function(x) sum(dnorm(x, mean = 0, sd = 1))
g <- function(x) sum(dnorm(x, mean = 5, sd = 2))
S1 <- 1

simulation <- function(n, h, g, S1){
  f <- function(x){0.7*h(x) + 0.3*g(x)}
  # 直接初始化x为初始值，省去首次循环特殊判断
  x <- c(3, 3)
  # 预先初始化alfa存储向量，避免循环动态扩容
  alfa <- numeric(n+1)
  for(i in 1:(n+1)){
    y <- rnorm(2, mean = x, sd = S1)
    alfa[i] <- f(y)/f(x)
    if(runif(1) < alfa[i]){
      x <- y
    }
  }
  return(alfa[n])
}

# 测试调用
simulation(10, h, g, S1)

向量初始化的可选方式 #

• 直接指定初始值：x <- c(3,3)，适合已知固定初始值的场景，无需后续分支判断调整
• 用rep()生成固定值向量：如果需要初始为其他统一值可以用x <- rep(0, 2)，和vector("numeric",2)效果一致但可读性更强
• 参数化初始化：如果初始值需要统一调整可以写成x <- rep(init_val, dim)，后续只需修改参数即可批量调整初始值

内容的提问来源于stack exchange，提问作者Sofía Contreras