Re: problems running DEseq2

Michael Love <[email protected]>
Newsgroups gmane.science.biology.informatics.conductor
Message-ID <CADqzidXrdumBcyVgStzSQxvTkimMx4NpveWDjvpfsfd2iFBDAA@mail.gmail.com>
hi Julian,

Here's some information from a previous thread:
https://stat.ethz.ch/pipermail/bioconductor/2014-July/060369.html

the brief version: not all GTF files contain sufficient/correct
information to create a txdb. the good news is that you can most
likely use the prebuilt txdb for the known genes track referenced in
that thread.

Mike

On Mon, Aug 11, 2014 at 3:35 PM, Julian Rozenberg [guest]
<[email protected]> wrote:
> rying to run deseq2.
>
> I am trying to use ref seq genes from UCSC to build a set of reference coordinates according to you manual.Refseq genes are downloaded from UCSC in .GTF format from gene and gene prediction group.
>
>>hse <- makeTranscriptDbFromGFF( "/proj/macklab/projects/JR/RNA/refGene.gtf", format="gtf" )
> Error in .parse_attrCol(attrCol, file, colnames) :
>   Some attributes do not conform to 'tag value' format
>
> Could you suggest how to proceed.
> Thank you very much in advance for you help.
>
>  -- output of sessionInfo():
>
> R version 3.0.1 (2013-05-16)
> Platform: x86_64-unknown-linux-gnu (64-bit)
>
> locale:
>  [1] LC_CTYPE=en_US.iso885915       LC_NUMERIC=C
>  [3] LC_TIME=en_US.iso885915        LC_COLLATE=en_US.iso885915
>  [5] LC_MONETARY=en_US.iso885915    LC_MESSAGES=en_US.iso885915
>  [7] LC_PAPER=C                     LC_NAME=C
>  [9] LC_ADDRESS=C                   LC_TELEPHONE=C
> [11] LC_MEASUREMENT=en_US.iso885915 LC_IDENTIFICATION=C
>
> attached base packages:
> [1] parallel  stats     graphics  grDevices utils     datasets  methods
> [8] base
>
> other attached packages:
> [1] GenomicFeatures_1.12.1 AnnotationDbi_1.22.5   Biobase_2.20.1
> [4] GenomicRanges_1.12.3   IRanges_1.20.7         BiocGenerics_0.8.0
> [7] knitr_1.5              BiocInstaller_1.12.1
>
> loaded via a namespace (and not attached):
>  [1] biomaRt_2.16.0     Biostrings_2.30.1  bitops_1.0-6       BSgenome_1.28.0
>  [5] DBI_0.2-7          evaluate_0.5.1     formatR_0.10       RCurl_1.95-4.1
>  [9] Rsamtools_1.12.3   RSQLite_0.11.4     rtracklayer_1.20.2 stats4_3.0.1
> [13] stringr_0.6.2      tools_3.0.1        XML_3.96-1.1       XVector_0.2.0
> [17] zlibbioc_1.6.0
>
>
> --
> Sent via the guest posting facility at bioconductor.org.

_______________________________________________
Bioconductor mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.