Tính lứa tuổi trong R

Question

Tôi có hai khung dữ liệu trong R. Một khung có một người năm sinh:

YEAR 
/1931 
/1924 

và sau đó cột khác cho thấy một thời gian gần đây.

RECENT 
09/08/2005 
11/08/2005 

Điều tôi muốn làm là trừ năm để tôi có thể tính tuổi của họ trong nhiều năm, tuy nhiên tôi không biết cách tiếp cận điều này. Bất kỳ giúp đỡ xin vui lòng?

Answer 1

Bạn có thể làm một số các định dạng:

as.numeric(format(as.Date("01/01/2010", format="%m/%d/%Y"), format="%Y")) - 1930 

Với dữ liệu của bạn:

> yr <- c(1931, 1924) 
> recent <- c("09/08/2005", "11/08/2005") 
> as.numeric(format(as.Date(recent, format="%m/%d/%Y"), format="%Y")) - yr 
[1] 74 81 

Kể từ khi bạn có dữ liệu của bạn trong một data.frame (Tôi sẽ cho rằng nó được gọi là df), nó sẽ giống như thế này:

as.numeric(format(as.Date(df$recent, format="%m/%d/%Y"), format="%Y")) - df$year 

Answer 2

Dựa trên câu trả lời trước, chuyển đổi cột của bạn thành đối tượng ngày s và trừ. Một số chuyển đổi các loại giữa ký tự và số là cần thiết:

> foo=data.frame(RECENT=c("09/08/2005","11/08/2005"),YEAR=c("/1931","/1924")) 
> foo 
     RECENT YEAR 
1 09/08/2005 /1931 
2 11/08/2005 /1924 
> foo$RECENTd = as.Date(foo$RECENT, format="%m/%d/%Y") 
> foo$YEARn = as.numeric(substr(foo$YEAR,2,999)) 
> foo$AGE = as.numeric(format(foo$RECENTd,"%Y")) - foo$YEARn 
> foo 
     RECENT YEAR RECENTd YEARn AGE 
1 09/08/2005 /1931 2005-09-08 1931 74 
2 11/08/2005 /1924 2005-11-08 1924 81 

Lưu ý: Tôi giả định bạn có dấu gạch chéo trong cột năm của bạn.

Ngoài ra, mẹo khi đặt câu hỏi về ngày là bao gồm ngày quá mười hai để chúng tôi biết bạn là người tháng/ngày/năm hay người ngày/tháng/năm.

Answer 3

Với dữ liệu trong ví dụ của bạn:

> m <- data.frame(YEAR=c("/1931", "/1924"),RECENT=c("09/08/2005","11/08/2005")) 
> m 
    YEAR  RECENT 
1 /1931 09/08/2005 
2 /1924 11/08/2005 

Extract năm với strptime chức năng:

> strptime(m[,2], format = "%m/%d/%Y")$year - strptime(m[,1], format = "/%Y")$year 
[1] 74 81 

Answer 4

Bạn có thể giải quyết điều này với các gói lubridate.

> library(lubridate) 

Tôi không nghĩ/1931 là lớp ngày chung. Vì vậy, tôi sẽ giả định tất cả các mục là chuỗi ký tự.

> RECENT <- data.frame(recent = c("09/08/2005", "11/08/2005")) 
> YEAR <- data.frame(year = c("/1931", "/1924")) 

Trước tiên, hãy thông báo cho R rằng ngày gần đây là ngày. Tôi sẽ giả định ngày tháng theo thứ tự tháng/ngày/năm, vì vậy tôi sử dụng mdy(). Nếu họ đặt hàng ngày/tháng/năm, chỉ cần sử dụng dmy().

> RECENT$recent <- mdy(RECENT$recent) 
     recent 
1 2005-09-08 
2 2005-11-08 

Bây giờ, hãy biến năm thành số để chúng tôi có thể thực hiện một số phép tính với chúng.

> YEAR$year <- as.numeric(substr(YEAR$year, 2, 5)) 

Bây giờ, chỉ cần thực hiện phép tính. year() chiết xuất giá trị năm của ngày RECENT.

> year(RECENT$recent) - YEAR 
    year 
1 74 
2 81 

p.s.nếu mục năm của bạn đang thực sự đầy đủ ngày, bạn có thể nhận được sự khác biệt trong năm với

> YEAR1 <- data.frame(year = mdy("01/08/1931","01/08/1924")) 
> as.period(RECENT$recent - YEAR1$year, units = "year") 
[1] 74 years and 8 months 81 years and 10 months 

Answer 5

Hàm ở phía dưới có một vectơ của các đối tượng ngày và tính toán các thời đại, chiếm một cách chính xác cho năm nhuận. Dường như là một giải pháp đơn giản hơn bất kỳ câu trả lời nào khác.

age = function(from, to) { 
    from_lt = as.POSIXlt(from) 
    to_lt = as.POSIXlt(to) 

    age = to_lt$year - from_lt$year 

    ifelse(to_lt$mon < from_lt$mon | 
     (to_lt$mon == from_lt$mon & to_lt$mday < from_lt$mday), 
     age - 1, age) 
} 

Answer 6

Tôi nghĩ rằng đây có thể là một chút trực quan hơn và không cần định dạng hoặc tước:

as.numeric(as.Date("2002-02-02") - as.Date("1924-08-03"))/365 

cho đầu ra:

77.55342 

Sau đó, bạn có thể sử dụng sàn(), tròn (), hoặc trần() để làm tròn thành một số nguyên.

Answer 7

cách thực sự vững chắc mà còn hỗ trợ vectơ sử dụng gói lubridate:

age <- function(date.birth, date.ref = Sys.Date()) { 
    if (length(date.birth) > 1 & length(date.ref) == 1) { 
    date.ref <- rep(date.ref, length(date.birth)) 
    } 

    date.birth.monthdays <- paste0(month(date.birth), day(date.birth)) %>% as.integer() 
    date.ref.monthdays <- paste0(month(date.ref), day(date.ref)) %>% as.integer() 

    age.calc <- 0 

    for (i in 1:length(date.birth)) { 
    if (date.birth.monthdays[i] <= date.ref.monthdays[i]) { 
     # didn't had birthday 
     age.calc[i] <- year(date.ref[i]) - year(date.birth[i]) 
    } else { 
     age.calc[i] <- year(date.ref[i]) - year(date.birth[i]) - 1 
    } 
    } 
    age.calc 
} 

này cũng chiếm tới năm nhuận. Tôi chỉ kiểm tra xem ai đó đã có một sinh nhật chưa.

Answer 8

Tôi sử dụng chức năng tùy chỉnh, xem mã bên dưới, thuận tiện để sử dụng trong biến đổi và khá linh hoạt (bạn cần gói lubridate).

Ví dụ

get_age("2000-01-01") 
# [1] 17 
get_age(lubridate::as_date("2000-01-01")) 
# [1] 17 
get_age("2000-01-01","2015-06-15") 
# [1] 15 
get_age("2000-01-01",dec = TRUE) 
# [1] 17.92175 
get_age(c("2000-01-01","2003-04-12")) 
# [1] 17 14 
get_age(c("2000-01-01","2003-04-12"),dec = TRUE) 
# [1] 17.92176 14.64231 

Chức năng

#' Get age 
#' 
#' Returns age, decimal or not, from single value or vector of strings 
#' or dates, compared to a reference date defaulting to now. Note that 
#' default is NOT the rounded value of decimal age. 
#' @param from_date vector or single value of dates or characters 
#' @param to_date date when age is to be computed 
#' @param dec return decimal age or not 
#' @examples 
#' get_age("2000-01-01") 
#' get_age(lubridate::as_date("2000-01-01")) 
#' get_age("2000-01-01","2015-06-15") 
#' get_age("2000-01-01",dec = TRUE) 
#' get_age(c("2000-01-01","2003-04-12")) 
#' get_age(c("2000-01-01","2003-04-12"),dec = TRUE) 
get_age <- function(from_date,to_date = lubridate::now(),dec = FALSE){ 
    if(is.character(from_date)) from_date <- lubridate::as_date(from_date) 
    if(is.character(to_date)) to_date <- lubridate::as_date(to_date) 
    if (dec) { age <- lubridate::interval(start = from_date, end = to_date)/(lubridate::days(365)+lubridate::hours(6)) 
    } else { age <- lubridate::year(lubridate::as.period(lubridate::interval(start = from_date, end = to_date)))} 
    age 
}