Gerrit Eichner
2014-Mar-18 12:28 UTC
[R] format of output of residual standard errors of manova()
Dear all,
I wasn't successful in finding any related "bug" or
"wish" report on
bugs.r-project.org, svn.R-project.org/R/trunk/doc/NEWS.Rd, or RSeek
regarding the following:
The code of stats:::print.aov contains the two commands
cat("Residual standard error: ", sapply(sqrt(ss/rdf), format),
"\n", sep =
"")
and
cat("Residual standard error: ", sapply(rs, format), "\n",
sep = "")
which provide the output of the RSEs. However, the RSE-values are glued
together because of sep = "" in the call to cat(), which makes them
hard
or even impossible to read off free of doubt. (See the respective output
of npk.aov2 in example( manova) or of fit in example( summary.manova).)
Is it intended to change those lines of code into something like
cat("Residual standard error: ", paste( sapply(sqrt(ss/rdf), format),
collapse = "; "), "\n", sep = "")
and
cat("Residual standard error: ", paste( sapply(rs, format), collapse =
";
"), "\n", sep = "")
respectively? (Of course, any other sep-value would be "accepted". :))
If not, I would suggest this modification as a "wish" (and was
wondering
if my next step should indeed be submitting a wish report on
bugs.r-project.org ...)
Best regards -- Gerrit
PS:
> sessionInfo()
R version 3.0.2 (2013-09-25)
Platform: x86_64-w64-mingw32/x64 (64-bit)
locale:
[1] LC_COLLATE=German_Germany.1252 LC_CTYPE=German_Germany.1252
[3] LC_MONETARY=German_Germany.1252 LC_NUMERIC=C
[5] LC_TIME=German_Germany.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] fortunes_1.5-2
loaded via a namespace (and not attached):
[1] tools_3.0.2
---------------------------------------------------------------------
Dr. Gerrit Eichner Mathematical Institute, Room 212
gerrit.eichner at math.uni-giessen.de Justus-Liebig-University Giessen
Tel: +49-(0)641-99-32104 Arndtstr. 2, 35392 Giessen, Germany
Fax: +49-(0)641-99-32109 http://www.uni-giessen.de/cms/eichner