She*_*hep 1 r dataframe tapply dplyr
我想“总结” R中的一个因子变量,这样对于每条记录,我都知道存在哪些因子水平。
这是一个简化的示例数据框:
df <- data.frame(record= c("a","a","b","c","c","c"),
species = c("COD", "SCE", "COD", "COD","SCE","QSC"))
record species
a COD
a SCE
b COD
c COD
c SCE
c QSC
Run Code Online (Sandbox Code Playgroud)
这就是我想要实现的目标:
data.frame(record= c(a,b,c), species = c("COD, SCE", "COD", "COD, SCE, QSC"))
record species
a COD, SCE
b COD
c COD, SCE, QSC
Run Code Online (Sandbox Code Playgroud)
这是我所能获得的最接近的结果,但它会将所有级别的因素记录在每个记录中,而不仅仅是每个记录应显示的水平。
summarise(group_by(df, record),
species = (paste(levels(species), collapse="")))
record species
<fctr> <chr>
a CODQSCSCE <- this should be CODSCE
b CODQSCSCE <- this should just be COD
c CODQSCSCE <- this is correct as CODQSCSCE as it has all levels
Run Code Online (Sandbox Code Playgroud)
轻触返回相同的问题
tapply(df$species, df$record, function(x) paste(levels(x), collapse=""))
a b c
"CODQSCSCE" "CODQSCSCE" "CODQSCSCE"
Run Code Online (Sandbox Code Playgroud)
我需要找到一种方法来判断每种记录中存在哪些物种因子组合。
谢谢你的帮助!
用途unique():
library(dplyr)
df %>%
group_by(site) %>%
summarise(species = paste(unique(species), collapse = ', '))
# A tibble: 3 x 2
site species
<fctr> <chr>
1 a COD, SCE
2 b COD
3 c COD, SCE, QSC
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
4190 次 |
| 最近记录: |