Skip to content
Prev 6421 / 21312 Next

[Bioc-devel] NAMESPACE question

hi Martin,

thanks for your recommendations regarding the conditional loading of 
packages. I think however, that they are not related to the problem I'm 
referring. Let me put here a reproducible example which works with 
qpgraph version 0.99.7 that I have just pushed to svn:

library(qtl)
library(qpgraph)

map <- sim.map(len=100, n.mar=10, anchor.tel=FALSE, eq.spacing=TRUE, 
include.x=FALSE)

eqtlcross <- eQTLcross(map)
eqtlcross <- addGenes(eqtlcross, 5)
eqtlcross <- addeQTL(eqtlcross, "g1", location=map[[1]][1])

sim.eqtl <- reQTLcross(eqtlcross, rho=0.5, a=1)
cross <- sim.cross(map, sim.eqtl, n.ind=100)

gstarts <-runif(5, min=range(map[[1]])[1], max=range(map[[1]])[2])

annot <- data.frame(chr=rep(names(map)[1], 5),
                     start=gstarts, end=gstarts+1,
                     strand=rep("+", 5),
                     row.names=sim.eqtl$model$Y,
                     stringsAsFactors=FALSE)

## the following is the method that triggers the
## unexpected behavior. Its last but one instruction
## in line 208 of file qpgraph/R/eQTLnetworkEstimationParam-methods.R
## is the following:
##
## geneAnnotation <- geneAnnotation[genes]
##
## and should be using the method "[" imported from GenomicRanges
## however it starts loading a number of packages to do the job

param <- eQTLnetworkEstimationParam(cross, geneAnnotation=annot, 
genome="simulatedGenome")
Loading required package: parallel

Attaching package: ?BiocGenerics?

The following objects are masked from ?package:parallel?:

[...etc...]

sessionInfo()
R version 3.1.0 (2014-04-10)
Platform: x86_64-unknown-linux-gnu (64-bit)

locale:
  [1] LC_CTYPE=en_US.UTF8       LC_NUMERIC=C 
LC_TIME=en_US.UTF8        LC_COLLATE=en_US.UTF8
  [5] LC_MONETARY=en_US.UTF8    LC_MESSAGES=en_US.UTF8 
LC_PAPER=en_US.UTF8       LC_NAME=C
  [9] LC_ADDRESS=C              LC_TELEPHONE=C 
LC_MEASUREMENT=en_US.UTF8 LC_IDENTIFICATION=C

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

other attached packages:
[1] IRanges_1.99.28     S4Vectors_0.2.4     BiocGenerics_0.11.5 
qpgraph_1.99.7      qtl_1.33-7
[6] vimcom_0.9-93       setwidth_1.0-3      colorout_1.0-2

loaded via a namespace (and not attached):
  [1] annotate_1.43.5          AnnotationDbi_1.27.16    base64enc_0.1-2 
          BatchJobs_1.4
  [5] BBmisc_1.7               Biobase_2.25.0 
BiocParallel_0.99.22     biomaRt_2.21.1
  [9] Biostrings_2.33.14       bitops_1.0-6             brew_1.0-6 
          checkmate_1.4
[13] codetools_0.2-9          DBI_0.3.1                digest_0.6.4 
         fail_1.2
[17] foreach_1.4.2            futile.logger_1.3.7 
futile.options_1.0.0     GenomeInfoDb_1.1.23
[21] GenomicAlignments_1.1.30 GenomicFeatures_1.17.19 
GenomicRanges_1.17.42    graph_1.43.0
[25] grid_3.1.0               iterators_1.0.7          lambda.r_1.1.6 
         lattice_0.20-29
[29] Matrix_1.1-4             mvtnorm_1.0-0            RCurl_1.95-4.3 
         Rgraphviz_2.9.1
[33] Rsamtools_1.17.34        RSQLite_0.11.4 
rtracklayer_1.25.17      sendmailR_1.2-1
[37] stringr_0.6.2            tools_3.1.0              XML_3.98-1.1 
         xtable_1.7-4
[41] XVector_0.5.8            zlibbioc_1.11.1


cheers,
robert.
On 10/07/2014 05:54 PM, Martin Morgan wrote: