On Wed, Apr 4, 2012 at 7:57 AM, Fred Foucault
<frederic...@gmail.com> wrote:
> Hello,
>
> I want to get allele specific copy number from SNP6 CL files
> I´m following the vignette for ACNE.
> I installed the annotation data in \test\annotationData\chipTypes
> \GenomeWideSNP6.0
The name of that subdirectory is not correct. Make sure it is
*exactly* as shown in the vignette. See also
http://aroma-project.org/setup/
What puzzles me though is how you manage to setup the
AffymetrixCdfFile object 'cdf' in the first place. You must have done
something different than:
chipType <- "GenomeWideSNP_6";
cdf <- AffymetrixCdfFile$byChipType(chipType, tags="Full");
print(cdf);
because that should also give a similar error.
/Henrik
>
> "GenomeWideSNP_6,Full,HB20080710.acs"
> "GenomeWideSNP_6,Full,na31,hg19,HB20110328.ufl"
> "GenomeWideSNP_6,Full,na31,hg19,HB20110328.ugp"
> "GenomeWideSNP_6,Full.cdf"
>
> GenomeWideSNP_6,Full.cdf is from Affymetrix while the three others are
> from aroma
>
> here is the error :
> gi <- getGenomeInformation(cdf);
> Error in getGenomeInformation.AffymetrixCdfFile(cdf) :
> [2012-04-04 16:48:48] Exception: Failed to retrieve genome information
> for this chip type: GenomeWideSNP_6
>
> at #02. getGenomeInformation.AffymetrixCdfFile(cdf)
> - getGenomeInformation.AffymetrixCdfFile() is in environment
> 'aroma.affymetrix'
>
> at #01. getGenomeInformation(cdf)
> - getGenomeInformation() is in environment
> 'aroma.affymetrix'
>
>
> Thank you for your help
> Frederic
>
> --
> When reporting problems on aroma.affymetrix, make sure 1) to run the latest version of the package, 2) to report the output of sessionInfo() and traceback(), and 3) to post a complete code example.
>
>
> You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group with website http://www.aroma-project.org/.
> To post to this group, send email to aroma-af...@googlegroups.com
> To unsubscribe and other options, go to http://www.aroma-project.org/forum/
Did you place all the annotation files in the
annotationData/chipTypes/GenomeWideSNP_6 folder ?
It looks that this one is missing :
"GenomeWideSNP_6,Full,na31,hg19,HB20110328.ugp"
Do these file appear in the output of
path <- getPath(cdf)
list.files(path)
If not, then add them and retry.
If yes then I'm missing something.
Cheers,
Pierre
Pierre
first of all, there is a very simple reason for this, we just have to
find out why. My best guess now is that one of the files are
incomplete. Start a fresh R session and verify that you get the
exactly same *file sizes* as below:
> library("aroma.affymetrix");
> cdf <- AffymetrixCdfFile$byChipType("GenomeWideSNP_6", tags="Full");
> ds <- GenericDataFileSet$byPath(getPath(cdf));
> print(as.list(ds));
[[1]]
GenericDataFile:
Name: GenomeWideSNP_6
Tags: Full,na31,hg19,HB20110328
Full name: GenomeWideSNP_6,Full,na31,hg19,HB20110328
Pathname: annotationData/chipTypes/GenomeWideSNP_6/GenomeWideSNP_6,Full,na31,hg19,HB20110328.ufl
File size: 7.18 MB (7526452 bytes)
RAM: 0.00 MB
[[2]]
GenericDataFile:
Name: GenomeWideSNP_6
Tags: Full,na31,hg19,HB20110328
Full name: GenomeWideSNP_6,Full,na31,hg19,HB20110328
Pathname: annotationData/chipTypes/GenomeWideSNP_6/GenomeWideSNP_6,Full,na31,hg19,HB20110328.ugp
File size: 8.97 MB (9407867 bytes)
RAM: 0.00 MB
[[3]]
GenericDataFile:
Name: GenomeWideSNP_6
Tags: Full
Full name: GenomeWideSNP_6,Full
Pathname: annotationData/chipTypes/GenomeWideSNP_6/GenomeWideSNP_6,Full.CDF
File size: 470.44 MB (493291745 bytes)
RAM: 0.00 MB
[[4]]
GenericDataFile:
Name: GenomeWideSNP_6
Tags: HB20080710
Full name: GenomeWideSNP_6,HB20080710
Pathname: annotationData/chipTypes/GenomeWideSNP_6/GenomeWideSNP_6,HB20080710.acs
File size: 170.92 MB (179217531 bytes)
RAM: 0.00 MB
/Henrik
> library("aroma.affymetrix");
> sessionInfo();
/Henrik
GSE14996 dataset. Other questions will come...
Here is my R session info