[Bioc-devel] Problem with "RGChannelSetExtended" class form minfi

2017-04-02 Thread Mohammad Tanvir Ahamed via Bioc-devel
Hi, 


In minfi (1.20.2), in current release version Bioconductor 3.4,  the 
"RGChannelSetExtended" class is defined as 

> getClass("RGChannelSetExtended") 

Class "RGChannelSetExtended" [package "minfi"] 

Slots: 

Name:   assayData  phenoData featureData   experimentData  
annotationprotocolData  .__classVersion__ 
Class:   AssayData AnnotatedDataFrame  AnnotatedDataFrameMIAxE  
 character  AnnotatedDataFrame  Versions 

Extends: 

Class "RGChannelSet", directly 
Class "eSet", by class "RGChannelSet", distance 2 
Class "VersionedBiobase", by class "RGChannelSet", distance 3 
Class "Versioned", by class "RGChannelSet", distance 4 
> 





But in minfi (1.21.5), in the development version Bioconductor 3.5, the 
"RGChannelSetExtended" class is defined as
> getClass("RGChannelSetExtended") 

Class "RGChannelSetExtended" [package "minfi"] 

Slots: 

Name:  annotation colData  assays NAMES  elementMetadata  
metadata 
Class:  character   DataFrame  Assays character_OR_NULL   DataFrame  
list 

Extends: 
Class "RGChannelSet", directly 
Class "SummarizedExperiment", by class "RGChannelSet", distance 2 
Class "Vector", by class "RGChannelSet", distance 3 

Class "Annotated", by class "RGChannelSet", distance 4 


Now to make a new package, which class definition I need to use ? 
And also , for dependency some other packages (wateRmelon, ChAMP), in their 
development version  Bioconductor 3.5, still using the old definition of 
"RGChannelSetExtended".
If so, how to solve the problem ? 


Tanvir Ahamed 

Göteborg, Sweden  |  mashra...@yahoo.com 

___
Bioc-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/bioc-devel

Re: [Bioc-devel] Problem with champ.load ( wateRmelon::beadcount function )

2017-03-29 Thread Mohammad Tanvir Ahamed via Bioc-devel
Dear Yuan Tian, 

Thank you very much for your mail . 

I have already add a session info with the previous mail. 

Problem arise in your script "champ.load.R" line  70 ( bc=beadcount(rgSet) ). 
"beadcount" belongs to package "wateRmelon". 
In "beadcount" function  from line 3 (nb <- assayDataElement(x, "NBeads")) 
problem arise !!

I am using : 
R version 3.4.0 alpha (2017-03-28 r72427) 
Platform: x86_64-w64-mingw32/x64 (64-bit) 
Biobase (2.35.1), minfi (1.21.5), wateRmelon (1.19.3) 

Yes you are right. minfi has change their class definition for 
"RGChannelSetExtended" 


In my current version of R, 

> getClass("RGChannelSetExtended") 
Class "RGChannelSetExtended" [package "minfi"] 

Slots: 

Name:  annotation  colData  assays NAMES   elementMetadata  
metadata 
Class:  characterDataFrame  Assays character_OR_NULL   DataFrame  
list 

Extends: 
Class "RGChannelSet", directly 
Class "SummarizedExperiment", by class "RGChannelSet", distance 2 
Class "Vector", by class "RGChannelSet", distance 3 
Class "Annotated", by class "RGChannelSet", distance 4 




but in 
R version 3.3.3 (2017-03-06) 
Platform: x86_64-w64-mingw32/x64 (64-bit) 
Running under: Windows 10 x64 (build 14393) 

minfi (1.20.2) 


> getClass("RGChannelSetExtended") 
Class "RGChannelSetExtended" [package "minfi"] 

Slots: 

Name:  assayData  phenoData featureData experimentData  
annotationprotocolData  .__classVersion__ 
Class:  AssayData AnnotatedDataFrame  AnnotatedDataFrame  MIAxE 
  character  AnnotatedDataFrame  Versions 

Extends: 
Class "RGChannelSet", directly 
Class "eSet", by class "RGChannelSet", distance 2 
Class "VersionedBiobase", by class "RGChannelSet", distance 3 
Class "Versioned", by class "RGChannelSet", distance 4 
> 



No "assayData" is present in present version of minfi !!


Can anybody suggest how to solve the problem ? 

 
Tanvir Ahamed 
Göteborg, Sweden  |  mashra...@yahoo.com 




From: Tian Yuan 

Sent: Thursday, 30 March 2017, 7:35
Subject: Re: Problem with champ.load



Below is my Test Code:

> testDir=system.file("extdata",package="ChAMPdata")
> myLoad <- champ.load(testDir,arraytype="450K")
[===]
[ ChAMP.LOAD START >]
-
Loading data from 
/home/tianyuan/R/x86_64-pc-linux-gnu-library/3.3/ChAMPdata/extdata
[read.metharray.sheet] Found the following CSV files:

[1] 
"/home/tianyuan/R/x86_64-pc-linux-gnu-library/3.3/ChAMPdata/extdata/lung_test_set.csv"
Loading required package: IlluminaHumanMethylation450kmanifest
<< Read DataSet Success. >>

The fraction of failed positions per sample

(You may need to delete samples with high proportion of failed 
probes
):
   Failed CpG Fraction.
C1 0.0013429122
C2 0.0022162171
C3 0.0003563249
C4 0.0002842360
T1 0.0003831007
T2 0.0011946152
T3 0.0014953286
T4 0.0015447610
Filtering probes with a detection p-value above 0.01 in one or more samples has 
removed 2728 probes from the analysis. If a large number of probes have been 
removed, ChAMP suggests you to identify potentially bad samples.
<< Filter DetP Done. >>

Filtering probes with a beadcount <3 in at least 5% of samples, has removed 
9291 from the analysis.
<< Filter Beads Done. >>

Filtering non-cg probes, has removed 2959 from the analysis.
<< Filter NoCG Done. >>

Using general 450K SNP list for filtering.
Filtering probes with SNPs as identified in Zhou's Nucleic Acids Research 
Paper, 2016, has removed 49231 from the analysis.
<< Filter SNP Done. >>

Filtering probes that align to multiple locations as identified in Nordlund et 
al, has removed 7003 from the analysis.
<< Filter MultiHit Done. >>

Filtering probes on the X or Y chromosome has removed 9917 from the analysis.
<< Filter XY chromosome Done. >>

[Beta value is selected as output.]

Zeros in your dataset have been replaced with 0.01

The analysis will be proceed with 404383 probes and 8 samples.

[< ChAMP.LOAD END >>]
[===]
[You may want to process champ.QC() next.]


Best
Yuan Tian

On Thu, 30 Mar 2017 at 12:02 Tian Yuan  wrote:

Hello Tanvir:
>
>
>The error looks wired because seems no one reported that before, I guess it's 
>a problem caused by version of Bioba, whihc is a bug in minfi, which is used 
>by ChAMP. The problem is the latest version minfi is conflict with old version 
>Biobase. The solution is replace your current Biobase package with the latest 
>version (2.34), then it should be OK.
>
>Please try it see if it works. You are using Demo data, which should be 
>correct as I tried now.
>
>
>Best
>Yuan Tian
>

wrote:
>
>Hi,
>>I am writing the following commend and get the following error.
>>Can you ple

[Bioc-devel] Problem with champ.load (Package ChAMP)

2017-03-29 Thread Mohammad Tanvir Ahamed via Bioc-devel
Hi, 


I am writing the following command and get the following error.
Can you please look into the issue ? 

> testDir=system.file("extdata",package="ChAMPdata") 
> myLoad <- champ.load(testDir,arraytype="450K") 
[===] 
[ ChAMP.LOAD START >] 
- 
Loading data from C:/Program Files/R/R-3.4.0alpha/library/ChAMPdata/extdata 
[read.metharray.sheet] Found the following CSV files: 

[1] "C:/Program 
Files/R/R-3.4.0alpha/library/ChAMPdata/extdata/lung_test_set.csv" 
Loading required package: IlluminaHumanMethylation450kmanifest 
<< Read DataSet Success. >> 

The fraction of failed positions per sample 

(You may need to delete samples with high proportion of failed probes 
): 
Failed CpG Fraction. 
C1        0.0013429122 
C2        0.0022162171 
C3        0.0003563249 
C4        0.0002842360 
T1        0.0003831007 
T2        0.0011946152 
T3        0.0014953286 
T4        0.0015447610 
Filtering probes with a detection p-value above 0.01 in one or more samples has 
removed 2728 probes from the analysis. If a large number of probes have been 
removed, ChAMP suggests you to identify potentially bad samples. 
<< Filter DetP Done. >> 

Error in (function (classes, fdef, mtable)  : 
unable to find an inherited method for function 'assayData' for signature 
'"RGChannelSetExtended"' 
> 

> sessionInfo()

R version 3.4.0 alpha (2017-03-28 r72427) 
Platform: x86_64-w64-mingw32/x64 (64-bit) 
Running under: Windows >= 8 x64 (build 9200) 

Matrix products: default 

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

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

other attached packages: 
[1] IlluminaHumanMethylation450kmanifest_0.4.0 ChAMP_2.6.0 
[3] IlluminaHumanMethylationEPICmanifest_0.3.0 
Illumina450ProbeVariants.db_1.11.0 
[5] DMRcate_1.11.2                            DMRcatedata_1.10.1 
[7] DSS_2.15.0                                bsseq_1.11.0 
[9] FEM_3.3.1                                  graph_1.53.0 
[11] org.Hs.eg.db_3.4.0                        impute_1.49.0 
[13] igraph_1.0.1                              corrplot_0.77 
[15] marray_1.53.0                              limma_3.31.19 
[17] Matrix_1.2-9                              AnnotationDbi_1.37.4 
[19] ChAMPdata_2.6.0                            minfi_1.21.5 
[21] bumphunter_1.15.0                          locfit_1.5-9.1 
[23] iterators_1.0.8                            foreach_1.4.3 
[25] Biostrings_2.43.6                          XVector_0.15.2 
[27] SummarizedExperiment_1.5.7                DelayedArray_0.1.7 
[29] matrixStats_0.51.0                        Biobase_2.35.1 
[31] GenomicRanges_1.27.23                      GenomeInfoDb_1.11.9 
[33] IRanges_2.9.19                            S4Vectors_0.13.15 
[35] BiocGenerics_0.21.3 

loaded via a namespace (and not attached): 
[1] R.utils_2.5.0                                      
IlluminaHumanMethylationEPICanno.ilm10b2.hg19_0.6.0 
[3] RSQLite_1.1-2                                      htmlwidgets_0.8 
[5] trimcluster_0.1-2                                  grid_3.4.0 
[7] BiocParallel_1.9.5                                  munsell_0.4.3 
[9] codetools_0.2-15                                    preprocessCore_1.37.0 
[11] statmod_1.4.29                                      colorspace_1.3-2 
[13] fastICA_1.2-0                                      BiocInstaller_1.25.3 
[15] knitr_1.15.1                                        robustbase_0.92-7 
[17] JADE_2.0-0                                          isva_1.9 
[19] GenomeInfoDbData_0.99.0                            biovizBase_1.23.2 
[21] diptest_0.75-7                                      R6_2.2.0 
[23] doParallel_1.0.10                                  illuminaio_0.17.0 
[25] clue_0.3-53                                        flexmix_2.3-13 
[27] bitops_1.0-6                                        reshape_0.8.6 
[29] assertthat_0.1                                      
IlluminaHumanMethylation450kanno.ilmn12.hg19_0.6.0 
[31] scales_0.4.1                                        nnet_7.3-12 
[33] gtable_0.2.0                                        methylumi_2.21.0 
[35] sva_3.23.0                                          ensembldb_1.99.12 
[37] genefilter_1.57.0                                  rtracklayer_1.35.9 
[39] lazyeval_0.2.0                                      acepack_1.4.1 
[41] GEOquery_2.41.0                                    dichromat_2.0-0 
[43] checkmate_1.8.2                                    yaml_2.1.14 
[45] reshape2_1.4.2                                      
GenomicFeatures_1.27.10 
[47] backports_1.0.5                                    httpuv_1.3.3 
[49] qvalue_2.7.0                                        Hmisc_4.0-2 
[51] tools_3.4.0                                   

[Bioc-devel] Problem to install bioconductor for R devel (3.5)

2017-03-28 Thread Mohammad Tanvir Ahamed via Bioc-devel
Hi, 

I have install development version of r from 

https://cran.r-project.org/bin/windows/base/rdevel.html

But when i have tried to install  Bioconductor, i got following error.
Can any one please explain why ?

It said "Bioconductor does not yet support R version 3.5.0"  but why ?? 

##
> ## try http:// if https:// URLs are not supported 
> source("https://bioconductor.org/biocLite.R";) 
Bioconductor version 3.5 (BiocInstaller 1.25.3), ?biocLite for help 
Bioconductor does not yet support R version 3.5.0 
> biocLite() 
BioC_mirror: https://bioconductor.org 
Using Bioconductor 3.5 (BiocInstaller 1.25.3), R Under development (unstable) 
(2017-03-27 r72423). 
Installing package(s) ‘Biobase’, ‘IRanges’, ‘AnnotationDbi’ 
also installing the dependency ‘S4Vectors’ 

Warning: unable to access index for repository 
https://bioconductor.org/packages/3.5/bioc/bin/windows/contrib/3.5: 
cannot open URL 
'https://bioconductor.org/packages/3.5/bioc/bin/windows/contrib/3.5/PACKAGES' 
Warning: unable to access index for repository 
https://bioconductor.org/packages/3.5/data/annotation/bin/windows/contrib/3.5: 
cannot open URL 
'https://bioconductor.org/packages/3.5/data/annotation/bin/windows/contrib/3.5/PACKAGES'
 
Warning: unable to access index for repository 
https://bioconductor.org/packages/3.5/data/experiment/bin/windows/contrib/3.5: 
cannot open URL 
'https://bioconductor.org/packages/3.5/data/experiment/bin/windows/contrib/3.5/PACKAGES'
 
Warning: unable to access index for repository 
https://bioconductor.org/packages/3.5/extra/bin/windows/contrib/3.5: 
cannot open URL 
'https://bioconductor.org/packages/3.5/extra/bin/windows/contrib/3.5/PACKAGES' 
Error in readRDS(dest) : unknown input format 
> 



Tanvir Ahamed 

___
Bioc-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/bioc-devel