[aroma.affymetrix] Re: extract \hat\sigma values for each mean in “regions.xls”

2009-09-03 Thread Max Moldovan
Thanks a lot Henrik. Your info is really helpful and I made everything running. Just as a suggestion, these \hat\sigma_{\delta} estimates could be either stored somewhere for each cbs sub-region, or simlply added as one more column to “regions.xls”. Best wishes Max On Sep 3, 2:27 pm, Henrik Be

[aroma.affymetrix] Re: extract \hat\sigma values for each mean in “regions.xls”

2009-09-02 Thread Henrik Bengtsson
Hi. On Wed, Sep 2, 2009 at 8:56 PM, Max Moldovan wrote: > > Thanks for your information Henrik. > > 1. I can extract sigma estimates using you first suggestion: > > rawCNs <- extractRawCopyNumbers(cbs, array=3, chromosome=9); > sigma <- estimateStandardDeviation(rawCNs); > > A colleague of mine r

[aroma.affymetrix] Re: extract \hat\sigma values for each mean in “regions.xls”

2009-09-02 Thread Max Moldovan
Thanks for your information Henrik. 1. I can extract sigma estimates using you first suggestion: rawCNs <- extractRawCopyNumbers(cbs, array=3, chromosome=9); sigma <- estimateStandardDeviation(rawCNs); A colleague of mine recommended an alternative way apparently leading to the same results (in

[aroma.affymetrix] Re: extract \hat\sigma values for each mean in “regions.xls”

2009-09-02 Thread Henrik Bengtsson
Hi. On Tue, Sep 1, 2009 at 10:04 PM, Max Moldovan wrote: > > Hi People, > > I have the following sequence of commands: > > chipType <- "GenomeWideSNP_5"; > cdfTags <- "Full,r2"; > cdf <- AffymetrixCdfFile$byChipType(chipType, tags=cdfTags) > print(cdf) > gi <- getGenomeInformation(cdf) > print(gi