Hello everyone,

Now I have got peakArd1(RangedData) with 1120 rows and peakArd2 with 4665 rows 
and peak1_with_peak2_overlaps 944 row..

Question 1:Now how can I  get peak1_without_overlaps?

OQuestion 2:when I want to >write.table(peakArd1,"peakArd1")
error:as.data.frame.default(x[[i]], optional = TRUE) : 
  cannot coerce class 'structure("RangedData", package = "IRanges")' into a 
data.frame
what can I do with this problem?

 I really need your help.Thank you!


> head(peak1_with_overlaps)
RangedData with 6 rows and 1 value column across 23 spaces
        space             ranges |      strand
     <factor>          <IRanges> | <character>
0002     chr1 [2389275, 2389694] |           1
0003     chr1 [3408545, 3408894] |           1
0004     chr1 [3928820, 3929099] |           1
0005     chr1 [4509890, 4510239] |           1
0006     chr1 [4961425, 4961704] |           1
0007     chr1 [5829985, 5830369] |           1

> head(peakArd1)
RangedData with 6 rows and 15 value columns across 23 spaces
     space             ranges |      rank   strand peaklenth       FDR
  <factor>          <IRanges> | <integer> <factor> <integer> <numeric>
1     chr1 [1819580, 1819929] |      1032        +       350  0.063628
2     chr1 [2389275, 2389694] |       319        +       420  0.003003
3     chr1 [3408545, 3408894] |      1070        +       350  0.073905
4     chr1 [3928820, 3929099] |       627        +       280  0.015823
5     chr1 [4509890, 4510239] |      1043        +       350  0.067733
6     chr1 [4961425, 4961704] |       979        +       280  0.050556


_______________________________________________
Bioc-sig-sequencing mailing list
Bioc-sig-sequencing@r-project.org
https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing

Reply via email to