Ejercicio 1.
Se realizó un estudio sobre factores estructurales relacionados con la permeabilidad de 30 tipos de tejidos. Los datos están disponibles en el archivo airpermeabilitywoven.csv y Clay et al 2007.pdf. Realice una lectura de la publicación para comprender que se efectuó en el trabajo.
La metadata de las variables en el archivo son las siguientes:
warp (warp yarn density, ends/cm) weft (weft " " , picks/cm) mass (mass per unit area, gr/m^2) airPermMn (average air permeability, cm2/s/cm2) airPermSD (Std Dev of air perm measurements) waterCont (Water Content after vacuum drying, %)
Pruebe los siguientes modelos de regresión lineal múltiple.
Model 1: E{airPermMn} = B0 + B1Warp + B2Weft + B3*Mass
Model 2: E{airPermMn} = C0 + C1Water + c2(Water^2)
Entregue los componentes necesarios para justificar cual modelo utilizaría para explicar la permeabilidad de los tejidos.
datos1 = read.csv("C:/Users/GOLLO/Desktop/airpermeabilitywoven.csv", sep = ",", dec = ".")
datos1## sampNum warp weft mass airPermMn airPermSD waterCont
## 1 1 54 20 101.0 45.74 1.80 53.65
## 2 2 54 25 110.3 27.02 2.44 45.65
## 3 3 54 30 119.5 15.68 0.68 40.34
## 4 4 54 35 127.3 8.76 0.39 38.21
## 5 5 54 40 136.7 4.37 0.15 35.66
## 6 6 54 45 144.1 2.90 0.40 32.35
## 7 7 57 20 101.2 47.94 1.83 50.53
## 8 8 57 25 111.1 27.52 0.78 44.91
## 9 9 57 30 120.5 14.84 0.97 40.93
## 10 10 57 35 130.7 8.99 0.48 37.49
## 11 11 57 40 142.0 3.98 0.24 33.62
## 12 12 57 45 146.3 2.68 0.12 32.33
## 13 13 60 20 107.9 33.98 1.16 50.19
## 14 14 60 25 118.7 17.01 0.79 45.68
## 15 15 60 30 129.5 9.75 0.41 39.31
## 16 16 60 35 141.4 4.56 0.41 36.00
## 17 17 60 40 151.2 2.12 0.07 33.80
## 18 18 60 45 153.7 1.70 0.05 32.04
## 19 19 63 20 109.0 27.68 0.55 47.89
## 20 20 63 25 119.1 15.24 0.64 44.52
## 21 21 63 30 129.1 7.23 0.13 40.20
## 22 22 63 35 148.5 3.22 0.42 35.72
## 23 23 63 40 150.3 1.89 0.08 33.40
## 24 24 63 42 154.4 1.61 0.05 32.31
## 25 25 66 20 111.2 27.28 1.68 47.77
## 26 26 66 25 123.0 14.42 1.39 43.59
## 27 27 66 30 134.0 6.90 0.21 37.66
## 28 28 66 35 144.4 3.19 0.17 37.09
## 29 29 66 40 155.0 1.65 0.06 35.69
## 30 30 66 42 160.2 1.38 0.05 29.49is.data.frame(datos1)## [1] TRUE#Modelo 1:
regre1 <- lm(formula = airPermMn ~ warp + weft + mass, data = datos1)
regre1##
## Call:
## lm(formula = airPermMn ~ warp + weft + mass, data = datos1)
##
## Coefficients:
## (Intercept) warp weft mass
## 105.035 -0.392 -0.680 -0.355summary(regre1)##
## Call:
## lm(formula = airPermMn ~ warp + weft + mass, data = datos1)
##
## Residuals:
## Min 1Q Median 3Q Max
## -6.888 -4.247 -0.656 2.405 14.770
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 105.035 16.355 6.42 8.4e-07 ***
## warp -0.392 0.569 -0.69 0.50
## weft -0.680 0.728 -0.93 0.36
## mass -0.355 0.366 -0.97 0.34
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 5.59 on 26 degrees of freedom
## Multiple R-squared: 0.84, Adjusted R-squared: 0.822
## F-statistic: 45.6 on 3 and 26 DF, p-value: 1.68e-10## Respuesta:
#Modelo 1: E{airPermMn} = 105.035 + -0.392 + -0.680 + -0.3
#Modelo2:
regre2 <- lm(formula = airPermMn ~ waterCont + (waterCont^2), data = datos1)
regre2 ##
## Call:
## lm(formula = airPermMn ~ waterCont + (waterCont^2), data = datos1)
##
## Coefficients:
## (Intercept) waterCont
## -63.04 1.92## Respuesta:
#Modelo 2: E{airPermMn} = -63.04 + 1.92 + 1.92^2Ejercicio 2.
A continuación se entrega un breve resumen del experimento explicado en Mimic 2004.pdf
“Experimental Results involving transesterification of rapeseed oil in 3 alcohols (methanol, ethanol, and 1-propanol). Three factors were varied: temperature (250, 300, 350C), pressure (8, 10, 12 MPa), and time (7,15,20,…,60,65 min). Yields of each alcohol type (%) and Energy consumption of each alcohol type (kW h/kg) were reported.”
Archivo disponible en biodiseltransect.csv.
Intente replicar los resultados de la Tabla 4.
datos2 = read.csv("C:/Users/GOLLO/Desktop/biodiesel_transest.csv", sep = ",", dec = ".")
datos2## run_id b_temp b_pres b_time meth_pct meth_ec eth_pct eth_ec prop1_pct
## 1 1 350 8 7 73.4 0.743 44.4 1.252 38.3
## 2 2 350 8 15 76.5 0.758 62.1 0.950 54.2
## 3 3 350 8 20 78.8 0.759 66.2 0.919 57.2
## 4 4 350 8 25 80.9 0.766 68.3 0.923 59.9
## 5 5 350 8 30 82.2 0.783 70.1 0.933 61.3
## 6 6 350 8 35 83.1 0.801 70.9 0.953 61.7
## 7 7 350 8 40 84.3 0.815 71.1 0.981 62.1
## 8 8 350 8 45 84.9 0.834 72.2 0.995 62.5
## 9 9 350 8 50 84.1 0.868 72.9 1.015 65.3
## 10 10 350 8 55 83.9 0.898 73.3 1.041 68.8
## 11 11 350 8 60 82.9 0.932 72.1 1.085 70.1
## 12 12 350 8 65 82.6 0.948 70.2 1.144 69.1
## 13 13 350 10 7 76.1 0.730 51.3 1.103 43.2
## 14 14 350 10 15 87.5 0.674 77.3 0.776 70.9
## 15 15 350 10 20 90.5 0.672 81.1 0.762 76.2
## 16 16 350 10 25 89.1 0.707 82.8 0.773 77.3
## 17 17 350 10 30 87.2 0.750 84.6 0.785 80.2
## 18 18 350 10 35 87.0 0.793 83.1 0.825 81.5
## 19 19 350 10 40 86.5 0.836 82.9 0.866 83.3
## 20 20 350 10 45 84.3 0.880 82.0 1.034 80.7
## 21 21 350 12 7 92.2 0.612 63.4 0.906 61.3
## 22 22 350 12 15 93.0 0.643 87.1 0.698 84.3
## 23 23 350 12 20 92.2 0.669 91.9 0.682 87.2
## 24 24 350 12 25 91.5 0.698 90.1 0.720 91.1
## 25 25 350 12 30 89.9 0.737 89.3 0.753 89.3
## 26 26 300 8 7 28.8 1.640 18.0 2.679 17.7
## 27 27 300 8 15 50.6 0.991 33.7 1.517 22.2
## 28 28 300 8 20 53.2 0.976 37.7 1.404 33.6
## 29 29 300 8 25 57.6 0.933 43.3 1.264 41.3
## 30 30 300 8 30 60.7 0.915 51.5 1.098 47.4
## 31 31 300 8 35 61.3 0.936 53.2 1.097 50.3
## 32 32 300 8 40 61.7 0.959 55.1 1.092 52.0
## 33 33 300 8 45 62.3 0.979 57.4 1.080 52.8
## 34 34 300 8 50 62.8 1.000 58.8 1.085 53.1
## 35 35 300 8 55 64.1 1.008 59.1 1.111 55.2
## 36 36 300 8 60 63.2 1.051 60.9 1.107 55.8
## 37 37 300 8 65 63.0 1.094 59.3 1.168 56.4
## 38 38 300 8 70 62.8 1.137 59.0 1.228 57.2
## 39 39 300 8 75 61.0 1.180 58.5 1.289 56.9
## 40 40 300 10 7 36.4 1.324 23.4 2.103 20.1
## 41 41 300 10 15 65.6 0.779 44.4 1.174 45.6
## 42 42 300 10 20 67.7 0.782 52.2 1.033 50.3
## 43 43 300 10 25 71.2 0.769 57.4 0.971 56.4
## 44 44 300 10 30 74.0 0.764 62.0 0.928 60.8
## 45 45 300 10 35 72.2 0.808 64.3 0.923 61.1
## 46 46 300 10 40 72.0 0.852 66.3 0.923 63.2
## 47 47 300 10 45 71.9 0.897 64.1 0.983 64.2
## 48 48 300 10 50 70.8 0.941 63.9 1.043 63.1
## 49 49 300 12 7 40.1 1.223 26.4 1.896 23.1
## 50 50 300 12 15 68.4 0.760 64.4 0.822 52.2
## 51 51 300 12 20 74.3 0.724 72.1 0.760 55.3
## 52 52 300 12 25 82.9 0.671 77.3 0.732 61.8
## 53 53 300 12 30 81.1 0.708 80.6 0.725 68.8
## 54 54 300 12 35 80.9 0.745 78.3 0.769 69.7
## 55 55 300 12 40 80.1 0.782 76.9 0.813 71.1
## 56 56 300 12 45 79.5 0.820 76.6 0.858 68.3
## 57 57 250 8 7 13.8 2.816 9.8 4.068 8.3
## 58 58 250 8 15 26.7 1.545 15.1 2.798 15.4
## 59 59 250 8 20 30.7 1.392 27.7 1.579 19.7
## 60 60 250 8 25 33.2 1.332 29.2 1.549 21.8
## 61 61 250 8 30 35.9 1.273 30.8 1.517 22.2
## 62 62 250 8 35 36.4 1.297 31.6 1.525 24.4
## 63 63 250 8 40 37.6 1.295 31.9 1.558 25.2
## 64 64 250 8 45 38.2 1.314 32.3 1.585 26.4
## 65 65 250 8 50 39.6 1.305 34.2 1.540 27.6
## 66 66 250 8 55 41.2 1.290 35.6 1.521 29.1
## 67 67 250 8 60 44.2 1.236 37.3 1.492 29.8
## 68 68 250 8 65 46.2 1.215 38.7 1.477 30.9
## 69 69 250 8 70 48.3 1.193 39.9 1.469 32.1
## 70 70 250 8 75 49.9 1.185 41.2 1.459 34.4
## 71 71 250 8 80 51.9 1.168 42.1 1.463 35.5
## 72 72 250 8 85 50.9 1.220 43.3 1.457 37.2
## 73 73 250 8 90 50.1 1.272 44.7 1.445 38.9
## 74 74 250 8 95 49.6 1.324 43.9 1.505 40.3
## 75 75 250 8 100 48.6 1.377 43.2 1.565 39.2
## 76 76 250 10 7 14.9 2.675 12.0 3.405 10.8
## 77 77 250 10 15 29.5 1.432 17.2 2.514 17.3
## 78 78 250 10 20 34.0 1.286 33.2 1.347 22.3
## 79 79 250 10 25 35.1 1.288 35.6 1.298 28.7
## 80 80 250 10 30 38.8 1.204 37.1 1.286 32.2
## 81 81 250 10 35 43.2 1.116 39.8 1.236 32.7
## 82 82 250 10 40 46.5 1.069 43.2 1.173 32.9
## 83 83 250 10 45 47.6 1.075 44.3 1.178 33.1
## 84 84 250 10 50 48.7 1.082 46.3 1.159 36.2
## 85 85 250 10 55 53.4 1.014 48.2 1.144 37.5
## 86 86 250 10 60 55.0 1.012 49.7 1.140 38.7
## 87 87 250 10 65 52.9 1.080 50.3 1.156 41.1
## 88 88 250 10 70 52.6 1.148 52.1 1.144 43.2
## 89 89 250 10 75 52.0 1.217 51.7 1.133 44.1
## 90 90 250 10 80 51.6 1.220 51.2 1.240 45.3
## 91 91 250 10 85 51.1 1.272 50.4 1.248 44.2
## 92 92 250 12 7 17.7 2.300 17.8 2.343 11.3
## 93 93 250 12 15 33.8 1.275 19.4 2.273 19.5
## 94 94 250 12 20 37.5 1.189 34.6 1.317 25.3
## 95 95 250 12 25 38.7 1.190 37.5 1.255 31.1
## 96 96 250 12 30 40.3 1.180 41.7 1.164 33.2
## 97 97 250 12 35 44.1 1.112 42.9 1.167 39.7
## 98 98 250 12 40 48.8 1.036 47.3 1.090 43.1
## 99 99 250 12 45 53.6 0.971 49.7 1.067 47.2
## 100 100 250 12 50 56.1 0.954 51.1 1.067 48.3
## 101 101 250 12 55 53.6 1.026 52.3 1.071 49.6
## 102 102 250 12 60 53.2 1.099 53.2 1.081 50.8
## 103 103 250 12 65 52.9 1.171 49.8 1.184 51.5
## 104 104 250 12 70 52.5 1.243 48.9 1.288 52.1
## 105 105 250 12 75 52.0 1.315 47.9 1.392 51.7
## prop1_ec
## 1 1.477
## 2 1.107
## 3 1.081
## 4 1.069
## 5 1.083
## 6 1.111
## 7 1.139
## 8 1.165
## 9 1.149
## 10 1.124
## 11 1.130
## 12 1.177
## 13 1.332
## 14 0.860
## 15 0.824
## 16 0.841
## 17 0.840
## 18 0.854
## 19 0.861
## 20 0.915
## 21 0.953
## 22 0.733
## 23 0.730
## 24 0.723
## 25 0.764
## 26 2.724
## 27 2.303
## 28 1.576
## 29 1.326
## 30 1.193
## 31 1.161
## 32 1.158
## 33 1.174
## 34 1.202
## 35 1.189
## 36 1.209
## 37 1.228
## 38 1.243
## 39 1.281
## 40 2.448
## 41 1.143
## 42 1.072
## 43 0.988
## 44 0.946
## 45 0.971
## 46 0.968
## 47 0.981
## 48 1.027
## 49 2.167
## 50 1.015
## 51 0.990
## 52 0.916
## 53 0.849
## 54 0.864
## 55 0.872
## 56 0.935
## 57 4.803
## 58 2.743
## 59 2.220
## 60 2.074
## 61 2.104
## 62 1.976
## 63 1.972
## 64 1.939
## 65 1.909
## 66 1.861
## 67 1.868
## 68 1.849
## 69 1.827
## 70 1.748
## 71 1.735
## 72 1.696
## 73 1.660
## 74 1.640
## 75 1.724
## 76 3.783
## 77 2.499
## 78 2.006
## 79 1.610
## 80 1.482
## 81 1.505
## 82 1.541
## 83 1.576
## 84 1.483
## 85 1.471
## 86 1.464
## 87 1.415
## 88 1.380
## 89 1.386
## 90 1.382
## 91 1.450
## 92 3.691
## 93 2.261
## 94 1.801
## 95 1.513
## 96 1.463
## 97 1.261
## 98 1.196
## 99 1.123
## 100 1.129
## 101 1.129
## 102 1.132
## 103 1.145
## 104 1.161
## 105 1.198as.factor(datos2$eth_pct)## [1] 44.4 62.1 66.2 68.3 70.1 70.9 71.1 72.2 72.9 73.3 72.1 70.2 51.3 77.3 81.1
## [16] 82.8 84.6 83.1 82.9 82 63.4 87.1 91.9 90.1 89.3 18 33.7 37.7 43.3 51.5
## [31] 53.2 55.1 57.4 58.8 59.1 60.9 59.3 59 58.5 23.4 44.4 52.2 57.4 62 64.3
## [46] 66.3 64.1 63.9 26.4 64.4 72.1 77.3 80.6 78.3 76.9 76.6 9.8 15.1 27.7 29.2
## [61] 30.8 31.6 31.9 32.3 34.2 35.6 37.3 38.7 39.9 41.2 42.1 43.3 44.7 43.9 43.2
## [76] 12 17.2 33.2 35.6 37.1 39.8 43.2 44.3 46.3 48.2 49.7 50.3 52.1 51.7 51.2
## [91] 50.4 17.8 19.4 34.6 37.5 41.7 42.9 47.3 49.7 51.1 52.3 53.2 49.8 48.9 47.9
## 96 Levels: 9.8 12 15.1 17.2 17.8 18 19.4 23.4 26.4 27.7 29.2 30.8 31.6 ... 91.9as.factor(datos2$meth_pct)## [1] 73.4 76.5 78.8 80.9 82.2 83.1 84.3 84.9 84.1 83.9 82.9 82.6 76.1 87.5 90.5
## [16] 89.1 87.2 87 86.5 84.3 92.2 93 92.2 91.5 89.9 28.8 50.6 53.2 57.6 60.7
## [31] 61.3 61.7 62.3 62.8 64.1 63.2 63 62.8 61 36.4 65.6 67.7 71.2 74 72.2
## [46] 72 71.9 70.8 40.1 68.4 74.3 82.9 81.1 80.9 80.1 79.5 13.8 26.7 30.7 33.2
## [61] 35.9 36.4 37.6 38.2 39.6 41.2 44.2 46.2 48.3 49.9 51.9 50.9 50.1 49.6 48.6
## [76] 14.9 29.5 34 35.1 38.8 43.2 46.5 47.6 48.7 53.4 55 52.9 52.6 52 51.6
## [91] 51.1 17.7 33.8 37.5 38.7 40.3 44.1 48.8 53.6 56.1 53.6 53.2 52.9 52.5 52
## 95 Levels: 13.8 14.9 17.7 26.7 28.8 29.5 30.7 33.2 33.8 34 35.1 35.9 ... 93as.factor(datos2$prop1_pct)## [1] 38.3 54.2 57.2 59.9 61.3 61.7 62.1 62.5 65.3 68.8 70.1 69.1 43.2 70.9 76.2
## [16] 77.3 80.2 81.5 83.3 80.7 61.3 84.3 87.2 91.1 89.3 17.7 22.2 33.6 41.3 47.4
## [31] 50.3 52 52.8 53.1 55.2 55.8 56.4 57.2 56.9 20.1 45.6 50.3 56.4 60.8 61.1
## [46] 63.2 64.2 63.1 23.1 52.2 55.3 61.8 68.8 69.7 71.1 68.3 8.3 15.4 19.7 21.8
## [61] 22.2 24.4 25.2 26.4 27.6 29.1 29.8 30.9 32.1 34.4 35.5 37.2 38.9 40.3 39.2
## [76] 10.8 17.3 22.3 28.7 32.2 32.7 32.9 33.1 36.2 37.5 38.7 41.1 43.2 44.1 45.3
## [91] 44.2 11.3 19.5 25.3 31.1 33.2 39.7 43.1 47.2 48.3 49.6 50.8 51.5 52.1 51.7
## 98 Levels: 8.3 10.8 11.3 15.4 17.3 17.7 19.5 19.7 20.1 21.8 22.2 22.3 ... 91.1#Etanol
regre3 <- lm(datos2$eth_pct ~ datos2$b_temp + datos2$b_pres + datos2$b_time, data=datos2)
regre3##
## Call:
## lm(formula = datos2$eth_pct ~ datos2$b_temp + datos2$b_pres +
## datos2$b_time, data = datos2)
##
## Coefficients:
## (Intercept) datos2$b_temp datos2$b_pres datos2$b_time
## -143.444 0.455 4.908 0.429#Metanol
regre4 <- lm(datos2$meth_pct ~ datos2$b_temp + datos2$b_pres + datos2$b_time, data=datos2)
regre4##
## Call:
## lm(formula = datos2$meth_pct ~ datos2$b_temp + datos2$b_pres +
## datos2$b_time, data = datos2)
##
## Coefficients:
## (Intercept) datos2$b_temp datos2$b_pres datos2$b_time
## -132.661 0.498 3.579 0.342#Propanol
regre5 <- lm(datos2$prop1_pct ~ datos2$b_temp + datos2$b_pres + datos2$b_time, data=datos2)
regre5##
## Call:
## lm(formula = datos2$prop1_pct ~ datos2$b_temp + datos2$b_pres +
## datos2$b_time, data = datos2)
##
## Coefficients:
## (Intercept) datos2$b_temp datos2$b_pres datos2$b_time
## -155.924 0.468 5.194 0.450Ejercicio 3.
La publicación relacionada a este ejercicio está disponible en el archivo Liu 2012. Se realizó un experimento en el cual se valoraron 35 variedades de cervezas. El cambio químico (Chemical Shift, no incluído en la publicación) corresponde a una variable de respuesta. Se valoraron 30 variables predictoras.
Elabore un análisis de correlación y considere los parámetros estadísticos necesarios para la reducción en la cantidad de variables.
Las varibles son las siguientes:
EA (Ethyl Acetate) propnol (propanol) isoBut (Iso-butanol) isoAmlAce (Iso-amyl acetate) isoAmlAlc (Isoamyl alcohol) AD (Alcohol Degree, vol %) OA (Oxalic acid) TA (Tartaric Acid) PA (Pyruvic Acid) MA (Malic Acid) KA (Ketoglutaric Acid) LA (Lactic Acid) AA (Acetic Acid) CA (Citric Acid) SA (Succinic Acid) Kp (K^+) Nap (Na^+) Mg2p (Mg^{2+}) Ca2p (CA^{2+}) Clm (Cl^-) SO42m (SO_4^{2-}) GA (Gallic acid) catechin VA (Vanillic acid) Caffeic (Caffeic acid) Syringic (Syringic acid) PCA (P-coumarilic acid) FA (Ferulic acid)
datos3 = read.csv("C:/Users/GOLLO/Desktop/hydrogen_beer.csv", sep = ",", dec = ".")
datos3## beerNum EA propnol isoBut isoAmlAce isoAmlAlc AD OA TA PA
## 1 1 14.33 8.22 11.06 2.32 69.91 4.22 0.66 121.60 12.69
## 2 2 8.87 10.26 18.02 1.13 80.42 3.43 0.35 37.11 16.81
## 3 3 10.30 8.40 10.23 1.09 57.44 3.57 0.87 76.33 14.19
## 4 4 9.35 5.71 5.67 1.67 37.83 3.43 0.71 113.78 12.44
## 5 5 9.86 9.93 7.01 1.76 52.22 4.00 0.67 133.04 9.51
## 6 6 9.90 8.28 10.16 1.41 61.99 3.85 0.31 126.45 16.43
## 7 7 11.18 8.00 10.59 2.31 51.80 3.43 0.45 94.07 6.78
## 8 8 12.14 6.62 8.91 2.20 51.80 3.78 0.43 124.52 5.47
## 9 9 15.59 9.67 11.89 3.41 63.23 3.85 1.09 127.55 8.47
## 10 10 16.37 9.88 15.51 3.82 71.96 4.29 0.66 101.32 8.16
## 11 11 11.87 6.84 7.12 1.80 47.77 3.86 0.96 126.99 9.44
## 12 12 18.33 10.02 13.08 3.95 82.97 4.94 0.65 137.19 14.85
## 13 13 11.19 10.91 7.87 2.30 59.57 3.92 0.69 140.07 9.82
## 14 14 11.56 11.08 7.91 2.16 61.20 3.92 0.73 130.05 9.16
## 15 15 12.73 11.93 8.30 2.54 66.50 3.92 0.78 147.14 8.95
## 16 16 11.06 11.43 11.61 1.69 58.43 3.71 0.25 110.72 11.79
## 17 17 17.51 7.94 15.16 3.89 84.74 3.00 0.29 106.39 6.31
## 18 18 15.19 8.81 8.38 2.09 55.46 4.07 0.76 110.61 11.77
## 19 19 7.06 13.34 20.66 1.26 111.44 4.36 0.17 65.25 19.64
## 20 20 9.10 8.21 13.64 1.08 69.40 3.85 0.13 82.99 13.87
## 21 21 14.60 9.77 15.33 2.68 80.76 4.58 0.36 116.66 15.13
## 22 22 5.96 13.39 19.26 0.90 108.43 4.14 0.84 61.14 21.77
## 23 23 13.28 6.86 9.16 1.70 64.21 3.85 0.56 123.28 9.79
## 24 24 18.25 7.48 15.17 3.97 81.49 4.36 0.52 118.97 16.12
## 25 25 13.89 10.55 10.37 2.32 74.54 4.00 0.96 156.97 8.90
## 26 26 13.56 13.20 14.32 2.06 71.45 4.36 0.25 125.87 11.07
## 27 27 22.86 15.53 21.12 5.99 124.48 4.94 1.09 189.40 13.95
## 28 28 11.82 7.80 10.17 2.08 66.14 3.15 0.32 81.24 9.80
## 29 29 12.93 20.54 17.25 1.78 100.80 5.16 0.58 195.84 19.11
## 30 30 10.77 8.87 11.03 1.61 73.08 2.93 0.65 54.51 8.85
## 31 31 5.43 7.39 9.45 0.64 57.86 3.29 0.13 88.47 11.16
## 32 32 10.00 9.44 13.90 1.73 66.44 3.85 0.29 108.95 14.16
## 33 33 6.78 6.62 6.06 0.71 41.35 2.51 0.33 103.14 10.57
## 34 34 6.89 7.15 5.05 0.55 37.88 2.72 0.33 93.49 9.06
## 35 35 6.99 6.23 5.54 0.68 35.84 2.65 0.92 134.27 8.86
## MA KA LA AA CA SA Kp Nap Mg2p Ca2p Clm
## 1 14.06 2.43 108.28 64.90 131.85 235.39 69.65 78.25 63.25 38.20 89.49
## 2 4.72 1.87 60.82 26.59 99.37 160.67 69.90 51.60 48.00 69.15 113.18
## 3 6.33 1.94 87.46 25.53 83.62 155.60 96.25 34.45 28.90 32.45 110.32
## 4 10.30 1.37 115.01 14.83 146.99 291.72 127.25 50.75 53.00 30.15 108.96
## 5 17.12 3.81 72.56 55.68 105.49 189.29 94.60 45.75 57.75 41.30 61.25
## 6 14.94 2.58 150.56 93.02 145.65 274.02 185.00 72.50 59.75 22.35 100.20
## 7 8.07 1.95 45.70 24.82 153.95 92.45 129.75 74.50 49.25 47.10 120.65
## 8 14.01 1.38 67.08 15.84 119.64 230.27 165.75 34.25 47.25 25.38 56.38
## 9 11.89 2.86 51.95 29.01 121.14 94.98 132.75 80.25 53.50 36.80 89.16
## 10 11.94 2.31 50.23 26.63 101.47 74.54 167.60 72.95 50.50 30.25 86.89
## 11 11.70 1.62 97.07 17.72 166.51 187.46 167.00 39.25 52.00 27.30 58.98
## 12 14.04 3.76 82.95 17.88 78.05 145.30 89.70 72.00 59.00 53.20 118.38
## 13 19.13 3.92 58.37 21.03 74.73 265.48 47.45 141.70 38.55 40.95 50.21
## 14 18.19 3.19 69.81 13.40 79.40 187.72 54.40 58.35 42.15 48.75 56.05
## 15 18.67 3.78 66.31 19.27 70.84 293.66 69.60 82.50 48.95 24.50 53.13
## 16 10.90 3.02 60.87 19.32 109.67 254.79 99.60 82.50 58.95 24.50 85.27
## 17 11.13 2.71 45.75 15.43 70.74 71.01 70.70 47.75 52.00 36.70 83.71
## 18 12.40 1.90 73.83 19.50 95.30 186.25 67.25 50.75 53.03 30.15 108.96
## 19 6.78 4.43 126.38 15.98 66.75 233.11 104.00 36.25 54.05 28.05 104.82
## 20 8.90 1.90 51.85 43.60 54.90 208.21 41.60 38.50 45.25 57.30 99.87
## 21 8.22 1.85 79.33 11.48 74.89 205.41 116.40 74.00 61.75 30.45 112.86
## 22 7.27 4.55 72.15 16.52 68.98 238.88 148.80 50.75 54.25 59.40 93.26
## 23 16.28 2.77 100.35 17.05 101.94 267.64 41.20 48.90 44.80 49.65 86.89
## 24 12.44 2.93 46.78 23.91 42.36 81.46 65.60 40.50 36.90 46.80 93.01
## 25 18.11 3.92 68.68 77.69 74.06 172.96 148.75 69.20 58.25 22.75 59.30
## 26 14.32 2.10 78.76 52.72 143.53 217.35 76.90 61.00 51.75 38.90 109.61
## 27 16.05 4.18 85.87 15.34 259.02 136.37 62.00 66.75 72.95 58.90 115.59
## 28 11.41 2.38 90.01 17.49 36.24 231.74 59.30 47.25 38.75 34.04 83.78
## 29 7.95 5.98 207.04 47.47 53.71 190.55 102.50 25.90 46.65 29.45 165.77
## 30 2.81 1.41 39.79 16.02 50.06 138.18 65.40 69.80 42.40 38.40 85.27
## 31 2.84 1.61 45.15 16.42 29.37 101.51 69.70 33.60 58.70 48.80 104.49
## 32 4.80 2.45 103.75 17.45 38.21 261.46 66.50 30.75 33.70 36.65 84.31
## 33 6.73 3.70 47.34 17.17 71.77 147.56 62.00 66.75 52.95 58.90 75.87
## 34 8.99 2.63 51.58 16.72 85.55 247.17 59.30 47.25 38.75 34.02 81.71
## 35 6.63 3.94 65.05 17.46 162.57 234.78 66.50 30.75 33.70 36.65 81.38
## SO42m GA catechin VA Caffeic Syringic PCA FA OG RG
## 1 266.55 5.59 0.39 0.46 0.22 0.38 0.65 2.23 10.00 3.43
## 2 162.16 3.96 0.36 0.34 0.16 0.42 0.54 1.94 8.18 2.80
## 3 189.50 4.12 0.38 0.39 0.20 0.45 0.62 2.05 8.12 2.50
## 4 56.52 4.21 0.32 0.35 0.21 0.43 0.59 2.13 8.06 2.67
## 5 172.10 6.25 0.76 0.51 0.08 0.19 0.65 1.79 9.49 3.25
## 6 266.55 6.10 0.59 0.46 0.23 0.22 0.46 1.58 9.17 3.14
## 7 153.46 3.48 0.40 0.39 0.20 0.28 0.42 1.45 8.05 2.66
## 8 153.46 4.21 0.46 0.31 0.12 0.21 0.46 1.23 8.54 2.61
## 9 83.86 3.96 0.45 0.38 0.15 0.32 0.41 1.25 9.19 3.17
## 10 90.08 3.87 0.43 0.42 0.17 0.29 0.38 1.34 9.48 2.78
## 11 163.40 4.14 0.35 0.37 0.25 0.41 0.52 1.95 8.48 2.86
## 12 158.43 4.37 0.13 0.39 0.23 0.39 0.56 1.54 11.40 3.52
## 13 168.37 6.09 0.48 0.49 0.41 0.49 0.77 1.83 9.57 3.32
## 14 163.40 5.57 0.46 0.51 0.23 0.14 0.58 1.59 9.53 3.40
## 15 159.67 5.89 0.50 0.56 0.37 0.28 0.63 1.64 9.43 3.30
## 16 100.02 3.86 0.47 0.42 0.21 0.28 0.40 1.49 8.44 2.61
## 17 106.23 4.02 0.58 0.49 0.24 0.31 0.58 1.53 7.77 3.03
## 18 56.52 4.53 0.40 0.38 0.22 0.45 0.60 2.10 9.88 2.84
## 19 114.93 3.03 0.22 0.14 0.09 0.24 0.32 1.21 9.72 3.13
## 20 180.80 4.03 0.37 0.35 0.18 0.44 0.57 1.99 8.12 2.07
## 21 169.62 4.75 0.47 0.38 0.13 0.14 0.51 1.37 10.49 3.37
## 22 113.69 3.14 0.25 0.19 0.10 0.28 0.37 1.24 9.56 3.08
## 23 229.27 5.87 0.66 0.50 0.18 0.14 0.63 1.55 9.05 3.02
## 24 81.38 3.96 0.57 0.46 0.21 0.22 0.54 1.49 9.59 2.78
## 25 169.62 5.17 0.30 0.48 0.18 0.23 0.52 1.59 9.29 3.03
## 26 95.05 4.02 0.35 0.32 0.17 0.39 0.58 1.95 9.71 2.90
## 27 143.52 3.45 0.29 0.28 0.13 0.36 0.47 1.35 11.12 3.45
## 28 209.38 4.30 0.36 0.42 0.20 0.21 0.52 1.54 7.51 2.55
## 29 152.22 5.18 0.29 0.33 0.24 0.06 0.49 1.48 11.24 3.11
## 30 201.93 3.82 0.35 0.32 0.15 0.40 0.50 1.89 7.16 2.52
## 31 198.20 3.74 0.38 0.34 0.17 0.43 0.54 1.94 8.26 3.09
## 32 154.70 4.09 0.20 0.36 0.18 0.29 0.51 1.52 9.08 3.06
## 33 150.97 3.06 0.25 0.26 0.10 0.32 0.43 1.30 5.88 1.90
## 34 131.09 3.34 0.28 0.28 0.12 0.35 0.45 1.33 6.27 1.95
## 35 104.99 2.96 0.26 0.24 0.11 0.28 0.40 1.31 6.27 2.05componentes_p <- prcomp(datos3, scale. = T)
names(componentes_p)## [1] "sdev" "rotation" "center" "scale" "x"componentes_p$center## beerNum EA propnol isoBut isoAmlAce isoAmlAlc AD OA
## 18.0000 11.9286 9.6086 11.5989 2.0937 68.0237 3.8254 0.5640
## TA PA MA KA LA AA CA SA
## 113.2963 11.8529 11.1449 2.8323 77.8420 27.4540 96.2377 191.5697
## Kp Nap Mg2p Ca2p Clm SO42m GA catechin
## 93.1614 57.3700 49.7523 39.0940 91.1131 150.6183 4.3466 0.3931
## VA Caffeic Syringic PCA FA OG RG
## 0.3791 0.1869 0.3063 0.5191 1.6203 8.8886 2.8843componentes_p$scale## beerNum EA propnol isoBut isoAmlAce isoAmlAlc AD OA
## 10.24695 3.90419 2.98077 4.34077 1.16264 20.50790 0.62836 0.27607
## TA PA MA KA LA AA CA SA
## 33.61683 3.89011 4.58526 1.08932 34.03843 19.51072 46.96325 64.79031
## Kp Nap Mg2p Ca2p Clm SO42m GA catechin
## 40.41614 22.48578 9.56239 12.13839 24.13179 51.26917 0.94648 0.13363
## VA Caffeic Syringic PCA FA OG RG
## 0.09500 0.06957 0.10647 0.09626 0.29569 1.31780 0.43122componentes_p$rotation## PC1 PC2 PC3 PC4 PC5 PC6 PC7
## beerNum 0.111800 -0.15880 0.023382 -0.080744 0.41767 -0.1059004 0.01032
## EA -0.257622 0.01323 -0.314405 0.064157 0.03705 0.0469272 0.23849
## propnol -0.232550 -0.21772 0.111858 0.012029 0.13779 -0.1202063 -0.04533
## isoBut -0.150489 -0.31415 -0.041321 0.178156 -0.01484 0.1867439 -0.02028
## isoAmlAce -0.247982 -0.02023 -0.374220 0.044106 0.05436 0.0625909 0.12373
## isoAmlAlc -0.211866 -0.27827 -0.012995 0.152491 0.08449 0.1146941 -0.09946
## AD -0.320451 -0.12927 0.040141 0.062049 -0.05235 -0.0159458 0.11451
## OA -0.134196 0.05596 -0.204985 -0.130462 -0.04632 -0.3546035 -0.01756
## TA -0.255021 0.05986 -0.037673 -0.237711 0.13464 -0.2425633 0.10453
## PA -0.088846 -0.27130 0.274590 0.169655 -0.07354 -0.0393653 -0.05586
## MA -0.247755 0.22582 -0.071153 -0.139874 0.09379 -0.0569992 -0.13697
## KA -0.167292 -0.15621 0.106920 -0.176124 0.32599 -0.1710756 -0.22353
## LA -0.166857 -0.10409 0.377722 -0.120424 -0.12264 -0.1381698 0.21868
## AA -0.135091 0.07933 0.228854 -0.093055 -0.24388 0.1125211 -0.14933
## CA -0.091722 0.01005 -0.213413 -0.166609 -0.33157 -0.2711761 -0.16111
## SA -0.003212 0.08768 0.360680 -0.141353 0.01885 -0.2805909 -0.14559
## Kp -0.065577 -0.04263 -0.023416 -0.330580 -0.43310 0.1184700 0.05496
## Nap -0.131150 0.16358 -0.136983 0.061143 0.02889 -0.1027922 -0.40488
## Mg2p -0.183688 -0.07437 -0.105974 -0.028042 -0.28341 0.0713745 -0.35393
## Ca2p 0.048108 -0.11165 -0.113380 0.341327 0.12411 0.0375728 -0.39458
## Clm -0.057759 -0.26100 0.095623 0.187181 -0.17029 -0.0237018 0.27360
## SO42m -0.037220 0.12557 0.276073 0.084017 -0.03897 0.2073609 -0.30422
## GA -0.219801 0.27108 0.232926 -0.001178 0.05276 0.0931644 -0.03595
## catechin -0.089167 0.24421 -0.008720 -0.047306 0.01862 0.4126910 0.02529
## VA -0.159763 0.32351 -0.003959 0.050329 0.11025 0.2137446 0.10798
## Caffeic -0.142019 0.22741 0.060185 0.164877 0.12287 -0.1998594 0.17073
## Syringic 0.113228 0.07090 -0.143392 0.365843 -0.22352 -0.3586570 -0.08670
## PCA -0.101625 0.28668 0.101278 0.318105 0.10536 -0.1095152 0.07604
## FA 0.006892 0.19628 0.139789 0.382527 -0.24350 -0.1762767 0.13097
## OG -0.339360 -0.08668 0.032890 0.094336 -0.03194 0.0001798 0.07261
## RG -0.312917 0.01105 0.019479 0.100072 -0.01313 0.0792094 -0.09000
## PC8 PC9 PC10 PC11 PC12 PC13
## beerNum -0.2352239 0.344374 -0.055060 -0.1486866 0.279317 -0.0503773
## EA -0.1479420 0.127289 -0.058642 -0.1180526 -0.076153 -0.2094586
## propnol 0.0914652 0.074273 0.145160 0.2150400 0.006274 0.0005329
## isoBut 0.1161031 0.010032 -0.018006 0.0181038 -0.079940 -0.3384094
## isoAmlAce -0.0483660 0.089068 -0.076826 -0.1263826 -0.108691 -0.1461975
## isoAmlAlc 0.0710849 -0.016804 0.009513 -0.1220678 0.023583 -0.2951564
## AD 0.0486897 -0.069796 -0.094703 0.0018138 0.155672 -0.0184852
## OA -0.0448635 -0.458007 0.402664 -0.2375834 0.067170 -0.0081947
## TA -0.2989276 0.152353 0.004188 0.0305481 0.028666 0.1752155
## PA 0.1017329 -0.178173 0.007344 0.0954587 -0.038898 -0.0662303
## MA -0.0190774 -0.060392 -0.127836 0.1442311 0.085565 -0.2472622
## KA -0.0002206 -0.082716 0.290884 0.2152190 -0.036557 0.0956520
## LA -0.1242036 -0.007593 -0.057735 -0.1148486 -0.227947 0.0485379
## AA -0.2374141 0.288034 0.426789 0.2324572 0.162799 -0.3693687
## CA -0.3138704 -0.082440 -0.227166 0.0851035 -0.346629 -0.1462388
## SA 0.1475419 -0.066516 -0.511972 -0.0077738 -0.059223 -0.2626285
## Kp 0.2879028 -0.029102 0.205675 -0.1090018 0.048855 0.0606804
## Nap 0.3604067 0.343218 0.059431 0.0754851 -0.276953 0.1091678
## Mg2p -0.1168172 0.266904 -0.243208 0.0556952 0.309223 0.3596538
## Ca2p -0.2807702 -0.306851 0.046846 0.1040601 -0.181199 0.1464182
## Clm -0.2222977 0.177467 0.044668 0.1846028 -0.342660 0.3104620
## SO42m -0.2559544 0.106611 0.111654 -0.6248708 -0.237220 -0.0521461
## GA -0.0340112 -0.074242 0.019285 0.0008333 0.014500 0.0789549
## catechin -0.0753222 -0.200009 -0.110617 0.3870947 -0.104238 -0.0207736
## VA 0.0028377 0.023062 0.118585 0.0417245 -0.064362 0.0398774
## Caffeic 0.3113438 0.221518 0.111704 -0.0335997 -0.313694 0.0545836
## Syringic 0.1033100 0.140965 0.069731 0.0211049 0.132554 -0.1345170
## PCA -0.1301657 -0.113150 -0.051583 0.0367373 0.123754 0.0027378
## FA -0.1377315 0.008127 0.059395 0.0710053 0.256512 -0.0447513
## OG 0.0780948 -0.079058 -0.102789 -0.0411572 0.182557 0.1065932
## RG 0.1427189 -0.126619 -0.111356 -0.2471265 0.169288 0.2957157
## PC14 PC15 PC16 PC17 PC18 PC19 PC20
## beerNum -0.08663 0.287619 0.106497 0.449557 -0.06428 -0.154608 -0.136370
## EA 0.10353 0.133549 0.138472 -0.178245 -0.09023 0.104165 -0.048428
## propnol -0.40304 -0.047497 -0.477217 -0.005623 -0.22668 0.098574 -0.004427
## isoBut -0.15465 -0.075673 -0.032411 0.075711 -0.01859 -0.124954 -0.063231
## isoAmlAce 0.08008 -0.053763 0.159185 -0.106932 0.08484 -0.042318 -0.080330
## isoAmlAlc -0.25529 -0.091503 0.131882 0.001959 -0.06668 0.102559 -0.202989
## AD 0.29067 0.130521 -0.292894 0.145313 0.03971 0.262198 0.264986
## OA -0.19728 0.364732 0.157048 -0.011968 -0.19717 0.024521 0.002547
## TA 0.09841 -0.120573 -0.037621 0.248318 0.16123 0.064137 0.010798
## PA 0.23400 0.198907 0.477789 0.204298 0.04022 0.033618 0.212455
## MA 0.24131 -0.166269 -0.052224 -0.084725 -0.13114 0.109851 -0.116364
## KA -0.02580 -0.233558 0.242227 -0.216602 0.21926 0.145009 0.109794
## LA -0.05590 -0.146818 0.153172 -0.149132 0.09534 -0.020347 -0.287636
## AA 0.17694 0.039243 0.003664 -0.039160 0.07579 -0.331095 0.016237
## CA -0.20821 -0.058969 -0.117355 0.189326 0.06723 -0.207369 0.186585
## SA 0.03877 0.152573 0.065153 -0.043794 -0.21026 -0.081952 -0.045474
## Kp 0.13666 -0.079312 -0.026147 0.373175 -0.10188 0.088313 -0.420553
## Nap 0.04111 0.469891 -0.061755 -0.094378 0.09887 -0.019564 -0.148499
## Mg2p -0.14274 -0.123998 0.327258 -0.081405 -0.27746 0.226483 0.009447
## Ca2p 0.34188 -0.142236 -0.102662 0.230440 -0.24447 -0.126694 -0.240049
## Clm 0.03369 0.274587 -0.001460 -0.122516 -0.07231 -0.020498 -0.105421
## SO42m -0.05305 -0.005649 -0.119715 0.049253 0.07857 0.308537 0.188775
## GA -0.01243 0.132064 -0.050637 -0.041421 -0.01485 0.005027 -0.181055
## catechin -0.35602 0.197633 0.151933 0.345461 0.25855 0.209213 0.012692
## VA 0.02021 -0.024817 0.150371 -0.025662 -0.36747 -0.209848 0.267562
## Caffeic -0.01699 -0.332801 0.151911 0.382529 -0.19397 -0.003840 0.114414
## Syringic -0.01721 -0.111872 0.072840 0.104049 0.37326 0.136279 0.028991
## PCA -0.01169 -0.009168 -0.018457 -0.032700 0.19469 0.103781 -0.457834
## FA -0.21440 0.001960 -0.003796 -0.012833 -0.24298 0.009524 0.061601
## OG 0.18585 0.131279 -0.203513 0.056470 0.10333 -0.004055 0.175426
## RG -0.16638 -0.090264 0.006278 -0.003035 0.26923 -0.609692 0.020179
## PC21 PC22 PC23 PC24 PC25 PC26 PC27
## beerNum 0.219150 0.0694027 0.032203 0.05884 0.260760 -0.011380 -0.017993
## EA -0.040051 0.0736038 0.099743 -0.02738 -0.019531 -0.544532 -0.411188
## propnol -0.122104 0.0616793 0.294812 0.39443 -0.105311 -0.038485 -0.144142
## isoBut -0.242941 -0.0804054 -0.104523 -0.26998 0.130362 -0.197642 0.370496
## isoAmlAce -0.176963 0.0621301 0.171085 0.17376 0.034283 0.660676 -0.171816
## isoAmlAlc 0.139787 -0.1207621 -0.173248 -0.08106 -0.129420 0.164440 0.092804
## AD 0.015805 0.0771795 -0.105420 -0.18400 0.226144 0.099226 0.094371
## OA 0.045513 -0.1134660 0.061328 -0.15414 -0.115437 0.025194 0.135648
## TA -0.433091 -0.0656582 -0.048902 -0.06836 -0.288662 0.061235 0.271445
## PA -0.114332 0.2005522 -0.020865 0.46244 -0.217198 -0.083530 0.029471
## MA 0.571422 0.0101499 -0.024686 0.24683 -0.089808 0.018491 0.293568
## KA 0.024546 -0.1049587 -0.174201 -0.08396 0.413868 0.004545 -0.275340
## LA 0.098487 0.4608049 0.258729 -0.14586 0.188706 -0.028673 0.203417
## AA -0.003854 -0.0816820 0.021194 -0.12634 -0.173621 0.040403 -0.128518
## CA 0.064267 0.0698320 -0.212638 0.14759 0.163633 -0.127147 -0.048529
## SA -0.201206 -0.3536302 0.111050 -0.12931 0.058952 0.132107 -0.210168
## Kp -0.063438 -0.0984869 -0.058705 0.16768 0.263469 0.003966 -0.173762
## Nap -0.061793 0.2754796 -0.156785 -0.04903 0.059401 0.031928 0.115052
## Mg2p -0.043139 -0.0422304 0.034359 -0.08522 -0.060742 -0.080458 0.045683
## Ca2p -0.049419 0.0862959 0.191041 -0.16172 0.038560 0.027798 -0.115965
## Clm 0.273119 -0.4618915 -0.049911 0.05204 0.009648 0.123069 0.064417
## SO42m 0.028030 -0.1003420 -0.020231 0.12649 0.024177 0.031996 -0.040811
## GA -0.037479 0.0215961 0.297549 -0.12350 -0.059523 -0.078288 0.002368
## catechin 0.096730 -0.0007225 0.136848 -0.12634 -0.001309 0.084089 -0.092040
## VA -0.180243 -0.1402766 0.107674 0.22103 0.461701 -0.044940 0.305426
## Caffeic 0.177637 -0.0447439 -0.085278 -0.21189 -0.228592 -0.036006 -0.155173
## Syringic 0.056828 -0.2183925 0.502624 0.03821 0.193154 -0.068384 0.135874
## PCA -0.198893 -0.1721231 -0.381989 0.25327 0.052115 -0.160347 0.003790
## FA -0.008825 0.3363171 -0.252885 -0.08868 0.138652 0.260165 -0.160870
## OG 0.131518 -0.0599857 0.032417 -0.14410 0.061041 -0.033863 -0.155954
## RG 0.136658 -0.0567369 -0.003767 0.09413 -0.067440 -0.038067 -0.082785
## PC28 PC29 PC30 PC31
## beerNum -0.061182 -0.0901302 -0.032874 -0.0279160
## EA 0.150009 0.1823454 0.090691 -0.1280397
## propnol 0.113796 -0.0883968 0.002011 0.0152795
## isoBut 0.003568 -0.1691031 -0.379085 -0.3213230
## isoAmlAce -0.095266 -0.2442690 -0.119450 -0.0472821
## isoAmlAlc -0.170560 0.4107079 0.361176 0.3429674
## AD 0.050737 -0.2480012 0.512514 -0.1464746
## OA 0.102997 -0.1821321 -0.039123 0.0173172
## TA 0.109266 0.3789118 -0.045191 -0.0757736
## PA -0.081886 0.0165470 -0.043135 -0.0584752
## MA 0.120598 0.0667268 -0.184673 -0.2193959
## KA -0.053704 0.0586501 -0.064071 -0.1512398
## LA 0.189552 -0.0867880 0.021420 0.2031402
## AA 0.087827 -0.1652226 0.088625 0.1313474
## CA -0.293394 -0.0470778 0.062824 0.1035815
## SA 0.180436 0.0538584 0.015214 -0.0422281
## Kp 0.012443 0.1451893 -0.042151 -0.0869520
## Nap 0.085927 0.1067051 0.004523 0.0444914
## Mg2p -0.005514 -0.2462196 -0.009643 0.0547609
## Ca2p 0.119955 0.0567301 0.044600 0.0320012
## Clm 0.013014 0.0543574 -0.018183 -0.1577005
## SO42m 0.117589 0.0003933 -0.123502 -0.0864615
## GA -0.745693 0.0365304 0.045418 -0.2607927
## catechin 0.246151 0.0115768 -0.020246 0.0005374
## VA 0.036996 0.1284404 0.101191 0.2270463
## Caffeic -0.038734 -0.2348697 0.004107 0.0301108
## Syringic 0.008435 0.1158046 0.058168 -0.0108112
## PCA 0.023493 -0.3331339 0.022969 0.1960167
## FA 0.048077 0.3075085 -0.157476 -0.2531713
## OG -0.104344 0.1109310 -0.545673 0.5109236
## RG 0.204441 0.0620920 0.157492 -0.2376392componentes_p$rotation[1:5,1:4]## PC1 PC2 PC3 PC4
## beerNum 0.1118 -0.15880 0.02338 -0.08074
## EA -0.2576 0.01323 -0.31440 0.06416
## propnol -0.2326 -0.21772 0.11186 0.01203
## isoBut -0.1505 -0.31415 -0.04132 0.17816
## isoAmlAce -0.2480 -0.02023 -0.37422 0.04411sd <- componentes_p$sdev
pv <- sd^2
pv[1:10]## [1] 7.1805 5.7986 3.1313 2.6411 2.2322 1.7801 1.3369 1.2260 0.9071 0.8425propv <- pv/sum(pv)
propv[1:20]## [1] 0.231630 0.187051 0.101009 0.085197 0.072005 0.057423 0.043124 0.039549
## [9] 0.029261 0.027176 0.022433 0.020261 0.017351 0.013543 0.010056 0.009658
## [17] 0.007273 0.005603 0.005085 0.003818#EA
EA1 <- datos3$EA
#propnol
propnol1 <- datos3$propnol
#isoBut
isoBut1 <- datos3$isoBut
#isoAmlAce
isoAm1 <- datos3$isoAmlAce
regre6 <- lm(formula= datos3$Syringic ~ EA1 + propnol1 + isoBut1 + isoAm1 , data=datos3)
regre6##
## Call:
## lm(formula = datos3$Syringic ~ EA1 + propnol1 + isoBut1 + isoAm1,
## data = datos3)
##
## Coefficients:
## (Intercept) EA1 propnol1 isoBut1 isoAm1
## 0.42000 -0.00352 -0.01157 0.00175 0.00912