Question

我是R的新手。我正在使用一个脚本，该脚本根据国家/地区将atlas1006微生物组数据分为3种疾病患病率（低，中，高）。我想将每组细分：低，中，高，每组又分为男性和女性，即低男性，低女性，中男性，中女性，高男性，高女性。我也想保持3个现有组的低，中，高。这是将数据分为低，中，高的代码：

# Create mnd variable
# grabs the nationality from phyloseq
test <- get_variable(pseq, "nationality" )

# check you've installed this library
library(forcats)

# this collapses the two variables into a new variable AB in this case. The command for multiple changes is fct_collapse(x, AB = c("A","B"), DE = c("D","E"))
# test <- fct_collapse(test, LOW = c("Scandinavia","EasternEurope") , MEDIUM = c("SouthEurope","CentralEurope", "UKIE"), HIGH = c("US") )

test <- fct_collapse(test, LOW = c("EasternEurope","Scandinavia") , MEDIUM = c("SouthEurope","CentralEurope", "UKIE" ), HIGH = c("US"))

# reorder
test <- factor(test, levels = (c("LOW", "MEDIUM", "HIGH")))
levels(test)

# creates a new variable in the phyloseq called mnd
sample_data(pseq)$mnd = test

# checks that it has worked.
get_variable(pseq, "mnd")

关于如何从此处拆分性别的任何想法？

Answer 1

如果要创建单独的phyloseq对象，其中包含较高的群体和女性，则可以使用subset_samples

subset_samples(pseq, mnd=="high" & gender_variable == "female")

如何在R-微生物组上将流行人群细分为性别

1 个答案: