# Re: [R] 2^2 problem revisited

From: Rolf Turner <r.turner_at_auckland.ac.nz>
Date: Fri, 14 Nov 2008 11:42:02 +1300

Rolf Turner

On 14/11/2008, at 10:44 AM, Edna Bell wrote:

> Dear R gurus:
>
> Here is the following from Montgomery's Design and Analysis of
> Experiments, 5th edition.
>
>> str(rout1.df)
> 'data.frame': 16 obs. of 3 variables:
> \$ resp: num 18.2 18.9 12.9 14.4 27.2 24 22.4 22.5 15.9 14.5 ...
> \$ A : Factor w/ 2 levels "-1","1": 1 1 1 1 2 2 2 2 1 1 ...
> \$ B : Factor w/ 2 levels "-1","1": 1 1 1 1 1 1 1 1 2 2 ...
>> rout1.df
> resp A B
> 1 18.2 -1 -1
> 2 18.9 -1 -1
> 3 12.9 -1 -1
> 4 14.4 -1 -1
> 5 27.2 1 -1
> 6 24.0 1 -1
> 7 22.4 1 -1
> 8 22.5 1 -1
> 9 15.9 -1 1
> 10 14.5 -1 1
> 11 15.1 -1 1
> 12 14.2 -1 1
> 13 41.0 1 1
> 14 43.9 1 1
> 15 36.3 1 1
> 16 39.9 1 1
>> rout1.aov <- aov(resp~A*B,data=rout1.df)
>> summary(rout1.aov)
> Df Sum Sq Mean Sq F value Pr(>F)
> A 1 1107.23 1107.23 185.252 1.175e-08 ***
> B 1 227.26 227.26 38.023 4.826e-05 ***
> A:B 1 303.63 303.63 50.801 1.201e-05 ***
> Residuals 12 71.72 5.98
> ---
> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>
>
> When you test on interaction, you reject (of course).
>
> Now, I thought that you could not test on the main effects, A and B.
>
> Thanks,
> Edna Bell
>
> ______________________________________________
> R-help_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> guide.html
> and provide commented, minimal, self-contained, reproducible code.

```######################################################################
```
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}

R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Thu 13 Nov 2008 - 22:43:56 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 13 Nov 2008 - 23:30:24 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.