Я пытаюсь написать функцию R, чтобы вычислить количество будних дней между двумя датами. Например, Nweekdays ('01/30/2011 ',' 02/04/2011 ') будет равен 5.
Как и этот вопрос. Спасибо!
/edit: @J. Ответ Винчестера велик, но мне было интересно, может ли кто-нибудь подумать о возможности его векторизации, чтобы он работал на двух столбцах дат. Благодарю! /edit 2: Еще раз спасибо!
Date1 <- as.Date("2011-01-30")
Date2 <- as.Date("2011-02-04")
sum(!weekdays(seq(Date1, Date2, "days")) %in% c("Saturday", "Sunday"))
EDIT: И Зак сказал, пусть будет Vectorize:)
Dates1 <- as.Date("2011-01-30") + rep(0, 10)
Dates2 <- as.Date("2011-02-04") + seq(0, 9)
Nweekdays <- Vectorize(function(a, b)
sum(!weekdays(seq(a, b, "days")) %in% c("Saturday", "Sunday")))
Nweekdays(Dates1, Dates2)
Эти модифицированные функции учитывают различия дат как положительных, так и отрицательных, тогда как принятое решение учитывает положительную разницу дат.
library("dplyr")
e2 <- structure(list(date.pr = structure(c(16524, 16524, 16524, 16524, 16524, 16524, 16524, 16524, 16524, 16524, 16545, 5974), class = "Date"),
date.po = structure(c(16524, 16525, 16526, 16527, 16528, 16529, 16530, 16531, 16538, 16545, 16524, 15974), class = "Date")),
.Names = c("date.1", "date.2"), class = c("tbl_df", "data.frame"), row.names = c(NA, -12L))
1. Локальное зависимое решение: Nweekdays() функция адаптирована из @J. Выиграл. решение. Он работает для locale = "English_United States.1252"
Nweekdays <- Vectorize(
function(a, b)
{
ifelse(a < b,
return(sum(!weekdays(seq(a, b, "days")) %in% c("Saturday", "Sunday")) - 1),
return(sum(!weekdays(seq(b, a, "days")) %in% c("Saturday", "Sunday")) - 1))
})
а. Английский язык
> Sys.setlocale(category="LC_ALL", locale = "English_United States.1252")
[1] "LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252"
> Sys.getlocale()
[1] "LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252"
> e2 %>%
mutate(wkd1 = format(date.1, "%A"),
wkd2 = format(date.2, "%A"),
ndays_with_wkends = ifelse((date.2 > date.1), (date.2 - date.1), (date.1 - date.2)),
ndays_no_wkends = Nweekdays(date.1, date.2))
Source: local data frame [12 x 6]
date.1 date.2 wkd1 wkd2 ndays_with_wkends ndays_no_wkends
(date) (date) (chr) (chr) (dbl) (dbl)
1 2015-03-30 2015-03-30 Monday Monday 0 0
2 2015-03-30 2015-03-31 Monday Tuesday 1 1
3 2015-03-30 2015-04-01 Monday Wednesday 2 2
4 2015-03-30 2015-04-02 Monday Thursday 3 3
5 2015-03-30 2015-04-03 Monday Friday 4 4
6 2015-03-30 2015-04-04 Monday Saturday 5 4
7 2015-03-30 2015-04-05 Monday Sunday 6 4
8 2015-03-30 2015-04-06 Monday Monday 7 5
9 2015-03-30 2015-04-13 Monday Monday 14 10
10 2015-03-30 2015-04-20 Monday Monday 21 15
11 2015-04-20 2015-03-30 Monday Monday 21 15
12 1986-05-11 2013-09-26 Sunday Thursday 10000 7143
б. Китайский язык
> Sys.setlocale(category="LC_ALL", locale = "chinese")
[1] "LC_COLLATE=Chinese (Simplified)_People Republic of China.936;LC_CTYPE=Chinese (Simplified)_People Republic of China.936;LC_MONETARY=Chinese (Simplified)_People Republic of China.936;LC_NUMERIC=C;LC_TIME=Chinese (Simplified)_People Republic of China.936"
> Sys.getlocale()
[1] "LC_COLLATE=Chinese (Simplified)_People Republic of China.936;LC_CTYPE=Chinese (Simplified)_People Republic of China.936;LC_MONETARY=Chinese (Simplified)_People Republic of China.936;LC_NUMERIC=C;LC_TIME=Chinese (Simplified)_People Republic of China.936"
> e2 %>%
mutate(wkd1 = format(date.1, "%A"),
wkd2 = format(date.2, "%A"),
ndays_with_wkends = ifelse((date.2 > date.1), (date.2 - date.1), (date.1 - date.2)),
ndays_no_wkends = Nweekdays(date.1, date.2))
Source: local data frame [12 x 6]
date.1 date.2 wkd1 wkd2 ndays_with_wkends ndays_no_wkends
(date) (date) (chr) (chr) (dbl) (dbl)
1 2015-03-30 2015-03-30 ÐÇÆÚÒ» ÐÇÆÚÒ» 0 0
2 2015-03-30 2015-03-31 ÐÇÆÚÒ» ÐÇÆÚ¶þ 1 1
3 2015-03-30 2015-04-01 ÐÇÆÚÒ» ÐÇÆÚÈý 2 2
4 2015-03-30 2015-04-02 ÐÇÆÚÒ» ÐÇÆÚËÄ 3 3
5 2015-03-30 2015-04-03 ÐÇÆÚÒ» ÐÇÆÚÎå 4 4
6 2015-03-30 2015-04-04 ÐÇÆÚÒ» ÐÇÆÚÁù 5 5
7 2015-03-30 2015-04-05 ÐÇÆÚÒ» ÐÇÆÚÈÕ 6 6
8 2015-03-30 2015-04-06 ÐÇÆÚÒ» ÐÇÆÚÒ» 7 7
9 2015-03-30 2015-04-13 ÐÇÆÚÒ» ÐÇÆÚÒ» 14 14
10 2015-03-30 2015-04-20 ÐÇÆÚÒ» ÐÇÆÚÒ» 21 21
11 2015-04-20 2015-03-30 ÐÇÆÚÒ» ÐÇÆÚÒ» 21 21
12 1986-05-11 2013-09-26 ÐÇÆÚÈÕ ÐÇÆÚËÄ 10000 10000
2. Локальное независимое решение: Nweekdays() функция адаптирована из решения @Sacha Epskamp. Он работает для всех локалей, однако @Sacha Epskamp использовал c(0,6) для отсеивания выходных дней, который отличается от этого решения, которое использует c(2,3) для извлечения выходных дней.
Nweekdays <- Vectorize(
function(a, b) {
return(sum(!(((as.numeric(b:a)) %% 7) %in% c(2,3))) - 1) # 2: Saturday and 3: Sunday
})
а. Английский язык
> Sys.setlocale(category="LC_ALL", locale = "English_United States.1252")
[1] "LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252"
> Sys.getlocale()
[1] "LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252"
> e2 %>%
mutate(wkd1 = format(date.1, "%A"),
wkd2 = format(date.2, "%A"),
ndays_with_wkends = ifelse((date.2 > date.1), (date.2 - date.1), (date.1 - date.2)),
ndays_no_wkends = Nweekdays(date.1, date.2))
Source: local data frame [12 x 6]
date.1 date.2 wkd1 wkd2 ndays_with_wkends ndays_no_wkends
(date) (date) (chr) (chr) (dbl) (dbl)
1 2015-03-30 2015-03-30 Monday Monday 0 0
2 2015-03-30 2015-03-31 Monday Tuesday 1 1
3 2015-03-30 2015-04-01 Monday Wednesday 2 2
4 2015-03-30 2015-04-02 Monday Thursday 3 3
5 2015-03-30 2015-04-03 Monday Friday 4 4
6 2015-03-30 2015-04-04 Monday Saturday 5 4
7 2015-03-30 2015-04-05 Monday Sunday 6 4
8 2015-03-30 2015-04-06 Monday Monday 7 5
9 2015-03-30 2015-04-13 Monday Monday 14 10
10 2015-03-30 2015-04-20 Monday Monday 21 15
11 2015-04-20 2015-03-30 Monday Monday 21 15
12 1986-05-11 2013-09-26 Sunday Thursday 10000 7143
б. Китайский язык
> Sys.setlocale(category="LC_ALL", locale = "chinese")
[1] "LC_COLLATE=Chinese (Simplified)_People Republic of China.936;LC_CTYPE=Chinese (Simplified)_People Republic of China.936;LC_MONETARY=Chinese (Simplified)_People Republic of China.936;LC_NUMERIC=C;LC_TIME=Chinese (Simplified)_People Republic of China.936"
> Sys.getlocale()
[1] "LC_COLLATE=Chinese (Simplified)_People Republic of China.936;LC_CTYPE=Chinese (Simplified)_People Republic of China.936;LC_MONETARY=Chinese (Simplified)_People Republic of China.936;LC_NUMERIC=C;LC_TIME=Chinese (Simplified)_People Republic of China.936"
> e2 %>%
mutate(wkd1 = format(date.1, "%A"),
wkd2 = format(date.2, "%A"),
ndays_with_wkends = ifelse((date.2 > date.1), (date.2 - date.1), (date.1 - date.2)),
ndays_no_wkends = Nweekdays(date.1, date.2))
Source: local data frame [12 x 6]
date.1 date.2 wkd1 wkd2 ndays_with_wkends ndays_no_wkends
(date) (date) (chr) (chr) (dbl) (dbl)
1 2015-03-30 2015-03-30 ÐÇÆÚÒ» ÐÇÆÚÒ» 0 0
2 2015-03-30 2015-03-31 ÐÇÆÚÒ» ÐÇÆÚ¶þ 1 1
3 2015-03-30 2015-04-01 ÐÇÆÚÒ» ÐÇÆÚÈý 2 2
4 2015-03-30 2015-04-02 ÐÇÆÚÒ» ÐÇÆÚËÄ 3 3
5 2015-03-30 2015-04-03 ÐÇÆÚÒ» ÐÇÆÚÎå 4 4
6 2015-03-30 2015-04-04 ÐÇÆÚÒ» ÐÇÆÚÁù 5 4
7 2015-03-30 2015-04-05 ÐÇÆÚÒ» ÐÇÆÚÈÕ 6 4
8 2015-03-30 2015-04-06 ÐÇÆÚÒ» ÐÇÆÚÒ» 7 5
9 2015-03-30 2015-04-13 ÐÇÆÚÒ» ÐÇÆÚÒ» 14 10
10 2015-03-30 2015-04-20 ÐÇÆÚÒ» ÐÇÆÚÒ» 21 15
11 2015-04-20 2015-03-30 ÐÇÆÚÒ» ÐÇÆÚÒ» 21 15
12 1986-05-11 2013-09-26 ÐÇÆÚÈÕ ÐÇÆÚËÄ 10000 7143
Я написал этот, но другой ответ лучше:)
Nweekdays <- function(a,b)
{
dates <- as.Date(as.Date(a,"%m/%d/%y",origin="1900-01-01"):as.Date(b,"%m/%d/%y",origin="1900-01-01"),origin="1900-01-01")
days <- format(dates,"%w")[c(-1,-length(dates))]
return(sum(!days%in%c(0,6)))
}
Nweekdays('01/30/2011','02/04/2011')
[1] 3
EDIT: подсчет количества дней недели между двумя указанными днями.
Взяв рекомендацию J. Winchesters, эта функция может быть упорядочена как:
Nweekdays <- function(a,b)
{
dates <- as.numeric((as.Date(a,"%m/%d/%y")):(as.Date(b,"%m/%d/%y")))
dates <- dates[- c(1,length(dates))]
return(sum(!dates%%7%in%c(0,6)))
}
Некоторые результаты:
> Nweekdays('01/30/2011','02/04/2011')
[1] 4
>
> Nweekdays('01/30/2011','01/30/2011')
[1] 0
>
> Nweekdays('01/30/2011','01/25/2011')
[1] 3
Обратите внимание, что это не зависит от языка. (По этой теме, как мне изменить язык в любом случае?)
Я использую следующий подход: сначала помощник:
weekDays <- function(UPPER = TRUE) {
days <- c('MONDAY', 'TUESDAY', 'WEDNESDAY',
'THURSDAY', 'FRIDAY', 'SATURDAY',
'SUNDAY')
if(!UPPER) return(.Internal(tolower(days)))
days
}
... и теперь основная функция:
NumWeekDays <- function(dd, Xdays = c('saturday', 'sunday')) {
# a function to count the number of non-Xdays in a month
# >
# first check if Xdays is of correct format
stopifnot( all(.Internal(tolower(Xdays)) %in% weekDays(UP = FALSE)))
# >
# a helper function to find the number of non-X days between two dates
NonXDays <- function(startDate, endDate, Xdays) {
sum(!(.Internal(tolower(weekdays(seq(startDate, endDate, 'day')))) %in%
.Internal(tolower(Xdays))))
}
startDate <- as.Date(as.yearmon(index(dd)), frac = 0)
endDate <- as.Date(as.yearmon(index(dd)), frac = 1)
vapply(1:nrow(dd),
FUN = function(i) NonXDays(startDate[i],
endDate[i],
Xdays = c('saturday', 'sunday')),
FUN.VALUE = numeric(1))
}
Пример:
set.seed(1)
dx <- apply.monthly(xts(rnorm(600), order.by = Sys.Date() + 1:600), mean)
R> NumWeekDays(dx)
[1] 23 21 22 23 20 23 22 20 22 22 21 22 23 21 22 22 21 23 21 21
Работая с использованием lubridate, вы можете создать такую функцию:
library(lubridate)
WorkingDays_function <- function(StartDate,EndDate){
startDate <- dmy(StartDate)
endDate <- dmy(EndDate)
#Now build a sequence between the dates:
myDates <-seq(from = startDate, to = endDate, by = "days")
#Week starts on Sunday (1) so to exclude Sun (1) and Sat (7)
#use > 1 & < 7
working_days <- sum(wday(myDates)>1 & wday(myDates)<7)
print(working_days)
}
WorkingDays_function("11/07/2019","20/07/2019")
Дж. Вин. ответ хороший, но может быть немного быстрее с lubridate.
require(lubridate)
count_weekdays<- Vectorize(function(from,to) sum(!wday(seq(from, to, "days")) %in% c(1,7)))
Вот результаты моего компьютера:
> v1<- seq(from = ymd(19000101), to = ymd(20000101), by='month')
> v2<- seq(from = ymd(20000101), to = ymd(21000101), by='month')
> require(tictoc)
> tic(); out<- Nweekdays(v1,v2); toc();
293.06 sec elapsed
> tic(); out<- count_weekdays(v1,v2); toc();
9.95 sec elapsed
Примерно в 30 раз быстрее. Имеет смысл, если вы делаете много периодов.