Next Article in Journal
The Effect of Zinc, Copper, and Silver Ions on Oat (Avena sativa L.) Androgenesis
Previous Article in Journal
The Ubiquitin Switch in Plant Stress Response
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 10
Issue 2
10.3390/plants10020247
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Machine Learning Applications and Optimization of Clustering Methods Improve the Selection of Descriptors in Blackberry Germplasm Banks

by
Juan Camilo Henao-Rojas
1,
María Gladis Rosero-Alpala
1,
Carolina Ortiz-Muñoz
1,
Carlos Enrique Velásquez-Arroyo
1,
William Alfonso Leon-Rueda
2 and
Joaquín Guillermo Ramírez-Gil
2,*
1
Corporación Colombiana de Investigación Agropecuaria—AGROSAVIA, Centro de Investigación La Selva- Km 7, 250047 Ríonegro, Colombia
2
Departamento de Agronomía, Facultad de Ciencias Agrarias, Universidad Nacional de Colombia, 111321 Sede Bogotá, Colombia
*
Author to whom correspondence should be addressed.
Plants 2021, 10(2), 247; https://doi.org/10.3390/plants10020247
Submission received: 6 October 2020 / Revised: 27 November 2020 / Accepted: 30 November 2020 / Published: 28 January 2021
Browse Figures
Review Reports Versions Notes

Abstract

:
Machine learning (ML) and its multiple applications have comparative advantages for improving the interpretation of knowledge on different agricultural processes. However, there are challenges that impede proper usage, as can be seen in phenotypic characterizations of germplasm banks. The objective of this research was to test and optimize different analysis methods based on ML for the prioritization and selection of morphological descriptors of Rubus spp. 55 descriptors were evaluated in 26 genotypes and the weight of each one and its ability to discriminating capacity was determined. ML methods as random forest (RF), support vector machines, in the linear and radial forms, and neural networks were optimized and compared. Subsequently, the results were validated with two discriminating methods and their variants: hierarchical agglomerative clustering and K-means. The results indicated that RF presented the highest accuracy (0.768) of the methods evaluated, selecting 11 descriptors based on the purity (Gini index), importance, number of connected trees, and significance (p value < 0.05). Additionally, K-means method with optimized descriptors based on RF had greater discriminating power on Rubus spp., accessions according to evaluated statistics. This study presents one application of ML for the optimization of specific morphological variables for plant germplasm bank characterization.

1. Introduction

Machine learning (ML) is a form of artificial intelligence (AI) that gives machines the ability to learn through the use of algorithms and a training process [1] and is used in tandem with big data technologies and high-performance computing [2,3], which, together with information and communication technologies (ICTs) and the Internet of Things (IoT), deep learning, among others tools, have created new opportunities for data-intensive science. These tools are being applied in multiple areas, including agriculture, as an emerging technology [2,4]. These tools helped create what is now known as digital agriculture, a new agricultural revolution [2,5].
The practical applications of ML in agriculture are broad and are applied in various fields, such as yield prediction [6], pest detection, classification, monitoring, and management [3,7], and species recognition [8], which is why ML is undoubtedly a powerful tool within digital agriculture, providing information for correct decision-making including the study and efficient use of tropical genetic resources. The ML approach in agriculture presents many challenges, because it is highest impacted by environmental factors that cannot be controlled, demanding rigorous process and extensive data for validation and testing [9]. This implies avoiding the overfitting, select very complex models and data without standardization, improve the statists used, and data security, among others [2,9,10]. The goal of good practices in ML could improve real applications in agriculture with biological, ecological, agronomic, economic, cultural, and social contexts.
The subcategories of ML algorithms present different approaches, with variations in mathematical and statistical principles, computational requirements, advantages and disadvantages, and current and potential applications [9,10]. The most common algorithms associated with ML such as random forests, decision trees, supper vector machine, Bayesian networks, neuronal networks, regression analysis, among others are named supervised learning [9]. On the other hand, when unlabeled datasets without prior knowledge of the input and output variable response are named unsupervised learning such as Artificial Neural Networks, genetic algorithm, Instance-based learning models, deep learning, and clustering [2,7,9,10].
Classification is one of the most common processes in which ML algorithms are used. For this purpose, decision trees are the most popular tools, highlighting the random forest algorithm [2,10]. On the other hand, supper vector machine and neuronal networks are algorithms that have been used for this objective [2,9]. In addition, in recent years deep learning has gained popularity in the classification processes of biological phenomena [3,7].
Colombia is one of the more biodiverse countries on the planet because it is rich in natural resources and has a variety of species and ecosystems [11]. In order to conserve and preserve high genetic diversity in species of agricultural interest, the System of Germplasm Banks of the Nation for Food and Agriculture (SBGNAA) was created by the Colombian Agricultural Research Corporation, AGROSAVIA. Banks include collections of species of high value for the nation as a result of their potential for agribusiness, pharmaceuticals, textiles, and food security and sovereignty [12,13]. The SBGNAA genetic pool includes an ex situ blackberry collection (Rubus spp.) that is made up of different species with high potential for the domestic and export markets given their nutritional and functional characteristics [14,15].
The value of conserved plant genetic resources is given by the information used to promote their use, where morphological and genetic descriptors are the basic parameters that contribute to the knowledge of variability and parental selection, providing useful information for the development of new cultivars and the protection of varieties [16]. Morphological characterization is the determination of a set of phenotypic traits measurement based on key botanical-taxonomic named as descriptors. The descriptors should have a highly polymorphic expression, stable in the population and expressible under different environmental conditions [17].
Currently, the Colombian Rubus spp. collection does not have validated descriptors. Therefore, a series of morphological, chemical, and physiological descriptors is needed in order to provide easy-to-measure characteristics that relate the shape, structure or functions of an accession, allowing compliance with distinguishability, uniformity, and stability requirements [18].
After a proper validation process for descriptors in germplasm banks, the phenotypic variability classification stage continues, which has traditionally been based on multivariate analysis, commonly using principal components, clustering methods such as cluster analysis, and linear and quadratic discriminants [19]. These tools have been used to recognize cultivated and wild species of Rubus spp. and commercial crops of Rubus glaucus Benth [20,21], without a prior optimization and validation process for descriptors or evaluating which discrimination method had the best performance.
Therefore, given the need to have morphological descriptors for the Rubus spp. germplasm bank along with the lack of knowledge on the ability of other discrimination methods to provide specific information on possible descriptors, machine learning tools, and discrimination method optimization have provided an opportunity for the creation of more efficient processes in terms of information capacity, computational performance and optimization of work time in the field and laboratory. The present study aimed to evaluate different machine learning techniques for the selection of morphological variables in 26 Rubus spp. accessions cultivated under tropical conditions and for the optimization of discrimination methods to increase the differentiation capacity of accessions.

2. Materials and Methods

2.1. Rubus spp. Accession Genotypes and Agronomic Management

A total of 26 blackberry accessions (Rubus spp.) were selected from the germplasm bank at the La Selva Research Center in Rionegro, Antioquia—Colombia located at an altitude of 2100 masl with an average temperature of 16 °C and 74.83% relative humidity (Longitude: 075°24′51.9 and Latitude: 06°07′52.7), a lower montane humid forest life zone (bh-MB).
The selected accessions had interspecific morphological variation at the species level, obtained as a representative sample of all variation of the Rubus spp. germplasm bank collection. The origin of each accession was as follows: Cultivated natives (CN) with 18 accessions; introduced (I) with 6 accessions; and wild native (WN) with 2 accessions (Table 1). The plant material was collected in the phenological phases of flowering; the primary floriferous stems were obtained after renewal pruning to guarantee their quality. For each accession, five plants were selected, which had three floriferous primary stems extracted from each one. The sanitary and renewal pruning were standard for all accessions, as well as the fertilization practices. When there was a water deficit, water was added with automatic drip irrigation.

2.2. Morphological Descriptors

The selected morphological descriptors were based on previous studies on morphological characterizations carried out by Evans and Weber [22], combined with the proposals of Ligarreto-Moreno, Espinosa, Barrero and Medina [21] for species from the Rubus L genus. Once selected, they were grouped using typology criteria, botanical terminology, and measurement scales. These descriptors were verified with taxonomic keys of the genus Rubus L. [20] and the characters for the distinction of varieties of blackberry [21]. This validation confirms that morphological traits are associated with stable phenotypic expressions into Rubus genus and low environmental effect.
The descriptors were recorded at the plant and organ levels, including qualitative variables such as growth habit, organ shape, and other vegetative structures. In addition, the organs were colored using the Royal Horticultural Society’s color chart for plants [23]. The quantitative variables focused on morphometric measurements such as length (linear dimension) and length and width (horizontal dimensions). Discrete variables were also used for the vegetative structures in the plant organs (Table 2).

2.3. Selection and Optimization of Morphological Descriptors from the Rubus spp. Germplasm Bank Using Machine Learning Tools

The Rubus spp. germplasm bank currently does not have validated morpho-agronomic descriptors, meaning much methodological effort is needed for characterization, which is why three ML tools were used to discriminate and determine the weight of the descriptors in terms of diversity for optimal clustering of the evaluated accessions. This study used an internal optimization process from ML algorithms random forest (RF), linear (SVMl) and radial (SVMr) support vector machines, and neural networks (NN).
The data matrix was associated with 55 columns (morphological descriptors) per 78 rows (26 accessions with 3 replicates by each one). In addition, one-hot encode data was realized through a numerical representation of categorical descriptors. The arrays data were randomly divided into two data sets: (i) training (75%) and (ii) testing (25%) to be evaluated under different ML methods.
The RF algorithm is widely used as a classifier given its simplicity in terms of the parameters required for decision making [24]. It was implemented using the libraries randomForest [25], caret [26], ranger [27], and h2o [28] available to free software R.
The optimization of the RF algorithm and determination of the importance of the morphological descriptors (numerical and non-numerical) were carried out through a multi-step analysis. The first step was evaluating the number of trees (set of combinations from 100 to 1000), the size through the number of nodes (1 to 4000), and the hyperparameter alpha (0 to 20) with of the error rate, Bayes error and out-of-bag root-mean-square error (OOB-RMSE), using the caret, ranger and h2o libraries, selecting 500 trees, 4000 nodes and an alpha of 4 as a balance between the robustness, stabilization of the error rate and computational performance. Subsequently, the importance of the morphological descriptors was determined by calculating the confusion matrices [29] and then their global accuracy [30] using the “caret” package. These parameters identified the weight and interactions within the decision tree for each variable. As a complement, the importance of the descriptors was evaluated by calculating the Mean Decreases in the Gini Index metric. This process was corroborated using the metrics: mean minimal depth, number of nodes, mse increase, node purity increase, number of trees times a root, root variable, interaction occurrences, uncond mean minimal depth and significance (p < 0.05), evaluated with the randomForestExplainer library [31]. Finally, these metrics were graphed. The type of assembly used in the RF algorithm for the selection of the models was Bagging.
The second evaluated algorithm was SVM, which belongs to the general category of kernel methods, widely used in classification and regression because of its high precision and capacity to handle high-dimensional data [32]. The linear and radial forms of the SVM algorithm (SVMl and SVMr) were evaluated using the e1071 library [33], optimized with two steps: (i) selection of parameters and (ii) final training [34]. To select the fitness parameters, the classification error, and the mean square error of the regression were determined using the e1071 library. In the case of SVMr, gammas of 0.5, 1, 2, 3, 4, and 5 were tested, and the computational cost was determined with the sum of the hyperparameters and an indirect measurement of the computational simplicity of the code.
The third algorithm was NN, commonly used in identifying and predicting patterns between multiple variables [35]. It was implemented with the nnet [36], NeuronalNetTools, and RSNNS packages [37]. The importance and sensitivity of the descriptors were evaluated using the Garson and Olden algorithms and the Lek profile method; then, the network was optimized through a step-by-step process [35]: (i) Normalize entries, standardize responses, and evaluate the influence of outliers. (ii) Network architecture, which includes the size or number of units in the hidden layer, the number of nodes in each layer, inclusion of bias layers, and weights or inputs. (iii) Decay by decreasing the specific weight of the regularization in the neural network and (iv) interactions by evaluating different amounts of interactions. Additionally, the correlation matrix between variables was calculated [35].
With the optimization of each algorithm (RF, SVMr, SVMl, and NN), the results obtained in the selection of descriptors and its ability to classify adequately the accessions of Rubus spp. were evaluated based on training/validation accuracy and training/validation missclass Error using area under receiver operating characteristic (ROC) curve (AUC) [38] implemented in the free software R [39] with own code. The accuracy quantified by AUC it is considered a good metric that has been used in the comparison of ML algorithms [40].

2.4. Rubus spp. Germplasm Bank k Genotype Clustering Methods

The tools selected for the analysis of discrimination of the accessions in the Rubus spp. germplasm bank included hierarchical agglomerative clustering, using the Ward D2 method as the clustering union strategy, and k-means [41]. The analyses were developed in the R software [39] through the creation of an own code with the help of the libraries vegan [42], pvclust [43], ape [44], and rgl [45]. The consolidated morphological descriptors in Table 2 were used as discriminant variables.
In order to optimize the discrimination of genotypes based on morphological descriptors, a multistage analysis was carried out: (i) Standardization of variables carried out by means of the Z score using the clusterSim library [46], guaranteeing equal or similar measurement scales especially in measurements of dissimilarity sensitive to magnitude such as the Euclidean distance [47,48]. (ii) For both discrimination methods, the optimal number of clusters was estimated with the gap statistic [49] using the libraries factorextra [50] and stat [51]. (iii) Since there are no validated morphological descriptors for the germplasm bank, the effect of the number of variables in the discrimination methods was determined using two data sets: (a) all variables (Table 2) and (b) those selected using the RF algorithm (V23, V44, V24, V49, V25, V27, V50, V29, V28, and V42).
The variations in the hierarchical agglomerative clustering and k-means methods were evaluated using the statistics normalized variation index (NVI), Adjusted Rand index (ARI), separation index (IS), Calinski-Harabasz index (CH), entropy (EN), and Pearson Gamma (PG) [52]. The parameters were determined with the fpc v2.2-5 [53] and aricode v0.1.2 libraries [54], implemented in the free R software [39].
All combinations of morphological descriptors used on discriminating methods included the results of ML algorithms were evaluated for their ability to discriminate each accession used as replicas in order to avoid the artifacts such as the environmental condition in differential expression of morphological markers. In addition, the results were taxonomically corroborated to detect anomalies in the algorithms evaluated.

3. Results

3.1. Selection and Optimization of Machine Learning Algorithms for the Prioritization and Selection of Morphological Descriptors in Rubus spp.

Table 3 shows the results of the ability of each of the machine learning methods to prioritize the importance of morphological variables and their discriminating ability in the studied blackberry accessions after an internal optimization process. It was determined that the RF algorithm had the best performance based on the test statistics. In decreasing order of ability to discriminate adequately based on the appropriate selection of descriptors, the algorithms were RF (descriptive and numerical variables), RF (numerical variables), neuronal networks (NN), support vector machine (SVM) radial (SVMr) and linear (SVMl) with area under curve (AUC) accuracy classification values of 0.76, 0.64, 0.31, 0.21, and 0.09, respectively.
The RF algorithm, after optimization based on the number of trees, and size and reduction of the hyperparameters, was able to determine that the quantitative and qualitative descriptors of greatest importance for use in the discrimination of blackberry genotypes were those with the highest value for the mean decrease in the Gini index (Figure 1a and Figure 2a), minimum depth within the decision forest (Figure 1b and Figure 2b), maximum number of connected nodes (Figure 1c and Figure 2c), increase in purity (Figure 1d and Figure 2d) and number of most frequent interactions in the decision trees (Figure 1e and Figure 2e). According to these criteria, the variables in decreasing order of importance for quantitative and quantitative-descriptive descriptors were V44, V24, V42, V49, V25, V29, V50, V27, V28, V26 and V23, V44, V9, V24, V53, V22, V25, V49, V28, V50, and V26 (Figure 2a–e).
The optimized SVMl, varying parameter C, as a balance between the massification of the algorithm margin and the error with a selection factor of 0.04, determined that the descriptors in increasing order of importance were V51, V29, V19, V7, V27, V32, V46, V31, and V23 (Figure 3a,b). The SVMr found that the selection factor 50 in the Kernel function and a gamma of 0.5 minimized the error at the lowest possible computational cost, which prioritized the descriptors in decreasing order as V19, V37, V46, V48, V23, V30, V42, V44, V51, V7, V32 and V31 (Figure 3b).
NN found that the best relationship as a function of the decrease in the specific weight of the regularization of the network with respect to the number of hidden layers quantified using the Bootstrap indicator determined that the highest value of occurrences was seen with a weight of decay of 0.1 and a number of layers of 5.0 (Figure 4b). The relationship between the variables and the network determined using Olsen’s connection weights algorithm showed that the descriptors with importance greater than 2.5 in absolute value were V30, V34, V32, V31, V49, V12, V50, V26, V33, V36, V48 and V19 (Figure 4a). Additionally, as a result of the correlation analysis, two groups were found: (i) directly proportional relationships and (ii) inversely proportional relationships (Figure 4c). The descriptors associated with the flower, such as V49 and V50 (inversely proportional), helped differentiate the accessions at the inter- and intraspecific levels. The leaf descriptors, such as V31 and V32 (directly proportional), were able to discriminate using a single descriptor. Grouping also occurred in five clusters associated with the importance value of the variables: High importance: V34 to V44, medium importance: V26 to V33, and low importance V12 to V25. The importance of variables V27, V37, V43, and V44 was notable in all groups (Figure 4d).

3.2. Genotype Discriminating Methods from the Rubus spp. Germplasm Bank

It was found that the massification of the Gap statistic was obtained when the number of clusters for the standardized hierarchical agglomerative clustering, non-standardized hierarchical agglomerative clustering, standardized k-means, and non-standardized k-means methods were 10, 6, 10, and 7, respectively (Figure 5), indicating a uniform, non-random distribution of the accessions within each group. Additionally, the effect of standardization on the two methods tested was notable (Table 4). Superior performance was found in the K-means and hierarchical agglomerative clustering grouping when the reduction of variables was carried out using the RF method, indicating that the selection process was highly informative (Figure 6 and Table 4).
The evaluation of the behavior and discrimination capacity of the hierarchical agglomerative clustering and k-means methods with their different evaluated variations based on the Normalized Variation Index (NVI); Adjusted Rand Index (ARI); Separation index (IS); Calinski-Harabasz Index (CH); Entropy (EN); Pearson Gamma (PG) indices showed that the method that presented the best performance was K-means with standardized and reduced variables as a function of the RF optimization process, followed by the K-means method with all standardized variables, standardized agglomerative hierarchical grouping with optimization of variables using RF, agglomerative hierarchical grouping with all standardized variables, K-means and non-standardized agglomerative hierarchical grouping (Table 4).
Figure 7 shows the relationship between the descriptors prioritized by the RF algorithm and the results of the discrimination using the K-means method, without contradictory variations in the discriminating morphological characteristics for each group or accession.

4. Discussion

4.1. Selection and Discrimination of Descriptors Applied to the Rubus spp. Germplasm Bank Using Machine Learning Tools

Of the analyzed machine learning tools, RF presented the best performance when compared to SVMl, SVMr, and NN for the ability to prioritize highly discriminating descriptors of accessions from the Rubus genus. This algorithm has been widely used in classification processes given its good behavior, simplicity in terms of requirements and parameters, computational optimization capacity, high precision, and robustness to noise, among other reasons [55]. The classification values (AUC) found in this study agree with similar processes where the RF algorithm was used to determine classes associated with different phenotypes and landscape uses with information from remote sensors, among other uses [56].
Machine learning is one of the most used tools today in various fields, including agriculture [2,57], which has established itself as one of the most effective methods for detecting and predicting patterns [58]. Generalized use poses a challenge in the user community for adequate applications; an optimization and analysis process is needed for each algorithm, such as the one developed in this study.
The biological interpretation of the optimization and selection of morphological descriptors in the Rubus spp germplasm bank based on RF algorithm were variables and informative into Rubus genus, maximizing the contrast in the phenotypic discrimination [59]. These results suggest that some numerical morphometric characteristics largely define the interspecific morphological variability of the evaluated Rubus accessions, especially those related to the size of the plant organs.
The numerical descriptors with the highest level of discrimination, such as V26 (number of stingers in the stem internode) and V24 (length of the base of the stinger in the stem), are characteristics related to V23 (shape of the stinger in the stem), which makes variations in this vegetative structure a very valuable indicator that easily discriminates accessions or materials. Descriptors related to the size and arrangement of the leaf on the stem were also very important. Since species are represented by accessions and intraspecific variability, our results generally confirmed that the descriptors associated with the leaf and stem (Figure 7) tend to be the most informative [60]. Studies on Rubus subgenus Rubus highlight these descriptors in the determination of qualitative and quantitative variation among accessions [60]. These results suggest that many of the numerical descriptors prioritized by the RF method show the possibility of generating scale ratios, which is useful for comparing intervals, differences, and derivatives in absolute or dimensionless values [60]. Therefore, prioritizing these descriptors would allow the generation of composite indicators and more informative comparisons between collections and plant germplasm banks from different regions [61].
Likewise, when descriptive variables were incorporated both interspecific and intraspecific characteristics helped discriminate the phenotype of the Rubus genus materials. With this combination, the descriptors V9 (primary color of the stem surface), V23 (shape of the stinger in the stem), and V53 (primary color of the petals) contributed greatly to the definition of the morphological descriptors for this genus, highlighting the importance of the stinger shape characteristic in the stem, which has proven to be a highly discriminating characteristic at the species level and between species [21].
The 11 phenotypic variables allowed the discrimination of accessions, considered as the minimum morphological characters that would facilitate the study of Rubus germplasm. This optimization will allow characterizing the phenotypic traits of the accessions and covering a large number in a short time, reducing the time for the characterization of the entire germplasm bank [17].

4.2. Selection of the Genotype Clustering Method of the Rubus spp. Germplasm Bank

In the hierarchical agglomerative clustering and K-means, the standardization process combined with the determination of the optimal number of clusters presented a better discriminant power for Rubus spp. accessions. If the initial population and its distribution have large distances between individuals, produced by the use of non-standardized values, the number of clusters produced by the Gap statistic tends to be low and, therefore, the discrimination power is lower [49]. This affects the K-means and hierarchical agglomerative clustering methods since the Gap calculation includes the logarithm of the Sum of Square Errors in modal distribution in all terms as a quotient [62].
The NVI, ARI, CH, and PG statistics indicated that the greatest effectiveness was seen with the standardized K-means variants with reduced variables, standardized K- means with all variables, and standardized hierarchical agglomerative clustering with all the variables, because said statistics related the number of classes, their internal normalized variations and the order of the scores between the different classes formed. Therefore, they were strongly influenced by the number of clusters, the distances between individuals, and the nature of the grouping methods [63].
The K-means method performed better than hierarchical agglomerative clustering. This result was due to the reassignment nature of the K-means method, which allows each permutation to have an individual assigned to a group, independent of the group it was assigned to in the immediately previous permutation, contrary to the hierarchical methods, where individuals are assigned to a cluster depending on the initial parameters and remain in that group until the end of the analysis, creating subgroups in lower hierarchies [64].
This study indicated that the K-means method and the reduction of variables using the RF algorithm are an excellent alternative for the descriptors optimization and discrimination of accessions from the Rubus spp. germplasm bank, with high potential for use in fast, efficient characterizations in other germplasm banks. Even with the promising results presented here, these methods require internal validation, proper selection of the combination of variables, and specific clustering models for replication in different plant matrices [65].
The knowledge of the morphological variability of germplasm improves the understanding of the relationship between structural morphology and their corresponding functional botany [66]. It is considered that in the case of financial or human resource limitations, less relevant characters can be eliminated with objective elements, such as the process performed in this work. In addition, morphological descriptors must be easily determined and have a constant phenotypic expression in all environments—that is, high heritability and low environmental influence. Optimization of descriptors improves the availability of information quickly and accurately inducing efficient management of conservation and maximizing the use of financial resources [13]. Based on previous assumptions, our work constitutes an important contribution in the evaluation of the morphological variation of the Rubus germplasm with statistical, botanical, and taxonomic validity.

5. Conclusions

The correct optimization process of the RF algorithm allowed stable morphological descriptors with high taxonomic concordance to be selected, thus eliminating redundant and obsolete descriptors that present a high cost–benefit ratio. The adequate combination of discriminant morphological descriptors in combination with the optimal parameters of the K-means clustering method showed a promising approach to discriminate different materials from a population with high phenotypic variability, such as the Rubus spp. germplasm bank. This is particularly valuable since it is the first report on the use of machine learning tools and optimization of discriminant methods for the prioritization of quantitative and qualitative morphological descriptors and the ability to differentiate genotypes from plant germplasm banks for the Rubus spp. genus in tropical environments.

Author Contributions

Conceptualization, J.C.H.-R., J.G.R.-G., M.G.R.-A., C.E.V.-A.; formal analysis, J.G.R.-G., W.A.L.-R.; investigation, J.G.R.-G., J.C.H.-R., M.G.R.-A, C.E.V.-A, C.O.-M.; writing—original draft preparation, J.C.H.-R., J.G.R.-G., M.G.R.-A.; writing—review and editing, J.C.H.-R., J.G.R.-G., M.G.R.-A., C.O.-M., C.E.V.-A.; supervision, J.G.R.-G., J.C.H.-R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Corporacion Colombiana de Investigacion Agropecuaria—Agrosavia through recurring transfer from the government of Colombia.

Acknowledgments

Thanks are expressed to the Corporacion Colombiana de Investigacion Agropecuaria—Agrosavia, and the Department of Agronomia of the Faculty of Agricultural Sciences of the Universidad Nacional de Colombia Sede Bogotá for the time spent by researchers and the use of the assets inherent in the preparation of this document.

Conflicts of Interest

The authors declare no conflict of interest.

Data Availability Statement

Data associated with morphological characterization of Rubus genus germplasms are part of the country’s genetic resource; the nation’s protection laws do not allow the publication of specific results without prior authorization. In addition, the R code generated during the current study is available from the corresponding author on reasonable request. For the future as a group, we are working on the development of an R package and a jupyter notebook for Python.

References

  1. Paul, A.; Ghosh, S.; Das, A.K.; Goswami, S.; Choudhury, S.D.; Sen, S. A review on agricultural advancement based on computer vision and machine learning. In Emerging Technology in Modelling and Graphics; Springer: Berlin/Heidelberg, Germany, 2020; pp. 567–581. [Google Scholar]
  2. Liakos, K.G.; Busato, P.; Moshou, D.; Pearson, S.; Bochtis, D. Machine learning in agriculture: A review. Sensors 2018, 18, 2674. [Google Scholar] [CrossRef] [Green Version]
  3. Saleem, M.H.; Potgieter, J.; Arif, K.M. Plant disease detection and classification by deep learning. Plants 2019, 8, 468. [Google Scholar] [CrossRef] [Green Version]
  4. Chlingaryan, A.; Sukkarieh, S.; Whelan, B. Machine learning approaches for crop yield prediction and nitrogen status estimation in precision agriculture: A review. Comput. Electron. Agric. 2018, 151, 61–69. [Google Scholar] [CrossRef]
  5. Rose, D.C.; Chilvers, J. Agriculture 4.0: Broadening responsible innovation in an era of smart farming. Front. Sustain. Food Syst. 2018, 2, 87. [Google Scholar] [CrossRef] [Green Version]
  6. Ramos, P.; Prieto, F.A.; Montoya, E.; Oliveros, C.E. Automatic fruit count on coffee branches using computer vision. Comput. Electron. Agric. 2017, 137, 9–22. [Google Scholar] [CrossRef]
  7. Saleem, M.H.; Potgieter, J.; Arif, K.M. Plant Disease Classification: A Comparative Evaluation of Convolutional Neural Networks and Deep Learning Optimizers. Plants 2020, 9, 1319. [Google Scholar] [CrossRef]
  8. Kumar, M.; Gupta, S.; Gao, X.-Z.; Singh, A. Plant Species Recognition Using Morphological Features and Adaptive Boosting Methodology. IEEE Access 2019, 7, 163912–163918. [Google Scholar] [CrossRef]
  9. Sharma, R.; Kamble, S.S.; Gunasekaran, A.; Kumar, V.; Kumar, A. A systematic literature review on machine learning applications for sustainable agriculture supply chain performance. Comput. Oper. Res. 2020, 104926. [Google Scholar] [CrossRef]
  10. Jordan, M.I.; Mitchell, T.M. Machine learning: Trends, perspectives, and prospects. Science 2015, 349, 255–260. [Google Scholar] [CrossRef]
  11. Aldana-Domínguez, J.; Montes, C.; Martínez, M.; Medina, N.; Hahn, J.; Duque, M. Biodiversity and ecosystem services knowledge in the Colombian Caribbean: Progress and challenges. Trop. Conserv. Sci. 2017, 10, 1940082917714229. [Google Scholar] [CrossRef]
  12. Galeano, C.H.; Tehelen, K.; Jiménez, H.R.; Gonzalez, C.; Ceron-Souza, I. Opportunities and challenges for improving a Colombian public research program in plant breeding and plant genetic resources lead by Agrosavia. bioRxiv 2020. [Google Scholar] [CrossRef]
  13. Valencia-Ramírez, R.A.; Lobo-Arias, M.; Ligarreto, G.A. Estado del arte de los recursos genéticos vegetales en Colombia: Sistema de Bancos de Germoplasma. Cienc. Tecnol. Agropecu. 2010, 11, 85–94. [Google Scholar] [CrossRef]
  14. Arias, M.L. Recursos genéticos y mejoramiento de frutales andinos: Una visión conceptual. Cienc. Tecnol. Agropecu. 2006, 7, 40–54. [Google Scholar] [CrossRef] [Green Version]
  15. Lobo-Arias, M.; Medina, C.I. Conservación de recursos genéticos de la agrobiodiversidad como apoyo al desarrollo de sistemas de producción sostenibles. Cienc. Tecnol. Agropecu. 2009, 10, 33–42. [Google Scholar] [CrossRef] [Green Version]
  16. Gotor, E.; Alercia, A.; Rao, V.R.; Watts, J.; Caracciolo, F. The scientific information activity of Bioversity International: The descriptor lists. Genet. Resour. Crop Evol. 2008, 55, 757–772. [Google Scholar] [CrossRef]
  17. Ghosh, K. Normas Para Bancos de Germoplasma de Recursos Fitogenéticos Para la Alimentación y la Agricultura; FAO: Rome, Italy, 2013. [Google Scholar]
  18. Iza, M.; Viteri, P.; Hinojosa, M.; Martínez, A.; Sotomayor, A.; Viera, W. Morphological, phenological and pomological differentiation of commercial cultivars of blackberry (Rubus glaucus Benth.). Enfoque Ute 2020, 11, 47–57. [Google Scholar] [CrossRef] [Green Version]
  19. Mekonnen, F.; Mekbib, F.; Kumar, S.; Ahmed, S.; Sharma, T. Phenotypic variability and characteristics of lentil (Lens culinaris Medik.) germplasm of Ethiopia by multivariate analysis. J. Agric. Crop Res. 2014, 2, 104–116. [Google Scholar]
  20. Cancino-Escalante, G.O.; Sánchez-Montaño, L.R.; Quevedo-García, E.; Díaz-Carvajal, C. Caracterización fenotípica de accesiones de especies de Rubus L. de los municipios de Pamplona y Chitagá, región Nororiental de Colombia. Univ. Sci. 2011, 16, 219–233. [Google Scholar] [CrossRef]
  21. Ligarreto-Moreno, G.A.; Espinosa, N.; Barrero, L.S.; Medina, C.I. Variabilidad morfológica de variedades nativas de mora (Rubus sp.) en los Andes de Colombia. Rev. Colomb. Cienc. Hortic. 2016, 10, 211–221. [Google Scholar] [CrossRef] [Green Version]
  22. Evans, K.J.; Weber, H.E. Rubus anglocandicans (Rosaceae) is the most widespread taxon of European blackberry in Australia. Aust. Syst. Bot. 2003, 16, 527–537. [Google Scholar] [CrossRef]
  23. Voss, D.H. The Royal Horticultural Society Colour Chart 2001; Flower Council Of Holland: London, UK, 2002. [Google Scholar]
  24. Pal, M. Random forest classifier for remote sensing classification. Int. J. Remote Sens. 2005, 26, 217–222. [Google Scholar] [CrossRef]
  25. Liaw, A.; Wiener, M. randomForest: Breiman and Cutler’s Random Forests for classification and regression (R package). version 4.6. 2018. [Google Scholar]
  26. Kuhn, M.; Wing, J.; Weston, S.; Williams, A.; Keefer, C.; Engelhardt, A.; Cooper, T.; Mayer, Z.; Kenkel, B.; Team, R.C. Package ‘caret’. R J. 2020. [Google Scholar]
  27. Wright, M.N.; Wager, S.; Probst, P.; Wright, M.M.N. Package ‘ranger’. 2018. [Google Scholar]
  28. LeDell, E.; Poirier, S. H2O automl: Scalable automatic machine learning. In Proceedings of the 7th ICML Workshop on Automated Machine Learning; 2020; Volume 1, pp. 1–16. [Google Scholar]
  29. Müller, A.C.; Guido, S. Introduction to Machine Learning with Python: A Guide for Data Scientists; O’Reilly Media Inc.: Sebastopol, CA, USA, 2016. [Google Scholar]
  30. Hasmadi, M.; Pakhriazad, H.; Shahrin, M. Evaluating supervised and unsupervised techniques for land cover mapping using remote sensing data. Geogr. Malays. J. Soc. Space 2009, 5, 1–10. [Google Scholar]
  31. Paluszynska, A.; Biecek, P.; Jiang, Y. randomForestExplainer: Explaining and Visualizing Random Forests in Terms of Variable Importance. R Package version 0.10. 0. 2019. [Google Scholar]
  32. Ben-Hur, A.; Weston, J. A user’s guide to support vector machines. In Data Mining Techniques for the Life Sciences; Springer: Berlin/Heidelberg, Germany, 2010; pp. 223–239. [Google Scholar]
  33. Meyer, D.; Dimitriadou, E.; Hornik, K.; Weingessel, A.; Leisch, F.; Chang, C.; Lin, C. e1071: Misc functions of the Department of Statistics (e1071), TU Wien. R Package Version. 2014; Volume 1. [Google Scholar]
  34. Fan, R.-E.; Chen, P.-H.; Lin, C.-J. Working set selection using second order information for training support vector machines. J. Mach. Learn. Res. 2005, 6, 1889–1918. [Google Scholar]
  35. Beck, M.W. NeuralNetTools: Visualization and analysis tools for neural networks. J. Stat. Softw. 2018, 85, 1. [Google Scholar] [CrossRef]
  36. Ripley, B.; Venables, W. nnet: Feed-forward neural networks and multinomial log-linear models. R Package Version 2016, 7. [Google Scholar]
  37. Bergmeir, C.N.; Benítez Sánchez, J.M. Neural Networks in R Using the Stuttgart Neural Network Simulator: RSNNS; American Statistical Association: Alexandria, VA, USA. [CrossRef] [Green Version]
  38. Bradley, A.P. The use of the area under the ROC curve in the evaluation of machine learning algorithms. Pattern Recognit. 1997, 30, 1145–1159. [Google Scholar] [CrossRef] [Green Version]
  39. R Development Core Team. R: The R Project for Statistical Computing. R Found. Stat. Comput. 2020. [Google Scholar]
  40. Huang, J.; Ling, C.X. Using AUC and accuracy in evaluating learning algorithms. IEEE Trans. Knowl. Data Eng. 2005, 17, 299–310. [Google Scholar] [CrossRef] [Green Version]
  41. Murtagh, F.; Legendre, P. Ward’s hierarchical agglomerative clustering method: Which algorithms implement Ward’s criterion? J. Classif. 2014, 31, 274–295. [Google Scholar] [CrossRef] [Green Version]
  42. Szoecs, E.; Wagner, H. Vegan: Community ecology package. R Package Version. 2019; Volume 2, 4–5. [Google Scholar]
  43. Suzuki, R.; Shimodaira, H.; Suzuki, M.R. Package ‘pvclust’. R Top. Doc. 2015, 14. [Google Scholar]
  44. Paradis, E.; Blomberg, S.; Bolker, B.; Brown, J.; Claude, J.; Cuong, H.S.; Desper, R. Package ‘ape’. Anal. Phylogenetics Evol. Version 2015, 2, 1–4. [Google Scholar]
  45. Adler, D.; Murdoch, M.D. Available online: https://cran.r-project.org/web/packages/rgl/ (accessed on 9 December 2020).
  46. Walesiak, M.; Dudek, A. mdsOpt–Searching for Optimal MDS Procedure for Metric and Interval-Valued Data, R package version 0.3-3. 2018. Available online: https://CRAN.R-project.org/package=mdsOpt (accessed on 9 December 2020).
  47. Aggarwal, S.; Phoghat, P.; Maitrey, S. Hierarchical clustering-an efficient technique of data mining for handling voluminous data. Int. J. Comput. Appl. 2015, 129, 31–36. [Google Scholar] [CrossRef]
  48. Milligan, G.W.; Cooper, M.C. A study of standardization of variables in cluster analysis. J. Classif. 1988, 5, 181–204. [Google Scholar] [CrossRef]
  49. Tibshirani, R.; Walther, G.; Hastie, T. Estimating the number of clusters in a data set via the gap statistic. J. R. Stat. Soc. Ser. B 2001, 63, 411–423. [Google Scholar] [CrossRef]
  50. Kassambara, A.; Mundt, F. Factoextra: Extract and visualize the results of multivariate data analyses. R Package Version 1.0.7. 2020. [Google Scholar]
  51. Bolar, K. STAT: Interactive document for Working with Basic Statistical Analysis. R Package Version 0.1. 0. 2019. [Google Scholar]
  52. Vinh, N.X.; Epps, J.; Bailey, J. Information theoretic measures for clusterings comparison: Variants, properties, normalization and correction for chance. J. Mach. Learn. Res. 2010, 11, 2837–2854. [Google Scholar]
  53. Hennig, C. fpc: Flexible Procedures for Clustering. R package version 2.2-5 (2020). 2020. [Google Scholar]
  54. Chiquet, J.; Rigaill, G.; Dervieux, V. aricode: Efficient Computations of Standard Clustering Comparison Measures. 2019. Available online: https://rdrr.io/cran/aricode/ (accessed on 9 December 2020).
  55. Breiman, L. Using iterated bagging to debias regressions. Mach. Learn. 2001, 45, 261–277. [Google Scholar] [CrossRef] [Green Version]
  56. Christodoulou, M.D.; Battey, N.H.; Culham, A. Can you make morphometrics work when you know the right answer? Pick and mix approaches for apple identification. PLoS ONE 2018, 13, e0205357. [Google Scholar] [CrossRef] [Green Version]
  57. Rehman, T.U.; Mahmud, M.S.; Chang, Y.K.; Jin, J.; Shin, J. Current and future applications of statistical machine learning algorithms for agricultural machine vision systems. Comput. Electron. Agric. 2019, 156, 585–605. [Google Scholar] [CrossRef]
  58. Wang, Z.; O’Boyle, M. Machine learning in compiler optimization. Proc. IEEE 2018, 106, 1879–1901. [Google Scholar] [CrossRef] [Green Version]
  59. Tabarki, S.; Mechergui, K.; Khaldi, S.; Toumi, L.; Raies, A. Analysis of Morphological Variability in Five Spontaneous Populations of Rubus ulmifolius Schott in Tunisia. Asian J. Biol. 2017, 1–9. [Google Scholar] [CrossRef] [Green Version]
  60. Guérin, J.; Gibaru, O.; Thiery, S.; Nyiri, E. Clustering for different scales of measurement-the gap-ratio weighted k-means algorithm. arXiv 2017, arXiv:1703.07625. [Google Scholar]
  61. Mereu, R.F. A note on the ratio of the moment magnitude scale to other magnitude scales: Theory and applications. Seismol. Res. Lett. 2017, 88, 193–205. [Google Scholar] [CrossRef]
  62. Oliva, C. Métodos para la Segmentación de Datos Longitudinales. Aplicación a Datos de Rendimientos de Cultivos en Argentina. Ph.D. Thesis, Universidad de Buenos Aires, Buenos Aires, Argentina, 2015; p. 72. [Google Scholar]
  63. Reichart, R.; Rappoport, A. The NVI clustering evaluation measure. In Proceedings of the Thirteenth Conference on Computational Natural Language Learning (CoNLL-2009), Boulder, CO, USA, 4–5 June 2009; pp. 165–173. [Google Scholar]
  64. Malik, A.; Tuckfield, B. Applied Unsupervised Learning with R: Uncover Hidden Relationships and Patterns with K-Means Clustering, Hierarchical Clustering, and PCA; Packt Publishing Ltd.: Birmingham, UK, 2019. [Google Scholar]
  65. Ochieng, J.A.; Gesimba, R.M.; Oyoo, M.E.; Korir, P.C.; Owuoche, J.O.; Miheso, M. Morphological characterization of blackberry (Rubus subgenus Rubus Watson) genetic resources in Kenya. Afr. J. Plant Sci. 2019, 13, 297–308. [Google Scholar]
  66. Lauri, P.-É.; Normand, F. Are leaves only involved in flowering? Bridging the gap between structural botany and functional morphology. Tree Physiol. 2017, 37, 1137–1139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Plants 10 00247 g001 550
Figure 1. Optimization of numerical morphological descriptors of blackberry (Rubus spp.) using machine learning Random Forest algorithm. (a) importance of variables based on purity of classification nodes using Gini Index; (b) importance of variable based on minimum depth locations between of tree and number of trees; (c) classification of top and not variables based on minimum average depth, number of nodes and times to root; (d) significance (p < 0.05) based on increase in the purity of the nodes and mean minimum depth; (e) mean minimum depth for the 30 most frequent interactions variables.
Figure 1. Optimization of numerical morphological descriptors of blackberry (Rubus spp.) using machine learning Random Forest algorithm. (a) importance of variables based on purity of classification nodes using Gini Index; (b) importance of variable based on minimum depth locations between of tree and number of trees; (c) classification of top and not variables based on minimum average depth, number of nodes and times to root; (d) significance (p < 0.05) based on increase in the purity of the nodes and mean minimum depth; (e) mean minimum depth for the 30 most frequent interactions variables.
Plants 10 00247 g001
Plants 10 00247 g002 550
Figure 2. Optimization of all morphological descriptors (numerical and categorical) of blackberry (Rubus spp.) using machine learning Random Forest algorithm. (a) importance of variables based on purity of classification nodes using Gini Index; (b) importance of variable based on minimum depth locations between of tree and number of trees; (c) classification of top and not variables based on minimum average depth, number of nodes and times to root; (d) significance (p < 0.05) based on increase in the purity of the nodes and mean minimum depth; (e) mean minimum depth for the 30 most frequent interactions variables.
Figure 2. Optimization of all morphological descriptors (numerical and categorical) of blackberry (Rubus spp.) using machine learning Random Forest algorithm. (a) importance of variables based on purity of classification nodes using Gini Index; (b) importance of variable based on minimum depth locations between of tree and number of trees; (c) classification of top and not variables based on minimum average depth, number of nodes and times to root; (d) significance (p < 0.05) based on increase in the purity of the nodes and mean minimum depth; (e) mean minimum depth for the 30 most frequent interactions variables.
Plants 10 00247 g002
Plants 10 00247 g003 550
Figure 3. Optimization of morphological descriptors of blackberry (Rubus spp.) using machine learning support vector machine algorithms. (a) variable relevance based on support vector machine form linear algorithm; (b) variable relevance based on support vector machine form radial; (c) Parameter optimization associate with support vector machine form radial.
Figure 3. Optimization of morphological descriptors of blackberry (Rubus spp.) using machine learning support vector machine algorithms. (a) variable relevance based on support vector machine form linear algorithm; (b) variable relevance based on support vector machine form radial; (c) Parameter optimization associate with support vector machine form radial.
Plants 10 00247 g003
Plants 10 00247 g004 550
Figure 4. Optimization of morphological descriptors of blackberry (Rubus spp.) using machine learning Neuronal Network algorithms. (a) Variable importance for models using Olden’s connection weights algorithm.; (b) calibration validation and prediction optimizing neural network parameters through Bootstrap; (c) Correlogram of interaction of predictive variables in the neuron network; (d) Cluster centers for each variable in each group.
Figure 4. Optimization of morphological descriptors of blackberry (Rubus spp.) using machine learning Neuronal Network algorithms. (a) Variable importance for models using Olden’s connection weights algorithm.; (b) calibration validation and prediction optimizing neural network parameters through Bootstrap; (c) Correlogram of interaction of predictive variables in the neuron network; (d) Cluster centers for each variable in each group.
Plants 10 00247 g004
Plants 10 00247 g005 550
Figure 5. Cluster number selection based on GAP statistics. (a) Hierarchical agglomerative clustering (standardized); (b) Hierarchical agglomerative clustering (not standardized); (c) K-means (standardized); (d) K-means (not standardized).
Figure 5. Cluster number selection based on GAP statistics. (a) Hierarchical agglomerative clustering (standardized); (b) Hierarchical agglomerative clustering (not standardized); (c) K-means (standardized); (d) K-means (not standardized).
Plants 10 00247 g005
Plants 10 00247 g006 550
Figure 6. Selection of the genotype clustering method of the blackberry (Rubus spp.) germplasm bank. (a) Hierarchical agglomerative clustering (not standardized); (b) K-means (not standardized); (c) Hierarchical agglomerative clustering (standardized); (d) K-means (standardized) (e) Hierarchical agglomerative clustering (standardized and with reduced variables based on Random Forest process); (f) K-means (standardized and with reduced variables based on Random Forest process). HC: Hierarchical agglomerative clustering.
Figure 6. Selection of the genotype clustering method of the blackberry (Rubus spp.) germplasm bank. (a) Hierarchical agglomerative clustering (not standardized); (b) K-means (not standardized); (c) Hierarchical agglomerative clustering (standardized); (d) K-means (standardized) (e) Hierarchical agglomerative clustering (standardized and with reduced variables based on Random Forest process); (f) K-means (standardized and with reduced variables based on Random Forest process). HC: Hierarchical agglomerative clustering.
Plants 10 00247 g006
Plants 10 00247 g007 550
Figure 7. Morphological characteristics of the ten clusters formed by standardized K-means and fed with the variables prioritized by the optimized RF method.
Figure 7. Morphological characteristics of the ten clusters formed by standardized K-means and fed with the variables prioritized by the optimized RF method.
Plants 10 00247 g007
Table
Table 1. Blackberry (Rubus spp.) accessions, selected from the blackberry germplasm bank for morphological characterization.
Table 1. Blackberry (Rubus spp.) accessions, selected from the blackberry germplasm bank for morphological characterization.
SpecieUnique Identification Code (CUI)SpecieUnique Identification Code (CUI)
Rubus sp. (I)72600104Rubus glaucus (CN)72600111
Rubus sp. (I)72600101Rubus glaucus (CN)72600087
Rubus glaucus (CN)72600121Rubus glaucus (CN)72600079
Rubus glaucus (CN)72600015Rubus sp. (I)72600086
Rubus glaucus (CN)72600027Rubus glaucus (WN)72600072
Rubus glaucus (CN)72600028Rubus sp. (CN)72600090
Rubus glaucus (CN)72600034Rubus glaucus (CN)72600068
Rubus glaucus (CN)72600042Rubus sp. (I)72600115
Rubus glaucus (CN)72600044Rubus glaucus (CN)72600106
Rubus glaucus (CN)72600045Rubus glaucus (CN)72600116
Rubus sp. (WN)72600108Rubus glaucus (CN)72600080
Rubus glaucus (CN)72600050Rubus sp. (I)72600081
Rubus glaucus (CN)72600094Rubus sp. (I)72600112
Cultivated natives (CN); introduced (I); wild native (WN).
Table
Table 2. Morphological variables used in the characterization of blackberry (Rubus spp.) accessions in this study.
Table 2. Morphological variables used in the characterization of blackberry (Rubus spp.) accessions in this study.
Plant and Organs Morphological Variables
Plant growth habit(V7) Plant size (types of plant growth)
Stem(V8, V9) Primary color of the stem surface
(V10, V11) Secondary color of the stem surface
(V12) Density of stem pubescence (tector hairs)
(V13) Presence of trichomes on the stem (glandular hairs)
(V14) Color of stem trichomes
(V15) Number of trichomes on the stem [Area 0.5 cm2]
(V16, V17, V18) Color of stem pubescence
(V19) Stem waxiness
(V20) Stem shape
(V21, V22) Color of the base of the stinger on the stem
(V23) Shape of the stinger on the stem
(V24) Length of the base of the stinger on the stem [mm]
(V25) Length of the stinger on the stem [mm]
(V26) Number of stingers in the stem internode
Leaf(V27) Length of the internode on the stem [mm]
(V28) Stipule length (from leaf base and petiole) [mm]
(V29) Stipule protrusion length on petiole [mm]
(V30) Shape of stipules on petiole
(V31) Blade shape
(V32) Number of leaflets on the leaf
(V33) Margin of terminal leaflet
(V34) Shape of the base of the terminal leaflet
(V35) Shape of terminal leaflet apex
(V36) Terminal leaflet shape
(V37) Pubescence of terminal leaflet (tector hairs)
(V38, V39) Color of the terminal leaflet bundle
(V40, V41) Color of the underside of the terminal leaflet
(V42) Length of terminal leaflet [mm]
(V43) Terminal leaflet width [mm]
(V44) Petiole length at terminal leaflet [mm]
(V45) Petiole length in terminal leaflet [mm]
Flower(V46) Number of stingers in the terminal leaflet
(V47) Type of inflorescence
(V48) Petal shape
(V49) Petal length [mm]
(V50) Petal width [mm]
(V51) Color distribution on the petals (pigmentation of certain areas)
Table
Table 3. Ability of machine learning methods to select the importance and discriminant blackberry (Rubus spp.) accessions based on morphology descriptors.
Table 3. Ability of machine learning methods to select the importance and discriminant blackberry (Rubus spp.) accessions based on morphology descriptors.
Algorithm TypeTrain Accuracy 1Train Missclass Error 1Validation Accuracy 1 Validation Missclass Error 1
Random Forest with numerical variables0.6120.3870.6430.356
Random Forest with all variables0.7340.2650.7680.231
Support Vector Machine linearNaNa0.09030.906
Support Vector Machine radialNaNa0.2180.781
Neuronal NetworksNaNa0.3120.687
1 Determined parameter using area under receiver operating characteristic (ROC) curve (AUC).
Table
Table 4. Statistics test associated with evaluation of clustering methods to discriminate blackberry (Rubus spp.) accessions based on morphology descriptors.
Table 4. Statistics test associated with evaluation of clustering methods to discriminate blackberry (Rubus spp.) accessions based on morphology descriptors.
Clustering Method type NVIARIISCHENPG
HAC (not standardized)0.831890.0659411.11800147.10611.457030.45055
K-means (not standardized)0.649810.194197.72195103.67781.924560.31012
HAC (standardized)0.640300.181118.0372765.750581.965500.30829
K-means (standardized)0.323660.454147.5572424.582372.186740.12498
HAC (standardized and with reduced variables based on random Forest process)0.650260.179331.3629869.086291.949820.40440
K-means (standardized and with reduced variables based on random Forest process)0.321220.467301.0558920.307522.194650.14359
Normalized Variation Index (NVI); Adjusted Rand Index (ARI); Separation index (IS); Calinski-Harabasz Index (CH); Entropy (EN); Pearson Gamma (PG). HC: Hierarchical agglomerative clustering.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Henao-Rojas, J.C.; Rosero-Alpala, M.G.; Ortiz-Muñoz, C.; Velásquez-Arroyo, C.E.; Leon-Rueda, W.A.; Ramírez-Gil, J.G. Machine Learning Applications and Optimization of Clustering Methods Improve the Selection of Descriptors in Blackberry Germplasm Banks. Plants 2021, 10, 247. https://doi.org/10.3390/plants10020247

AMA Style

Henao-Rojas JC, Rosero-Alpala MG, Ortiz-Muñoz C, Velásquez-Arroyo CE, Leon-Rueda WA, Ramírez-Gil JG. Machine Learning Applications and Optimization of Clustering Methods Improve the Selection of Descriptors in Blackberry Germplasm Banks. Plants. 2021; 10(2):247. https://doi.org/10.3390/plants10020247

Chicago/Turabian Style

Henao-Rojas, Juan Camilo, María Gladis Rosero-Alpala, Carolina Ortiz-Muñoz, Carlos Enrique Velásquez-Arroyo, William Alfonso Leon-Rueda, and Joaquín Guillermo Ramírez-Gil. 2021. "Machine Learning Applications and Optimization of Clustering Methods Improve the Selection of Descriptors in Blackberry Germplasm Banks" Plants 10, no. 2: 247. https://doi.org/10.3390/plants10020247

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop