我正在使用 phyloseq 包。
test <- function( ...){
bar <- unique(sampleData[,'pH'])
foo <- subset_samples(phyloseqObject, pH == as.numeric(bar[1]@.Data))
print(foo)
}
test(pH)
我想将pH作为参数传递给test()
但unique()
不接受它是有效的。我可以将“pH”传递给test()
但subset_samples()
不接受它是有效的。我曾尝试将论点强制转换为几种不同的类型,但没有成功。
子集样本的来源:
subset_samples <- function(physeq, ...){
if( is.null(sample_data(physeq)) ){
cat("Nothing subset. No sample_data in physeq.\n")
return(physeq)
} else {
oldDF <- as(sample_data(physeq), "data.frame")
newDF <- subset(oldDF, ...)
if( class(physeq) == "sample_data" ){
return(sample_data(newDF))
} else {
sample_data(physeq) <- sample_data(newDF)
return(physeq)
}
}
}