2013-06-05 26 views
5

我正在使用GenomicRange R軟件包。我有一個這樣的輸入的文件:顯示GenomicRange軟件包輸出中的所有行

> tabla <- read.table(file="~/Desktop/ncRNA/Data/Inputs/rRNAs/60/Infernal/test3.txt") 
> colnames(tabla) <- c('scal', 'strand', 'start', 'end') 
> all <- with(tabla, GRanges(scal, IRanges(start, end), strand)) 
> reduce(all) 

ř產生輸出:

GRanges with 24 ranges and 0 metadata columns: 
    seqnames  ranges strand 
     <Rle> <IRanges> <Rle> 
    [1] dvex108056 [ 87, 226]  + 
    [2] dvex108056 [ 59, 240]  - 
    [3] dvex114041 [6255, 6383]  - 
    [4] dvex144086 [2557, 2678]  + 
    [5] dvex144086 [2549, 2678]  - 
    ...  ...   ... ... 
    [20] dvex758265 [ 1, 172]  - 
    [21] dvex777870 [ 80, 200]  + 
    [22] dvex777870 [ 72, 200]  - 
    [23] dvex792767 [2083, 2211]  - 
    [24] dvex833399 [ 434, 545]  - 
    --- 
    seqlengths: 
    dvex108056 dvex114041 dvex144086 ... dvex777870 dvex792767 dvex833399 
      NA   NA   NA ...   NA   NA   NA 

dvex108056 + 87 206 
dvex108056 + 87 226 
dvex108056 - 101 240 
dvex108056 - 104 240 
dvex108056 - 59 188 
dvex108056 - 68 197 
dvex108056 - 70 208 
dvex108056 - 75 211 
dvex108056 - 78 217 
dvex108056 - 79 218 
dvex108056 - 84 223 
dvex108056 - 85 220 
dvex108056 - 87 226 
dvex108056 - 88 226 
dvex108056 - 88 227 
dvex108056 - 91 210 
dvex108056 - 91 230 
dvex114041 - 6255 6383 
dvex144086 + 2557 2678 
dvex144086 + 2561 2678 
dvex144086 + 2562 2678 
dvex144086 + 2564 2678 
dvex144086 - 2549 2678 
dvex186339 + 971 1108 
dvex186339 + 971 1112 
dvex186339 + 972 1112 
dvex186339 + 979 1115 
dvex209845 - 428 553 
dvex209845 - 436 553 
dvex209845 - 445 553 
dvex238049 + 435 572 
dvex238049 + 435 575 
dvex238049 + 435 576 
dvex238049 - 400 541 
dvex490175 - 5 144 
dvex564118 - 476 604 
dvex567944 + 1 121 
dvex567944 + 1 124 
dvex567944 + 1 125 
dvex567944 + 1 129 
dvex567944 + 6 125 
dvex567944 - 1 130 
dvex567944 - 6 125 
dvex600495 - 706 844 
dvex619713 - 26 165 
dvex619713 - 28 167 
dvex619713 - 34 173 
dvex619713 - 49 185 
dvex667132 + 1 119 
dvex667132 + 1 99 
dvex667132 - 1 120 
dvex667132 - 1 121 
dvex685093 - 113 242 
dvex685093 - 127 243 
dvex685093 - 144 243 
dvex700401 - 1 120 
dvex700401 - 1 130 
dvex700401 - 1 132 
dvex758265 - 1 139 
dvex758265 - 15 154 
dvex758265 - 16 155 
dvex758265 - 21 160 
dvex758265 - 36 172 
dvex758265 - 4 140 
dvex758265 - 7 146 
dvex777870 + 100 200 
dvex777870 + 104 200 
dvex777870 + 80 200 
dvex777870 + 84 200 
dvex777870 + 85 200 
dvex777870 + 87 200 
dvex777870 + 90 200 
dvex777870 + 94 200 
dvex777870 - 72 200 
dvex792767 - 2083 2211 
dvex833399 - 434 545 
dvex833399 - 435 545 
dvex833399 - 437 545 

當我使用GRANGES所有這些行減少到overlaping與這些代碼座標

我想打印由GenomicRanges生成的所有輸出。避免... ...線路。 這可能嗎?我相信這是一個data.frame結構,但我無法訪問所有記錄。

在此先感謝。

回答

5

在當前版本中有此指令

## The number of lines displayed in the 'show' method 
## are controlled with two global options. 
longGR <- c(gr[,"score"], gr2[,"score"], gr3) 
longGR 
options("showHeadLines"=7) 
options("showTailLines"=2) 
longGR 

## Revert to default values 
options("showHeadLines"=NULL) 
options("showTailLines"=NULL) 

所以例如

options(showTailLines=Inf) 

給你屏幕充滿了一切。

+0

非常感謝!那就是解決方案。來自哥倫比亞波哥大的問候。 –

0

嘗試打印它像

write.table(YourData ,"Howyouwannanameit.tsv", quote=F, coln.names=T, row.names=F)

我肯定不會有行....我認爲這正是[R並試圖向你展示了第一個和最後時行。

+0

> write.table(牀, 「Test.tsv」,報價= F,col.names = T,row.names = F) 其產生這樣的錯誤: 錯誤as.data.frame.default( x [[i]],可選= TRUE): 不能強制類「structure(」GRanges「,package =」GenomicRanges「)」data.frame –

+0

屬,el mismo錯誤: as.data.frame中的錯誤.default(x [[i]],可選= TRUE): 不能強制類「structure(」GRanges「,package =」GenomicRanges「)」data.frame –

+0

y si usas solo:> all para ver si sale toda lainformación(las 24 lineas) – Abdocia

1

將其轉換爲data.frame,它「以通常的方式」寫出來的:(訴1.12.4)上?GenomicRanges

red.all <- reduce(all) 
write.table(as.data.frame(red.all), ...)