[Rd] bug in summary.aovlist() with split= and (PR#13579)
rmh at temple.edu
rmh at temple.edu
Fri Mar 6 16:55:04 CET 2009
-------62a8e378fd5c9332aae960888fd28459
Content-Type: text/plain; charset=us-ascii
Content-Transfer-Encoding: 7bit
# R for Windows will not send your bug report automatically.
# Please copy the bug report (after finishing it) to
# your favorite email program and send it to
#
# r-bugs at r-project.org
#
######################################################
summary.aovlist() with split= and expand.split=TRUE
gives two different types of nonsensical results for a:b
in the Within stratum in the two different expansions of tmp3.aov.
S-Plus gives appropriate results and I attach them for comparison.
There are three attached files.
split.r source
split.rt R transcript showing nonsense results
split.st S-Plus transcript showing appropriate results
Rich
--please do not edit the information below--
Version:
platform = i386-pc-mingw32
arch = i386
os = mingw32
system = i386, mingw32
status =
major = 2
minor = 8.1
year = 2008
month = 12
day = 22
svn rev = 47281
language = R
version.string = R version 2.8.1 (2008-12-22)
Windows XP (build 2600) Service Pack 3
Locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252
Search Path:
.GlobalEnv, package:RcmdrPlugin.HH, package:Rcmdr, package:car, package:tcltk, package:fortunes,
package:VGAM, package:stats4, package:splines, package:HH, package:leaps, package:multcomp,
package:mvtnorm, package:grid, package:lattice, package:stats, package:graphics, package:datasets,
package:grDevices, package:rcom, package:rscproxy, package:utils, package:methods, RExcelEnv,
RcmdrEnv, Autoloads, package:base
-------62a8e378fd5c9332aae960888fd28459
Content-Type: text/plain;
name="split.r"
Content-Disposition: inline;
filename="split.r"
Content-Transfer-Encoding: quoted-printable
tmp <- data.frame(y=3Drnorm(48),
a=3Drep(letters[1:3]=
, 16),
b=3Drep(rep(LETTERS[1:4], each=3D3),4),
=
block=3Drep(LETTERS[5:6], each=3D24)
)
=
tmp.aov <- aov(y ~ a*b, data=3Dtmp)
summary(tmp.aov,
split=3D=
list(a=3Dlist(t=3D1,u=3D2), b=3Dlist(v=3D1,w=3D2,x=3D3)),
expan=
d.split=3DTRUE)
tmp2.aov <- aov(y ~ Error(block) + a*b, data=3Dtmp)
=
summary(tmp2.aov,
split=3Dlist(a=3Dlist(t=3D1,u=3D2), b=3Dlist(=
v=3D1,w=3D2,x=3D3)),
expand.split=3DTRUE)
summary(tmp2.aov,
=
split=3Dlist(a=3Dlist(t=3D1,u=3D2)),
expand.split=3DTRUE=
)
tmp3.aov <- aov(y ~ Error(block/a) + a*b, data=3Dtmp)
summary(tmp3=
.aov,
split=3Dlist(a=3Dlist(t=3D1,u=3D2), b=3Dlist(v=3D1,w=3D2,=
x=3D3)),
expand.split=3DTRUE)
summary(tmp3.aov,
split=
=3Dlist(a=3Dlist(t=3D1,u=3D2)),
expand.split=3DTRUE)
-------62a8e378fd5c9332aae960888fd28459
Content-Type: text/plain;
name="split.rt"
Content-Disposition: inline;
filename="split.rt"
Content-Transfer-Encoding: quoted-printable
> tmp <- data.frame(y=3Drnorm(48),
+ a=3Drep(letters[=
1:3], 16),
+ b=3Drep(rep(LETTERS[1:4], each=3D3),4),
+ block=3Drep(LETTERS[5:6], each=3D24)
+ =
)
> =
> tmp.aov <- aov(y ~ a*b, data=3Dtmp)
> summary(tmp.aov,
+ spl=
it=3Dlist(a=3Dlist(t=3D1,u=3D2), b=3Dlist(v=3D1,w=3D2,x=3D3)),
+ =
expand.split=3DTRUE)
Df Sum Sq Mean Sq F value Pr(>F)
a =
2 2.060 1.030 1.0528 0.3595
a: t 1 1.411 1.41=
1 1.4416 0.2377
a: u 1 0.650 0.650 0.6639 0.4205
b =
3 0.839 0.280 0.2859 0.8353
b: v 1 0.264 0.264 0.2=
702 0.6063
b: w 1 0.001 0.001 0.0013 0.9711
b: x 1=
0.573 0.573 0.5860 0.4489
a:b 6 2.300 0.383 0.3918 0.=
8794
a:b: t.v 1 0.556 0.556 0.5685 0.4558
a:b: u.v 1 0.99=
8 0.998 1.0203 0.3192
a:b: t.w 1 0.171 0.171 0.1747 0.6785
=
a:b: u.w 1 0.092 0.092 0.0942 0.7607
a:b: t.x 1 0.361 0.=
361 0.3685 0.5476
a:b: u.x 1 0.122 0.122 0.1246 0.7261
Residu=
als 36 35.226 0.978 =
> =
> tmp2.aov <- aov(y ~ Error(block) + a*b, data=3Dtmp)
> summary(tmp2.ao=
v,
+ split=3Dlist(a=3Dlist(t=3D1,u=3D2), b=3Dlist(v=3D1,w=3D2,x=
=3D3)),
+ expand.split=3DTRUE)
Error: block
Df S=
um Sq Mean Sq F value Pr(>F)
Residuals 1 0.57849 0.57849 =
=
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
a =
2 2.060 1.030 1.0406 0.3639
a: t 1 1.411 1.411 1.4250 =
0.2406
a: u 1 0.650 0.650 0.6563 0.4234
b 3 0.83=
9 0.280 0.2826 0.8376
b: v 1 0.264 0.264 0.2671 0.6085
=
b: w 1 0.001 0.001 0.0013 0.9713
b: x 1 0.573 0.573=
0.5793 0.4517
a:b 6 2.300 0.383 0.3873 0.8822
a:b: t.v=
1 0.556 0.556 0.5619 0.4585
a:b: u.v 1 0.998 0.998 1.0085 =
0.3222
a:b: t.w 1 0.171 0.171 0.1727 0.6803
a:b: u.w 1 0.09=
2 0.092 0.0931 0.7621
a:b: t.x 1 0.361 0.361 0.3642 0.5500
=
a:b: u.x 1 0.122 0.122 0.1232 0.7277
Residuals 35 34.647 0.990=
=
> summary(tmp2.aov,
+ split=3Dlist(a=3Dlist(t=3D1,u=3D2)),
+ =
expand.split=3DTRUE)
Error: block
Df Sum Sq Mean S=
q F value Pr(>F)
Residuals 1 0.57849 0.57849 =
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
a =
2 2.060 1.030 1.0406 0.3639
a: t 1 1.411 1.411 1.4250 0.2=
406
a: u 1 0.650 0.650 0.6563 0.4234
b 3 0.839 0=
.280 0.2826 0.8376
a:b 6 2.300 0.383 0.3873 0.8822
a:b: =
t 3 1.088 0.363 0.3663 0.7777
a:b: u 3 1.212 0.404 0.4082=
0.7480
Residuals 35 34.647 0.990 =
> =
> tmp3.aov <- aov(y ~ Error(block/a) + a*b, data=3Dtmp)
> summary(tmp3.=
aov,
+ split=3Dlist(a=3Dlist(t=3D1,u=3D2), b=3Dlist(v=3D1,w=3D2=
,x=3D3)),
+ expand.split=3DTRUE)
Error: block
Df =
Sum Sq Mean Sq F value Pr(>F)
Residuals 1 0.57849 0.57849 =
=
Error: block:a
Df Sum Sq Mean Sq F value Pr(>F)
a =
2 2.0603 1.0301 0.4883 0.6719
a: t 1 1.4106 1.4106 0.6686 0.=
4994
a: u 1 0.6497 0.6497 0.3079 0.6347
Residuals 2 4.2194 2=
.1097 =
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
b =
3 0.8392 0.2797 0.3034 0.8227
b: v 1 0.2644 0.2644 0.28=
68 0.5959
b: w 1 0.0013 0.0013 0.0014 0.9703
b: x 1 =
0.5734 0.5734 0.6219 0.4360
a:b 6 2.3002 0.3834 0.4158 0.8=
632
a:b: t.v 1 =
a:b: u.v 1 =
a:b: t.w 1 =
a:b: u.w 1 =
a:b: t.x 1 =
a:b: u.x 1 =
Residuals 33 30.4279 0.9221 =
> summary(tmp3.aov,
+ split=3Dlist(a=3Dlist(t=3D1,u=3D2)),
+ =
expand.split=3DTRUE)
Error: block
Df Sum Sq Mean S=
q F value Pr(>F)
Residuals 1 0.57849 0.57849 =
Error: block:a
Df Sum Sq Mean Sq F value Pr(>F)
a =
2 2.0603 1.0301 0.4883 0.6719
a: t 1 1.4106 1.4106 0.6686 0.=
4994
a: u 1 0.6497 0.6497 0.3079 0.6347
Residuals 2 4.2194 2=
.1097 =
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
b =
3 0.8392 0.2797 0.3034 0.8227
a:b 6 2.3002 0.3834 0.4158 =
0.8632
a:b: t 0 0.0000 =
a:b: u 0 0.0000 =
Residuals 33 30.4279 0.9221 =
> =
> version
_ =
platform i386-pc-mingw32 =
arch i386 =
os mingw32 =
system i386, mingw32 =
status =
major 2 =
minor 8.1 =
year 2008 =
month 12 =
day 22 =
svn rev 47281 =
language R =
version.string R version 2.8.1 (2008-12-22)
> =
-------62a8e378fd5c9332aae960888fd28459
Content-Type: text/plain;
name="split.st"
Content-Disposition: inline;
filename="split.st"
Content-Transfer-Encoding: quoted-printable
> tmp <- data.frame(y=3Drnorm(48),
+ a=3Drep(letters[=
1:3], 16),
+ b=3Drep(rep(LETTERS[1:4], each=3D3),4),
+ block=3Drep(LETTERS[5:6], each=3D24)
+ =
)
> =
> tmp.aov <- aov(y ~ a*b, data=3Dtmp)
> summary(tmp.aov,
+ spl=
it=3Dlist(a=3Dlist(t=3D1,u=3D2), b=3Dlist(v=3D1,w=3D2,x=3D3)),
+ =
expand.split=3DTRUE)
Df Sum of Sq Mean Sq F Value Pr=
(F) =
a 2 0.36016 0.180078 0.141799 0.8682780
a: t 1 0.0=
3839 0.038393 0.030232 0.8629384
a: u 1 0.32176 0.321762 0.253=
366 0.6177813
b 3 4.96430 1.654765 1.303014 0.2883921
=
b: v 1 4.84859 4.848586 3.817929 0.0585102
b: w 1 0.10289=
0.102887 0.081017 0.7775555
b: x 1 0.01282 0.012823 0.010097 =
0.9205175
a:b 6 3.78177 0.630295 0.496314 0.8068076
a:b: t=
.v 1 0.08038 0.080378 0.063292 0.8027964
a:b: u.v 1 0.68080 0.6=
80800 0.536083 0.4688007
a:b: t.w 1 2.29276 2.292763 1.805394 0.18=
74701
a:b: u.w 1 0.47283 0.472829 0.372321 0.5455749
a:b: t.x =
1 0.10382 0.103815 0.081747 0.7765822
a:b: u.x 1 0.15119 0.15118=
5 0.119048 0.7320772
Residuals 36 45.71826 1.269952 =
=
> =
> tmp2.aov <- aov(y ~ Error(block) + a*b, data=3Dtmp)
> summary(tmp2.ao=
v,
+ split=3Dlist(a=3Dlist(t=3D1,u=3D2), b=3Dlist(v=3D1,w=3D2,x=
=3D3)),
+ expand.split=3DTRUE)
Error: block =
Df Sum of Sq Mean Sq F Value Pr(F) =
Residuals 1 0.01256877 0.01256877 =
Error: Within =
Df Sum of Sq Mean Sq F Value Pr(F) =
a 2 0.36016 0.180078 0.137898 0.8716586
a: t 1 0.0=
3839 0.038393 0.029401 0.8648450
a: u 1 0.32176 0.321762 0.246=
395 0.6227280
b 3 4.96430 1.654765 1.267168 0.3006705
=
b: v 1 4.84859 4.848586 3.712896 0.0621432
b: w 1 0.10289=
0.102887 0.078788 0.7805996
b: x 1 0.01282 0.012823 0.009819 =
0.9216304
a:b 6 3.78177 0.630295 0.482660 0.8166555
a:b: t=
.v 1 0.08038 0.080378 0.061551 0.8055105
a:b: u.v 1 0.68080 0.6=
80800 0.521335 0.4750693
a:b: t.w 1 2.29276 2.292763 1.755727 0.19=
37442
a:b: u.w 1 0.47283 0.472829 0.362078 0.5512326
a:b: t.x =
1 0.10382 0.103815 0.079498 0.7796389
a:b: u.x 1 0.15119 0.15118=
5 0.115773 0.7356989
Residuals 35 45.70569 1.305877 =
=
> summary(tmp2.aov,
+ split=3Dlist(a=3Dlist(t=3D1,u=3D2)),
+=
expand.split=3DTRUE)
Error: block =
Df Sum of Sq Mean Sq F Value Pr(F) =
Residuals 1 0.01256877 0.01256877 =
Error: Within =
Df Sum of Sq Mean Sq F Value Pr(F) =
a 2 0.36016 0.180078 0.137898 0.8716586
a: t 1 0.038=
39 0.038393 0.029401 0.8648450
a: u 1 0.32176 0.321762 0.246395=
0.6227280
b 3 4.96430 1.654765 1.267168 0.3006705
a:=
b 6 3.78177 0.630295 0.482660 0.8166555
a:b: t 3 2.47696 0.825=
652 0.632259 0.5991524
a:b: u 3 1.30481 0.434938 0.333062 0.80147=
81
Residuals 35 45.70569 1.305877 =
> =
> tmp3.aov <- aov(y ~ Error(block/a) + a*b, data=3Dtmp)
> summary(tmp3.=
aov,
+ split=3Dlist(a=3Dlist(t=3D1,u=3D2), b=3Dlist(v=3D1,w=3D2=
,x=3D3)),
+ expand.split=3DTRUE)
Error: block =
Df Sum of Sq Mean Sq F Value Pr(F) =
Residuals 1 0.01256877 0.01256877 =
Error: a %in% block =
Df Sum of Sq Mean Sq F Value Pr(F) =
a 2 0.360156 0.180078 0.0884982 0.9186970
a: t 1 0.038=
393 0.038393 0.0188683 0.9033255
a: u 1 0.321762 0.321762 0.1581=
281 0.7293139
Residuals 2 4.069638 2.034819 =
Error: Within =
Df Sum of Sq Mean Sq F Value Pr(F) =
b 3 4.96430 1.654765 1.311538 0.2871289
b: v 1 4.8=
4859 4.848586 3.842903 0.0584446
b: w 1 0.10289 0.102887 0.081=
546 0.7769973
b: x 1 0.01282 0.012823 0.010163 0.9203096
=
a:b 6 3.78177 0.630295 0.499561 0.8040295
a:b: t.v 1 0.08038=
0.080378 0.063706 0.8022963
a:b: u.v 1 0.68080 0.680800 0.539590 =
0.4677908
a:b: t.w 1 2.29276 2.292763 1.817204 0.1868277
a:b: u=
.w 1 0.47283 0.472829 0.374756 0.5446185
a:b: t.x 1 0.10382 0.1=
03815 0.082282 0.7760218
a:b: u.x 1 0.15119 0.151185 0.119827 0.73=
14199
Residuals 33 41.63605 1.261699 =
> summary(tmp3.aov,
+ split=3Dlist(a=3Dlist(t=3D1,u=3D2)),
+=
expand.split=3DTRUE)
Error: block =
Df Sum of Sq Mean Sq F Value Pr(F) =
Residuals 1 0.01256877 0.01256877 =
Error: a %in% block =
Df Sum of Sq Mean Sq F Value Pr(F) =
a 2 0.360156 0.180078 0.0884982 0.9186970
a: t 1 0.038=
393 0.038393 0.0188683 0.9033255
a: u 1 0.321762 0.321762 0.1581=
281 0.7293139
Residuals 2 4.069638 2.034819 =
Error: Within =
Df Sum of Sq Mean Sq F Value Pr(F) =
b 3 4.96430 1.654765 1.311538 0.2871289
a:b 6 3.781=
77 0.630295 0.499561 0.8040295
a:b: t 3 2.47696 0.825652 0.654397=
0.5859203
a:b: u 3 1.30481 0.434938 0.344724 0.7931477
Residual=
s 33 41.63605 1.261699 =
> version
TIBCO Spotfire S+ Version 8.1.1 for Microsoft Windows : 200=
8 =
> =
-------62a8e378fd5c9332aae960888fd28459--
More information about the R-devel
mailing list