我有一个包含很多列的 data.frame,其中一列包含样本区域的代码,另一列包含样本的编号。我想从每个样本区域的倒数第二个样本中提取信息子集。我尝试了很多不同的事情......最后这是我最好的猜测......但它仍然不起作用。
site <- sample (1:3, 10, replace= T)
d2 <- sample (1:5, 10, replace= T)
d3 <- sample (1:5, 10, replace= T)
samplet <- sample (1:4, 10, replace= T)
mydata <- data.frame (cbind(site, d2, d3, samplet))
penultimate <- matrix(NA,,) # here I dont know how the return will be, as I dont know how the dataframe will change
si <- matrix (NA, , )
pl <- unique (site)
for (i in 1:(length (pl))) {
si <- mydata[which (samplet==pl[i]),] # I tried to create a temporary matrix, so I can calculate each site at a time
penultimate <- si[which (si$samplet!=(max(si$samplet[si$samplet!=max(si$samplet)]))),]
}
Cheers!
一个简单的方法是使用data.table
及其内置的.N
value
# assuming `d1` is the column from which you want to find the penultimate
mydata <- data.frame(d1=strsplit("AAABBCCCCCDD", "")[[1]], d2=rnorm(12), d3=LETTERS[1:12], d4=c(101:103, 201:202, 301:305, 401:402))
DT <- data.table(mydata)
DT[, .SD[.N-1], by=d1]
d1 d2 d3 d4
1: A 1.6906714 B 102
2: B -0.1239458 D 201
3: C -0.2976339 I 304
4: D 0.6858120 K 401
与我的数据比较
> mydata
d1 d2 d3 d4
1 A 0.5986002 A 101
2 A 1.6906714 B 102 <~~~~ \
3 A -0.3253657 C 103
4 B -0.1239458 D 201 <~~~~ -\
5 B 0.8261401 E 202
6 C 0.0601318 F 301 Penultimate Values by d1
7 C -0.9766622 G 302
8 C 0.1028259 H 303
9 C -0.2976339 I 304 <~~~~~ -/
10 C -1.1467000 J 305
11 D 0.6858120 K 401 <~~~~~ /
12 D -0.6160335 L 402
编辑,用新的示例数据更新。
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)