Venu, Thanks for your kind words and feedback! Which version of ChIPpeakAnno did you use? This function has been modified extensively recently. Could you please try the new one at http://www.bioconductor.org/packages/2.12/bioc/html/ChIPpeakAnno.html? Thanks! Best regards, Julie On 4/10/13 12:19 PM, "Venu Pullabhatla" <venu.pullabhatla at="" kcl.ac.uk=""> wrote: > Dear Lihua > > Firstly many thanks for providing a wonderful tool for the scientific > community. I tried using chippeakanno software for identifying the > distribution of enriched regions with "assignChromosomeRegion". The only > difference in the way I used is to make pie charts rather boxplots. It > was great to see the output. It worked amazingly well for 2 samples > (example attached), but ended up with an error to the other two. This is > the error encountered. Google search really did not help me in this > case. Can you please help me in finding the problem and fixing it. > > Error in ann.utr5.peaks[ann.utr5.peaks$insideFeature == "upstream", ] : > selecting rows: subscript contains NAs > Calls: assignChromosomeRegion -> [ -> [ > Execution halted > > > Many thanks in advance for all the help. > > Best wishes > > Venu > > >