M^{3}G: Maximum Margin Microarray Gridding
 Dimitris Bariamis^{1}Email author,
 Dimitris K Iakovidis^{2} and
 Dimitris Maroulis^{1}
https://doi.org/10.1186/147121051149
© Bariamis et al; licensee BioMed Central Ltd. 2010
Received: 31 July 2009
Accepted: 25 January 2010
Published: 25 January 2010
Abstract
Background
Complementary DNA (cDNA) microarrays are a well established technology for studying gene expression. A microarray image is obtained by laser scanning a hybridized cDNA microarray, which consists of thousands of spots representing chains of cDNA sequences, arranged in a twodimensional array. The separation of the spots into distinct cells is widely known as microarray image gridding.
Methods
In this paper we propose M^{3}G, a novel method for automatic gridding of cDNA microarray images based on the maximization of the margin between the rows and the columns of the spots. Initially the microarray image rotation is estimated and then a preprocessing algorithm is applied for a rough spot detection. In order to diminish the effect of artefacts, only a subset of the detected spots is selected by matching the distribution of the spot sizes to the normal distribution. Then, a set of grid lines is placed on the image in order to separate each pair of consecutive rows and columns of the selected spots. The optimal positioning of the lines is determined by maximizing the margin between these rows and columns by using a maximum margin linear classifier, effectively facilitating the localization of the spots.
Results
The experimental evaluation was based on a reference set of microarray images containing more than two million spots in total. The results show that M^{3}G outperforms state of the art methods, demonstrating robustness in the presence of noise and artefacts. More than 98% of the spots reside completely inside their respective grid cells, whereas the mean distance between the spot center and the grid cell center is 1.2 pixels.
Conclusions
The proposed method performs highly accurate gridding in the presence of noise and artefacts, while taking into account the input image rotation. Thus, it provides the potential of achieving perfect gridding for the vast majority of the spots.
Background
The process of protein synthesis inside the cells begins with the transcription of a gene sequence from DNA to messenger RNA (mRNA) in the cell nucleus. The mRNA is then transported outside the nucleus and the sequence encoded in the mRNA chain is translated into amino acids that form the corresponding protein for that particular sequence. Since proteins are translated directly from mRNA chains, the quantity of each mRNA chain that is present in a cell is indicative of the corresponding protein synthesis, i.e. the gene expression. The goal of a microarray experiment is the quantification of the amount of mRNA present in a test sample compared to that of a reference sample.
The first step of such an experiment is the isolation of the test and reference mRNA samples. These two samples are reversetranscribed into complementary DNA (cDNA), amplified using polymerase chain reaction and labelled, usually by means of two distinct fluorescent dyes such as the red Cy5 and the green Cy3. The labelled cDNA is hybridized on a microarray device that consists of a solid substrate and a large number of spots, where singlestranded chains of known DNA sequences are attached. Each of these sequences corresponds to a part of a specific gene. The sample cDNA can only be hybridized with its complementary sequence. The hybridized microarray is then scanned at the wavelength of each dye and the output of the experiment is a high resolution greyscale digital image for each wavelength. Such an image consists of a matrix of blocks, each of which contains a number of rows and columns of spots. The grey level intensity of each spot signifies the degree of hybridization of the labelled cDNA sample to the known DNA sequences, thereby indicating the expression levels of the respective genes.
The gene expression levels are extracted from microarray images in three steps. The first step of this process is the separation of the blocks present in the image. The next step is gridding, i.e. the construction of a grid covering each block so that it isolates each spot into a distinct cell, enabling the localization of each spot. The last step involves the segmentation of the spots from the background of the image and the quantification of the intensity of each spot, which corresponds to the expression level of the respective gene.
The distance between the blocks of each image is significantly larger than the distance between the spots of each block, thus the blocks can easily be separated. A variety of approaches have been proposed for block separation and have achieved accurate results. These include the analysis of the distances between neighbouring spots [1] and the use of projections of the image pixels to the x and y axes [2, 3]. In contrast to the block separation step, the process of gridding poses several challenges and has a significant impact on the accuracy of a microarray experiment [4]. A gridding algorithm should be able to grid images that include spots of various shapes, sizes and intensities, while being robust to noise and artefacts introduced at a microarray preparation stage, as well as rotation due to slight misalignments of the scanning robot coordinate system to the image coordinate system [5]. Furthermore, it is desirable that the gridding be automatically performed, without any user intervention that would possibly affect the microarray experiment, as well as limit the processing throughput of large amounts of microarray images.
Several methods have been proposed for microarray image gridding; they can be viewed in terms of automation as manual, semiautomated and fully automated [6]. However, most of the proposed methods are not fully automated and require manual tuning of parameters or other user intervention. For example, the state of the art method implemented in ImaGene [7] is semiautomated, requiring the tuning of a multitude of parameters, whereas in the manual gridding methods implemented in ScanAlyze [8] and SpotFinder [9], the process of gridding is performed interactively by the user. The method proposed by Brändle et al. [10] is parametric, requiring estimated values for several parameters.
Only a few state of the art methods have been proposed as providing automatic gridding, but most of them do not address all requirements of fully automatic gridding, i.e. handling of irregular spots and robustness against noise, artefacts and image rotation. The state of the art method proposed by Angulo et al. [11] is based on mathematical morphology and requires that grid rows and columns are strictly aligned with the x and y axes of the microarray image. The same requirement is imposed by the hillclimbing approach proposed by Rueda et al. [12]. A fully automatic region segmentation approach based on Markov random fields was proposed by Katzer et al. [13] but the results showed that its performance is diminishing in the presence of weakly expressed spots. The Bayesian grid matching method proposed by Hartelius et al. [14] employs an iterative algorithm to solve a complex deformable model for accurate microarray gridding, whereas methods producing simpler linear grids such as [13] and [15] have been proved highly accurate as well. Blekas et al. [16] proposed a method based on Gaussian mixture model, whereas later a methodology combining a stochastic search approach for the grid positioning and a Markov Chain Monte Carlo method was proposed by Antoniol et al. [17] to account for local deformations of the microarray image. However, this approach cannot be considered as fully automatic since it requires prior knowledge about the number of rows and columns of the spots in the microarray image. Another method based on Voronoi diagrams was proposed by Giannakeas et al. [3], however it requires that artificial spots are introduced in place of the spots that are very weakly expressed. Recently, a heuristic gridding approach based on a genetic algorithm was proposed by Zacharia et al. [15]. This algorithm provides a near optimal gridding outperforming the method proposed by Blekas et al. [16], while being robust to both noise and rotation. However, it is well known that the genetic optimization processes tend to require long processing times to converge, since a multitude of possible solutions has to be created and evaluated.
In this paper we propose a novel methodology for automatic cDNA microarray gridding based on a computationally efficient optimization approach. The proposed methodology is based on the maximization of the margin between the consecutive rows and columns of the microarray spots, which is implemented by training a linear maximum margin classifier with an automatically detected subset of spots on the microarray image. The classifier determines the optimal positioning of each grid line, whereas the use of the softmargin variant provides robustness to outliers. This methodology, named M^{3}G (Maximum Margin Microarray Gridding) is supported by a nonparametric Radonbased rotation estimator of general applicability for cDNA microarray images. The novel contributions of the proposed method include:

Optimal grid line determination through margin maximization

Fully automatic gridding without user intervention

Robustness to a wide range of image imperfections, including rotation, artefacts, noise and weakly expressed spots

Higher accuracy than the state of the art methods
All algorithms have been implemented under a GNU/Linux environment. The M^{3}G software is publicly available [18] and also provided along with the manuscript [Additional file 1]. A preliminary version of the proposed methodology, which involved an arbitrary threshold setting in the spot detection process and a crude rotation estimator, has been presented in [19] along with a limited experimental evaluation.
Methods
Rotation estimation
In the above equation, the grey level intensity of the image in pixel (x, y) is denoted by I(x, y). In the transformed image illustrated in Fig. 2b, the intensity of each pixel with coordinates (a, r) is equal to the integral of the image brightness over a straight line with an angle a to the x axis and a distance r to the origin. The rotation angle θ of the microarray image is estimated by locating the column with the highest mean brightness in the transformed image, which is denoted by the arrow. The image is subsequently counterrotated by angle θ as illustrated by Fig. 2c.
Preprocessing
The next step involves the preprocessing of the microarray image by linearly normalizing it so as to fit the intensity histogram to the full dynamic range of the 16bit image. The image is then quantized to 256 grey levels in order to reduce the computational complexity of the next steps. Subsequently, the edges of the spots are detected by the application of the Sobel operator [20] on the normalized image. A threshold is determined automatically using the Otsu method [21] in order to binarize the image and isolate the sharpest edges that correspond to spots, as illustrated in Fig. 2d.
Since the image has been normalized to the full dynamic range and the result of the preprocessing step is a binarized image, the quantization error is expected to affect a very small number of pixels. Indeed, in the data set used for the evaluation of the proposed method (see Results), the percentage of the affected pixels in each binarized image was less than 0.01%, having no effect on the subsequent grid placement.
Spot detection
The binarized image is then analyzed so as to locate all the groups of consecutive white pixels that reside on the same spot edge. Each of these groups is characterized by the location of the center pixel and the size of the group. Fig. 2e illustrates the detected groups by representing them as rectangles that circumscribe the pixels of each group. Ideally, each rectangle should contain the edge pixels of a single microarray spot, however depending on the threshold used, it might also include artefacts or multiple merged spots, due to the noise present in the image and the interspot proximity. A spot selection step is introduced in the next subsection, in order to refine the spots based on their shape and size characteristics resulting in the rectangles shown in Fig. 2f.
Spot selection
The spot selection process aims to the removal of false spots introduced by noise and artefacts. In the spot selection step, the aspect ratios of the detected pixel groups are first evaluated. Considering that the ideal spot shape is circular, the rectangles (Fig. 2e) should not deviate much from being square, so that each rectangle contains only one microarray spot. Therefore, the aspect ratio of each spot must be close to unity. Then, a lower bound s_{ min }and an upper bound s_{ max }of the spot sizes are calculated so as to maximize the similarity of the spot size distribution to the normal distribution. The spots that have sizes out of the calculated bounds are considered false and discarded.
Distance estimation between consecutive rows and columns
A wide range of d_{ r }values is tested in order to find the optimal value, ensuring successful estimation without any user intervention. The standard deviation of the average stripes is calculated forall values of d_{ r }within that range, using a small real valued step. From all the tested values of d_{ r }, those that result in local maxima of the standard deviation are selected. These local maxima are most often located on multiples of the optimal d_{ r }, since such an estimation also results in highly overlapping spots. For each of the selected d_{ r }values, the mean of the resulting standard deviation in its neighbourhood is calculated. The neighbourhood for the calculation of each is equal to the range between its adjacent local maxima. The value of d_{ r }that results in the highest value of the ratio is selected as optimal.
Another method for estimating the distances between rows and columns of spots has been proposed by Ceccarelli et al. [22]. It employs the Orientation Matching (OM) and Radon transforms in order to extract the spot positions and grid rotation respectively. Subsequently, the spot locations are projected on the axes of the grid and the distance between rows and columns is estimated. This method requires prior knowledge about the radii of the spots and uses a filter for noise reduction. In contrast to this approach, the proposed method performs distance estimation without any parameters or arbitrarily selected filters, based on the maximization of the standard deviation of the average stripe. This maximization over a wide range of d_{ r }values allows successful estimation without any userdefined parameters, whereas the use of the average stripe acts as a low pass filter, allowing high tolerance to noise. An evaluation of the distance estimation for noisy images has been included in the Results section.
Maximum margin gridding
In order to determine the optimal grid lines, each spot that has been selected by the spot selection step (Fig. 2e) is represented by a vector , i = 1... N, where N is the total number of selected spots in the image and the components of each vector are the coordinates of the spot centre. These vectors are assigned into distinct rows and columns, based on the distances d_{ r }and d_{ c }, which were previously estimated. Each pair of consecutive rows or columns of spots can now be separated by a single separating line. The optimal separating line is positioned so as to maximize the margin between the rows or columns of the spots. For a pair of rows numbered k and k + 1, the vectors that belong to row k or to any row above it are assigned a class label c_{ i }= +1 and the vectors that belong to row k + 1 or to any row below it are assigned a class label c_{ i }= 1. These vectors , along with their respective class labels c_{ i }are provided as a training set to a linear Support Vector Machine (SVM) classifier [23], which produces the maximum margin grid line.
In this particular application, the SVM classifier is provided with the aforementioned training set . By solving a quadratic programming optimization problem, it produces the normal vector and the parameter b of the separating line , which maximizes the margin between vectors of different classes, i.e. the margin between spots of distinct rows or columns. The width of the margin is equal to 2/w, therefore the widest margin is found by minimizing w under the constraints c_{ i }(  b) ≥ 1, i.e. requiring that all the spots lie on the correct side of the resulting grid line.
where C is a cost parameter that determines the effect of outliers on the positioning of the resulting grid line. Large values of C result in a grid line that is mostly determined by any outliers, while on the other hand, smaller values of C result in a grid line that follows the general trend of the spot locations given to the classifier, virtually ignoring any outliers. The hardmargin classifier is equivalent to a softmargin classifier with an infinitely large C[24].
The softmargin SVM is employed, in order to diminish the effects of misdetected spots that result from artefacts or noise. A small fraction of these outliers might have a shape and size similar to valid spots and could therefore pass through the selection step without being discarded. The softmargin SVM ensures that such outliers will not have an impact on the produced grid lines. For ideal microarray images, where all spots could be successfully detected and no outliers are present, a hardmargin SVM could be used as well, but a gridding application for real microarray images requires robustness against outliers. Furthermore, in ideal noiseless images, the training set for the SVM classifier would consist only of the necessary spots, i.e. those residing on rows k and k + 1. However, in real microarray images, there are cases where several consecutive spots might be very weakly expressed and therefore not detected. In order to cope with this problem, spots from rows above k and below k + 1 are included in the training set, providing redundant data to the classifier to ensure successful gridding. Using an algorithm based on the Sequential Minimal Optimization (SMO) to solve the SVM optimization problem [25], the additional data introduces only a small computational overhead, since such algorithms evaluate vectors that are far from the separating line in only the first few iterations of their outer loops [26, 27]. The SVM has been selected over similar methods for the determination of the grid lines, such as a least squares fit, because the soft margin SVM is adjustable with regards to its tolerance to outliers through the cost parameter C[28].
In the case that row k contains less than two detected spots, the two grid lines that separate this row from rows k  1 and k + 1 cannot be determined by the use of the SVM classifier. This is a rather rare case considering that the image is normalized during the preprocessing step. To cope with this limitation, the endpoints of the two grid lines are positioned equidistantly between the endpoints of the first neighboring grid lines above and below them. In the case where the top or bottom rows of spots contain less than two spots, the endpoints of the grid lines that cannot be determined are positioned d_{ r }pixels further from the nearest grid lines.
Results and Discussion
Data sets and evaluation method
The second data set consists of 10 microarray blocks selected from distinct microarray images that were artificially created or obtained from public microarray databases. The blocks are stored in TIFF files with 16bit grey level depth. This data set has been used for the evaluation of the method proposed by Blekas et al. [16] and has been obtained from the authors of [16].
In order to produce directly comparable results with the aforementioned methods, the same statistical analysis is performed. Each spot was evaluated as being perfectly, marginally or incorrectly gridded when the percentage of its pixels contained within its grid cell is 100%, more than 80%, or less than 80% respectively.
Comparative evaluation
In comparison to the preliminary version [19], several improvements allow a significant enhancement of gridding accuracy:

Automatic selection of spots based on the distribution of their sizes

Automatic determination of the operating parameters, such as the edge detection threshold and the upper and lower bounds of the sizes of the spots

Inclusion of all selected spots into the training set of each SVM classifier instead of only pairs of rows or columns of selected spots

More accurate rotation detection using the Radon transform
Euclidean distances between spot centers and grid cell centers
Mean  Standard Deviation  

Giannakeas et al. (Red)  2.73 pixels  1.10 pixels 
Giannakeas et al. (Green)  2.34 pixels  1.07 pixels 
M^{3}G (all spots)  1.17 pixels  0.71 pixels 
M^{3}G (50% of spots)  1.21 pixels  1.04 pixels 
Time performance
The proposed method was also evaluated with regards to its computational time requirements. The evaluation platform is based on an Athlon 64 X2 3800+ processor and includes 3GB of RAM, whereas the blocks used belong to the first data set and their dimensions are roughly 450 × 450 pixels. For the first block in each microarray image, our method requires 18 seconds of processing time, most of which is used for estimating the distance between consecutive rows and columns of spots. For any subsequent blocks of the same image, the range of possible values for d_{ r }and d_{ c }is reduced, resulting in 10 seconds of processing time. It is worth noting that the processing times mentioned above can be significantly reduced by optimizing the implementation and by using multiple cores of the processor, as the most time consuming parts of the gridding process can be efficiently parallelized. The genetic algorithm approach [15] requires a processing time of 92 seconds for each microarray image block, which is nearly one order of magnitude larger than the time required by M^{3}G.
Evaluation of the rotation detection
In order to assess the performance of the rotation estimation step, the images in the first data set were randomly rotated by angle θ_{ real }ranging from 25° to +25°. The proposed rotation detection method was then used to compute an estimate θ_{ est }of the rotation for each image. Based on that estimate, the images were counterrotated and subsequently gridded. The evaluation was performed based on the mean and standard deviation of the difference between the real and estimated rotation of the images Δθ = θ_{ real } θ_{ est }. Both the mean difference μ_{Δθ}and the standard deviation σ_{Δθ}were below 0.1°. The accuracy achieved when gridding the counterrotated images was within 0.3% of the accuracy obtained using the original images, thus the introduction of rotation results in only negligible variation of the gridding accuracy.
Evaluation of the distance estimation
The distance estimation step was evaluated with regards to its tolerance to noise. Additive Gaussian noise was introduced to fifty randomly selected blocks from distinct real microarray images of the first data set. The standard deviations used were σ = 250, 500 and 1000, resulting in mean signal to noise ratios of 9 dB, 5.5 dB and 1 dB respectively, whereas in several cases the introduced noise was stronger than the signal in the image. For all the images and noise levels tested, the distance estimation step displayed a negligible variance of up to 0.02 pixels, illustrating that the use of the average stripe for distance estimation provides very high tolerance to noise.
Conclusions
In this paper we presented M^{3}G, a novel method for gridding cDNA microarray images without user intervention, based on the maximization of the margin between consecutive rows and columns of spots. The proposed method involves several preprocessing steps, including a Radonbased rotation estimation for the microarray image, as well as spot detection and selection. The distance between rows and columns of spots is then estimated and the positions of the selected spots are used to train a set of linear softmargin Support Vector Machine classifiers. The use of softmargin SVMs allows high tolerance to outliers that result from artefacts and noise, whereas the use of redundant vectors in the SVM training set and the automatic determination of the operating parameters facilitate a substantial increase in gridding accuracy. Overall, the proposed method achieves successful automatic gridding of cDNA microarray images in the presence of irregular spots, noise and artefacts, as well as image rotation.
The experimental results on reference DNA microarray images containing more than two million spots showed that the proposed method outperforms the most accurate state of the art methods, providing the potential of achieving perfect gridding for the vast majority of the spots.
Declarations
Acknowledgements
The authors would like to thank K. Blekas, N. P. Galatsanos, A. Likas and I. E. Lagaris for providing one of the data sets used for the evaluation of the proposed gridding method. The authors would also like to thank the anonymous reviewers for their constructive and insightful comments.
This work was realized under the framework of the Reinforcement Program of Human Research Manpower ("PENED 2003", 03ED324), cofunded 25% by the General Secretariat for Research and Technology, Greece, and 75% by the European Social Fund.
This work was partially funded by the Special Account of Research Grants of the National and Kapodestrian University of Athens ("Kapodistrias", 70/4/6445).
Authors’ Affiliations
References
 Jung HY, Cho HG: An automatic block and spot indexing with knearest neighbors graph for microarray image analysis. Bioinformatics 2002, 18(Suppl 2):141–151.View ArticleGoogle Scholar
 Hirata R Jr, Barrera J, Hashimoto RF, Dantas DO, Esteves GH: Segmentation of Microarray Images by Mathematical Morphology. RealTime Imaging 2002, 8(6):491–505. 10.1006/rtim.2002.0291View ArticleGoogle Scholar
 Giannakeas N, Fotiadis DI: An automated method for gridding and clusteringbased segmentation of cDNA microarray images. Computerized Medical Imaging and Graphics 2009, 33: 40–49. 10.1016/j.compmedimag.2008.10.003View ArticlePubMedGoogle Scholar
 Lawrence ND, Milo M, Niranjan M, Rashbass P, Soullier S: Reducing the variability in cDNA microarray image processing by Bayesian inference. Bioinformatics 2004, 20(4):518–26. 10.1093/bioinformatics/btg438View ArticlePubMedGoogle Scholar
 Katzer M, Kummert F, Sagerer G: Methods for automatic microarray image segmentation. IEEE Transactions on Nanobioscience 2003, 2(4):202–214. 10.1109/TNB.2003.817023View ArticlePubMedGoogle Scholar
 Bajcsy P: An Overview of DNA Microarray Grid Alignment and Foreground Separation Approaches. EURASIP Journal on Applied Signal Processing 2006, 2006: 1–13.View ArticlePubMedPubMed CentralGoogle Scholar
 Biodiscovery Inc: ImaGene.2007. [http://www.biodiscovery.com/imagene.asp]Google Scholar
 Eisen MB: ScanAlyze.2002. [http://rana.lbl.gov/EisenSoftware.htm]Google Scholar
 Hegde P, Qi R, Abernathy K, Gay C, Dharap S, Gaspard R, Hughes JE, Snesrud E, Lee N, Quackenbush J: A concise guide to cDNA microarray analysis. Biotechniques 2000, 29(3):548–550.PubMedGoogle Scholar
 Brändle N, Bischof H, Lapp H: Robust DNA microarray image analysis. Machine Vision and Applications 2003, 15: 11–28. 10.1007/s001380020114xView ArticleGoogle Scholar
 Angulo J, Serra J: Automatic analysis of DNA microarray images using mathematical morphology. Bioinformatics 2003, 19(5):553–562. 10.1093/bioinformatics/btg057View ArticlePubMedGoogle Scholar
 Rueda L, Vidyadharan V: A HillClimbing Approach for Automatic Gridding of cDNA Microarray Images. IEEE/ACM Transactions on Computational Biology and Bioinformatics 2006, 3: 72–83. 10.1109/TCBB.2006.3View ArticlePubMedGoogle Scholar
 Katzer M, Kummert F, Sagerer G: A Markov Random Field Model of Microarray Gridding. In Proceedings of the 2003 ACM symposium on Applied Computing. ACM; 2003:72–77. full_textView ArticleGoogle Scholar
 Hartelius K, Carstensen JM: Bayesian Grid Matching. IEEE Transactions on Pattern Analysis and Machine Intelligence 2003, 25(2):162–173. 10.1109/TPAMI.2003.1177149View ArticleGoogle Scholar
 Zacharia E, Maroulis D: An Original Genetic Approach to the Fully Automatic Gridding of Microarray Images. IEEE Transactions on Medical Imaging 2008, 27(6):805–813. 10.1109/TMI.2008.915561View ArticlePubMedGoogle Scholar
 Blekas K, Galatsanos NP, Likas A, Lagaris IE: Mixture model analysis of DNA microarray images. IEEE Transactions on Medical Imaging 2005, 24(7):901–909. 10.1109/TMI.2005.848358View ArticlePubMedGoogle Scholar
 Antoniol G, Ceccarelli M: Microarray image gridding with stochastic search based approaches. Image and Vision Computing 2007, 25(2):155–163. 10.1016/j.imavis.2006.01.023View ArticleGoogle Scholar
 M3G: Maximum Margin Microarray Gridding[http://rtsimage.di.uoa.gr/]
 Bariamis DG, Maroulis D, Iakovidis DK: Automatic DNA microarray gridding based on Support Vector Machines. Proceedings of the 8th IEEE International Conference on Bioinformatics and Bioengineering (BIBE 2008), IEEE 2008, 1–5. full_textView ArticleGoogle Scholar
 Gonzalez RC, Woods RE: Digital Image Processing. 3rd edition. Upper Saddle River, NJ, USA: PrenticeHall, Inc; 2006.Google Scholar
 Otsu N: A threshold selection method from graylevel histograms. IEEE Transactions on Systems, Man and Cybernetics 1979, 9: 62–66. 10.1109/TSMC.1979.4310076View ArticleGoogle Scholar
 Ceccarelli M, Antoniol G: A Deformable GridMatching Approach for Microarray Images. IEEE Transactions on Image Processing 2006, 15(10):3178–3188. 10.1109/TIP.2006.877488View ArticlePubMedGoogle Scholar
 Cortes C, Vapnik V: SupportVector Networks. Machine Learning 1995, 20(3):273–297.Google Scholar
 Theodoridis S, Koutroumbas K: Pattern Recognition. 4th edition. Academic Press; 2008.Google Scholar
 Chang CC, Lin CJ: LIBSVM: a library for support vector machines.2001. [http://www.csie.ntu.edu.tw/~cjlin/libsvm]Google Scholar
 Platt J: Sequential minimal optimization: A fast algorithm for training support vector machines. In Tech. rep. Microsoft, Inc; 1998.Google Scholar
 Fan RE, Chen PH, Lin CJ: Working Set Selection Using Second Order Information for Training Support Vector Machines. Journal of Machine Learning Research 2005, 6: 1889–1918.Google Scholar
 Burges CJC: A Tutorial on Support Vector Machines for Pattern Recognition. Data Min Knowl Discov 1998, 2(2):121–167. 10.1023/A:1009715923555View ArticleGoogle Scholar
 Stanford Microarray Database2007. [http://smd.stanford.edu/]
 Juric D, Lacayo NJ, Ramsey MC, Racevskis J, Wiernik PH, Rowe JM, Goldstone AH, O'Dwyer PJ, Paietta E, Sikic BI: Differential gene expression patterns and interaction networks in BCRABLpositive and negative adult acute lymphoblastic leukemias. Journal of Clinical Oncology 2007, 25(11):1341–1349. 10.1200/JCO.2006.09.3534View ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.