如何创建R函数遍历指定列行，生成基于特定公式的新列？

我明白你的需求啦——就是要根据ctV1和ctV2列的数值，用给定的公式计算出对应的拷贝数，分别存入samples2to9和samples11to18这两个新列。先帮你修正一下原来代码里的语法问题，再给你更高效的实现方式：

第一步：加载你的数据集 #

先把你提供的数据加载到R环境里：

dfqpcr <- structure(list(X = c(4000, 40000, 4e+05, 4e+06, 4e+07, 4e+08, NA, NA), 
                         X.1 = c(3.602059991, 4.602059991, 5.602059991, 6.602059991, 7.602059991, 8.602059991, NA, NA), 
                         X.2 = c(24L, 48L, 72L, 96L, 24L, 48L, 72L, 96L), 
                         ctV1 = c(22.92, 24.29, 24.29, 27.91, 21.44, 16.5, 17.79, 17.34), 
                         ctV2 = c(24.21, 22.16, 22.16, 23.76, 22.81, 18.31, 16.56, 19.09), 
                         ctSTD = c(33L, 28L, 24L, 20L, 15L, 12L, NA, NA), 
                         X.3 = c(NA, NA, NA, NA, NA, NA, NA, NA)), 
                    row.names = c(NA, 8L), class = "data.frame")

第二步：定义正确的RNAcopy函数 #

你原来的函数写法有语法错误，而且完全没必要手动循环——R是向量型语言，直接对整列（向量）做运算效率更高。下面是修正后的函数：

RNAcopy <- function(ct_values, na.rm = TRUE) {
  # 直接用向量运算计算拷贝数，自动处理每一行
  copy_numbers <- 10^((ct_values - 45.61196) / 4.229)
  
  # 如果需要处理NA值（比如替换成0），可以取消下面的注释
  # if (na.rm) {
  #   copy_numbers[is.na(copy_numbers)] <- 0
  # }
  
  return(copy_numbers)
}

这里的逻辑很简单：传入一个ct值的向量，用你给定的公式计算每个值对应的拷贝数，默认保留NA（因为原数据的NA计算后也应该是NA，更符合数据逻辑）。

第三步：生成新列 #

直接调用函数处理指定列，赋值给新列即可：

# 基于ctV1生成samples2to9
dfqpcr$samples2to9 <- RNAcopy(dfqpcr$ctV1)
# 基于ctV2生成samples11to18
dfqpcr$samples11to18 <- RNAcopy(dfqpcr$ctV2)

查看结果 #

运行完上面的代码后，你可以用print(dfqpcr)查看处理后的数据集，就能看到新增的两列计算结果啦。

为什么原来的代码有问题？ #

• 函数定义语法错误：你写的RNAcopy<- function(for i in x {10^((i - 45.61196)/(4.229))})不符合R的函数定义规则，正确的写法是function(参数) {函数体}。
• 没必要手动循环：R的向量运算会自动遍历每个元素，比手动写for循环简洁高效得多。
• na.rm参数的处理：原代码里加了这个参数但没实际用到，上面的函数里我添加了可选的NA处理逻辑，你可以根据需求调整。

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