Optimasi Ekstraksi Larva Hermetia illucens dengan Response Surface Modelling dan Kualitas Asam Amino serta Aktivitas Antibakteri

Rancangan BBD; 2 center point

# rancangan.bbd
a <- bbd(3, n0 = 2, coding = list(x1 ~ (NaOH - 2)/1, x2 ~ (HCl - 2)/1, x3 ~ (Presipitasi - 20)/10))
a <- decode.data(a)

# data.desain.bbd
kable(a) %>%
  kable_styling(bootstrap_options = c("striped", "hover", "condensed"))

run.order	std.order	NaOH	HCl	Presipitasi
1	1	1	1	20
2	11	2	1	30
3	8	3	2	30
4	14	2	2	20
5	10	2	3	10
6	5	1	2	10
7	12	2	3	30
8	9	2	1	10
9	2	3	1	20
10	4	3	3	20
11	13	2	2	20
12	3	1	3	20
13	7	1	2	30
14	6	3	2	10

Data

# input.data.awal
data.awal <- read_excel("Data_S2_2019_Produksi_Rekayasa Pengolahan BPM [HCl] dan [NaOH].xlsx", sheet = "InputYield")

# transformasi.data.ke.bentuk.(-1, 0, 1)
data.awal <- coded.data(data.awal, X1 ~ (NaOH - 2)/1, X2 ~ (HCl - 2)/1, X3 ~ (Precipitation - 20)/10)

# tabel.data.awal
kable(data.awal) %>%
  kable_styling(bootstrap_options = c("striped", "hover", "condensed"))

X1	X2	X3	Yield
1	0	-1	3.115500
0	1	-1	18.360972
0	-1	1	18.781400
-1	0	-1	3.447328
-1	1	0	5.231945
1	0	1	18.128745
0	-1	-1	7.761600
1	1	0	17.272564
-1	-1	0	5.268217
0	1	1	19.101212
1	-1	0	4.124200
-1	0	1	4.215000
0	0	0	20.156798
0	0	0	20.746737

Metode RSM

##    NaOH HCl Precipitation     Yield
## 1     3   2            10  3.115500
## 2     2   3            10 18.360972
## 3     2   1            30 18.781400
## 4     1   2            10  3.447328
## 5     1   3            20  5.231945
## 6     3   2            30 18.128745
## 7     2   1            10  7.761600
## 8     3   3            20 17.272564
## 9     1   1            20  5.268217
## 10    2   3            30 19.101212
## 11    3   1            20  4.124200
## 12    1   2            30  4.215000
## 13    2   2            20 20.156798
## 14    2   2            20 20.746737
## 
## Data are stored in coded form using these coding formulas ...
## X1 ~ (NaOH - 2)/1
## X2 ~ (HCl - 2)/1
## X3 ~ (Precipitation - 20)/10

## New names:
## * Temeprature -> Temeprature...4
## * Temeprature -> Temeprature...6

##    NaOH HCl Precipitation Temeprature...4       pH Temeprature...6
## 1     3   2            10        28.94300 11.10340        28.66773
## 2     2   3            10        28.15200 11.02920        27.63600
## 3     2   1            30        28.76400 11.00900        27.24400
## 4     1   2            10        28.10000 10.78000        28.10000
## 5     1   3            20        28.20000 10.99000        27.80000
## 6     3   2            30        28.94300 11.36090        28.87177
## 7     2   1            10        28.66200 11.12010        27.93000
## 8     3   3            20        29.04600 11.31970        28.36167
## 9     1   1            20        28.90000 11.89000        27.80000
## 10    2   3            30        28.15200 11.03930        27.63600
## 11    3   1            20        29.76700 12.24670        28.36167
## 12    1   2            30        28.10000 11.03000        28.30000
## 13    2   2            20        28.35600 11.34230        27.24400
## 14    2   2            20        27.42584 10.94435        28.39014
##          pI
## 1  3.242445
## 2  6.825750
## 3  0.636500
## 4  3.670000
## 5  0.720000
## 6  0.538935
## 7  3.914000
## 8  0.741600
## 9  2.060000
## 10 0.522500
## 11 1.820010
## 12 0.610000
## 13 1.273000
## 14 3.210570
## 
## Data are stored in coded form using these coding formulas ...
## X1 ~ (NaOH - 2)/1
## X2 ~ (HCl - 2)/1
## X3 ~ (Precipitation - 20)/10

## 
## Call:
## rsm(formula = Yield ~ SO(X1, X2, X3), data = Yield.data)
## 
##              Estimate Std. Error  t value  Pr(>|t|)    
## (Intercept)  20.45177    0.67984  30.0834 7.272e-06 ***
## X1            3.05981    0.33992   9.0016 0.0008432 ***
## X2            3.00391    0.33992   8.8372 0.0009051 ***
## X3            3.44262    0.33992  10.1278 0.0005350 ***
## X1:X2         3.29616    0.48072   6.8568 0.0023685 ** 
## X1:X3         3.56139    0.48072   7.4085 0.0017711 ** 
## X2:X3        -2.56989    0.48072  -5.3460 0.0059016 ** 
## X1^2        -10.62609    0.53746 -19.7710 3.861e-05 ***
## X2^2         -1.85144    0.53746  -3.4448 0.0261782 *  
## X3^2         -2.59903    0.53746  -4.8358 0.0084258 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Multiple R-squared:  0.9949, Adjusted R-squared:  0.9835 
## F-statistic: 87.15 on 9 and 4 DF,  p-value: 0.0003148
## 
## Analysis of Variance Table
## 
## Response: Yield
##                 Df Sum Sq Mean Sq  F value    Pr(>F)
## FO(X1, X2, X3)   3 241.90  80.634  87.2323 0.0004228
## TWI(X1, X2, X3)  3 120.61  40.203  43.4935 0.0016423
## PQ(X1, X2, X3)   3 362.50 120.833 130.7217 0.0001905
## Residuals        4   3.70   0.924                   
## Lack of fit      3   3.52   1.174   6.7493 0.2740366
## Pure error       1   0.17   0.174                   
## 
## Stationary point of response surface:
##        X1        X2        X3 
## 0.3497367 0.7560169 0.5281373 
## 
## Stationary point in original units:
##          NaOH           HCl Precipitation 
##      2.349737      2.756017     25.281373 
## 
## Eigenanalysis:
## eigen() decomposition
## $values
## [1]  -0.8788308  -2.8202727 -11.3774626
## 
## $vectors
##           [,1]       [,2]       [,3]
## X1 -0.03285412 -0.2940768  0.9552169
## X2 -0.81673783 -0.5429696 -0.1952520
## X3  0.57607285 -0.7865766 -0.2223449

## 
## Call:
## rsm(formula = pI ~ SO(X1, X2, X3), data = pH.data)
## 
##              Estimate Std. Error t value Pr(>|t|)  
## (Intercept)  2.241785   0.942693  2.3781  0.07615 .
## X1          -0.089626   0.471346 -0.1901  0.85845  
## X2           0.047417   0.471346  0.1006  0.92471  
## X3          -1.918032   0.471346 -4.0693  0.01523 *
## X1:X2        0.065398   0.666585  0.0981  0.92657  
## X1:X3        0.089122   0.666585  0.1337  0.90010  
## X2:X3       -0.756437   0.666585 -1.1348  0.31985  
## X1^2        -0.932862   0.745264 -1.2517  0.27887  
## X2^2         0.026480   0.745264  0.0355  0.97336  
## X3^2         0.706423   0.745264  0.9479  0.39686  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Multiple R-squared:  0.8399, Adjusted R-squared:  0.4798 
## F-statistic: 2.332 on 9 and 4 DF,  p-value: 0.2153
## 
## Analysis of Variance Table
## 
## Response: pI
##                 Df  Sum Sq Mean Sq F value  Pr(>F)
## FO(X1, X2, X3)   3 29.5130  9.8377  5.5351 0.06591
## TWI(X1, X2, X3)  3  2.3377  0.7792  0.4384 0.73796
## PQ(X1, X2, X3)   3  5.4570  1.8190  1.0234 0.47099
## Residuals        4  7.1094  1.7773                
## Lack of fit      3  5.2323  1.7441  0.9291 0.62419
## Pure error       1  1.8771  1.8771                
## 
## Stationary point of response surface:
##         X1         X2         X3 
## -0.1544199 -2.8319895 -0.1489394 
## 
## Stationary point in original units:
##          NaOH           HCl Precipitation 
##     1.8455801    -0.8319895    18.5106061 
## 
## Eigenanalysis:
## eigen() decomposition
## $values
## [1]  0.8754233 -0.1392041 -0.9361791
## 
## $vectors
##           [,1]        [,2]        [,3]
## X1  0.01516345 -0.06032313  0.99806372
## X2 -0.40642234 -0.91237213 -0.04896920
## X3  0.91355950 -0.40489285 -0.03835138

## 
## Call:
## rsm(formula = pH ~ SO(X1, X2, X3), data = pH.data)
## 
##              Estimate Std. Error t value  Pr(>|t|)    
## (Intercept) 11.143323   0.253875 43.8930 1.611e-06 ***
## X1           0.167588   0.126937  1.3202    0.2572    
## X2          -0.235950   0.126937 -1.8588    0.1366    
## X3           0.050812   0.126937  0.4003    0.7094    
## X1:X2       -0.006750   0.179516 -0.0376    0.9718    
## X1:X3        0.001875   0.179516  0.0104    0.9922    
## X2:X3        0.030300   0.179516  0.1688    0.8742    
## X1^2         0.243726   0.200705  1.2143    0.2914    
## X2^2         0.224551   0.200705  1.1188    0.3259    
## X3^2        -0.318474   0.200705 -1.5868    0.1878    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Multiple R-squared:  0.7475, Adjusted R-squared:  0.1795 
## F-statistic: 1.316 on 9 and 4 DF,  p-value: 0.4233
## 
## Analysis of Variance Table
## 
## Response: pH
##                 Df  Sum Sq  Mean Sq F value Pr(>F)
## FO(X1, X2, X3)   3 0.69072 0.230240  1.7861 0.2890
## TWI(X1, X2, X3)  3 0.00387 0.001290  0.0100 0.9984
## PQ(X1, X2, X3)   3 0.83218 0.277394  2.1519 0.2364
## Residuals        4 0.51562 0.128905               
## Lack of fit      3 0.43643 0.145478  1.8372 0.4859
## Pure error       1 0.07918 0.079184               
## 
## Stationary point of response surface:
##         X1         X2         X3 
## -0.3370911  0.5133518  0.1032032 
## 
## Stationary point in original units:
##          NaOH           HCl Precipitation 
##      1.662909      2.513352     21.032032 
## 
## Eigenanalysis:
## eigen() decomposition
## $values
## [1]  0.2443077  0.2243940 -0.3188980
## 
## $vectors
##            [,1]        [,2]         [,3]
## X1  0.985338030 -0.17060366 -0.001832879
## X2 -0.170587981 -0.98494800 -0.027877990
## X3 -0.002950796 -0.02778191  0.999609653

NaOH	HCl	Precipitation	C.NaOH	C.HCl	C.Precipitation	YieldActual	YieldPredict
3	2	10	1	0	-1	3.115500	3.283
2	3	10	0	1	-1	18.360972	18.133
2	1	30	0	-1	1	18.781400	19.011
1	2	10	-1	0	-1	3.447328	4.285
1	3	20	-1	1	0	5.231945	4.623
3	2	30	1	0	1	18.128745	17.291
2	1	10	0	-1	-1	7.761600	6.985
3	3	20	1	1	0	17.272564	17.335
1	1	20	-1	-1	0	5.268217	5.207
2	3	30	0	1	1	19.101212	19.879
3	1	20	1	-1	0	4.124200	4.735
1	2	30	-1	0	1	4.215000	4.049
2	2	20	0	0	0	20.156798	20.452
2	2	20	0	0	0	20.746737	20.452

## [1] 0.9974598

Optimasi Ekstraksi Larva Hermetia illucens dengan Response Surface Modelling dan Kualitas Asam Amino serta Aktivitas Antibakteri

Mohammad Miftakhus Sholikin

4/30/2019

Rancangan BBD; 2 center point

Data

Metode RSM

Pola Interaksi HCL; NaOH; Presipitasi >< pI

run.order	std.order	NaOH	HCl	Presipitasi
1	1	1	1	20
2	11	2	1	30
3	8	3	2	30
4	14	2	2	20
5	10	2	3	10
6	5	1	2	10
7	12	2	3	30
8	9	2	1	10
9	2	3	1	20
10	4	3	3	20
11	13	2	2	20
12	3	1	3	20
13	7	1	2	30
14	6	3	2	10

run.order	std.order	NaOH	HCl	Presipitasi
1	1	1	1	20
2	11	2	1	30
3	8	3	2	30
4	14	2	2	20
5	10	2	3	10
6	5	1	2	10
7	12	2	3	30
8	9	2	1	10
9	2	3	1	20
10	4	3	3	20
11	13	2	2	20
12	3	1	3	20
13	7	1	2	30
14	6	3	2	10

run.order	std.order	NaOH	HCl	Presipitasi
1	1	1	1	20
2	11	2	1	30
3	8	3	2	30
4	14	2	2	20
5	10	2	3	10
6	5	1	2	10
7	12	2	3	30
8	9	2	1	10
9	2	3	1	20
10	4	3	3	20
11	13	2	2	20
12	3	1	3	20
13	7	1	2	30
14	6	3	2	10