Window Length to Genomic Coordinates

0

Hello Guys,

I am trying to convert window length into genomic coordinates using below mentioned commands, following a post
Finding gene density from reference genome using R

library(Homo.sapiens)

human.genes = genes(TxDb.Hsapiens.UCSC.hg19.knownGene)

human.genes

GRanges object with 23056 ranges and 1 metadata column:

        seqnames                 ranges strand   |     gene_id
           <Rle>              <IRanges>  <Rle>   | <character>
      1    chr19 [ 58858172,  58874214]      -   |           1
     10     chr8 [ 18248755,  18258723]      +   |          10
    100    chr20 [ 43248163,  43280376]      -   |         100
   1000    chr18 [ 25530930,  25757445]      -   |        1000
  10000     chr1 [243651535, 244006886]      -   |       10000
    ...      ...                    ...    ... ...         ...
   9991     chr9 [114979995, 115095944]      -   |        9991
   9992    chr21 [ 35736323,  35743440]      +   |        9992
   9993    chr22 [ 19023795,  19109967]      -   |        9993
   9994     chr6 [ 90539619,  90584155]      +   |        9994
   9997    chr22 [ 50961997,  50964905]      -   |        9997

But I am not confirm how to integrate my data in .csv format containing fields as chr, start, end and copy number variations
Please have a look below:

chr.   start.  end.    logR_Copy_Number

Please suggest.


coordinates


Gene

• 134 views



Source link