[BioC] Design Matrix Question
Gordon K Smyth
smyth at wehi.EDU.AU
Wed Oct 22 21:14:29 MEST 2003
This isn't a design matrix, it's a data matrix (from which a design matrix
might be computed).
Assuming the intention is to set up an design matrix for use in limma, the
recommended command is
designMatrix(targets, ref)
where 'targets' has a columns 'cy3' giving the channel1 treatment and
'cy5' giving the channel2 treatment and 'ref' is just the name of the
common reference.
In limma 1.3.0 (available next week) there is a function 'readTargets' for
reading in the target info. In the meantime, just set up the info in
Excel, write to disk as a tab-delimited file, read into R using
'read.table' with 'as.is=TRUE'.
Gordon
> this is what I would put your experiment in a design matrix
>
> Slide Treatment Dye
> 1 DIF1 -1
> 2 DIF1 1
> .
> .
> 6 DIF1 -1
> 7 DIF2 -1
> .
> .
> 12 DIF2 -1
> 13 DIF3 -1
> 14 DIF3 1
> .
> .
> 18 DIF3 -1
>
> Then you can use ANOVA to compare the different drugs.
>
> Kenny
>
>
> Kenny Ye
> Assistant Professor
> Department of Applied Math and Statistics
> SUNY at Stony Brook
> Stony Brook, New York 11794-3600
> Phone (631)632-9344
> Fax (631)632-8490
>
> On Wed, 22 Oct 2003, Jason Skelton wrote:
>
>> Kenny Ye wrote:
>>
>> >what is your treatments? and what is the level 0? would you please
>> give more details on your experiment?
>> >
>> >Kenny
>> >
>> Hi kenny
>>
>> Apologies my understanding of how a matrix works isn't necessarily
>> correct The 0 indicates that they are NOT present in a particular
>> treatment
>>
>>
>> I'm currently using R: Version 1.7.1
>>
>> DIF1,2 and 3 are different but similar drugs...............
>>
>> Slides 1-6 are treatment 1 (DIF1) Vs No treatment
>> Slide1 Cy5/Cy3 (DIF1/no treatment)
>> Slide2 Cy3/Cy5 (DIF1/no treatment)
>> Slide3 Cy5/Cy3 (DIF1/no treatment)
>> Slide4 Cy3/Cy5 (DIF1/no treatment)
>> Slide5 Cy3/Cy5 (DIF1/no treatment)
>> Slide6 Cy5/Cy3 (DIF1/no treatment)
>>
>> Slides 7-12 are treatment 2 (DIF2) Vs No treatment
>> Slide7 Cy5/Cy3 (DIF2/no treatment)
>> Slide8 Cy3/Cy5 (DIF2/no treatment)
>> Slide9 Cy5/Cy3 (DIF2/no treatment)
>> Slide10 Cy3/Cy5 (DIF2/no treatment)
>> Slide11 Cy3/Cy5 (DIF2/no treatment)
>> Slide12 Cy5/Cy3 (DIF2/no treatment)
>>
>> Slides 13-18 are treatment 3(DIF3) Vs No treatment
>> Slide13 Cy5/Cy3 (DIF3/no treatment)
>> Slide14 Cy3/Cy5 (DIF3/no treatment)
>> Slide15 Cy5/Cy3 (DIF3/no treatment)
>> Slide16 Cy3/Cy5 (DIF3/no treatment)
>> Slide17 Cy3/Cy5 (DIF3/no treatment)
>> Slide18 Cy5/Cy3 (DIF3/no treatment)
>>
>> I'd obviously like to compare across the different treatments DIF1,2
>> and 3
>>
>>
>>
>> Thanks
>>
>> Jason
>>
>>
>> >>Hi
>> >>
>> >>I'd like to construct a design matrix similar to the following but
>> have no idea how to do it
>> >>
>> >>I'm trying the model.matrix function
>> >>
>> >>Slide Treatment1 Treatment2 Treatment3
>> >>1 -1 0 0
>> >>2 1 0 0
>> >>3 -1 0 0
>> >>4 1 0 0
>> >>5 1 0 0
>> >>6 -1 0 0
>> >>7 0 -1 0
>> >>8 0 1 0
>> >>9 0 -1 0
>> >>10 0 1 0
>> >>11 0 1 0
>> >>12 0 -1 0
>> >>13 0 0 -1
>> >>14 0 0 1
>> >>15 0 0 -1
>> >>16 0 0 1
>> >>17 0 0 1
>> >>18 0 0 -1
>> >>
>> >>-1 represents Cy5/cy3
>> >>1 represents Cy3/Cy5
>> >>
>> >>
>> >>
>>
>>
>> --
>> --------------------------------
>> Jason Skelton
>> Pathogen Microarrays
>> Wellcome Trust Sanger Institute
>> Hinxton
>> Cambridge
>> CB10 1SA
>>
>> Tel +44(0)1223 834244 Ext 7123
>> Fax +44(0)1223 494919
>> --------------------------------
>>
>>
>>
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://www.stat.math.ethz.ch/mailman/listinfo/bioconductor
More information about the Bioconductor
mailing list