Landini Massimiliano
2005-Jan-22 21:49 UTC
[R] statistical test improvement of readability (was average disjunction)
Dear all ReadeRs I was finding a quick method to improve test readability adding or constructing (with your help....) one or more function that allow what follow. Please consider Trt<-c(1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4,5,5,5,5) Block<-c(1,2,3,4,1,2,3,4,1,2,3,4,1,2,3,4,1,2,3,4) AD2DAT<-c(1.44,2.32,1.68,1.28,0.12,0.08,0.24,0.52,0.12,0.08,0.16,0.28,0.16,0.08,0.12,0.12,0.16,0.08,0.32,0.76) Tab<-cbind(Trt,Block,AD2DAT);Tab AD2DAT.aov.lm<-aov(lm(AD2DAT~as.factor(Block)+as.factor(Trt),method="qr")) TukeyHSD(AD2DAT.aov.lm,"as.factor(Trt)",ordered=TRUE,conf.level=.95) TukeyHSD(AD2DAT.aov.lm,"as.factor(Trt)",ordered=TRUE,conf.level=.99) TabMEANS<-aggregate(Tab,list(Trt),FUN=mean);TabMEANS Group.1 Trt Block AD2DAT 1 1 2.5 1.68 2 2 2.5 0.24 3 3 2.5 0.16 4 4 2.5 0.12 5 5 2.5 0.33 So I'm thinking that was more readable something like one function that permit me to juxtapose tiny letters when we are talking about .95 confidence level and capital letters when we are talking about .99 confidence level so results will be approx this Group.1 Trt Block AD2DAT .95 .99 1 1 2.5 1.68 a A 2 2 2.5 0.24 b B 3 3 2.5 0.16 b B 4 4 2.5 0.12 b B 5 5 2.5 0.33 b B Does anyone have suggestion??? Thanks in advance!! ------------------------------------------------------------------------------------------------------------------------- Landini dr. Massimiliano Tel. mob. (+39) 347 140 11 94 Tel./Fax. (+39) 051 762 196 e-mail: numero (dot) primo (at) tele2 (dot) it ------------------------------------------------------------------------------------------------------------------------- Legge di Hanggi: Pi? stupida ? la tua ricerca, pi? verr? letta e approvata. Corollario alla Legge di Hanggi: Pi? importante ? la tua ricerca, meno verr? capita.