2016-01-05 40 views
0

我有一個data.table說,dat。以下是其dput保留x行並刪除2列的休息

structure(list(Sepal.Length = c(5.1, 4.9, 4.7, 4.6, 5, 5.4, 4.6, 
5, 4.4, 4.9), Sepal.Width = c(3.5, 3, 3.2, 3.1, 3.6, 3.9, 3.4, 
3.4, 2.9, 3.1), Petal.Length = c(1.4, 1.4, 1.3, 1.5, 1.4, 1.7, 
1.4, 1.5, 1.4, 1.5), Petal.Width = c(0.2, 0.2, 0.2, 0.2, 0.2, 
0.4, 0.4, 0.2, 0.1, 0.1), Species = c("a", "a", "a", "a", "b", 
"b", "b", "b", "b", "b")), .Names = c("Sepal.Length", "Sepal.Width", 
"Petal.Length", "Petal.Width", "Species"), row.names = c(NA, 
-10L), class = c("data.table", "data.frame"), .internal.selfref = <pointer: 0x0000000000330788>) 


dat 

    Sepal.Length Sepal.Width Petal.Length Petal.Width Species 
1:   5.1   3.5   1.4   0.2  a 
2:   4.9   3.0   1.4   0.2  a 
3:   4.7   3.2   1.3   0.2  a 
4:   4.6   3.1   1.5   0.2  a 
5:   5.0   3.6   1.4   0.2  b 
6:   5.4   3.9   1.7   0.4  b 
7:   4.6   3.4   1.4   0.4  b 
8:   5.0   3.4   1.5   0.2  b 
9:   4.4   2.9   1.4   0.1  b 
10:   4.9   3.1   1.5   0.1  b 

我想要做以下,

1)計算的列的數據Petal.WidthSpecies的任何組合(其爲2以上的例子的最小頻率)。我們稱這個號碼爲x。以下是這2列的數據組合

Species Petal.Width Frequency 
a   0.2   4 
b   0.2   2 
b   0.4   2 
b   0.1   2  

2)只保留x(2在我們的例子),在隨機觀察了Petal.Width列的每個組合和Species

當我們只有保持1例在這種情況下,我們隨機做

dat <- unique(dat, by = c("Petal.Width", "Species") 

但是我們什麼做的時候,我們要保持x案件而不只是1(即太隨機)?

輸出可能看起來像以下

Sepal.Length Sepal.Width Petal.Length Petal.Width Species 
2:   4.9   3.0   1.4   0.2  a 
4:   4.6   3.1   1.5   0.2  a 
5:   5.0   3.6   1.4   0.2  b 
6:   5.4   3.9   1.7   0.4  b 
7:   4.6   3.4   1.4   0.4  b 
8:   5.0   3.4   1.5   0.2  b 
9:   4.4   2.9   1.4   0.1  b 
10:   4.9   3.1   1.5   0.1  b 

在這個例子中,由於其他的組合滿足最小頻率標準僅2例組合的(A,0.2)將被刪除。

回答

2

第一部分:

x <- dat[, .N, by = list(Species, Petal.Width)][, min(N)] 

第二部分:

dat[, .SD[sample(1:.N, x)], by = list(Species, Petal.Width)] 
+0

'樣品(.N,X)'可以代替樣品的'被使用(1:.N,x)的'。另外,你可以像'dat [dat [,.I [sample(.N,2)],by = c(「Petal.Width」,「Species」)] $ V1]那樣使用'.I''作爲eddi顯示的這裏:http://stackoverflow.com/a/16574176/1191259 – Frank