Diseño de Bloques Completamente al azar

Problema a analizar: Se desea probar el efecto de la marca de la enzima y el pH sobre la actividad enzimatica, teniendo como tratamiento los niveles de pH y como factor bloque las cuatro diferentes marcas (A, B, C, D). Variable de respuesta: actividad enzimatica en UI obtenidas al final de su accion en diversos pHs segun la marca.

## Warning: package 'readxl' was built under R version 4.4.2
##     Marca                PH                  AE       
##  Length:16          Length:16          Min.   :0.470  
##  Class :character   Class :character   1st Qu.:0.650  
##  Mode  :character   Mode  :character   Median :0.695  
##                                        Mean   :0.685  
##                                        3rd Qu.:0.760  
##                                        Max.   :0.840

transformacion de variables

PH=factor(PH)
Marca=factor(Marca)

modelo lineal y calculo anova

## 
## Call:
## lm(formula = AE ~ (PH + Marca))
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.04250 -0.01812 -0.00875  0.01937  0.04500 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  0.80000    0.02385  33.545 9.16e-11 ***
## PHPH6       -0.14500    0.02550  -5.687 0.000299 ***
## PHPH7        0.03250    0.02550   1.275 0.234321    
## PHPH8        0.05250    0.02550   2.059 0.069568 .  
## MarcaB      -0.16750    0.02550  -6.570 0.000103 ***
## MarcaC      -0.13500    0.02550  -5.295 0.000497 ***
## MarcaD      -0.09750    0.02550  -3.824 0.004063 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.03606 on 9 degrees of freedom
## Multiple R-squared:  0.9314, Adjusted R-squared:  0.8857 
## F-statistic: 20.37 on 6 and 9 DF,  p-value: 9.227e-05

Supuestos del Diseño Experimental

1.Normalidad

## Warning: package 'car' was built under R version 4.4.3
## Cargando paquete requerido: carData
## Warning: package 'carData' was built under R version 4.4.2
## 
##  Shapiro-Wilk normality test
## 
## data:  resid(modelo)
## W = 0.9316, p-value = 0.2584

## [1]  2 16

2.Homocedasticidad e independencia de los residuos

## 
##  Bartlett test of homogeneity of variances
## 
## data:  resid(modelo) by PH
## Bartlett's K-squared = 2.4076, df = 3, p-value = 0.4922
## 
##  Bartlett test of homogeneity of variances
## 
## data:  resid(modelo) by PH
## Bartlett's K-squared = 2.4076, df = 3, p-value = 0.4922

Pruebas Aposteriori

## Warning: package 'agricolae' was built under R version 4.4.3
## 
## Study: anova ~ "PH"
## 
## LSD t Test for AE 
## 
## Mean Square Error:  0.0013 
## 
## PH,  means and individual ( 95 %) CI
## 
##         AE        std r         se       LCL       UCL  Min  Max    Q25   Q50
## PH5 0.7000 0.06218253 4 0.01802776 0.6592184 0.7407816 0.65 0.79 0.6650 0.680
## PH6 0.5550 0.10214369 4 0.01802776 0.5142184 0.5957816 0.47 0.70 0.4925 0.525
## PH7 0.7325 0.07274384 4 0.01802776 0.6917184 0.7732816 0.68 0.84 0.6950 0.705
## PH8 0.7525 0.07320064 4 0.01802776 0.7117184 0.7932816 0.65 0.81 0.7250 0.775
##        Q75
## PH5 0.7150
## PH6 0.5875
## PH7 0.7425
## PH8 0.8025
## 
## Alpha: 0.05 ; DF Error: 9
## Critical Value of t: 2.262157 
## 
## least Significant Difference: 0.05767392 
## 
## Treatments with the same letter are not significantly different.
## 
##         AE groups
## PH8 0.7525      a
## PH7 0.7325      a
## PH5 0.7000      a
## PH6 0.5550      b