version 1.1-5

cran · Nov 4, 2013 · 43f53bd · 43f53bd
1 parent d7f6615
commit 43f53bd
Show file tree

Hide file tree

Showing 54 changed files with 1,857 additions and 350 deletions.
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,27 +1,15 @@
 Package: agricolae
 Type: Package
 Title: Statistical Procedures for Agricultural Research
-Version: 1.1-4
-Date: 2013-4-8
+Version: 1.1-5
+Date: 2013-11-04
 Author: Felipe de Mendiburu
 Maintainer: Felipe de Mendiburu <f.mendiburu@cgiar.org>
 Suggests: akima, klaR, MASS, lme4, cluster
-Description: Agricolae was presented in the thesis "A statistical
-        analysis tool for agricultural research" to obtain the degree
-        of Master on science, mention Systems Engineering, Faculty of
-        Industrial and Systems Engineering, National Engineering
-        University (UNI), Lima-Peru Agricolae offers extensive
-        functionality on experimental design especially for
-        agricultural and plant breeding experiments, which can also be
-        useful for other purposes. It supports planning of lattice,
-        Alpha, Cyclic, Complete Block, Latin Square, Graeco-Latin
-        Squares, augmented block, split and strip plot designs. There
-        are also various analysis facilities for experimental data,
-        e.g. treatment comparison procedures and several non-parametric
-        tests comparison, biodiversity indexes and consensus cluster.
+Description: Agricolae was presented in the thesis "A statistical analysis tool for agricultural research" to obtain the degree of Master on science, mention Systems Engineering, Faculty of Industrial and Systems Engineering, National Engineering University (UNI), Lima-Peru Agricolae offers extensive functionality on experimental design especially for agricultural and plant breeding experiments, which can also be useful for other purposes. It supports planning of lattice, Alpha, Cyclic, Complete Block, Latin Square, Graeco-Latin Squares, augmented block, split and strip plot designs. There are also various analysis facilities for experimental data, e.g. treatment comparison procedures and several non-parametric tests comparison, biodiversity indexes and consensus cluster.
 License: GPL
 URL: http://tarwi.lamolina.edu.pe/~fmendiburu
-Packaged: 2013-04-09 14:25:17 UTC; fdemendiburu
+Packaged: 2013-11-04 13:09:54 UTC; fdemendiburu
 NeedsCompilation: no
 Repository: CRAN
-Date/Publication: 2013-04-09 17:19:01
+Date/Publication: 2013-11-04 14:53:40
diff --git a/MD5 b/MD5
@@ -1,15 +1,15 @@
-abdbd7ed265b9e4a83bf044e9cad185c *DESCRIPTION
+cea2c9d879e61fb253ebcae11c7c8691 *DESCRIPTION
 e79b1262fe29225e12e495d0c93a8cfd *NAMESPACE
 ac915672b4823dd46f1f67401b2e9aae *R/AMMI.R
 63bc79669decc3e830d4bc48db2f9b87 *R/AMMI.contour.R
-2b960905c916bd9564d7787516593a94 *R/BIB.test.R
-a71a16fc1f7603104d084b85435bc5e5 *R/DAU.test.R
-676b35dcf80103399c801f50e86e37bd *R/HSD.test.R
-2a9ebe175d495af2fdaf86f3f5084056 *R/LSD.test.R
-5e8f692c7d6af0542fe0b38f856b9d1d *R/PBIB.test.R
-13710a0d6b6f774c5769436af75dc468 *R/SNK.test.R
+b9c3c5f6813e4a4cf286838e888ca0d8 *R/BIB.test.R
+7e240b61c2a28f1cbe51df8b4a1f3090 *R/DAU.test.R
+c098195677cd71433eeebd51440c9321 *R/HSD.test.R
+327dd6d0acbcccb773820e0ab65c2659 *R/LSD.test.R
+38327fe3947c2c3e15b85c4d4adf3437 *R/PBIB.test.R
+de7a7868d7301e7de15636e24e461ff5 *R/SNK.test.R
 8394b67492266932ed75c813173b23d8 *R/audpc.R
-dba13680665d6c6f63f540e258238818 *R/bar.err.R
+879499bd1ecc934458b5e7d2475ebed6 *R/bar.err.R
 99a20c4d8a54834b071110ed1ad26768 *R/bar.group.R
 8a5727dca117688179b2999008027ef7 *R/carolina.R
 fb7b420d454122819738ebe35faba042 *R/consensus.R
@@ -30,37 +30,38 @@ fd0138319106167813b91f7fd80ee1f3 *R/design.lsd.R
 27e36d5280cd40704d68e382e7a4997d *R/design.rcbd.R
 be0c3e013c67afe325bb028d02999cdb *R/design.split.R
 ba1664d8911031906315890547073d57 *R/design.strip.R
-85aef6722d0f10923b98731a72563034 *R/duncan.test.r
-03b3659d4cb5056534eb02b4f1c30269 *R/durbin.test.R
+a387557878e142d7a82f580301f8de1a *R/duncan.test.r
+cf58d05146820ca2943ecfa18bf53f89 *R/durbin.test.R
 00d964a077169c76247cc389bd1d1423 *R/fact.nk.R
-2a77e5a049a94fd6b00c65263b31595a *R/friedman.R
+3825659a484fd9a062e273315d148e96 *R/friedman.R
 47e5c328d08f1a0b685fd8123f2789df *R/graph.freq.R
 c78890976271a2feb91e9a2b3eddd40c *R/grid3p.R
 e2362a1e694e883a567dd1f173d19525 *R/hcut.R
 3cf607a30269cf57b6a939c6d2d7e299 *R/hgroups.R
-b0d2a83fbdc99d129b4f71456e9a63ef *R/index.bio.R
+a1d4da2c1b53f70adede38b0d317cfb5 *R/index.bio.R
 b6aa6c034d4993c6fe74da578ec2f474 *R/index.smith.R
 67f2717076b8238a7555abeba36e6cd2 *R/intervals.freq.R
 8421821ac99f4e2e46064dea978c6c32 *R/join.freq.R
 b6203baf64a53238a80450b7f087ad79 *R/kendall.R
-b7d76faffb84fef1c113b1335b5497a2 *R/kruskal.R
+cb36090c45cc042ae7b985ea2e8473fe *R/kruskal.R
 4b424b6237f7930ee72f411deaac1876 *R/kurtosis.R
 6eb2a5ff933d46fa46cc9868e1b4fcf0 *R/lastC.R
+d40cdc0c3cec3be934b51f4c64ac7c76 *R/lateblight.R
 4bb9f4bd7404d8f9fcf7e52c64fdc90c *R/lineXtester.R
 fe9fa7a7dbb3983f6916dfec7548b9e7 *R/montecarlo.R
 97ad21836bfbd79e6b01155367e791c6 *R/nonadditivity.R
 3f30ff05f70772335fcc5a3e7711c1d2 *R/normal.freq.R
 2be0a2014c08b0da000ba1a1413281a3 *R/ojiva.freq.R
-ac326ed12c60a60cf956eb5659a8cc80 *R/order.group.R
-d48563090411ea98323e2a36a1f5d421 *R/order.stat.R
+d3965427438791b7ab400c2e456399d8 *R/order.group.R
+87fd0291ed06c27510dfb28bf5f54de2 *R/order.stat.R
 7d2166dd46e32dccda43fd1124a06acf *R/path.analysis.R
 a76382b84367fa44a352db3013017640 *R/plot.graph.freq.R
 33510c05a613f6bf250e98203f1079a5 *R/polygon.freq.R
 35a23572c169a49a5197a2c55234ac7c *R/random.ab.R
-4d60fa15a02edcd1249eca207e196349 *R/reg.homog.R
+8e136196ddb74e2e0d9dcfabcbc54c49 *R/reg.homog.R
 79f5d456dabed74f3aadbb9a787d7efc *R/resampling.cv.R
 bd0612f814df5a663e3f7e5c299595f0 *R/resampling.model.R
-46a1c93569fbd6f68f2e3ff6ffae998a *R/scheffe.test.R
+773f3e044b36f0aba026ea18da515440 *R/scheffe.test.R
 d6e3b2794952971fc5f427a55198f2fe *R/similarity.R
 a98e09f60ed734ef5924cebd5dbff750 *R/simulation.model.R
 cc4e5133b27731530e7fcda9f5748bdf *R/skewness.R
@@ -75,9 +76,10 @@ a11b0289caeb05954bd49d5266788499 *R/summary.graph.freq.R
 b7d9a984cdd216cfed7fd22ce919fbd0 *R/table.freq.R
 fccb15117eb001a7a0d5f4c2914ae287 *R/tapply.stat.R
 b07b11123ddd0420ec6ec95a6e2e7245 *R/vark.R
-51dc030121837ca8b0cc7ef8bcdc7be8 *R/waerden.test.R
+5fe98750f765d62b97c89318cbe2671d *R/waerden.test.R
 52f800d6c62727f3d1ec4dc2daedf5a6 *R/waller.R
-e4a4b255a52f4526e600aca918986268 *R/waller.test.R
+15ef67850f26d7e3260ff10c1d75f7ba *R/waller.test.R
+995b006d6389bd2f9730040f140738f4 *R/weatherSeverity.r
 1845e24b5bc2de19bbfe2dd1ec255f1c *R/wxyz.R
 3edde456b8a76e65f36c62888d86c78c *data/CIC.rda
 a5ad58e54034bbdbc9918b6e0b6aae4e *data/Chz2006.rda
@@ -114,29 +116,31 @@ ec644fe2f1fb85a8d096cc5b61888ddd *data/sinRepAmmi.rda
 f89bebe0b89dcc945e5d4850104f5570 *data/sweetpotato.rda
 93feb78eface8720dbdeff6d343a39e2 *data/wilt.rda
 06d77e44fb004716afa58983e3c0ff08 *data/yacon.rda
-8df48053c2d8e6519b752e4003beddf9 *inst/doc/ENagricolae.pdf
-e779d0ea3500a53d285211b7afed78cd *inst/doc/SPagricolae.pdf
-7206ba43e8d2e7e4e5896bb03667cff4 *man/AMMI.Rd
+3d239e241f8a5511976c134668f7ff50 *inst/doc/ENagricolae.pdf
+a7bdd20166b12fe95b9451a9635b574f *inst/doc/SPagricolae.pdf
+8860533aaab5520cf3529100f55f6187 *inst/external/severity.csv
+8ca1b77e68c92959eb509d6bc331fb0d *inst/external/weather.csv
+ad3a8bbc5a117fe1ede0a9b849b9bd0a *man/AMMI.Rd
 a996c0a1566a9620c631175230988a38 *man/AMMI.contour.Rd
-ad066e57cd843805932f970484fb5e01 *man/BIB.test.Rd
+26b40210c762bf41ca5caf934a5a6612 *man/BIB.test.Rd
 1134504b5a2f4c2619d3f7bc3a3dca8c *man/CIC.Rd
 eaed94c9f44ffd5226bfa1b6aaecfb8c *man/Chz2006.Rd
 e09abb490eb8e60ce05f99657941bb22 *man/ComasOxapampa.Rd
-d94bd3a55d2e4a9190b88f3e4e50c357 *man/DAU.test.Rd
+e3eeed465d019b026a2097ff7a498fa8 *man/DAU.test.Rd
 6744257e94b685d95f7a3e30a2cc8cf0 *man/DC.Rd
 4adceee2ea5542373b071e50cd273003 *man/Glycoalkaloids.Rd
-9c6c9792eddda384f1f02b3b0a9c4966 *man/HSD.test.Rd
+0a1ab97631948a33a72a0c200809bab8 *man/HSD.test.Rd
 872c2db9c0da0e98657fbd1439c919b8 *man/Hco2006.Rd
-cde3497722f040dcd8717be796d37ab5 *man/LSD.test.Rd
+0aafec176ec4db23784b5a281e64fff2 *man/LSD.test.Rd
 81c57d5edf1d7e2aee447026aead30da *man/LxT.Rd
-42e61db3f88713a811e097963bcc2a2e *man/PBIB.test.Rd
+f138acbd14447d57e26a209201276bcd *man/PBIB.test.Rd
 ab9296959478e556bc8f16c88ca3498c *man/RioChillon.Rd
-fb29a382767dab6dcc5ce8efdcd9f550 *man/SNK.test.Rd
-c800c773ee254efe906f5c40c9d33af6 *man/VanderWarden.Rd
-840ef103397c94ef6f394573a79f5cf6 *man/agricolae-package.Rd
+36cc5654c57818ec22e22a2a7d928ae2 *man/SNK.test.Rd
+2085a7e0ae45af87124ea2516e8a42da *man/VanderWarden.Rd
+c1b7e52ea858f71f2d65a4d090cb7084 *man/agricolae-package.Rd
 f4371b480b19a062d0366349142113db *man/audpc.Rd
-2ab99fe2399d5c6d3bfbfde1a53dba83 *man/bar.err.Rd
-6037ffdc1f5a2276b1973b71dbdf6198 *man/bar.group.Rd
+135ff25faa6e9aa661b0d1e9dd30ece0 *man/bar.err.Rd
+8d8ee493075ae261ab28c2a67f230aa5 *man/bar.group.Rd
 0fa231584f9418e8bacbde6380024cfb *man/carolina.Rd
 2ba514a973d1917a4a11c6c2b14412da *man/clay.Rd
 d9ba004a5d5c5d299eafcbf73b831ea7 *man/consensus.Rd
@@ -157,16 +161,16 @@ c498d37e3a28d7222e98fa27dae2a158 *man/design.dau.Rd
 74a6fc1f9716356cec0f77bc5d0b2299 *man/design.lattice.Rd
 19301ea4ebb7fa12de348441ece8a62d *man/design.lsd.Rd
 cd8edf98708f8e13102bf59efce7ab2b *man/design.rcbd.Rd
-cc4b7f74138f165920121cacf1939b0f *man/design.split.Rd
+4d29d5dab00c2ce0099d8ed125b89120 *man/design.split.Rd
 d60e3bb174d10db79be3c0c511663edc *man/design.strip.Rd
 2cea03b130dcc79ac02f7bd16507d65e *man/disease.Rd
-e62fe78ad7b495a58b84a9298994892d *man/duncan.test.Rd
-94ec2fc9664f1392cde8d9937b8d06a4 *man/durbin.test.Rd
+102db3f9712b73e1d8db3b2ad272d09d *man/duncan.test.Rd
+4c15d0566ce8a74aa4099dba4b4c6b6d *man/durbin.test.Rd
 82baf0121d82aad2b3d6bc4a9f009d28 *man/fact.nk.Rd
-42cf6be5b7cc8199a28294f80b5a684c *man/friedman.Rd
+4ac1f4d59abcc43acd486f85b111a763 *man/friedman.Rd
 4f201c0ad1a51f46826285abce1ef034 *man/frijol.Rd
 d110f25a7c11c055be0b2a08e7fc9e64 *man/genxenv.Rd
-c819330c5c5f3939e2519fff36009dd9 *man/graph.freq.Rd
+fd6f0a569ab3a35acb0678c329f59037 *man/graph.freq.Rd
 fc95477bcd5c333af28a7faa183da01c *man/grass.Rd
 f78032b7ef60b888495b7271674130a6 *man/greenhouse.Rd
 cb7e95125c036d7fddac34744515d491 *man/grid3p.Rd
@@ -176,14 +180,15 @@ be17faa0d1cb03a1afc14bac215cad32 *man/hcut.Rd
 419b8821f36672af21b529bb2f0bcb13 *man/heterosis.Rd
 8d42122c6e61a4f7deb2e691e5fe85dd *man/hgroups.Rd
 a6db0456a5141fc7d987dd1e451ac3ac *man/huasahuasi.Rd
-967f6c7033c52486841d1cd0440eb923 *man/index.bio.Rd
+96e33863b6646dcbdebf86714c075447 *man/index.bio.Rd
 340b27a8c896a918d415ae4ff554aa97 *man/index.smith.Rd
 8367b0c33e8cc01a9b0cca069c176246 *man/intervals.freq.Rd
 0039703e9b1f39e10c4eab3b5002c4bd *man/join.freq.Rd
 0457a41c85ad47a640a705cb2e086435 *man/kendall.Rd
-9aea6739239bc2a30aebf837c2b17cc1 *man/kruskal.Rd
+28a089d69a835b54772c03f6fac23f09 *man/kruskal.Rd
 73c4ad7dc61baa98a9ca82bef84b4942 *man/kurtosis.Rd
 ecf8b6aae6657343f66a6cac483823b0 *man/lastC.Rd
+68cda5c3b137e7e2ca5f97072edab1e3 *man/lateblight.Rd
 f71d1f504e0ff0b080e578e298c3eb8a *man/lineXtester.Rd
 cc61330de1f1e93bf2df61c890082609 *man/markers.Rd
 6da30d9aa36968804925b980b51087ce *man/melon.Rd
@@ -192,42 +197,43 @@ c9e24bc8c39bc7b9ddfc0300ad11f9db *man/natives.Rd
 b172d8394908d54d6acadf6c114c3acd *man/nonadditivity.Rd
 3fa02bad360444b789b114915f6e42d3 *man/normal.freq.Rd
 de71bc928be6c1e2b97e0c7c601b8396 *man/ojiva.freq.Rd
-9de0c0f7a4f92347ddbee07dc2f964a5 *man/order.group.Rd
-43be974a85148e83a6e56f81e4f7204a *man/order.stat.Rd
+426830028718bd738b5995630a85e729 *man/order.group.Rd
+f453f0dfbbc1e140b5931abcc13c4f9f *man/order.stat.Rd
 a20cc2ee6a10d07a551160e500e60163 *man/pamCIP.Rd
 2fa51c63ea06feb12882badb0ae7cf70 *man/paracsho.Rd
 93fd2b44325180cbd06c2872a0eb3f6b *man/path.analysis.Rd
-62e1f4b49bb1d863897ac4cbe12b1d8b *man/plot.graph.freq.Rd
+c5d66ae32b2b19c51f13978022e656a0 *man/plot.graph.freq.Rd
 666a614d408be31639980a6006166325 *man/plots.Rd
 50309c3d12fd5dff1ea549a218996f4a *man/plrv.Rd
 c8c99926cf8dee7cca599e0aa4ff3d72 *man/polygon.freq.Rd
 2ed8511b3248581a682f9898a5e3213f *man/potato.Rd
 baba6e3489bacdfd8fa7dea92e696326 *man/ralstonia.Rd
 a42fe6d9423622017fcfed1d1a8cfbea *man/random.ab.Rd
-e653fad20dd815ced3a69693684502c2 *man/reg.homog.Rd
+64817c98eb6980273df093c607a3096d *man/reg.homog.Rd
 63d57b62e51d1c067e3441bfd469f68e *man/resampling.cv.Rd
 10638a25d89d34c4bda2e2aece9b0e7b *man/resampling.model.Rd
 c45b903ab60dd151c3b86d42d833b13a *man/rice.Rd
-3294a664e4df7a80c0e36e603cfdebd3 *man/scheffe.test.Rd
+4ded27758281858fc92ba28ad305b4dc *man/scheffe.test.Rd
 a64a453df6181e073d209ffe74ea4079 *man/similarity.Rd
 9ea03515105547a094a366b3faf1fcba *man/simulation.model.Rd
 f5e5c1695de3d0c6253c62def909e4dd *man/sinRepAmmi.Rd
 058c7bdd1d7764ee231132cdb3072121 *man/skewness.Rd
 7be86c2db63d64198feb90505b32430b *man/soil.Rd
 394491ffd885da894d0e331c1a129e37 *man/sp.plot.Rd
-36f8a4ee139ab12e36bdc4dc77553b7a *man/ssp.plot.Rd
+ff30516915d31591ce7ef0cbaacb4e6e *man/ssp.plot.Rd
 ebdca7deb6522d9f103e457a34c9bf1d *man/stability.nonpar.Rd
 703d482881d10d88942e225c692e27d3 *man/stability.par.Rd
 60e1b8b24c38c241780eed798fc40c7c *man/stat.freq.Rd
-30b758e3205d916d6582d107621e6718 *man/strip.plot.Rd
+7bf3173e8cc582a8b37ea1c5fa5fa587 *man/strip.plot.Rd
 769625973fc6e95d65c5d1a6a318c0e3 *man/sturges.freq.Rd
 25fe8d237e3abeff87aa94ce9fc1e271 *man/summary.graph.freq.Rd
 7789d52e47110ede1a808524079300c6 *man/sweetpotato.Rd
 eb4447c25712d2954b0ae40f93b67c68 *man/table.freq.Rd
 a972a88fecf2fc7d71dc3ed9866ad37b *man/tapply.stat.Rd
 6d1141894359d35f36629d0a691c3474 *man/vark.Rd
 03d34f9bb1d06627fb2ce899a5a98e43 *man/waller.Rd
-3db7c6be5557baa5f19446d8642215c3 *man/waller.test.Rd
+161d9e45c18d15eb5adc04cfb840f878 *man/waller.test.Rd
+587c6d0f33a20c383fe2a551755cba43 *man/weatherSeverity.rd
 60509e2db335f9901840d28514a19168 *man/wilt.Rd
 90baa893e045dd26748c9f41c69b6fe5 *man/wxyz.Rd
 010f8deecefba432ecc5a85132a947bf *man/yacon.Rd
diff --git a/R/BIB.test.R b/R/BIB.test.R
@@ -1,5 +1,5 @@
 `BIB.test` <-
-		function (block, trt, y, test = c("lsd","tukey","duncan","waller","snk"), alpha = 0.05, group = TRUE)
+		function (block, trt, y, test = c("lsd","tukey","duncan","waller","snk"), alpha = 0.05, group = TRUE,console=FALSE)
 {
 	test<-match.arg(test)
 	block.unadj <- as.factor(block)
@@ -34,97 +34,104 @@
 	mean.adj <- mean(y,na.rm=TRUE) + Q * k/(lambda * ntr)
 	StdError.adj <- sqrt(MSerror * (1 + k * r * (ntr - 1)/(lambda *
 							ntr))/(r * ntr))
+	CV<-cv.model(model)
+	Mean<-mean(y,na.rm=TRUE)
+	if(console){
 	cat("\nANALYSIS BIB: ", name.y, "\nClass level information\n")
 	cat("\nBlock: ", unique(as.character(block)))
 	cat("\nTrt  : ", unique(as.character(trt)))
 	cat("\n\nNumber of observations: ", length(y), "\n\n")
 	print(anova(model))
-	CV<-cv.model(model)
-	Mean<-mean(y,na.rm=TRUE)
 	cat("\ncoefficient of variation:", round(CV, 1),
 			"%\n")
 	cat(name.y, "Means:", Mean, "\n\n")
 	cat(paste(name.t,",",sep="")," statistics\n\n")
+    }
 	nameTrt<-row.names(Y)
 	std.err <- sds/sqrt(r)
 	means <-data.frame( means=Y/r,mean.adj, SE=StdError.adj,r,std.err,Min.=mi,Max.=ma)
 	rownames(means)<-nameTrt
 	names(means)[1]<-name.y
-	print(means)
-
+	if(console)print(means)
 	parameter <- k/(lambda * ntr)
 	snk<-0
 	if (test == "lsd") {
 		Tprob <- qt(1 - alpha/2, DFerror)
+		if(console){
 		cat("\nLSD test")
 		cat("\nStd.diff   :", sdtdif)
 		cat("\nAlpha      :", alpha)
-		cat("\nLSD        :", Tprob * sdtdif)
+		cat("\nLSD        :", Tprob * sdtdif)}
 	}
 	if (test == "tukey") {
 		Tprob <- qtukey(1 - alpha, ntr, DFerror)
 		sdtdif<-sdtdif/sqrt(2)
+		if(console){
 		cat("\nTukey")
 		cat("\nAlpha      :", alpha)
 		cat("\nStd.err    :", sdtdif)
-		cat("\nHSD        :", Tprob * sdtdif)
+		cat("\nHSD        :", Tprob * sdtdif)}
 		parameter <- parameter/2
 	}
 	if (test == "waller") {
 		K <- 650 - 16000 * alpha + 1e+05 * alpha^2
 		Tprob <- waller(K, ntr - 1, DFerror, Fvalue)
+		if(console){
 		cat("\nWaller-Duncan K-ratio")
 		cat("\nThis test minimizes the Bayes risk under additive")
 		cat("\nloss and certain other assumptions.\n")
 		cat("\nk Ratio    : ", K)
-		cat("\nMSD        :", Tprob * sdtdif)
+		cat("\nMSD        :", Tprob * sdtdif)}
 	}
 	if (test == "snk") {
 		snk<-1
 		sdtdif<-sdtdif/sqrt(2)
-		cat("\nStudent Newman Keuls")
-		cat("\nAlpha     :", alpha)
-		cat("\nStd.err   :", sdtdif)
 		Tprob <- qtukey(1-alpha,2:ntr, DFerror)
 		SNK <- Tprob * sdtdif
 		names(SNK)<-2:ntr
+		if(console){cat("\nStudent Newman Keuls")
+		cat("\nAlpha     :", alpha)
+		cat("\nStd.err   :", sdtdif)
 		cat("\nCritical Range\n")
-		print(SNK)
+		print(SNK)}
 	}
 	if (test == "duncan") {
 		snk<-2
 		sdtdif<-sdtdif/sqrt(2)
-		cat("\nDuncan's new multiple range test")
-		cat("\nAlpha     :", alpha)
-		cat("\nStd.err   :", sdtdif)
 		Tprob <- qtukey((1-alpha)^(1:(ntr-1)),2:ntr, DFerror)
 		duncan <- Tprob * sdtdif
 		names(duncan)<-2:ntr
+		if(console){
+		cat("\nDuncan's new multiple range test")
+		cat("\nAlpha     :", alpha)
+		cat("\nStd.err   :", sdtdif)
 		cat("\n\nCritical Range\n")
-		print(duncan)
+		print(duncan)}
 	}
 	E <- lambda * ntr/(r * k)
+	if(console){
 	cat("\nParameters BIB")
 	cat("\nLambda     :", lambda)
 	cat("\ntreatmeans :", ntr)
 	cat("\nBlock size :", k)
 	cat("\nBlocks     :", b)
 	cat("\nReplication:", r, "\n")
-	cat("\nEfficiency factor", E, "\n\n<<< Book >>>\n")
+	cat("\nEfficiency factor", E, "\n\n<<< Book >>>\n")}
 	parameters<-data.frame(lambda= lambda,treatmeans=ntr,blockSize=k,blocks=b,r=r)
 	statistics<-data.frame(Mean=Mean,Efficiency=E,CV=CV)
 	rownames(parameters)<-" "
 	rownames(statistics)<-" "
 		if (group) {
+		if(console){
 		cat("\nMeans with the same letter are not significantly different.")
-		cat("\n\nComparison of treatments\n\nGroups, Treatments and means\n")
+		cat("\n\nComparison of treatments\n\nGroups, Treatments and means\n")}
 		if (snk==0)
 			groups <- order.group(names(mean.adj), as.numeric(mean.adj),
 					rep(1, ntr), MSerror, Tprob, std.err = StdError.adj,
-					parameter,sdtdif=sdtdif)
+					parameter,sdtdif=sdtdif,console=console)
 		else groups <- order.group(names(mean.adj), as.numeric(mean.adj),
 					rep(1, ntr), MSerror, Tprob, std.err = StdError.adj,
-					parameter, snk=snk,DFerror,alpha,sdtdif=sdtdif)
+					parameter, snk=snk,DFerror,alpha,sdtdif=sdtdif,console=console)
 		groups<-groups[,1:3]
 		names(groups)[2]<-"mean.adj"
 		comparison<-NULL
@@ -172,11 +179,11 @@
 			significant = abs(dif) > Tprob * sdtdif
 		tr.i <- nameTrt[comb[1, ]]
 		tr.j <- nameTrt[comb[2, ]]
-		cat("\nComparison between treatments means\n")
+		if(console)cat("\nComparison between treatments means\n")
 		if (test == "waller") comparison<-data.frame("Difference" = dif, significant)
 		else  comparison<-data.frame("Difference" = dif, pvalue=pvalue,"sig."=sig)
 		rownames(comparison)<-paste(tr.i,tr.j,sep=" - ")
-		print(comparison)
+		if(console)print(comparison)
 		groups <- NULL
 	}
 	invisible(list(parameters=parameters,statistics=statistics,comparison=comparison,means=means,groups=groups))