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: