根据列值从列表中选择数据帧

时间:2013-12-07 15:28:44

标签: r list

使用R,我试图通过仅在其中一列中存在特定值组合时从现有列表中选择数据帧来创建新的数据帧列表。让我解释一下工作正常的第一步。这是我在名为df:

的数据框中的原始数据
                                  Taxon      C     N    func.group  trophic.grp
1  Chrysomelidae.Phylotreta.exclamationis -30.23  5.06     grazer   herbivore
2        Chrysomelidae.Neocrepidodera.sp. -27.29  5.55     grazer   herbivore
3        Chrysomelidae.Neocrepidodera.sp. -27.84  5.54     grazer   herbivore
4        Chrysomelidae.Neocrepidodera.sp. -27.69  4.59     grazer   herbivore
5              Mitidulidae.Meligethes.sp. -26.99  5.30     grazer   herbivore
6           Chrysomelidae.Phylotreta.sp.2 -28.50  2.40     grazer   herbivore
7           Chrysomelidae.Phylotreta.sp.2 -28.36  4.17     grazer   herbivore
8           Chrysomelidae.Phylotreta.sp.2 -29.50  3.15     grazer   herbivore
9           Chrysomelidae.Phylotreta.sp.2 -27.69  3.72     grazer   herbivore
10          Chrysomelidae.Phylotreta.sp.2 -28.22  3.26     grazer   herbivore
11                  Gastropoda.snail.sp.1 -26.21  3.54     grazer   herbivore
12                  Gastropoda.snail.sp.1 -27.59  2.61     grazer   herbivore
13                  Gastropoda.snail.sp.1 -25.10  2.66     grazer   herbivore
14                  Gastropoda.snail.sp.2 -26.49  2.55     grazer   herbivore
15                  Gastropoda.snail.sp.4 -27.46 -0.38     grazer   herbivore
16       Lepidoptera.Arctidae.Ermine.moth -28.51  2.44     grazer   herbivore
17       Curculionidae.Ischapterapion.sp. -29.06  2.19     weevil   herbivore
18       Curculionidae.Ischapterapion.sp. -29.27  1.60     weevil   herbivore
19       Curculionidae.Ischapterapion.sp. -29.94  2.08     weevil   herbivore
20       Curculionidae.Ischapterapion.sp. -29.71  2.16     weevil   herbivore
21            Curculionidae.Protapion.sp. -28.45  1.91     weevil   herbivore
22            Curculionidae.Protapion.sp. -25.99  0.55     weevil   herbivore
23            Curculionidae.Protapion.sp. -28.27  1.52     weevil   herbivore
24            Curculionidae.Protapion.sp. -28.01  1.74     weevil   herbivore
25            Curculionidae.Protapion.sp. -27.06  0.54     weevil   herbivore
26             Curculionidae.Hypera.meles -25.41  3.38     weevil   herbivore
27               Curculionidae.Sitona.sp. -27.05  2.01     weevil   herbivore
28               Curculionidae.Sitona.sp. -26.70  3.07     weevil   herbivore
29               Curculionidae.Sitona.sp. -27.64  2.13     weevil   herbivore
30               Curculionidae.Sitona.sp. -27.50  1.47     weevil   herbivore
31            Curculionidae.Phylobius.sp. -28.27  2.66     weevil   herbivore
32      Curculionidae.Hypera.nigrorostris -25.52  2.43     weevil   herbivore

此数据框(df)包含14个不同的“Taxon”,其中一些具有多个样本,因此总共有32个样本。每个分类群也被“func.group”列分类为“grazer”或“象鼻虫”。

首先,我想从我的14中随机选择6个分类,对于所有可能的6个组合。因此,有6个分类单元的3003种组合,可以从14个(随机抽样而无需替换,顺序不重要) )。对于所选的每个分类,我想要包括该分类的所有样本。我使用这个代码,效果很好:

combos<-combn(unique(as.character(df$Taxon)), 6) 

接下来我想要包含所有其他信息列,因此我使用这个额外的代码行,对于每个选择的Taxon,它会添加其他数据列,并且运行良好:

mysamples <- apply(combos, 2, function(vec) df[ df$Taxon %in% vec, ] )

那么我们就解决了我的问题。从“mysamples”(现在应该是3003个数据帧的列表),我想选择包含至少一个“grazer”的Taxon和一个“象鼻虫”的Taxon的所有数据帧,并将这些数据帧存储在一个新的清单。

因此,我希望这个新列表仅包含包含象鼻虫和食草动物的分类的数据框。 (每个数据框中的6个分类中有多少是象鼻虫或食草动物,只要一个是grazer,一个是象鼻虫。)

非常感谢, 中号

2 个答案:

答案 0 :(得分:0)

试试此代码

mysamples[unlist(lapply(mysamples,
                        function(x) !any(is.na(match(levels(df$func.group),
                                                     x$func.group)))))]

如果丢失了grazer或象鼻虫,match将返回NA,any因此将返回TRUE,它将被反转(!),因此此数据帧不会在最后一个中使用< / p>

答案 1 :(得分:0)

试试这个。

df.list <- lapply(mysamples,
                  function(x){if(any(x$func.group=="grazer")&
                                 any(x$func.group=="weevil"))
                                 return(x)})

both <- Filter(Negate(is.null),df.list)

匿名函数以df为参数,如果df中有任何象鼻虫and,则返回df$func.group;否则返回NULLlapply(...)将此功能“应用”到列表中的每个df,返回df.list,其中包含您需要的所有数据框,以及大量NULLs。第二个语句仅返回列表中不是NULL的成员,例如df$func.group中至少有一个象鼻虫和至少一个grazer的数据框列表。