Hi Joshy,
In order to narrow down the problem it would be nice to have some
information about the output of your script.
First, can you add an 'str' instruction in your last loop, as below:
for( chr in 1:24)
{
units <- getUnitsOnChromosome(gi,chromosome=chr);
pos <- getPositions(gi,units = units)
unitNames <- getUnitNames(cdf,units = units)
theta <- extractMatrix(cesN,units = units)
#theta <- extractMatrix(ces,units = units)
str(theta)
rownames(theta) <- unitNames
all.data <- rbind(all.data,theta)
}
and paste the R output in your reply?
Cheers,
Pierre
On Sat, Sep 15, 2012 at 9:51 PM, joshy <[email protected]> wrote:
> Hi
> I have used the attached R script to normalize a bunch of SNP6 CEL files
> before. Recently I updated the aroama.affymetrix package and the output of
> the script does not make sense. Can anyone help me.
>
> Jg
>
> --
> 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 [email protected]
> To unsubscribe and other options, go to http://www.aroma-project.org/forum/
--
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 [email protected]
To unsubscribe and other options, go to http://www.aroma-project.org/forum/
