R：创建添加水年列的函数

问题描述

我试图在我的数据框中添加一列，列出每个样本的水文年份。水年定义为上一日历年的 10 月 1 日至 9 月 31 日。例如，WY 2014 从 10/1/2013 开始到 9/31/2014 结束。

我在这里找到了一个类似的问题，并给出了以下代码作为解决方案：

wtr_yr <- function(dates,start_month=9) {
  # Convert dates into POSIXlt
  dates.posix = as.POSIXlt(dates)
  # Year offset
  offset = ifelse(dates.posix$mon >= start_month - 1,1,0)
  # Water year
  adj.year = dates.posix$year + 1900 + offset
  # Return the water year
  adj.year
}

当我尝试将其应用于我自己的代码时，我不断收到错误消息：

rep(value,length.out = nrows) 中的错误：尝试复制“闭包”类型的对象

我将包括我的数据框示例以及我对上面给出的解决方案的改编版本。非常感谢任何帮助！

NN_Loads_Calculation <- structure(list(Year = c(2008,2008,2008),Month = c(1,1),Date.x = structure(c(13879,13879,13879),class = "Date"),Day.x = c(1L,1L,1L),SAC = c(26912331388.8,26912331388.8,26912331388.8),YOLO = c(183493168.56,183493168.56,183493168.56),CSMR = c(146794534.848,146794534.848,146794534.848),MOKE = c(433043877.8016,433043877.8016,433043877.8016),MISC = c(210405499.9488,210405499.9488,210405499.9488),SJR = c(3425205813.12,3425205813.12,3425205813.12),CCC = c(107649325.5552,107649325.5552,107649325.5552),SWP = c(3393400330.5696,3393400330.5696,3393400330.5696),CVP = c(2544438604.032,2544438604.032,2544438604.032),WEST = c(2532205726.128,2532205726.128,2532205726.128),RIO = c(19797689599.8336,19797689599.8336,19797689599.8336),CD = c(3038646871.3536,3038646871.3536,3038646871.3536),Day.y = c(31,31,31),Drainage = c(4799740774.51718,4799740774.51718,4799740774.51718),diversions = c(130231217.046031,130231217.046031,130231217.046031),C10A = c(2.36,2.36,2.36),C3A = c(0.46,0.46,0.46),D12 = c(NA_real_,NA_real_,NA_real_),D19 = c(0.91,0.91,0.91),D22 = c(NA_real_,D26 = c(0.98,0.98,0.98),D28A = c(1,D4 = c(0.72,0.72,0.72),MD10A = c(1.95,1.95,1.95),P8 = c(2.3,2.3,2.3),Island_Div = c(1.32,1.32,1.32),Island_Drain = c(1.94033333333333,1.94033333333333,1.94033333333333),Date.x.x = structure(c(1199145600,1199145600,1199145600),class = c("POSIXct","POSIXt"),tzone = "UTC"),Day.x.x = c(1,Stockton_Flow = c(50989496.946,50989496.946,50989496.946),Stockton_NN = c(18.6241666666667,18.6241666666667,18.6241666666667),Date.y = structure(c(1199145600,1199232000,1199318400,1199404800),Day.y.y = c(1,2,3,4),Tracy_Flow = c(35847849.746,33992997.964,33538748.548,34485101.498),Tracy_NN = c(4.9,4.9,4.9),Banks_NN = c(1.71,1.71,1.71),Jones_NN = c(1.71,SAC_kgd = c(12379.672438848,12379.672438848,12379.672438848),YOLO_kgd = c(84.4068575376,84.4068575376,84.4068575376),SACYOLO_kgd = c(12464.0792963856,12464.0792963856,12464.0792963856),CSMR_kgd = c(286.2493429536,286.2493429536,286.2493429536),MOKE_kgd = c(844.43556171312,844.43556171312,844.43556171312),MISC_kgd = c(410.29072490016,410.29072490016,410.29072490016),EAST_kgd = c(1540.97562956688,1540.97562956688,1540.97562956688),SJR_kgd = c(8083.4857189632,8083.4857189632,8083.4857189632),STOCKTON_kgd = c(949.636889371797,949.636889371797,949.636889371797),TRACY_kgd = c(175.6544637554,166.5656900236,164.3398678852,168.9769973402),Islan_Drain_kgd = c(9313.09701615484,9313.09701615484,9313.09701615484),Total_Input_kgd = c(32526.9290141977,32517.8402404659,32515.6144183275,32520.2515477825),WEST_kgd = c(1823.18812281216,1823.18812281216,1823.18812281216),RIO_kgd = c(14254.3365118802,14254.3365118802,14254.3365118802),WESTRIO_kgd = c(16077.5246346923,16077.5246346923,16077.5246346923),CCC_kgd = c(107.6493255552,107.6493255552,107.6493255552),SWP_kgd = c(5802.71456527402,5802.71456527402,5802.71456527402),CVP_kgd = c(4350.99001289472,4350.99001289472,4350.99001289472),Water_Exp_kg = c(10261.3539037239,10261.3539037239,10261.3539037239),Island_Div_kgd = c(171.90520650076,171.90520650076,171.90520650076),Total_Output_kgd = c(26510.783744917,26510.783744917,26510.783744917),Delta_Uptake_kgd = c(6016.14526928067,6007.05649554887,6004.83067341047,6009.46780286547),Percent_Uptake = c(18.4958907945311,18.4731103023059,18.4675294649386,18.4791553473554)),row.names = c(NA,4L),class = "data.frame")

这是我试图用来添加水年列的代码：

NN_Loads_Calculation$wtr_yr <- function(Date.x,start_month=9) {
  # Convert dates into POSIXlt
  dates.posix = as.POSIXlt(Date.x)
  # Year offset
  offset = ifelse(dates.posix$mon >= start_month - 1,0)
  # Water year
  adj.year = dates.posix$year + 1900 + offset
  # Return the water year
  adj.year
}

解决方法

您正在将函数（闭包）写入 NN_Loads_Calculation$wtr_yr，而不是函数的矢量化输出。

因为函数的长度为 1（即只有一个函数），R 尝试重复它，使其成为与数据框中的列长度相同的函数向量。您不能在 rep() 上运行 closure，所以这就是代码崩溃的地方。

您需要定义您的函数，然后使用它来计算新值：

## Define the function
get_water_year <- function(date_x = NULL,start_month = 9) {
  # Convert dates into POSIXlt
  dates.posix = as.POSIXlt(date_x)
  # Year offset
  offset = ifelse(dates.posix$mon >= start_month - 1,1,0)
  # Water year
  adj.year = dates.posix$year + 1900 + offset
  # Return the water year
  adj.year
}

### Write the function output to the df,rather than the function.
NN_Loads_Calculation$wtr_yr <- get_water_year(
date_x = NN_Loads_Calculations$Date.x,start_month = 9)

posixlt r r