[BioC] M vs A plot

Richard Friedman friedman at cancercenter.columbia.edu
Mon Feb 9 22:45:15 MET 2004


Dear Wolfgang and Everybody,

	I tried the span as high as .8 but it didn't cure the dip in the 
curve. I believe that higher
than .8 is not recommended (I saw this in either the microarray or the 
Loess literature, I am
not sure which

Thanks and best wishes,
Rich

On Jan 30, 2004, at 2:48 PM, Wolfgang Huber wrote:

> Hi Richard
>
> loess normalization has a parameter "span", which determines the 
> so-called bandwidth of the smoothing window. Apparently the default is 
> too small for you (leading to a too flexible regression curve), so you 
> have to make it larger.
>
> There is a lot of literature on choosing bandwidths (most of it seems 
> to involve some kind of cross-validation), but I am actually not aware 
> on recommendations for microarrays other than "by eye".
>
> Best wishes
>  Wolfgang
>
> -- 
> -------------------------------------
> Wolfgang Huber
> Division of Molecular Genome Analysis
> German Cancer Research Center
> Heidelberg, Germany
> Phone: +49 6221 424709
> Fax:   +49 6221 42524709
> Http:  www.dkfz.de/abt0840/whuber
> -------------------------------------
>
>
> Richard Friedman wrote:
>> Mick,
>> 	Thanks for the help. What concerns  me however is not a single
>> point being an outlier, but the whole loess fit to all the points 
>> leading
>> the lowess curve for a few printips to deviate significantly from 
>> being
>> a straight line practically colinear with the x-axis (abcissa). The 
>> two
>> test cases on which I learned to use marray - the apoE data that comes
>> with spot, and the swirl data that comes with marray, all had
>> significantly expressed genes - however they also had flat normalized
>> lowess curves. Significant curvature in the lowess curve leads me
>> to be concerned that the spots associated with that region of
>> the curve are improperly normalized.
>> 	Can anyone out there give me:
>> 1. Guidelines as to how flat the lowess curve should be for the
>>    data to be considered normalized.
>> 2. Advice as to what to do if the printtip normalization option
>>    in marray did not remove intensity dependence.
>> If anyone is willing to look at the M vs A curve, I would be grateful.
>> Thanks and best wishes,
>> Rich
>



More information about the Bioconductor mailing list