[R] statistical test improvement of readability (was average disjunction)
Landini Massimiliano
numero.primo at tele2.it
Sat Jan 22 22:49:53 CET 2005
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.
More information about the R-help
mailing list