Cách thêm cột đếm dựa trên giá trị trong các cột khác trong R

Question

Tôi có một tập dữ liệu tương đối lớn (16.000 + x ~ 31). Nói cách khác, nó đủ lớn đến mức tôi không muốn thao tác nó theo từng dòng trong Excel. Dữ liệu có dạng này:

block site  day X1 X2 
1  1  1 0.4 5.1 
1  1  2 0.8 1.1 
1  1  3 1.1 4.2 
1  2  1 ... ... 
1  2  2 
1  2  3 
2  3  1 
2  3  2 
2  3  3 
2  4  1 
2  4  2 
2  4  3 

Như bạn có thể thấy, số lượng trang web liên tục nhưng tôi muốn một cột mà số trang web đặt lại với mỗi khối. Ví dụ, tôi muốn một cái gì đó như thế này dưới đây:

block site  day X1 X2 site2 
1  1  1 0.4 5.1 1 
1  1  2 0.8 1.1 1 
1  1  3 1.1 4.2 1 
1  2  1 ... ... 2 
1  2  2    2 
1  2  3    2 
2  3  1    1 
2  3  2    1 
2  3  3    1 
2  4  1    2 
2  4  2    2 
2  4  3    2 

Tôi đã suy nghĩ về việc sử dụng hàm R rle nhưng không chắc chắn nếu nó sẽ làm việc vì biến chứng với ngày. Nếu không, tôi sẽ thử một cái gì đó như:

Data$site2 <- sequence(rle(Data$block)$lengths) 

Có ai có đề xuất thêm số lượng cột (chuỗi) số trang trong mỗi khối không? Nếu nó giúp, có số ngày tương tự (263) được ghi lại cho mỗi trang nhưng có số lượng trang web khác nhau trên mỗi khối.

Answer 1

Dưới đây là một giải pháp hơi vụng về sử dụng plyr và ddply:

ddply(df,.(block),transform, 
        site1 = rep(1:length(unique(site)), 
          times = rle(site)$lengths)) 

Hoặc một phiên bản hơi tinh quái:

ddply(df,.(block),transform,site1 = as.integer(as.factor(site))) 

Có thể có một cách thông minh để làm điều này trực tiếp, tuy nhiên, sử dụng các chức năng khác nhau seq, sequence và rle, nhưng bộ não của tôi hơi bị mờ một chút. Nếu bạn để nó mở cho một chút, một người nào đó có thể sẽ đi cùng với một giải pháp không trơn tru plyr.

Answer 2

Sử dụng tapply có thể làm việc

# Make some fake data 
dat <- data.frame(block = rep(1:3, each = 4), site = rep(1:6, each = 2), val = rnorm(12)) 
# For each block reset the count 
dat$site2 <- unlist(tapply(dat$site, dat$block, function(x){x - min(x) + 1})) 

Answer 3

Via ave:

df1 <- structure(list(block = c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2), 
    site = c(1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4), day = c(1, 
    2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3)), .Names = c("block", "site", 
"day"), row.names = c("2", "3", "4", "5", "6", "7", "8", "9", 
"10", "11", "12", "13"), class = "data.frame") 

df1$site2 <- ave(df1$site,df1$block,FUN=function(x) match(x,sort(unique(x))))