我有一个假设公司过去五年的季度股息支付历史。这是可重现的代码：
Date<-as.Date(c("2013-11-01", "2014-02-01", "2014-05-01", "2014-08-01", "2014-11-01", "2015-02-01", "2015-05-01", "2015-08-01", "2015-11-01", "2016-02-01", "2016-05-01", "2016-08-01", "2016-11-01", "2017-02-01", "2017-05-01", "2017-08-01", "2017-11-01", "2018-02-01", "2018-05-01", "2018-08-01"))
Dividend<-c(0.08, 0.10, 0.10, 0.10, 0.10, 0.11, 0.00, 0.11, 0.11, 0.13, 0.13, 0.13, 0.13, 0.14, 0.14, 0.00, 0.16, 0.15, 0.15, 0.15)
data.frame(Date,Dividend)
具有以下输出：
Date Dividend
1 2013-11-01 0.08
2 2014-02-01 0.10
3 2014-05-01 0.10
4 2014-08-01 0.10
5 2014-11-01 0.10
6 2015-02-01 0.11
7 2015-05-01 0.00
8 2015-08-01 0.11
9 2015-11-01 0.11
10 2016-02-01 0.13
11 2016-05-01 0.13
12 2016-08-01 0.13
13 2016-11-01 0.13
14 2017-02-01 0.14
15 2017-05-01 0.14
16 2017-08-01 0.00
17 2017-11-01 0.16
18 2018-02-01 0.15
19 2018-05-01 0.15
20 2018-08-01 0.15
我的问题是如何将其转换为显示每个完整年份支付的股息的产出，如果它们不完整(如本例，2013 年和 2018 年(，则忽略第一年和最后一年，并且不假设年度股息总是等于季度股息 x 4(在我的示例中，2015 年和 2017 年没有相同数量的季度股息(。
所以输出看起来像这样：
Date Dividend
2014 0.40
2015 0.33
2016 0.52
2017 0.44
你如何定义不完整？您要么知道该公司每年有 4 次付款，要么有 2 次或 1 次(或 12 次(。根据您的推理，其他答案都不是正确的，因为它们只是假设不应该考虑第一年和最后一年，但是当第四次付款时，2018年11月会发生什么？
由于您使用的是 quantmod，因此数据应采用 xts 形式。使用apply.yearly将使用数据集中每个年份的最后可用日期将数据汇总到年度行。我使用函数的 FUN 部分返回 2 列，1 列包含汇总的股息，1 列包含当年支付的股息数量。由于公司倾向于有结构化的股息支付(1、2、4 或 12(，您可以使用最大div 支付数量来过滤掉不符合此规定的年份。
当您有每年未支付的特殊股息之类的东西时，这可能会失败。此外，支付股息的开始并不总是符合规则。Microsoft 于 2013 年开始支付股息，在季度支付之后，只有 2 次支付。
df1 <- data.frame(Date = as.Date(c("2013-11-01", "2014-02-01", "2014-05-01", "2014-08-01", "2014-11-01", "2015-02-01", "2015-05-01", "2015-08-01", "2015-11-01", "2016-02-01", "2016-05-01", "2016-08-01", "2016-11-01", "2017-02-01", "2017-05-01", "2017-08-01", "2017-11-01", "2018-02-01", "2018-05-01", "2018-08-01")),
Dividend = c(0.08, 0.10, 0.10, 0.10, 0.10, 0.11, 0.00, 0.11, 0.11, 0.13, 0.13, 0.13, 0.13, 0.14, 0.14, 0.00, 0.16, 0.15, 0.15, 0.15))
# data in xts form since quantmod is being used.
my_xts <- xts(df1$Dividend, order.by = df1$Date)
annual_data <- apply.yearly(my_xts, function(x) as.matrix(data.frame(sum(x), length(x))) )
names(annual_data) <- c("total_divs", "no_divs")
# filter data to include only maximum dividens
annual_data[annual_data$no_divs == max(annual_data$no_divs)]
total_divs no_divs
2014-11-01 0.40 4
2015-11-01 0.33 4
2016-11-01 0.52 4
2017-11-01 0.44 4
我们可以从数据帧中删除最大和最小year，然后按year分组并进行求和。
library(dplyr)
library(lubridate)
df %>%
filter(year(Date) != min(year(Date)) & year(Date) != max(year(Date))) %>%
group_by(year = year(Date)) %>%
summarise(Dividend = sum(Dividend))
# year Dividend
# <dbl> <dbl>
#1 2014 0.4
#2 2015 0.33
#3 2016 0.52
#4 2017 0.44
它的基本 R 等效物将是
df$Year <- as.numeric(format(df$Date, "%Y"))
aggregate(Dividend~Year, df[with(df, Year != min(Year) & Year != max(Year)), ],sum)
# Year Dividend
#1 2014 0.40
#2 2015 0.33
#3 2016 0.52
#4 2017 0.44
有了data.table你可以试试
df<-data.frame(Date,Dividend)
library(data.table)
setDT(df)[,.(TotDiv=sum(Dividend)),
by=year(Date), ][-c(which.min(year),which.max(year))]
year TotDiv
1: 2014 0.40
2: 2015 0.33
3: 2016 0.52
4: 2017 0.44