有條件計數R

Question

我有一個問題，我希望你們中的一些人可以幫助我。我正在撰寫關於製藥和平行進口的影響的論文。我正在處理這個R，有一個面板數據集

我需要一個變量，計算給定的原始產品 - 在給定的時間段內有多少個並行進口商。

Product_ID  PI  t  
    1   0  1 
    1   1  1 
    1   1  1 
    1   0  2 
    1   1  2 
    1   1  2 
    1   1  2 
    1   1  2 
    2   0  1 
    2   1  1 
    2   0  2 
    2   1  2 
    2   0  3 
    2   1  3 
    2   1  3 
    2   1  3 

理想我想在這裏什麼是新的一列，等的PI-產品（PI = 1），用於在時間的原始（PI = 0），噸數。所以輸出會是這樣：

Product_ID  PI  t  nPIcomp 
    1   0  1  2 
    1   1  1 
    1   1  1 
    1   0  2  4 
    1   1  2 
    1   1  2 
    1   1  2 
    1   1  2 
    2   0  1  1 
    2   1  1 
    2   0  2  1 
    2   1  2 
    2   0  3  3 
    2   1  3 
    2   1  3 
    2   1  3 

我希望我已經作出了明確的問題:)

由於提前， 亨裏克

Answer 1

像這樣的事情？

x <- read.table(text = "Product_ID  PI  t  
    1   0  1 
    1   1  1 
    1   1  1 
    1   0  2 
    1   1  2 
    1   1  2 
    1   1  2 
    1   1  2 
    2   0  1 
    2   1  1 
    2   0  2 
    2   1  2 
    2   0  3 
    2   1  3 
    2   1  3 
    2   1  3", header = TRUE) 

find.count <- rle(x$PI) 
count <- find.count$lengths[find.count$values == 1] 
x[x$PI == 0, "nPIcomp"] <- count 

    Product_ID PI t nPIcomp 
1   1 0 1  2 
2   1 1 1  NA 
3   1 1 1  NA 
4   1 0 2  4 
5   1 1 2  NA 
6   1 1 2  NA 
7   1 1 2  NA 
8   1 1 2  NA 
9   2 0 1  1 
10   2 1 1  NA 
11   2 0 2  1 
12   2 1 2  NA 
13   2 0 3  3 
14   2 1 3  NA 
15   2 1 3  NA 
16   2 1 3  NA 

Answer 2

羅馬的回答讓正是你想要什麼。如果你想總結這將是得心應手的數據，使用plyr pacakge（df是我叫你data.frame）...

ddply(df , .(Product_ID , t) , summarise , nPIcomp = sum(PI)) 
# Product_ID t nPIcomp 
#1   1 1  2 
#2   1 2  4 
#3   2 1  1 
#4   2 2  1 
#5   2 3  3 

Answer 3

我會用ave和你的兩列Product_ID和t爲分組變量。然後，在每個組內，應用返回PI的隨後的NA S的適當數目之和的函數：

dat <- transform(dat, nPIcomp = ave(PI, Product_ID, t, 
            FUN = function(z) { 
             n <- sum(z) 
             c(n, rep(NA, n)) 
            })) 

同樣的想法可以用data.table包被使用，如果你的數據是大和速度是關心。