[BioC] a question about how to use ChIPseqR

wang peter wng.peter at gmail.com
Sun Dec 11 22:31:08 CET 2011


hello all

i think the strandPileup function support both of data.frame and
readAligned object
as import data
so that is the coding

library(ShortRead)
library(ChIPseqR)
bamfile <- "nuclosome_MG.bam"
aln1 <- readAligned(bamfile, type="BAM")
counts <- strandPileup(aln1, chrLen=2e6, extend=1, coords="leftmost",
plot=FALSE)
Error in is.null(arg) : 'arg' is missing



> sessionInfo()
R version 2.13.2 (2011-09-30)
Platform: x86_64-unknown-linux-gnu (64-bit)


-- 
shan gao
Room 231(Dr.Fei lab)
Boyce Thompson Institute
Cornell University
Tower Road, Ithaca, NY 14853-1801
Office phone: 1-607-254-1267(day)
Official email:sg839 at cornell.edu
Facebook:http://www.facebook.com/profile.php?id=100001986532253



More information about the Bioconductor mailing list