Hello

I'm trying to make my own package for Ovis Aries. I would like to annotate some DEG obtain with DESeq2 or EdgeR package.

But I obtained the following error :

Error in file(description = tmp, open = "r") : cannot open the connection
In addition: Warning message:
In file(description = tmp, open = "r") :
  cannot open file 'C:/Users/cagenet/Documents/AMHAROC/3_Analyses_bioinfo/Analyse_Stat/EdgeR/gene2pubmed.gz': No such file or directory

Any help will be greatly appreciate. Thank in advance

Carine

> rm(list=ls())
> workDir <- "C:/Users/cagenet/Documents/AMHAROC/3_Analyses_bioinfo/Analyse_Stat/EdgeR"
> setwd(workDir)
> makeOrgPackageFromNCBI(version = "0.1",
+                        author = "Carine GENET <carine.genet@..fr",
+                        maintainer = "Carine GENET <carine.genet@..fr",
+                        outputDir = ".",
+                        tax_id = "9940",
+                        genus = "Ovis",
+                        species = "aries",
+                        rebuildCache = FALSE)
preparing data from NCBI ...
starting download for 6 data files
getting data for gene2pubmed.gz
Error in file(description = tmp, open = "r") : cannot open the connection
In addition: Warning message:
In file(description = tmp, open = "r") :
  cannot open file 'C:/Users/cagenet/Documents/AMHAROC/3_Analyses_bioinfo/Analyse_Stat/EdgeR/gene2pubmed.gz': No such file or directory
> sessionInfo()
R version 3.3.1 RC (2016-06-17 r70798)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows >= 8 x64 (build 9200)

locale:
[1] LC_COLLATE=French_France.1252  LC_CTYPE=French_France.1252    LC_MONETARY=French_France.1252
[4] LC_NUMERIC=C                   LC_TIME=French_France.1252    

attached base packages:
[1] parallel  stats4    stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] AnnotationHub_2.4.2    edgeR_3.14.0           limma_3.28.14          AnnotationForge_1.14.2
 [5] org.Hs.eg.db_3.3.0     AnnotationDbi_1.34.3   IRanges_2.6.1          S4Vectors_0.10.1      
 [9] Biobase_2.32.0         BiocGenerics_0.18.0    BiocInstaller_1.22.3  

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.5                   splines_3.3.1                 xtable_1.8-2                 
 [4] lattice_0.20-33               R6_2.1.2                      httr_1.2.1                   
 [7] tools_3.3.1                   grid_3.3.1                    DBI_0.4-1                    
[10] htmltools_0.3.5               digest_0.6.9                  interactiveDisplayBase_1.10.3
[13] shiny_0.13.2                  rsconnect_0.4.3               mime_0.4                     
[16] RSQLite_1.0.0                 XML_3.98-1.4                  locfit_1.5-9.1               
[19] httpuv_1.3.3                 

You might also consider using the already existing OrgDb that you can get from AnnotationHub:

> hub <- AnnotationHub()
snapshotDate(): 2016-06-06
> query(hub, c("ovis aries","orgdb"))
AnnotationHub with 1 record
# snapshotDate(): 2016-06-06
# names(): AH48021
# $dataprovider: ftp://ftp.ncbi.nlm.nih.gov/gene/DATA/
# $species: Ovis aries
# $rdataclass: OrgDb
# $title: org.Ovis_aries.eg.sqlite
# $description: NCBI gene ID based annotations about Ovis_aries
# $taxonomyid: 9940
# $genome: NCBI genomes
# $sourcetype: NCBI/UniProt
# $sourceurl: ftp://ftp.ncbi.nlm.nih.gov/gene/DATA/, ftp://ftp.uniprot.org/p...
# $sourcelastmodifieddate: NA
# $sourcesize: NA
# $tags: NCBI, Gene, Annotation
# retrieve record with 'object[["AH48021"]]'

> org.Oa.eg.db <- hub[["AH48021"]]

> org.Oa.eg.db
OrgDb object:
| DBSCHEMAVERSION: 2.1
| DBSCHEMA: NOSCHEMA_DB
| ORGANISM: Ovis aries
| SPECIES: Ovis aries
| CENTRALID: GID
| Taxonomy ID: 9940
| Db type: OrgDb
| Supporting package: AnnotationDbi

Please see: help('select') for usage information
>

Try rebuildCache = TRUE. When rebuildCache = FALSE it's assumed you've already downloaded the files necessary to build the Org package. 

Valerie