[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