Rocket Propellant Example Extended

Description

In continuation of Example 13.4.7 of the Rocket Propellant data, we now have the added blocking factor in test assemblies.

Usage

data(rocket_Graeco)

Format

A data frame with 25 observations on the following 5 variables.

y

burning rate

batch

raw materials batch

op

experience of the operator

treat

formulation type of the propellant A B C D E

assembly

test assemblies a b c d e

References

Montgomery, D. C. (1976-2012). Design and Analysis of Experiments, 8e. J.Wiley.

Examples

data(rocket_Graeco)
plot(y~op+batch+treat+assembly,rocket_Graeco)
rocket.glsd.aov <- aov(y~factor(op)+factor(batch)+treat+assembly,rocket_Graeco)
summary(rocket.glsd.aov)

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(ACSWR)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/ACSWR/rocket_Graeco.Rd_%03d_medium.png", width=480, height=480)
> ### Name: rocket_Graeco
> ### Title: Rocket Propellant Example Extended
> ### Aliases: rocket_Graeco
> ### Keywords: Graeco Latin square design
> 
> ### ** Examples
> 
> data(rocket_Graeco)
> plot(y~op+batch+treat+assembly,rocket_Graeco)
> rocket.glsd.aov <- aov(y~factor(op)+factor(batch)+treat+assembly,rocket_Graeco)
> summary(rocket.glsd.aov)
              Df Sum Sq Mean Sq F value  Pr(>F)   
factor(op)     4    150   37.50   4.545 0.03293 * 
factor(batch)  4     68   17.00   2.061 0.17831   
treat          4    330   82.50  10.000 0.00334 **
assembly       4     62   15.50   1.879 0.20764   
Residuals      8     66    8.25                   
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>