Abstract
Background: Carpobrotus edulis is widely used in traditional medicine to treat several health problems in South Africa. Phytochemical screening of medicinal plants involves several assays and requires a reasonable amount of starting material. Several factors, including plant species and variety, growth conditions, plant age and harvesting time, and plant part, can influence the extraction yield and quality of medicinal plant extracts.
Aim: This research aimed to optimise extraction yields from C. edulis using an analytical 23 full factorial design.
Setting: The plant material for this study was collected from the dunes along the Bloubergstrand coast (latitude: −33°47’49.92” S and longitude: 18°27’43.20” E) in the Western Cape Province, South Africa.
Methods: The factorial design involved two extractants (aqueous and methanolic), two pH levels (7 and 9) and two extraction temperatures (25°C and 40°C) and two extraction periods of 72 and 168 h. Regression models and factorial analysis provided a robust framework for evaluating these variables, including interaction terms that allowed for a better understanding of how combined conditions affect yields differently from individual conditions.
Results: The most effective conditions were a pH of 9 and a temperature of 40°C, with extraction times of 72 h for the aqueous extracts and 168 h for the methanolic extracts.
Conclusion: A factorial design can systematically optimise the exact extraction parameters to produce the highest yield of bioactive compounds from plants.
Contribution: Obtaining the best yield in medicinal plant research can serve several purposes, including conservation of the plant and the need to drive innovation in herbal medicine research.
Keywords: Carpobrotus edulis; extraction yield; 23 factorial design; phytochemical analysis; regression models.
Introduction
Carpobrotus edulis L. Bolus (family Aizoaceae) is a succulent, drought-resistant, spreading plant with triangular and fleshy leaves. The plant is an element of the Materia medica of South Africa and is used traditionally for treating mouth and throat infections, diphtheria, dysentery, diarrhoea and tuberculosis (De Beer & Van Wyk 2011; Van Wyk 2008). The plant is also used to treat sores, scalds and skin diseases (De Beer & Van Wyk 2011). Omoruyi, Bradley and Afolayan (2012) reported that the plant is used to treat diabetes, hypertension, intestinal parasites and infrequent bowel movements in some regions of South Africa.
Phytochemical screening of medicinal plants involves several assays and requires a practical amount of starting material. Obtaining the best yield can serve several purposes, including economics, conservation of the plant, supply security, quality of the product and the need to drive innovation in herbal medicine research. Several factors, including plant species and variety, growth conditions, plant age and harvesting time, and plant part, can influence the extraction yield and quality of medicinal plant extracts.
The efficiency of any extraction technique is influenced by independent variables, including potential of hydrogen (pH), extraction time, temperature, solvent and solid-to-liquid ratio. These values cannot be generalised for all medicinal plants because of the diversity in their composition and bioactive compounds. Knowledge of the independent variables and their interactions in the extraction process is required to understand the optimal method for the mining of phytochemicals. This information dictates the selection of variables and their values to obtain the maximum benefit from the technique, ensuring maximum efficiency (extraction yield). Traditionally, optimisation is achieved by examining the effect of a single variable at a time. Factorial designs for the extraction of plant compounds have been reported in our previous studies (Laloo et al. 2024; Terblanche et al. 2017) by several authors, including Kim et al. (2022), Suksaeree and Monton (2024), Cvitković et al. (2024) and Assaggaf et al. (2024).
This study aimed to investigate the effects of pH, extraction time and temperature on the yield of crude aqueous and methanolic extracts from C. edulis using two 23 full factorial designs. Furthermore, the extracts were screened to establish the presence of diverse phytochemicals obtained using the listed independent variables.
Research methods and design
Reagents used in this study
Many of the reagents were obtained from Sigma-Aldrich Corporation (St. Louis, MO, United States), including aluminum chloride, cholesterol, copper(II) acetate monohydrate, copper(II) sulphate, 3,5-dinitrobenzoic acid, Dragendorff reagent, gallic acid, iron(III) chloride, magnesium ribbon, methanol, 1-naphtol, ninhydrin, nitric acid, picric acid, pyridine, quercetin and sulphuric acid. Glacial acetic acid, chloroform, glucose, Lugol’s iodine solution and sodium chloride were purchased from Rochelle, and blood agar plates (Selecta-Media Columbia Agar, ThermoScientific) were donated by the Biomedical Technology division of the Department of Health Sciences (VUT, South Africa). Acetic anhydride, aspartic acid, ethanol, gelatin, glycine, hydrochloric acid, lead acetate 3-hydrate, potassium acetate, potassium iodide, sodium carbonate anhydrous, sodium citrate dihydrate, sodium hydroxide, sucrose and tyrosine were obtained from Merck (Germany). Iodine and mercury(II) chloride were obtained from Thomas Baker Chemicals (Mumbai, India) and Radchem (Pty) Ltd (Gauteng, South Africa), respectively. All chemicals were analytical grade. We used Ultrapure water that was obtained from a Millipore Direct-Q water purifier system (Merck Millipore, Darmstadt, Germany).
Collection, authentication and preparation of the plant material
The plant material for this study was collected from the dunes along the Bloubergstrand coast (latitude: −33°47’49.92” S and longitude: 18°27’43.20” E) in the Western Cape Province, South Africa. The plant was authenticated by Professor Stefan Siebert, curator of the AP Goossens Herbarium (Northwest University, Potchefstroom), where two voucher specimens (numbers 148 300 and 148 301) are housed. The fresh leaves of the plant were prepared as described by Laloo et al. (2024). The dried pulverised plant material was stored at −20°C in an airtight container for further analysis.
Extraction and experimental design
The variables affecting the two solvents used to extract phytocompounds from C. edulis were screened using a factorial design comprising three independent variables, pH, extraction time and temperature, each at two levels. Box 1 shows the independent variables, their designated letters and the levels used in this study.
| BOX 1: Independent variables and their levels used in the 23 factorial design. |
The solid-to-solvent ratio (1:20) was kept constant for all experimental runs. The experimental runs were randomised to minimise the effects of unexpected variability in the observed response. The extraction yield (%) was the dependent variable or response.
The independent variables (pH, extraction time and temperature) for each experimental run were chosen as dictated by the experimental design matrix (Table 1).
| TABLE 1: Design matrix, response and predicted values for the aqueous extraction conditions. |
A shaking incubator (Labotec, Model 355, Johannesburg, South Africa) was used for the extractions at 150 rpm. After filtering the extracts through Whatman No. 1 filter paper, the aqueous filtrates were lyophilised, and the methanolic extracts were concentrated using a rotary evaporator (HB 10 basic, IKA®-Werke GmbH & Co. KG, Staufen, Germany). The concentrated methanolic extracts were then dried further until the weight remained constant. The extracts were stored at −20°C in dark containers until further analysis. The extraction yield (%) was calculated using Equation 1:

Regression and statistical analysis
The experiments were performed in triplicate to estimate experimental error, reduce noise and minimise bias in the observed response values (Vaux, Fidler & Cumming 2012). The Design Expert version 11 (Stat-Ease Inc., MN, United States) statistical software was used to analyse the response variables, as described in Laloo et al. (2023) The following equation expressed the regression model in this study (Equation 2):

In Equation 2, the predicted response is represented by y, and βn represents the regression coefficient correlated with variable n. The regression coefficients were obtained from the analysis of the experimental data. The main variables are represented by A, B and C, the two-way interactions are by AB, AC and BC, the three-way interaction by ABC and the experimental error by e. Analysis of variance (ANOVA) was used to determine the statistical significance (p < 0.05) of the model, and the accuracy was determined using the coefficients of determination (R2). The Fisher F-test was used to determine whether the constructed models were adequate to describe the observed response values. The relationship between the independent variables (main effects), their intervariable effects (interaction effects) and the response variable was demonstrated through statistical plots.
Qualitative phytochemical screening
Qualitative phytochemical screening was performed on all reconstituted (50 mg/mL) aqueous (Terblanche et al. 2017) and methanolic crude extracts using standard methods (Godghate, Sawant & Sutar 2012; Nayak et al. 2011; Obouayeba et al. 2015; Samejo et al. 2013). The following phytochemicals were analysed: carbohydrates, proteins, amino acids and secondary metabolites, including alkaloids, anthocyanins, diterpenes, flavonoids, glycosides, phenols, phytosterols, saponins, tannins and terpenoids.
Ethical considerations
Ethical clearance to conduct this study was obtained from the Vaal University of Technology Faculty Research Ethics Committee (No. FACSREC-16102020-A00015).
Results and discussion
Extraction optimisation
This study determined how varying the pH, extraction time and temperature can affect the extraction yield of phytochemicals from the leaves of C. edulis using maceration. A factorial design shows the influence of the independent variables and their interactions. Table 1 shows the design matrix with the values of the independent variables, the experimental response values (extraction yield as a percent) and predicted values (based on the regression analysis) for each of the 16 experimental runs (8 aqueous extracts and 8 methanolic extracts).
Conspicuously higher extraction yields (41.70% – 64.21%) were obtained for the methanolic extracts than the aqueous extracts (21.55% – 31.11%). The amphiphilic nature of methanol can form intermolecular forces between polar and nonpolar functional groups, resulting in the solubilisation and extraction of a larger range of phytochemicals (Tiwari et al. 2011). Lapornik, Prošek and Wondra (2005) stated that methanol is more efficient in the degradation of cell walls, which have nonpolar characteristics, releasing more phytochemicals from the cells. Methanol was also found to be a better extraction medium by Magar et al. (2023) and Nortjie et al. (2024). The observation from this study underscores the importance of extracting solvents in medicinal plant research.
The effects of the independent variables on the extraction yield
Statistical plots (i.e. Pareto plots, main effect and interaction plots) for the water and methanol extraction conditions were generated to: (1) identify the statistically significant variables; and (2) describe the effect of each of the independent variables and their interactions on the percent extraction yield.
The weight and statistical significance of the independent variables and the intervariable interactions are displayed on Pareto charts (Figure 1), main effect plots (Figure 2) and interaction plots (Figure 3).
 |
FIGURE 1: Pareto chart for the (a) aqueous extraction conditions and (b) methanolic extraction conditions. |
|
 |
FIGURE 2: Main effect plots illustrating the influence of each independent variable on the extraction yield (%) for the aqueous (a–c) and methanolic (d–f) extraction conditions. |
|
 |
FIGURE 3: Interaction plots for the aqueous (a–c) and methanolic (d–f) extraction conditions. The main variables are represented by A, B and C, the two-way interactions are by AB, AC and BC, the three-way interaction by ABC. |
|
According to Ozbay and Yargic (2015), Pareto charts show the absolute values of the standardised effects in decreasing order and establish two limit lines, namely, the Bonferroni limit line (t-value of effect = 3.08209) and the t limit line (t-value of effect = 2.11991), which defines the minimum statistically significant effect magnitude for the 95% confidence level. Regression coefficients with t-value effects above the Bonferroni line are labelled as significant, while those with t-value effects below the t-limit line are considered statistically nonsignificant (Shah & Pathak 2010). If the t-value of the effect of a regression coefficient lies between the Bonferroni and t-limit lines, the coefficient is likely to be significant. Main effect and interaction plots refer to simple line graphs obtained from connecting the mean values of each treatment. An effect with a zero-slope horizontal line (main effect plots) and an interaction displaying parallel lines (interaction plots) were interpreted as not significant effects, whereas increasing (positive effect – the response is higher at the high level) or decreasing (negative effect – the response is lower at the high level) effect lines were considered significant.
Figure 1 for the aqueous extraction conditions revealed that the extraction time (B) had the most significant effect in a negative mode (Figure 2) on the extraction yield (%), whereas it showed a positive significant effect (Figure 2) on the percent extraction yield (%) for the methanolic extraction conditions. The temperature (°C) had the most significant effect in a positive mode (Figure 2) on the extraction yield (%) for the methanolic extraction conditions (Figure 1) and exhibited a significantly positive effect (Figure 2) for the aqueous extraction conditions. The interactive effects, AB (pH*extraction time) and BC (extraction time*temperature), demonstrated both negative significant effects (Figure 3) for the aqueous extraction conditions but displayed no significant effects (Figure 3) on the extraction yield (%) for the methanolic extraction.
The contribution of pH (A) was positively significant (Figure 2) for the methanolic extraction conditions but was not significant (Figure 2) for the aqueous extraction conditions. The interaction between pH and temperature (AC) exhibited a weak, positively significant effect (Figure 3) on the extraction yield for the aqueous extraction conditions but had no significant role (Figure 3) during the methanolic extraction conditions. The three-way interaction ABC (pH*extraction time*temperature) was above the Bonferroni limit and displayed a significantly negative effect on the extraction yield (%) for both the aqueous and methanolic extraction conditions.
Factorial regression model development
After fitting the obtained response values for each extraction condition tested to the regression model equation (Equation 2) and performing ANOVA to determine the accuracy of the suggested models, the p-values for the three independent variables and the intervariable interactions were obtained and are presented in Table 2 (aqueous), Table 3 (aqueous), Table 4 (methanolic) and Table 5 (methanolic), respectively. The variables and interactions with p-values greater than 0.05 were considered not significant and excluded from the model development. These mathematical models, based on the regression equation (Equation 2), were built and are presented in Equation 3 (aqueous conditions) and Equation 4 (methanolic conditions):


where ŷ represents the predicted extraction yield, ET the extraction time and T the temperature.
| TABLE 2: Analysis of variance report for the full model for the aqueous extraction conditions. |
| TABLE 3: Analysis of variance report for the reduced model for the aqueous extraction conditions. |
| TABLE 4: Analysis of variance report for the full model for the methanolic extraction conditions. |
The regression coefficients (βn) were obtained from the ANOVA and describe the potential influence of each of the independent variables and their interactions on the response value (Montgomery 2013). The higher the value of the regression coefficient, the more significant the corresponding effect (Alcheikh Hamdon, Darwish & Hilal 2015). The positive sign of the regression coefficients for the aqueous extraction condition (β0, β3 and β13) and the methanolic extraction condition (β1, β2 and β3) is indicative of a synergistic effect, whereas the negative sign of the regression coefficients β2, β12, β23 and β123 (aqueous) and β123 (methanolic) denotes an antagonistic effect.
Significance and validation of the regression models
The Fisher F-test was used to evaluate the statistical significance of the constructed regression models. The F-values of 24.06 for the full model (Table 2) and 29.46 for the reduced model (Table 3), both with p-values of < 0.0001, suggested that the models were highly significant in describing the response values for the aqueous extraction. The improved F-value obtained for the reduced model was because of the elimination of the nonsignificant model terms.
The coefficient of determination (R2) and adjusted coefficient of determination ( ) were used to validate the accuracy of the model. The coefficient of determination (R2) determines the proportion of variability in the response value that can be explained by the independent variables, indicating the proximity of the predicted values to the fitted regression line (Ogee et al. 2013). The R2 for the full (R2 = 0.9132) and reduced (R2 = 0.9123) models (aqueous extraction conditions) indicated that the predicted values from the models are a good fit.
Note that the accuracy of the model should not solely be based on R2, as R2 always increases with an increase in the number of model variables, even if these variables are not significant (Kukreja et al. 2011). The presence of multiple variables and intervariable effects warrants the evaluation of the adjusted R2 values, which compensates for the introduction of additional variables on the normal R2 calculation by only increasing the value, should the additional variables significantly contribute to the model (Grant & Kenton 2019). The value of for the full ( = 0.8753) and reduced ( = 0.8813) models, even though lower than R2, still affirms an acceptable model fit.
Similar to the validation criteria used for the aqueous extraction conditions, F-tests, R2 and were used to ascertain the significance and validity of the regression models for the methanolic extraction conditions. The F-values of 58.96 for the full model (Table 4) and 107.17 for the reduced model (Table 5), both with p-values of < 0.0001, affirm the statistical significance in describing the response values for the methanolic extraction conditions. The R2 for the full (R2 = 0.9627) and reduced (R2 = 0.9576) models (methanolic extraction conditions) points to a better model prediction than that of the aqueous extraction conditions. This is underscored by the improved calculations ( = 0.9464) and reduced ( = 0.9486) for the methanolic extraction models.
| TABLE 5: Analysis of variance report for the reduced model for the methanolic extraction conditions. |
Numerical optimisation of the extraction conditions
The regression models for maximum extraction yield were optimised by using the maximum values of the aqueous and methanolic extraction conditions. The optimal yield for the aqueous model was obtained at a pH of 9, an extraction time of 72 h and a temperature of 40°C. The methanolic model produced the best yield at a pH of 9, extraction time of 168 h and temperature of 40°C.
Qualitative phytochemical analysis
Table 6 shows the results of the qualitative phytochemical screening of the extracts, showing the impact of the various extraction conditions. The extraction variables did not affect the presence of carbohydrates (Molisch and Benedict’s test), phenols, tannins, flavonoids, cardiac glycosides, diterpenes, phytosteroids and triterpenoids because they were present in all extracts under all experimental conditions. Many of the therapeutic properties of C. edulis are attributed to the presence of these metabolites. Phytochemical tests for the xanthoproteic reaction, ninhydrin alkaloids, anthraquinone glycosides, leucoanthocyanins, haemolysin test for saponins and phlobatannins were negative for both extracts. On the contrary, anthocyanins, starch, phytosteroids, triterpenoids, the Biuret test for proteins and saponins were present in some of the experimental runs. The sulphur test showed the presence of proteins only in the methanolic extracts. The presence or absence of these phytochemicals provides evidence that a factorial design is essential in obtaining a substantial quantity of a desired active compound from plants. The results also show that certain solvents are required to isolate specific compounds, as shown in the case of sulphur in this study. The negative results for starch, aromatic amino acids, free amino acids, alkaloids, anthocyanins, leucoanthocyanins, anthraquinone glycosides, saponins and phlobatannins may be because of the choice of extraction solvent or the absence of these compounds in detectable quantities in the extracts.
| TABLE 6: Qualitative phytochemical screening of aqueous and methanolic crude extracts of Carpobrotus edulis. |
Conclusion
Two 23 full factorial designs were used to determine the contribution of pH, extraction time and temperature to the extraction yield (%) of phytochemicals from the leaves of C. edulis. Statistical analysis of data for the aqueous extraction of C. edulis showed that the percentage extraction yield was significantly affected by the extraction time and temperature and their intervariable interactions. Although the effect of pH was not statistically significant, the pH*extraction time and pH*temperature interactions strongly affected the extraction yield. The methanolic extraction of C. edulis strongly depended on all three investigated extraction parameters (pH, extraction time and temperature) as well as the three-way interaction between them. The presence or absence of some of the phytochemicals observed in this study suggests the power of a factorial design to isolate specific compounds from C. edulis.
Acknowledgements
The authors are grateful to the Research Directorate of the Vaal University of Technology, South Africa, for their financial support for this study.
This article is based on research originally conducted as part of Unisa Terblanche’s Doctoral thesis titled ‘The effect of Carpobrotus edulis and Cotyledon orbiculata on the migration, proliferation and differentiation of keratinocytes in vitro’, submitted to the Faculty of Applied and Computer Sciences, Department of Natural Sciences (previously Department of Biotechnology), Vaal University of Technology, in 2020. The thesis was supervised by Prof. Michael Pillay, and co-supervised by Prof. Cornelius Cano Ssemakalu and Prof. Fanyana Mtunzi. The manuscript has since been revised and adapted for journal publication. The original thesis is available at: https://hdl.handle.net/10352/809.
Competing interests
The authors reported that they received funding from the Vaal University of Technology, which may be affected by the research reported in the enclosed publication. The authors have disclosed those interests fully and have implemented an approved plan for managing any potential conflicts arising from their involvement. The terms of these funding arrangements have been reviewed and approved by the affiliated university in accordance with its policy on objectivity in research.
Authors’ contributions
U.T. and C.C.S. designed the research, analysed the data and wrote the manuscript. N.L. contributed to the research. M.P. supervised the work.
Funding information
Funding was provided by the Vaal University of Technology.
Data availability
The data that support the findings of this study are available from the corresponding author, M.P., upon reasonable request.
Disclaimer
The views and opinions expressed in this article are those of the authors and are the product of professional research. They do not necessarily reflect the official policy or position of any affiliated institution, funder, agency or that of the publisher. The authors are responsible for this article’s results, findings and content.
References
Alcheikh Hamdon, A.A., Darwish, N.A. & Hilal, N., 2015, ‘The use of factorial design in the analysis of air-gap membrane distillation data’, Desalination 367, 90–102. https://doi.org/10.1016/j.desal.2015.03.014
Assaggaf, H., Jeddi, M., Mrabti, H.N., Ez-zoubi, A., Qasem, A., Attar, A. et al., 2024, ‘Design of three-component essential oil extract mixture from Cymbopogon flexuosus, Carum carvi, and Acorus calamus with enhanced antioxidant activity’, Scientific Reports 14, 9195. https://doi.org/10.1038/s41598-024-59708-x
Cvitković, D., Škarica, I., Dragović-Uzelac, V. & Balbino, S., 2024, ‘Supercritical CO2 extraction of fatty acids, phytosterols, and volatiles from myrtle (Myrtus communis L.) fruit’, Molecules 29, 1755. https://doi.org/10.3390/molecules29081755
De Beer, J.J.J. & Van Wyk, B.-E., 2011, ‘An ethnobotanical survey of the Agter–Hantam, Northern Cape Province, South Africa’, South African Journal of Botany 77, 741–754. https://doi.org/10.1016/j.sajb.2011.03.013
Godghate, A., Sawant, R. & Sutar, A., 2012, ‘Phytochemical analysis of ethanolic extracts of roots of Carrisa carandus Linn’, Rasayan Journal of Chemistry 5, 456–459.
Grant, M. & Kenton, W., 2019, How the coefficient of determination works, Investopedia, viewed 25 June 2019, from https://www.investopedia.com/terms/c/coefficient-of-determination.asp.
Kim, M.-K., Park, G., Ji, Y., Lee, Y.-G., Choi, M., Go, S.-H. et al., 2022, ‘Design of experiments-based optimization of flavonoids extraction from Daphne genkwa flower buds and flavonoids contents at different blooming stages’, Plants 11(7), 925. https://doi.org/10.3390/plants11070925
Kukreja, A., Chopra, P., Aggarwal, A. & Khanna, P., 2011, ‘Application of full factorial design for optimization of feed rate of stationary hook hopper’, International Journal of Modeling and Optimization 1(3), 205–209. https://doi.org/10.7763/IJMO.2011.V1.36
Laloo, N., Terblanche, U., Ssemakalu, C. & Pillay, M., 2024, ‘The effect of solvent, pH, extraction time and temperature on the extraction of phenolic compounds and antioxidant activity of Carpobrotus edulis’, Journal of Phytology 16, 1–7. https://doi.org/10.25081/jp.2024.v16.8393
Lapornik, B., Prošek, M. & Wondra, A.G., 2005, ‘Comparison of extracts prepared from plant by-products using different solvents and extraction time’, Journal of Food Engineering 71(2), 214–222. https://doi.org/10.1016/j.jfoodeng.2004.10.036
Magar, A.B., Shrestha, D., Pakka, S. & Sharma, K.R., 2023, ‘Phytochemistry, biological, and toxicity study on aqueous and methanol extracts of Chromolaena odorata’, Scientific World Journal 2023, 6689271. https://doi.org/10.1155/2023/6689271
Montgomery, D.C., 2013, Design and analysis of experiments, 8th edn., John Wiley & Sons, Inc, Hoboken, NJ.
Nayak, B.S., Kanhai, J., Milne, D.M., Pereira, L.P. & Swanston, W.H., 2011, ‘Experimental evaluation of ethanolic extract of Carapa guianensis L. leaf for its wound healing activity using three wound models’, Evidence-Based Complementary and Alternative Medicine 2011, 419612. https://doi.org/10.1093/ecam/nep160
Nortjie, E., Basitere, M., Moyo, D. & Nyamukamba, P., 2024, ‘Assessing the efficiency of antimicrobial plant extracts from Artemisia afra and Eucalyptus globulus as coatings for textiles’, Plants 13(4), 514. https://doi.org/10.3390/plants13040514
Obouayeba, A.P., Diarrassouba, M., Soumahin, E.F. & Kouakou, H., 2015, ‘Phytochemical analysis, purification and identification of Hibiscus anthocyanins’, Journal of Pharmaceutical, Chemical and Biological Sciences 3, 156–168.
Ogee, A., Ellis, M., Scibilia, B., Pammer, C. & Steele, C., 2013, ‘Regression analysis: How do I interpret R-squared and assess the goodness-of-fit?’, Minitab Blog, viewed 24 February 2019, from http://blog.minitab.com/blog/adventures-in-statistics-2/regression-analysis-how-do-i-interpret-r-squared-and-assess-the-goodness-of-fit.
Omoruyi, B.E., Bradley, G. & Afolayan, A.J., 2012, ‘Antioxidant and phytochemical properties of Carpobrotus edulis (L.) Bolus leaf used for the management of common infections in HIV/AIDS patients in Eastern Cape Province’, BMC Complementary & Alternative Medicine 125, 1186. https://doi.org/10.1186/1472-6882-12-215
Ozbay, N. & Yargic, A.S., 2015, ‘Factorial experimental design for Remazol yellow dye sorption using apple pulp/apple pulp carbon–titanium dioxide co-sorbent’, Journal of Cleaner Production 100, 333–343. https://doi.org/10.1016/j.jclepro.2015.03.050
Samejo, M.Q., Sumbul, A., Shah, S., Memon, S.B. & Chundrigar, S., 2013, ‘Phytochemical screening of Tamarix dioica Roxb. ex Roch’, Journal of Pharmacy Research 7, 181–183. https://doi.org/10.1016/j.jopr.2013.02.017
Shah, M. & Pathak, K., 2010, ‘Development and statistical optimization of solid lipid nanoparticles of simvastatin by using 23 full-factorial design’, AAPS PharmSciTech 11, 489–496. https://doi.org/10.1208/s12249-010-9414-z
Suksaeree, J. & Monton, C., 2024, ‘Maximizing curcuminoid extraction from Curcuma aromatica Salisb. rhizomes via environmentally friendly microwave-assisted extraction technique using full factorial design’, International Journal of Food Science 2024, 4566123. https://doi.org/10.1155/2024/4566123
Terblanche, U., Ssemakalu, C., Mtunzi, F. & Pillay, M., 2017, ‘Screening of variables influencing the extraction yield of Cotyledon orbiculata: 23 full factorial design’, International Journal of Pharmacognosy and Phytochemical Research 9, 303–312.
Tiwari, P., Kumar, B., Kaur, M., Kaur, G. & Kaur, H., 2011, ‘Phytochemical screening and extraction: A review’, Internationale Pharmaceutica Sciencia 1, 98–106.
Van Wyk, B.-E., 2008, ‘A review of Khoi-San and Cape Dutch medical ethnobotany’, Journal of Ethnopharmacology 119(3), 331–341. https://doi.org/10.1016/j.jep.2008.07.021
Vaux, D.L., Fidler, F. & Cumming, G., 2012, ‘Replicates and repeats – What is the difference and is it significant? A brief discussion of statistics and experimental design’, EMBO Reports 13(4), 291–296. https://doi.org/10.1038/embor.2012.36
|