R语言 数据为2001-01-01到2099-12-31的全部日期，格式为YYYY-MM-DD, 使用substr函数提取“年”的后两位，计算“月*日 = 年”的概率

在R语言中，如果你有一个从2001年1月1日到2099年12月31日的日期序列，可以首先将日期转换成"YYYYMMDD"格式，然后利用`substr()`函数提取出年份的最后两位。接着，你需要创建一个新的变量，计算每个月有多少天，并将其乘以该月的日期，形成"月*日"的结果。为了得到"年"的概率，你可以将这个结果除以总天数（365 * (2099 - 2001 + 1)），因为一般我们假设每年有365天。

以下是实现这一过程的一个步骤：

# 首先加载必要的包（如 lubridate 或者 base R 的 as.Date）
library(lubridate) # 如果未安装，运行 install.packages("lubridate")

# 假设你的日期数据存储在一个叫做'dates'的向量里
dates <- seq(from = "2001-01-01", to = "2099-12-31", by = "day")

# 将日期转为 "YYYYMMDD" 格式
formatted_dates <- format(dates, "%Y%m%d")

# 提取年份后两位
years_two_digits <- substr(formatted_dates, nchar(formatted_dates), nchar(formatted_dates))

# 创建一个包含每个月天数的向量
days_in_month <- c(31, rep(28, 7), rep(31, 4), rep(30, 4)) # 考虑闰年规则

# 计算每个组合"月*日"
monthly_combinations <- days_in_month * 100 # 因为我们要提取的是年份的最后两位，所以扩大100倍

# 现在，我们需要一个映射表，将两个数字组合对应到相应的年份
year_map <- tapply(monthly_combinations, years_two_digits, sum)

# 计算每个组合出现的概率
probability <- year_map / sum(year_map)

# 检查概率分布
probability