R语言嵌套循环迭代值丢失问题：多替换长度下的模型测试

问题原因 #

当前代码的核心问题出在数据绑定的逻辑上：

• 当i==1时，每次j循环都会执行ALL_modData = modData，直接用当前j对应的单行数据覆盖之前的结果。比如i=1的循环中，j=1生成的第一条数据会被j=2的覆盖，直到j=4循环结束，ALL_modData里只剩i=1、j=4的那一条数据，而非预期的4条。
• 后续i=2、i=3的循环用rbind追加数据，最终结果只有1+4+4=9条，而非12条。

另外原代码中nls模型的初始值start = list(b1 = 0)可能引发收敛问题：初始值为0时，模型预测值全为0，拟合过程容易报错或无法得到有效结果。

修复方法 #

1. 提前初始化一个空数据框，指定列名和数据类型，避免rbind自动转换类型导致的问题。
2. 去掉i==1的判断分支，每次循环都用rbind追加数据。
3. 调整nls模型的初始值为合理值（比如1），确保模型能正常收敛。

修正后的完整代码 #

data("iris")
head(iris)

Species <- as.character(unique(iris$Species))
expval <- c(0.25,0.50,0.75,0.80)

# 初始化空数据框，明确列类型
ALL_modData <- data.frame(
  Species = character(),
  expval = numeric(),
  est = numeric(),
  stderr = numeric(),
  tval = numeric(),
  pval = numeric(),
  stringsAsFactors = FALSE
)

for (i in 1:length(Species)) {
  for(j in 1:length(expval)) {
    
    currdata <- iris[iris$Species == Species[i],] 
    modData <- data.frame(Species = Species[i], stringsAsFactors = FALSE)
    
    # 调整初始值为1，避免收敛失败
    mod <- nls(Sepal.Width ~ b1*Petal.Width^expval[j],
               data = currdata,
               start = list(b1 = 1)) 
    
    modsum <- summary(mod)
    
    modData$expval <- expval[j]
    modData$est <- modsum[["parameters"]][1] 
    modData$stderr <- modsum[["parameters"]][2] 
    modData$tval <- modsum[["parameters"]][3] 
    modData$pval <- modsum[["parameters"]][4] 
    
    # 每次循环追加数据，不再覆盖
    ALL_modData <- rbind(ALL_modData, modData)
  }
}

# 验证结果行数（应为12）
nrow(ALL_modData)

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