Làm cách nào để chọn hàng đầu tiên trong khung dữ liệu R đáp ứng các tiêu chí nhất định?

Question

Làm cách nào để chọn hàng đầu tiên của khung dữ liệu R đáp ứng các tiêu chí nhất định?

Đây là bối cảnh:

Tôi có một khung dữ liệu với năm cột:

"pixel", "year","propvar", "component", "cumsum." 

Có 1225 kết hợp của pixel và year, vì dữ liệu đã được tính toán từ các chuỗi thời gian hàng năm của 49 điểm ảnh địa lý cho mỗi 25 năm nghiên cứu. Trong mỗi năm điểm ảnh, tôi đã tính toán propvar, tỷ lệ phương sai tổng được giải thích bởi một thành phần nhất định của phép biến đổi Fourier nhanh cho chuỗi thời gian của một năm pixel nhất định. Sau đó tôi tính được cumsum, là tổng tích luỹ của propvar cho từng thành phần tần suất trong một năm pixel. Cột component chỉ cung cấp cho bạn chỉ mục cho thành phần chuỗi Fourier (cộng 1) mà từ đó propvar được tính toán.

Tôi muốn xác định số lượng thành phần cần thiết để giải thích lớn hơn 99% phương sai. Tôi tìm một cách để thực hiện việc này là tìm hàng đầu tiên trong mỗi năm pixel trong đó cumsum> 0,99 và tạo khung dữ liệu từ đó với ba cột, pixel, year và numbercomps, trong đó numbercomps là số lượng thành phần bắt buộc trong một năm pixel nhất định để giải thích lớn hơn 99% phương sai. Tôi không biết làm thế nào để làm điều này trong R. Có ai có một giải pháp?

Answer 1

Chắc chắn. Một cái gì đó như thế này nên làm như lừa:

# CREATE A REPRODUCIBLE EXAMPLE! 
df <- data.frame(year = c("2001", "2003", "2001", "2003", "2003"), 
       pixel = c("a", "b", "a", "b", "a"), 
       cumsum = c(99, 99, 98, 99, 99), 
       numbercomps=1:5) 
df 
# year pixel cumsum numbercomps 
# 1 2001  a  99   1 
# 2 2003  b  99   2 
# 3 2001  a  98   3 
# 4 2003  b  99   4 
# 5 2003  a  99   5 

# EXTRACT THE SUBSET YOU'D LIKE. 
res <- subset(df, cumsum>=99) 
res <- subset(res, 
       subset = !duplicated(res[c("year", "pixel")]), 
       select = c("pixel", "year", "numbercomps")) 
# pixel year numbercomps 
# 1  a 2001   1 
# 2  b 2003   2 
# 5  a 2003   5 

EDIT Ngoài ra, đối với những quan tâm data.table, có là thế này:

library(data.table) 
dt <- data.table(df, key="pixel, year")  
dt[cumsum>=99, .SD[1], by=key(dt)]