manta.path <- system.file("extdata","PapaGO-BLAST.results-diatoms.Rdata", package="manta")
load(manta.path)
x <- calcNormFactors(x)
conditions <- factor(x=1:2, labels=c('ambient','plusFe')); names(conditions) <- c('ref','obs')
nr(x, meta.lev='genus_sp', pair=conditions)
Run the code above in your browser using DataLab