Re: funnction in LIMMA
0
0
Entering edit mode
@gordon-smyth
Last seen 19 hours ago
WEHI, Melbourne, Australia
Dear Ping, I do not know the paper by Calinksi & Harabasz (1974), but the function 'bwss' simply does one-way analysis of variance as explained in any introductory statistics text. I did not actually intend the 'bwss' function to be used for differential expression analyses as you describe. It is there rather as a utility function for the 'anova' function for MAList objects. However you could use it for your data by setting 'x' to be a vector of log-ratios or log-intensities and 'group' to be a vector which takes value '1' for group1, '2' for group2 etc etc. I cannot tell though whether this would be a sensible analysis for your data. Best wishes Gordon At 03:07 PM 10/06/2003, you wrote: >Dear Dr. Smyth, > >I am interested in your R package LIMMA. >I want to use the function BWSS in this package, but >I am not sure the meanings of the arguments of this function > >bwss <- function(x,group). > >Your docu. explains they are: > >x: a numeric vector giving the responses. >group: a vector or factor giving the grouping variable. > >For example, I have a microaray dataset. >gene1, .., gene 10 are grouped to group 1 >gene11,...,gene 20 are grouped to group 2 >gene21,.., gene 31 are grouped to group 3. > >The tumors are t1, t2, .., t10. > >For this dataset, how can I decide these two arguments: x and group???? > > >By the way, the meanings of sums of squares between the group means (bss) >and sums of squares within the groups (wss) are the same as those >in Calinski & Harabasz (1974) CH(K)=(B(k)/(k-1)) / (W(k)/(n-k))????? > >Your kind help is highly appreciated!!!! > >Best regards, > >ping ---------------------------------------------------------------------- ----------------- Dr Gordon K Smyth, Senior Research Scientist, Bioinformatics, Walter and Eliza Hall Institute of Medical Research, 1G Royal Parade, Parkville, Vic 3050, Australia Tel: (03) 9345 2326, Fax (03) 9347 0852, Email: smyth@wehi.edu.au, www: http://www.statsci.org
limma PING limma PING • 491 views
ADD COMMENT

Login before adding your answer.

Traffic: 237 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6