对R中矩阵的每一列的最后n个非NA值求和

# 创建一个示例矩阵
set.seed(123)
mat <- matrix(sample(c(NA, 1:10), 20, replace = TRUE), nrow = 4)

# 定义一个函数来计算每一列的最后n个非NA值的和
sum_last_n_non_na <- function(mat, n) {
  # 获取矩阵的列数
  num_cols <- ncol(mat)
  
  # 初始化结果向量
  result <- numeric(num_cols)
  
  # 遍历每一列
  for (i in 1:num_cols) {
    # 获取当前列的非NA值
    non_na_values <- na.omit(mat[, i])
    
    # 如果非NA值的数量大于等于n，则计算最后n个非NA值的和
    if (length(non_na_values) >= n) {
      result[i] <- sum(tail(non_na_values, n))
    } else {
      # 如果非NA值的数量小于n，则计算所有非NA值的和
      result[i] <- sum(non_na_values)
    }
  }
  
  return(result)
}

# 示例调用
n <- 2
sum_last_n_non_na(mat, n)

if (length(non_na_values) == 0) {
  result[i] <- 0
}

if (length(non_na_values) >= n) {
  result[i] <- sum(tail(non_na_values, n))
} else {
  result[i] <- sum(non_na_values)
}