Propelente

Ejemplo Propelente

(Montgomery et al. (2001))

Se considera un motor de cohete estudiando el propelente de encendido dentro de un depósito de metal. La fuerza para separar la unión entre los componentes del combustible es la respuesta, que depende de la edad del propelente.

Variable	Descrpción
Fuerza	Fuerza para separar la union entre los componentes del combustible
Edad	edad del propelente

Lectura de datos y características de los mismos:

dat <- propellant.data <- read.csv("./propellant.dat", sep="")
dim(propellant.data)

## [1] 20  3

str(propellant.data)

## 'data.frame':    20 obs. of  3 variables:
##  $ obs   : int  1 2 3 4 5 6 7 8 9 10 ...
##  $ fuerza: num  2159 1678 2316 2061 2208 ...
##  $ edad  : num  15.5 23.8 8 17 5.5 ...

head(propellant.data)

##   obs  fuerza  edad
## 1   1 2158.70 15.50
## 2   2 1678.15 23.75
## 3   3 2316.00  8.00
## 4   4 2061.30 17.00
## 5   5 2207.50  5.50
## 6   6 1708.30 19.00

summary(propellant.data)

##       obs            fuerza          edad       
##  Min.   : 1.00   Min.   :1678   Min.   : 2.000  
##  1st Qu.: 5.75   1st Qu.:1783   1st Qu.: 7.125  
##  Median :10.50   Median :2183   Median :12.750  
##  Mean   :10.50   Mean   :2131   Mean   :13.363  
##  3rd Qu.:15.25   3rd Qu.:2342   3rd Qu.:19.625  
##  Max.   :20.00   Max.   :2654   Max.   :25.000

## 
## Call:
## lm(formula = fuerza ~ edad, data = dat)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -215.98  -50.68   28.74   66.61  106.76 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2627.822     44.184   59.48  < 2e-16
## edad         -37.154      2.889  -12.86 1.64e-10
## 
## Residual standard error: 96.11 on 18 degrees of freedom
## Multiple R-squared:  0.9018, Adjusted R-squared:  0.8964 
## F-statistic: 165.4 on 1 and 18 DF,  p-value: 1.643e-10

## Analysis of Variance Table
## 
## Response: fuerza
##           Df  Sum Sq Mean Sq F value    Pr(>F)
## edad       1 1527483 1527483  165.38 1.643e-10
## Residuals 18  166255    9236

## 
## Respuesta observada, ajustada y residuales:

##       yobs   yhat    res
## 1  2158.70 2051.9  106.8
## 2  1678.15 1745.4  -67.3
## 3  2316.00 2330.6  -14.6
## 4  2061.30 1996.2   65.1
## 5  2207.50 2423.5 -216.0
## 6  1708.30 1921.9 -213.6
## 7  1784.70 1736.1   48.6
## 8  2575.00 2534.9   40.1
## 9  2357.90 2349.2    8.7
## 10 2256.70 2219.1   37.6
## 11 2165.20 2144.8   20.4
## 12 2399.55 2488.5  -88.9
## 13 1779.80 1699.0   80.8
## 14 2336.75 2265.6   71.2
## 15 1765.30 1810.4  -45.1
## 16 2053.50 1959.1   94.4
## 17 2414.40 2404.9    9.5
## 18 2200.50 2163.4   37.1
## 19 2654.20 2553.5  100.7
## 20 1753.70 1829.0  -75.3