Entering edit mode
Hi there! I need to use tximport.
I have some poblems on intallation of this package:
https://bioconductor.org/install/#troubleshoot-biocinstaller > biocLite("tximportData") BioC_mirror: https://bioconductor.org Using Bioconductor 3.2 (BiocInstaller 1.20.3), R 3.3.0 (2016-05-03). Installing package(s) ‘tximportData’ Warning message: package ‘tximportData’ is not available (for R version 3.3.0) I was able to install tximport but miss some commands..
rm(list=ls()) suppressMessages(library( "DESeq" )) suppressMessages(library("tximport")) setwd("~/Desktop/Prova_strand") dir<-"~/Desktop/Prova_strand/" metadata<-read.table("tipo.csv",header=F) files=file.path(dir,metadata$V1) txi.rsem<-tximport(files,type="rsem",reader=read.delim) sample<-data.frame(condition=factor(c("A","B"))) dds<-DESeqDataSetTxiport(txi.rsem,sample,~1) Error: could not find function "DESeqDataSetTxiport" R version 3.3.0 (2016-05-03) Platform: x86_64-pc-linux-gnu (64-bit) Running under: Ubuntu 14.04.4 LTS locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C LC_TIME=it_IT.UTF-8 [4] LC_COLLATE=en_US.UTF-8 LC_MONETARY=it_IT.UTF-8 LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=it_IT.UTF-8 LC_NAME=C LC_ADDRESS=C [10] LC_TELEPHONE=C LC_MEASUREMENT=it_IT.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] parallel stats4 stats graphics grDevices utils datasets methods base other attached packages: [1] edgeR_3.12.1 limma_3.26.9 tximport_1.1.3 [4] DESeq2_1.10.1 RcppArmadillo_0.7.100.3.1 Rcpp_0.12.5 [7] SummarizedExperiment_1.0.2 Biobase_2.30.0 GenomicRanges_1.22.4 [10] GenomeInfoDb_1.6.3 IRanges_2.4.8 S4Vectors_0.8.11 [13] BiocGenerics_0.16.1 loaded via a namespace (and not attached): [1] RColorBrewer_1.1-2 futile.logger_1.4.1 plyr_1.8.4 XVector_0.10.0 [5] futile.options_1.0.0 tools_3.3.0 zlibbioc_1.16.0 rpart_4.1-10 [9] RSQLite_1.0.0 annotate_1.48.0 gtable_0.2.0 lattice_0.20-33 [13] Matrix_1.2-6 DBI_0.4-1 gridExtra_2.2.1 genefilter_1.52.1 [17] cluster_2.0.4 locfit_1.5-9.1 grid_3.3.0 nnet_7.3-12 [21] data.table_1.9.6 AnnotationDbi_1.32.3 XML_3.98-1.4 survival_2.39-5 [25] BiocParallel_1.4.3 foreign_0.8-66 latticeExtra_0.6-28 Formula_1.2-1 [29] geneplotter_1.48.0 ggplot2_2.1.0 lambda.r_1.1.7 Hmisc_3.17-4 [33] scales_0.4.0 splines_3.3.0 xtable_1.8-2 colorspace_1.2-6 [37] acepack_1.3-3.3 munsell_0.4.3 chron_2.3-47
Someone could please help me?
hi,
This error is because the modeling cannot have transcripts of length zero included in the offset calculation. What does it mean to have a gene of length 0 after all?
I don't know the best solution here, the zeros are being produced by the upstream software, but I've recommended to others to just insert a length of 1 instead for these 0's, before moving on to DESeq2: