Parcial 3

1)

a.

##    c_agua molusco cons_o
## 1     100       A   7.16
## 2     100       A   8.26
## 3     100       A   6.78
## 4     100       A  14.00
## 5     100       A  13.60
## 6     100       A  11.10
## 7     100       A   8.93
## 8     100       A   9.66
## 9     100       B   6.14
## 10    100       B   6.14
## 11    100       B   3.68
## 12    100       B  10.00
## 13    100       B  10.40
## 14    100       B  11.60
## 15    100       B   5.49
## 16    100       B   5.80
## 17     75       A   5.20
## 18     75       A  13.20
## 19     75       A   5.20
## 20     75       A   8.39
## 21     75       A   7.18
## 22     75       A  10.40
## 23     75       A   6.37
## 24     75       A   7.18
## 25     75       B   4.47
## 26     75       B   4.95
## 27     75       B   9.96
## 28     75       B   6.49
## 29     75       B   5.75
## 30     75       B   5.44
## 31     75       B   1.80
## 32     75       B   9.90
## 33     50       A  11.11
## 34     50       A  10.50
## 35     50       A   9.74
## 36     50       A  14.60
## 37     50       A  18.80
## 38     50       A  11.11
## 39     50       A   9.74
## 40     50       A  11.80
## 41     50       B   9.63
## 42     50       B  14.50
## 43     50       B   6.38
## 44     50       B  10.20
## 45     50       B  13.40
## 46     50       B  17.70
## 47     50       B  14.50
## 48     50       B  12.30

Se observan 3 variables, concentración de agua, molusco y concentración de oxígeno, permitiendo medir el consumo de oxigeno midiendo la concentracion del mismo por unidad de peso del molusco.

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   1.800   6.312   9.700   9.305  11.232  18.800
## [1] 3.682652

Se observan los valores de desviacion estandar, que indican que los datos no estan dispersos, y el summary de los datos contenidos en la base de datos.

Se observa que el consumo de oxígeno en relacion a la concentracion del mismo tiene valores que se repiten mas en el intervalo 5 y 8, en las concentraciones de 75% y 100%. Una concentración del 50% interfiere en el consumo del oxigeno

Se observa en este gráfico la relación entre el consumo de O2 y el tipo de molusco.

En este gráfico se observa como los moluscos tipo A presentan menor variación en el consumo de oxigeno, pues tienen una variacion estandar mas grande. Los moluscos tipo B se dispersan en un rango mayor.

Se observa ahora la relación entre el consumo de O2 y el tipo de molusco, dependiendo de la concentración del agua. El molusco tipo B consume la mayor la cantidad de oxigeno en una concentracion de agua al 50%. Se observa tambien que, aunque suele tener mayor dispersion, el molusco tiene A tiende a consumir mas O2 en los otros dos valores de concentracion del agua. Ahora bien, la concentracion del agua al 75% genera un menor consumo de O2, y al 50% genera un mayor consumo de O2.

b.

## Analysis of Variance Table
## 
## Response: cons_o
##                Df Sum Sq Mean Sq F value   Pr(>F)   
## c_agua          1 102.50 102.495  9.0684 0.004298 **
## molusco         1  23.23  23.227  2.0550 0.158773   
## c_agua:molusco  1  14.38  14.378  1.2721 0.265481   
## Residuals      44 497.31  11.302                    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Call:
## lm(formula = cons_o ~ c_agua * molusco, data = BD_moluscos)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -6.8092 -2.2067 -0.1795  2.2469  7.6802 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)     13.35854    2.61316   5.112 6.65e-06 ***
## c_agua          -0.04477    0.03362  -1.332    0.190    
## moluscoB         2.63063    3.69556   0.712    0.480    
## c_agua:moluscoB -0.05363    0.04754  -1.128    0.265    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.362 on 44 degrees of freedom
## Multiple R-squared:  0.2198, Adjusted R-squared:  0.1666 
## F-statistic: 4.132 on 3 and 44 DF,  p-value: 0.01151

El valor P=0.01 indicaque la prueba tiene un parametro de confianza dentro del rango aceptable, es decir que estos datos son lógicos.

## $statistics
##   MSerror Df     Mean       CV  t.value      LSD
##   11.3025 44 9.304792 36.13105 2.015368 2.395502
## 
## $parameters
##         test p.ajusted name.t ntr alpha
##   Fisher-LSD      none c_agua   3  0.05
## 
## $means
##       cons_o      std  r       LCL       UCL  Min  Max    Q25    Q50     Q75
## 50  12.25062 3.199643 16 10.556750 13.944500 6.38 18.8 10.085 11.455 14.5000
## 75   6.99250 2.804093 16  5.298625  8.686375 1.80 13.2  5.200  6.430  8.7675
## 100  8.67125 3.000940 16  6.977375 10.365125 3.68 14.0  6.140  8.595 10.5750
## 
## $comparison
## NULL
## 
## $groups
##       cons_o groups
## 50  12.25062      a
## 100  8.67125      b
## 75   6.99250      b
## 
## attr(,"class")
## [1] "group"

Se evalúa ahora el efecto que tienen los tipos de moluscos y la concentracion del agua sobre el consumo de O2, y se concluye que el consumo de oxígeno mas efectivo entre los moluscos son aquellos pertenecientes al grupo A y que se encuentran en una concentración de agua del 50%.

Los resultados anteriores son mostrados de manera gráfica.

2)

load("/Users/sabm/Downloads/Salinidad.RData")
Salinidad
##     Biomasa   pH Salinidad    Zinc Potasio
## 1   765.280 5.00        33 16.4524 1441.67
## 2   954.017 4.70        35 13.9852 1299.19
## 3   827.686 4.20        32 15.3276 1154.27
## 4   755.072 4.40        30 17.3128 1045.15
## 5   896.176 5.55        33 22.3312  521.62
## 6  1422.836 5.50        33 12.2778 1273.02
## 7   821.069 4.25        36 17.8225 1346.35
## 8  1008.804 4.45        30 14.3516 1253.88
## 9  1306.494 4.75        38 13.6826 1242.65
## 10 1039.637 4.60        30 11.7566 1282.95
## 11 1193.223 4.10        30  9.8820  553.69
## 12  777.474 3.45        37 16.6752  494.74
## 13  818.127 3.45        33 12.3730  526.97
## 14 1203.568 4.10        36  9.4058  571.14
## 15  977.515 3.50        30 14.9302  408.64
## 16  369.823 3.25        30 31.2865  646.65
## 17  509.872 3.25        27 30.1652  514.03
## 18  448.315 3.20        29 28.5901  350.73
## 19  615.091 3.35        34 17.8795  496.29
## 20  545.538 3.30        36 18.5056  580.92
## 21  436.552 3.25        30 22.1344  535.82
## 22  465.907 3.25        28 28.6101  490.34
## 23  664.601 3.20        31 23.1908  552.39
## 24  502.466 3.20        31 24.6917  661.32
## 25  496.797 3.35        35 22.6758  672.12
## 26 2270.294 7.10        29  0.3729  525.65
## 27 2332.220 7.35        35  0.2703  563.13
## 28 2162.531 7.45        35  0.3205  497.96
## 29 2222.588 7.45        30  0.2648  458.38
## 30 2337.326 7.40        30  0.2105  498.25
## 31 1349.192 4.85        26 18.9875  936.26
## 32 1058.976 4.60        29 20.9687  894.79
## 33 1408.206 5.20        25 23.9841  941.36
## 34 1491.276 4.75        26 19.9727 1038.79
## 35 1254.872 5.20        26 21.3864  898.05
## 36 1152.341 4.55        25 23.7063  989.87
## 37  568.455 3.95        26 30.5589  951.28
## 38  612.447 3.70        26 26.8415  929.83
## 39  654.825 3.75        27 27.7292  925.42
## 40  991.829 4.15        27 21.5699  954.11
## 41 1895.942 5.60        24 19.6531  720.72
## 42 1346.880 5.35        27 20.3295  782.09
## 43 1482.793 5.50        26 19.5880  773.30
## 44 1145.643 5.50        28 20.1328  829.26
## 45 1137.193 5.40        28 19.2420  856.96

La base de datos tiene 5 variables, entre ellas la biomasa, el pH, la salinidad, el ZINC y el potasio. En el estudio a continuación, se analizaran las variables para conocer su correlación y realizar predicciones, mediante una regresión lineal múltiple.

require(psych)
## Loading required package: psych
## 
## Attaching package: 'psych'
## The following objects are masked from 'package:ggplot2':
## 
##     %+%, alpha
require(ggplot2)
tapply(Salinidad$Biomasa, Salinidad$pH, mean)
##       3.2      3.25       3.3      3.35      3.45       3.5       3.7      3.75 
##  538.4607  445.5385  545.5380  555.9440  797.8005  977.5150  612.4470  654.8250 
##      3.95       4.1      4.15       4.2      4.25       4.4      4.45      4.55 
##  568.4550 1198.3955  991.8290  827.6860  821.0690  755.0720 1008.8040 1152.3410 
##       4.6       4.7      4.75      4.85         5       5.2      5.35       5.4 
## 1049.3065  954.0170 1398.8850 1349.1920  765.2800 1331.5390 1346.8800 1137.1930 
##       5.5      5.55       5.6       7.1      7.35       7.4      7.45 
## 1350.4240  896.1760 1895.9420 2270.2940 2332.2200 2337.3260 2192.5595
plot(Salinidad$Biomasa~Salinidad$pH, col="blue", xlab= "Niveles pH", ylab="Biomasa")

En el gráfico podemos observar que a mayor pH hay mayor producción de biomasa, como es el caso de pH>7. Sin embargo, en los datos con menor pH son más comunes, por lo tanto, los niveles de biomasa se mantienen comunmento por debajo de 1500, y raramente sube por encima de 2000.

Ahora bien, para la medicion del impacto de los niveles de Zinc y producción de biomasa:

tapply(Salinidad$Biomasa, Salinidad$Zinc, mean)
##   0.2105   0.2648   0.2703   0.3205   0.3729   9.4058    9.882  11.7566 
## 2337.326 2222.588 2332.220 2162.531 2270.294 1203.568 1193.223 1039.637 
##  12.2778   12.373  13.6826  13.9852  14.3516  14.9302  15.3276  16.4524 
## 1422.836  818.127 1306.494  954.017 1008.804  977.515  827.686  765.280 
##  16.6752  17.3128  17.8225  17.8795  18.5056  18.9875   19.242   19.588 
##  777.474  755.072  821.069  615.091  545.538 1349.192 1137.193 1482.793 
##  19.6531  19.9727  20.1328  20.3295  20.9687  21.3864  21.5699  22.1344 
## 1895.942 1491.276 1145.643 1346.880 1058.976 1254.872  991.829  436.552 
##  22.3312  22.6758  23.1908  23.7063  23.9841  24.6917  26.8415  27.7292 
##  896.176  496.797  664.601 1152.341 1408.206  502.466  612.447  654.825 
##  28.5901  28.6101  30.1652  30.5589  31.2865 
##  448.315  465.907  509.872  568.455  369.823
  barplot(Salinidad$Biomasa~Salinidad$Zinc, col="pink", xlab= "Niveles zinc", ylab="Biomasa")

En el caso del Zinc, el gráfico nos muestra un aumento notorio en la producción de biomasa cuando los niveles de este elemento son más bien bajos, cercanos a cero. Sin embargo, existen niveles de zinc alrededor de 18 y 22 donde la planta produce biomasa a niveles cercanos a los que produce cuando el zinc es cercano a cero.

Al haber trabajado anteriormente con esta misma base de datos, se conoce que los niveles de salinidad en suelo no interfieren con el nivel de biomasa.

Se realizará ahora un analisis utilizando metodos de regresion lineal multiple, con el fin de conocer qué tan significativas son las variables Zinc y pH:

mod=lm(Biomasa~pH+Zinc, data=Salinidad)
summary(mod)
## 
## Call:
## lm(formula = Biomasa ~ pH + Zinc, data = Salinidad)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -467.50  -98.86   -0.22  113.45  514.41 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  -165.01     223.03  -0.740  0.46350    
## pH            330.51      32.45  10.187 6.44e-13 ***
## Zinc          -15.48       4.92  -3.147  0.00303 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 187.3 on 42 degrees of freedom
## Multiple R-squared:  0.8878, Adjusted R-squared:  0.8825 
## F-statistic: 166.2 on 2 and 42 DF,  p-value: < 2.2e-16
anova(mod)
## Analysis of Variance Table
## 
## Response: Biomasa
##           Df   Sum Sq  Mean Sq  F value    Pr(>F)    
## pH         1 11310631 11310631 322.5189 < 2.2e-16 ***
## Zinc       1   347361   347361   9.9049  0.003029 ** 
## Residuals 42  1472926    35070                       
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Con los resultados del modelo lineal se puede observar como el pH y el Zinc son variables significativas, sin embargo, el pH tiene una mayor influencia en la biomasa que el Zinc. Esto confirma los resultados obtenidos en los procedimientos anteriores. El coeficiente de correlación de estas variables (<2.2e-16) indica que es poco probable que estos datos se den al azar, es decir, efectivamente las variables nombradas anteriormente tienen influencia en el desarrollo de la biomasa. El valor R^2 (0.8825) tiene un valor cercano a 1, indicando que el modelo es bastante acertado para el caso.

Ahora se analizarán los datos teniendo en cuenta la biomasa en función de las covariables:

mod1=lm(Biomasa~pH+Zinc:pH,data=Salinidad)
anova(mod1)
## Analysis of Variance Table
## 
## Response: Biomasa
##           Df   Sum Sq  Mean Sq  F value    Pr(>F)    
## pH         1 11310631 11310631 307.7324 < 2.2e-16 ***
## pH:Zinc    1   276587   276587   7.5252  0.008907 ** 
## Residuals 42  1543700    36755                       
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
require(agricolae)
compara2=LSD.test(mod1,c("pH","Zinc"))
compara2
## $statistics
##    MSerror Df     Mean       CV  t.value      LSD
##   36754.76 42 1082.173 17.71578 2.018082 547.1552
## 
## $parameters
##         test p.ajusted  name.t ntr alpha
##   Fisher-LSD      none pH:Zinc  45  0.05
## 
## $means
##               Biomasa std r        LCL       UCL      Min      Max      Q25
## 3.2:23.1908   664.601  NA 1  277.70382 1051.4982  664.601  664.601  664.601
## 3.2:24.6917   502.466  NA 1  115.56882  889.3632  502.466  502.466  502.466
## 3.2:28.5901   448.315  NA 1   61.41782  835.2122  448.315  448.315  448.315
## 3.25:22.1344  436.552  NA 1   49.65482  823.4492  436.552  436.552  436.552
## 3.25:28.6101  465.907  NA 1   79.00982  852.8042  465.907  465.907  465.907
## 3.25:30.1652  509.872  NA 1  122.97482  896.7692  509.872  509.872  509.872
## 3.25:31.2865  369.823  NA 1  -17.07418  756.7202  369.823  369.823  369.823
## 3.3:18.5056   545.538  NA 1  158.64082  932.4352  545.538  545.538  545.538
## 3.35:17.8795  615.091  NA 1  228.19382 1001.9882  615.091  615.091  615.091
## 3.35:22.6758  496.797  NA 1  109.89982  883.6942  496.797  496.797  496.797
## 3.45:12.373   818.127  NA 1  431.22982 1205.0242  818.127  818.127  818.127
## 3.45:16.6752  777.474  NA 1  390.57682 1164.3712  777.474  777.474  777.474
## 3.5:14.9302   977.515  NA 1  590.61782 1364.4122  977.515  977.515  977.515
## 3.7:26.8415   612.447  NA 1  225.54982  999.3442  612.447  612.447  612.447
## 3.75:27.7292  654.825  NA 1  267.92782 1041.7222  654.825  654.825  654.825
## 3.95:30.5589  568.455  NA 1  181.55782  955.3522  568.455  568.455  568.455
## 4.1:9.4058   1203.568  NA 1  816.67082 1590.4652 1203.568 1203.568 1203.568
## 4.1:9.882    1193.223  NA 1  806.32582 1580.1202 1193.223 1193.223 1193.223
## 4.15:21.5699  991.829  NA 1  604.93182 1378.7262  991.829  991.829  991.829
## 4.2:15.3276   827.686  NA 1  440.78882 1214.5832  827.686  827.686  827.686
## 4.25:17.8225  821.069  NA 1  434.17182 1207.9662  821.069  821.069  821.069
## 4.4:17.3128   755.072  NA 1  368.17482 1141.9692  755.072  755.072  755.072
## 4.45:14.3516 1008.804  NA 1  621.90682 1395.7012 1008.804 1008.804 1008.804
## 4.55:23.7063 1152.341  NA 1  765.44382 1539.2382 1152.341 1152.341 1152.341
## 4.6:11.7566  1039.637  NA 1  652.73982 1426.5342 1039.637 1039.637 1039.637
## 4.6:20.9687  1058.976  NA 1  672.07882 1445.8732 1058.976 1058.976 1058.976
## 4.7:13.9852   954.017  NA 1  567.11982 1340.9142  954.017  954.017  954.017
## 4.75:13.6826 1306.494  NA 1  919.59682 1693.3912 1306.494 1306.494 1306.494
## 4.75:19.9727 1491.276  NA 1 1104.37882 1878.1732 1491.276 1491.276 1491.276
## 4.85:18.9875 1349.192  NA 1  962.29482 1736.0892 1349.192 1349.192 1349.192
## 5:16.4524     765.280  NA 1  378.38282 1152.1772  765.280  765.280  765.280
## 5.2:21.3864  1254.872  NA 1  867.97482 1641.7692 1254.872 1254.872 1254.872
## 5.2:23.9841  1408.206  NA 1 1021.30882 1795.1032 1408.206 1408.206 1408.206
## 5.35:20.3295 1346.880  NA 1  959.98282 1733.7772 1346.880 1346.880 1346.880
## 5.4:19.242   1137.193  NA 1  750.29582 1524.0902 1137.193 1137.193 1137.193
## 5.5:12.2778  1422.836  NA 1 1035.93882 1809.7332 1422.836 1422.836 1422.836
## 5.5:19.588   1482.793  NA 1 1095.89582 1869.6902 1482.793 1482.793 1482.793
## 5.5:20.1328  1145.643  NA 1  758.74582 1532.5402 1145.643 1145.643 1145.643
## 5.55:22.3312  896.176  NA 1  509.27882 1283.0732  896.176  896.176  896.176
## 5.6:19.6531  1895.942  NA 1 1509.04482 2282.8392 1895.942 1895.942 1895.942
## 7.1:0.3729   2270.294  NA 1 1883.39682 2657.1912 2270.294 2270.294 2270.294
## 7.35:0.2703  2332.220  NA 1 1945.32282 2719.1172 2332.220 2332.220 2332.220
## 7.4:0.2105   2337.326  NA 1 1950.42882 2724.2232 2337.326 2337.326 2337.326
## 7.45:0.2648  2222.588  NA 1 1835.69082 2609.4852 2222.588 2222.588 2222.588
## 7.45:0.3205  2162.531  NA 1 1775.63382 2549.4282 2162.531 2162.531 2162.531
##                   Q50      Q75
## 3.2:23.1908   664.601  664.601
## 3.2:24.6917   502.466  502.466
## 3.2:28.5901   448.315  448.315
## 3.25:22.1344  436.552  436.552
## 3.25:28.6101  465.907  465.907
## 3.25:30.1652  509.872  509.872
## 3.25:31.2865  369.823  369.823
## 3.3:18.5056   545.538  545.538
## 3.35:17.8795  615.091  615.091
## 3.35:22.6758  496.797  496.797
## 3.45:12.373   818.127  818.127
## 3.45:16.6752  777.474  777.474
## 3.5:14.9302   977.515  977.515
## 3.7:26.8415   612.447  612.447
## 3.75:27.7292  654.825  654.825
## 3.95:30.5589  568.455  568.455
## 4.1:9.4058   1203.568 1203.568
## 4.1:9.882    1193.223 1193.223
## 4.15:21.5699  991.829  991.829
## 4.2:15.3276   827.686  827.686
## 4.25:17.8225  821.069  821.069
## 4.4:17.3128   755.072  755.072
## 4.45:14.3516 1008.804 1008.804
## 4.55:23.7063 1152.341 1152.341
## 4.6:11.7566  1039.637 1039.637
## 4.6:20.9687  1058.976 1058.976
## 4.7:13.9852   954.017  954.017
## 4.75:13.6826 1306.494 1306.494
## 4.75:19.9727 1491.276 1491.276
## 4.85:18.9875 1349.192 1349.192
## 5:16.4524     765.280  765.280
## 5.2:21.3864  1254.872 1254.872
## 5.2:23.9841  1408.206 1408.206
## 5.35:20.3295 1346.880 1346.880
## 5.4:19.242   1137.193 1137.193
## 5.5:12.2778  1422.836 1422.836
## 5.5:19.588   1482.793 1482.793
## 5.5:20.1328  1145.643 1145.643
## 5.55:22.3312  896.176  896.176
## 5.6:19.6531  1895.942 1895.942
## 7.1:0.3729   2270.294 2270.294
## 7.35:0.2703  2332.220 2332.220
## 7.4:0.2105   2337.326 2337.326
## 7.45:0.2648  2222.588 2222.588
## 7.45:0.3205  2162.531 2162.531
## 
## $comparison
## NULL
## 
## $groups
##               Biomasa        groups
## 7.4:0.2105   2337.326             a
## 7.35:0.2703  2332.220             a
## 7.1:0.3729   2270.294             a
## 7.45:0.2648  2222.588             a
## 7.45:0.3205  2162.531             a
## 5.6:19.6531  1895.942            ab
## 4.75:19.9727 1491.276            bc
## 5.5:19.588   1482.793            bc
## 5.5:12.2778  1422.836           bcd
## 5.2:23.9841  1408.206           bcd
## 4.85:18.9875 1349.192          bcde
## 5.35:20.3295 1346.880           cde
## 4.75:13.6826 1306.494          cdef
## 5.2:21.3864  1254.872         cdefg
## 4.1:9.4058   1203.568        cdefgh
## 4.1:9.882    1193.223       cdefghi
## 4.55:23.7063 1152.341      cdefghij
## 5.5:20.1328  1145.643      cdefghij
## 5.4:19.242   1137.193      cdefghij
## 4.6:20.9687  1058.976     cdefghijk
## 4.6:11.7566  1039.637    cdefghijkl
## 4.45:14.3516 1008.804   cdefghijklm
## 4.15:21.5699  991.829  cdefghijklmn
## 3.5:14.9302   977.515 cdefghijklmno
## 4.7:13.9852   954.017 cdefghijklmno
## 5.55:22.3312  896.176 defghijklmnop
## 4.2:15.3276   827.686  efghijklmnop
## 4.25:17.8225  821.069  efghijklmnop
## 3.45:12.373   818.127  efghijklmnop
## 3.45:16.6752  777.474   fghijklmnop
## 5:16.4524     765.280   fghijklmnop
## 4.4:17.3128   755.072    ghijklmnop
## 3.2:23.1908   664.601     hijklmnop
## 3.75:27.7292  654.825      ijklmnop
## 3.35:17.8795  615.091       jklmnop
## 3.7:26.8415   612.447       jklmnop
## 3.95:30.5589  568.455        klmnop
## 3.3:18.5056   545.538        klmnop
## 3.25:30.1652  509.872         lmnop
## 3.2:24.6917   502.466         lmnop
## 3.35:22.6758  496.797         lmnop
## 3.25:28.6101  465.907          mnop
## 3.2:28.5901   448.315           nop
## 3.25:22.1344  436.552            op
## 3.25:31.2865  369.823             p
## 
## attr(,"class")
## [1] "group"
compara2$groups
##               Biomasa        groups
## 7.4:0.2105   2337.326             a
## 7.35:0.2703  2332.220             a
## 7.1:0.3729   2270.294             a
## 7.45:0.2648  2222.588             a
## 7.45:0.3205  2162.531             a
## 5.6:19.6531  1895.942            ab
## 4.75:19.9727 1491.276            bc
## 5.5:19.588   1482.793            bc
## 5.5:12.2778  1422.836           bcd
## 5.2:23.9841  1408.206           bcd
## 4.85:18.9875 1349.192          bcde
## 5.35:20.3295 1346.880           cde
## 4.75:13.6826 1306.494          cdef
## 5.2:21.3864  1254.872         cdefg
## 4.1:9.4058   1203.568        cdefgh
## 4.1:9.882    1193.223       cdefghi
## 4.55:23.7063 1152.341      cdefghij
## 5.5:20.1328  1145.643      cdefghij
## 5.4:19.242   1137.193      cdefghij
## 4.6:20.9687  1058.976     cdefghijk
## 4.6:11.7566  1039.637    cdefghijkl
## 4.45:14.3516 1008.804   cdefghijklm
## 4.15:21.5699  991.829  cdefghijklmn
## 3.5:14.9302   977.515 cdefghijklmno
## 4.7:13.9852   954.017 cdefghijklmno
## 5.55:22.3312  896.176 defghijklmnop
## 4.2:15.3276   827.686  efghijklmnop
## 4.25:17.8225  821.069  efghijklmnop
## 3.45:12.373   818.127  efghijklmnop
## 3.45:16.6752  777.474   fghijklmnop
## 5:16.4524     765.280   fghijklmnop
## 4.4:17.3128   755.072    ghijklmnop
## 3.2:23.1908   664.601     hijklmnop
## 3.75:27.7292  654.825      ijklmnop
## 3.35:17.8795  615.091       jklmnop
## 3.7:26.8415   612.447       jklmnop
## 3.95:30.5589  568.455        klmnop
## 3.3:18.5056   545.538        klmnop
## 3.25:30.1652  509.872         lmnop
## 3.2:24.6917   502.466         lmnop
## 3.35:22.6758  496.797         lmnop
## 3.25:28.6101  465.907          mnop
## 3.2:28.5901   448.315           nop
## 3.25:22.1344  436.552            op
## 3.25:31.2865  369.823             p

Ahora bien, finalmente los resultados del modelo indican que la mejor combinacion para obtener la mayor biomasa posible, será un pH de 7.4 y un nivel de Zinc de 0.2105.