On Fri, Oct 31, 2008 at 2:47 PM, Loren Engrav <engrav@u.washington.edu>wrote:

> Thank you, got it
>
>  exprs(obj) is missing from the output but shows in ?ExpressionSet and puts
> up the data so I can see it
>
> And now perhaps one final question to understanding this old data and
> progress
>
> Long long ago I did
> AA_ReadAffy <- ReadAffy()
> ABeset_gcrma < gcrma (AA_ReadAffy)
>
> Now I do
> AA_ReadAffy081031 <- ReadAffy() #on the same .cel files
> ABeset_gcrma081031 <- gcrma (AA_ReadAffy081031)
>
> Then I do
> Show(ABeset_gcrma) and exprs(ABeset_gcrma081031)
>
> And they do not match, for example
> Probe set Ssc.10026.1.A1_at
> ABeset_gcrma is 4.124689
> ABeset_gcrma081031 is 3.344443
>
>  ReadAffy changed?
>  gcrma changed?
>

See:

http://thread.gmane.org/gmane.science.biology.informatics.conductor/16664/focus=16664

The bioconductor archive is your friend.

Sean

	[[alternative HTML version deleted]]