1

我正在使用 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)
        }
    }
}
4

2 回答 2

1

试试这个:

test=function(x,...){
   bar=unique(mtcars[,x])
   foo=subset(mtcars,mtcars[,x]==bar[1])
   return(foo)
}
于 2016-01-28T00:32:49.103 回答
0

在@desc 所说的基础上,我设法解决了这个问题:

test <- function(...){

    bar <- unique(sampleData[,...])

    foo <- subset_samples(phyloseqObject, eval(parse(bar@names)) == as.numeric(bar[1]))
    print(foo)

}

test('pH')
于 2016-01-28T01:24:24.340 回答