Last data update: 2014.03.03

R: Rocket Propellant
rocketR Documentation

Rocket Propellant

Description

Five different formulations of a rocket propellant x1 may be used in an aircrew escape systems on the observed burning rate Y. Here, each of the formulation is prepared by mixing from a batch of raw materials x2 which can support only five formulations required for the purpose of testing.

Usage

data(rocket)

Format

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

y

burning rate

batch

raw materials batch

op

experience of the operator

treat

formulation type of the propellant A B C D E

References

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

Examples

data(rocket)
matrix(rocket$treat,nrow=5)
par(mfrow=c(1,3))
plot(y~factor(op)+factor(batch)+treat,rocket)
rocket_aov <- aov(y~factor(op)+factor(batch)+treat,rocket)

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.Rd_%03d_medium.png", width=480, height=480)
> ### Name: rocket
> ### Title: Rocket Propellant
> ### Aliases: rocket
> ### Keywords: latin square design
> 
> ### ** Examples
> 
> data(rocket)
> matrix(rocket$treat,nrow=5)
     [,1] [,2] [,3] [,4] [,5]
[1,] "A"  "B"  "C"  "D"  "E" 
[2,] "B"  "C"  "D"  "E"  "A" 
[3,] "C"  "D"  "E"  "A"  "B" 
[4,] "D"  "E"  "A"  "B"  "C" 
[5,] "E"  "A"  "B"  "C"  "D" 
> par(mfrow=c(1,3))
> plot(y~factor(op)+factor(batch)+treat,rocket)
> rocket_aov <- aov(y~factor(op)+factor(batch)+treat,rocket)
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>