【发布时间】:2016-05-05 01:13:29
【问题描述】:
我正在使用 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 的来源:
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)
}
}
}
【问题讨论】:
-
你能用吗?:
bar <- unique(sampleData$pH) -
@desc 这适用于
unique(),但我仍然无法传递给subset()。
标签: r parameter-passing subset phyloseq