h-Profile plots for the discovery and exploration of patterns in gene expression data with an application to time course data
© Pittelkow and Wilson; licensee BioMed Central Ltd. 2007
Received: 11 July 2007
Accepted: 20 December 2007
Published: 20 December 2007
An ever increasing number of techniques are being used to find genes with similar profiles from microarray studies. Visualization of gene expression profiles can aid this process, potentially contributing to the identification of co-regulated genes and gene function as well as network development.
We introduce the h-Profile plot to display gene expression profiles. Thumbnail versions of plots of gene expression profiles are plotted at coordinates such that profiles of similar shape are located in the same sector, with decreasing variance towards the origin. Negatively correlated profiles can easily be identified. A new method for selecting genes with fixed periodicity, but different phase and amplitude is described and used to demonstrate the use of the plots on cell cycle data.
Visualization tools for gene expression data are important and h-profile plots provide a timely contribution to the field. They allow the simultaneous visualization of many gene expression profiles and can be used for the identification of genes with similar or reversed profiles, the foundation step in many analyses.
In gene expression (or transcript profiling) experiments, gene expression is often compared under different experimental perturbations, such as between different genotypes, cell lines, diseased and healthy subjects or under different treatment conditions. Another, complementary approach, which is useful for understanding the dynamics of altered gene expression, is via time course data, where gene expression is measured at different time intervals, usually under the same experimental conditions.
Genes with similar expression profiles, i.e. co-expressed genes are potentially co-regulated and finding these genes is often a first step in inferring gene function and in developing gene networks.
An ever increasing number of techniques are being used to find genes with similar profiles from microarray studies. Many of these methods are adapted from statistics and/or machine learning. Visualizing the gene expression would aid this process but is a challenge because of the number of genes involved and the complexity of organizing the patterns.
In this paper we describe methods which can be used for visualizing gene expression profiles. The term, gene expression profile, is used differently by different authors but essentially it refers to gene expression values (or their transformations) on all arrays for a given gene, or a summary of the gene's expression values, such as the mean, in different groups or classes of arrays. Line plots of the gene expression measurements, or summaries of gene expression in subsets of the arrays, are arguably the most common plot used for visualization gene expression data. However any plots, icons, or glyphs such as Chernoff faces, castles, trees, survival curves which capture some aspect of a gene's expression profile can be used. See for example Chambers et al.  and references therein for further ideas. We use the term 'profile plot' to describe such plots.
The approach taken in this paper is straightforward. 'Thumbnail' versions of the plot for the genes selected to be displayed are plotted at the coordinates of the h-plot [2, 3]. The word thumbnail is used in the sense of a reduced image of a graphic, utilized in order to view multiple images on a visual device simultaneously. They make it easier to visually scan and recognize broad groupings of profiles and unexpected profiles. Details are of course sacrificed because of the small size but often it is vital to obtain a broad view of the profiles, followed by more detailed plots.
There are a number of decisions which the analyst needs to make. They include the choice of plot to describe the gene profile, which genes are to be displayed, and the placement of the thumbnails in the display. Preprocessing of the gene expression measurements, such as scaling, normalization and transformation can have a profound effect on the visual effect of the profile and these decisions also need to be considered carefully.
In this paper we focus on the time series graph where lines join the gene expression measurements (transformed or otherwise) in order of time. Line graphs, where lines connect mean gene expression in different classes or time intervals are also useful.
Other plots to use as thumbnails for time series data, include the autocorrelation function (acf), periodogram, spectral plots and fitted values from a statistical model. The acf plot, periodogram and spectral plots exhibit different characteristics for data which follows a pure random process to data where there is a periodic signature. Wichert et al.  proposed a version of the periodogram as a graphical exploratory device for microarray data.
There have been a number of interesting approaches to the identification of periodically expressed genes, but for reasons of space we cannot describe them all here. An early application was described in Spellman et al. . They combined Fourier analysis with SVD. A recent example is given by Luan and Li  who use linear combinations of B-splines basis to identify genes which follow the same periodic expression pattern, subject to horizontal shift and/or different amplitudes to a set of guide genes. A number of authors have used regression models. For example Liu et al.  used quadratic regression and Venezia et al.  used orthogonal polynomial regression models and Liu et al.  use a circular regression. Some selection methods have been based on the periodogram. For example Wichert et al  use an exact test applied to periodograms, Earl et al.  used the Lomb-Scargle periodogram.
In 2005 de Lichtenberg et al.  compared a number of different methods used to identify cell cycle-regulated genes. They show that some new and more sophisticated methods for selecting cell cycle-regulated genes performed worse than the analysis published with the original data sets. Also they showed that the genes selected visually by Cho et al. on two of their benchmark data sets 'was remarkable'. They conclude that methods need to accommodate both the periodicity and the magnitude of regulation to be biologically relevant.
In this paper we introduce a new selection method which provides estimates of amplitude and phase and show that it performs well. We demonstrate the use of the h-profile plot on genes selected by our method and compare the selections with genes previously identified by Cho et al.  and de Lichtenberg et al. .
Let Z be the gene centered matrix, with N g genes in the columns and N c microarrays in the rows. Further, let the SVD of Z, be
Z = U Λ V T ,
where U, size N c by N c , and V T , size N g by N g , are orthogonal matrices such that U T U = I and V T V = I. The notation I is used to denote a conformable identity matrix.
where U k and V k are matrices comprising the first k columns of U and V respectively and Λ k is a sub matrix of Λ formed from the first k columns and rows of Λ. For two dimensional representation k = 2.
In Figure 1 gene a has the highest standard deviation and lies most distant from the origin, (0,0). Gene d which is not changing over the 3 microarrays is placed close to the origin (it has zero standard deviation).
The angular separations between the gene points in the plot are approximations to cosines of correlation between the gene expression profiles. In Figure 1, α, the angular separation between genes f and g, is quite small, and therefore the cosine of α will be high, indicating, as one might expect from the shape of their profiles, high correlation. Gene h which has a profile of opposite 'shape' to these genes i.e. is highly negatively correlated with these two genes, is placed on the other side of the plot.
Genes a, b and c are also highly positively correlated and lie in the same sector, lying approximately on a line passing through the origin. As can be seen, their profiles are of similar 'shape', but with decreasing 'steepness', i.e. decreasing variance from a to c. Gene e has the reversed profile to gene a and its profile is placed on the opposite side of the plot also lying on the line passing through genes a, b and c and the origin.
In general, genes which lie on a line passing through the origin will be correlated, positive if on the same side and negative if on the other side. The distance between two genes approximates the standard deviation of the difference between the two genes.
It should be noted that when interpreting and comparing h-plots one uses the configuration of the points (or icons) without particular reference to the coordinate axes. The axes are used only for constructing the plots, not for interpreting them as is commonly done in PCA and factor analyses.
Sometimes N c is used instead of N c - 1 in the equations for the gene and microarray coordinates. The interpretation for the angular distance between the genes in this plot is the same as in the h-plot but additional information is possible using a geometric interpretation of inner products. This biplot is one of the biplots introduced by Gabriel  with variants called the covariance biplot  and the GE-Biplot .
It is easy to show that Z k = C k , where Z k is a rank k approximation to Z. Thus in biplots, where the microarrays, C k , and the genes, G k , are plotted on the same graph, the scalar product between the i th row point and j th column point with respect to the origin is approximately equal to the (i, j) th element, z ij , of Z, the transformed, mean corrected gene expression matrix used in the biplot.
In the particular biplot used here, called the GE-biplot (Gene Expression Biplot) to differentiate it from the original covariance biplot, the juxtapositions of the gene points to the microarray points provides an approximation of the value of the (transformed) gene expression values on the microarrays. Geometrically, the inner product can be viewed as the product of the signed length of one of the vectors and the length of the projection of the other vector onto it. Thus if a gene point is close to a microarray point then it can be deduced that the gene is relatively up regulated in the microarray represented by the point and if the gene point is on the opposite side of the plot to the microarray point then the gene is relatively down regulated on the microarray. How accurate these deductions are depends on how good the approximation is in the lower ranked space. Two measures of fit, I1 and I2, ranging from 0 to 1, are described in the Methods section.
With the biplot, using the factorization of the gene centered matrix described here, the distance between two microarray points approximates the Mahalanobis distance standardized in the plane of fit  between the microarrays. Mahalanobis distance, also called the standardized or generalized distance, is arguably more appropriate than Euclidean distance for microarray data, as it takes correlation into account.
Gene Selection for cyclic data
The profiles, and theoretical arguments, suggest that finding genes with a log profile described by
z = |Asin(2πBt + L)|,
is of interest. In this equation z is the gene expression (possibly transformed), A is the amplitude, B corresponds to the period and L to the part of the cycle at time zero i.e. related to phase. This model allows the selection of genes with different phases and/or amplitudes.
If B = 1 and time t is scaled to run from 0 to 1, there are exactly two cycles and the aim is to choose genes whose profiles complete two cycles and have approximately equal amplitude in both cycles. The model
z = |Asin(2πt + L)|,
To assess fit, using RSS, knowledge of the distribution of RSS under the null model is required. To obtain this we simulated 1000 data points, y, where y t = |Asin(2πt + L)| + ε and without loss of generality, A = 1 and L = 0 with ε ~ N(0, ). The average within phase variance across all genes was used as an estimate of σ ε in the simulation.
Different quantiles of the distribution of RSS for the simulated data, denoted by , can be used to select genes which fit the model. Estimates of RSS g for all genes, g, were calculated and compared to . All genes, where were selected as being compatible with the two cycle model. Different values of select genes with more or less compatability with the model.
Estimates of A from these fits provide information on the extent of gene expression change (amplitude), and the value of t when = A is an estimate of the time of maximal expression. The intercept parameter L determines the expression at the beginning of the cycle. Displays of this function, with A = 1 and B = 1 for values of 0 ≤ L ≤ 3 are shown in additional file 1.
Mitotic Cell Cycle Data
The time-course data set described in Cho et al.  and available from the website  consisted of 6601 gene expression measurements from 17 Affymetrix pairs of chips. The aim of the study  was to characterize mRNA transcript levels during the cell cycle of the budding yeast S. cerevisiae. Synchronous yeast cultures were arrested in late G1 and the cell cycle re-initiated with cells collected at 10 minute intervals, covering nearly two full cell cycles. The time course was divided into early G1, late G1, S, G2 and M phases based on the size of the buds, the cellular position of the nucleus, and standardization to known transcripts.
In the initial phases of their analysis, Cho et al. selected 1348 out of 6565 genes which had a 2 fold difference more than 40 minutes past the point of release from arrest. They visually inspected the profiles for these genes for evidence of periodicity of expression and chose 416 genes with identifiable periodic changes in normalized (but not, apparently, logged) gene expression.
The normalization method used by Cho et al. was not described in their paper. However the distributions of gene expression measurements from each microarray were similar, although positively skewed, as is common with untransformed gene expression values. There were however a large number of negative measurements, especially for samples taken at 90 and 100 minutes. In our analyses the negative values were truncated to .01, the data logged using base 2, and finally standardized so that the microarray means are zero and their variance one. This later transformation effectively 'normalizes' the distributions so that the first two moments agree. Control genes were removed leaving a total of 6565 genes. Although not technically correct we will refer to the transformed gene expression as simply gene expression to avoid cumbersome expressions. But it should be kept in mind that the preprocessing can have a large impact on analyses. The sample at time zero, that is immediately after arrest, was eliminated from the following analyses, leaving two microarrays in each phase.
The difference in gene expression between the two cycles in the early phases has a number of potential causes such as problems with the arrest method, loss of synchronization  or temperature induced effects . The profiles shown in Figure 4 also support this interpretation for the selected genes.
We now focus on finding genes with similar profiles to CLB1, the other landmark gene, for the M phase, and used by Cho et al. Here, we used the median of the simulated values of fit to define the criteria for selection. 157 genes satisfied this criteria and of these, 79 were also in the CHO list and 57 in the L300 list.
The profiles of 9 genes in a sector around CLB1 are shown in the middle and the right of the figure. 6. All of these genes were in the CHO list and 4 in the L300 list. They all peak between 60 and 70 minutes and have a correlation coefficient of .7 or higher with CLB1. Interestingly, CLB1 is a landmark gene for the M phase, but in this data it has peak expression on the arrays taken at 60 minutes and 140 minutes past arrest, which was identified by Cho et al to be in the G2 phase. This is not a product of the transformation and scaling of the chips used in this analysis, as this profile is consistent with Figure 4 of Cho et al. . In practice one would want to check the profiles of all the landmark genes for phase M and see if this were a consistent pattern for these data. The ability to visualize the profiles allows the researcher to question meaningfully the criteria used for searching and selection.
In almost all microarray applications, the analysis starts with performing some transformation on a data matrix, depending on the nature of the data. Examples of transformations are centring with respect to the overall mean, centring with respect to variable (gene) or sample (microarray) means, normalization of variables, scaling by row or column means, and finally the square root and logarithmic transforms. The effect of these transformations radically affects what is plotted and the utility of the resultant plot for a given application.
Many authors standardize or 'norm' the gene expression values for each gene. If this transformation is used for the type of plots described here then the gene points will tend to lie on the circumference of the circle. This is because the variances, which are approximated by the square of the distance of the gene point to the origin, have been set to one, for all genes. Gene points which lie off the circumference have poor fit in lower dimensions.
The standard h-plot centers the gene profile values so their means are zero and scales the singular vectors by the square root of one less than the number of samples. This ensures that the variances and covariances between the genes is reflected in the geometry of the plot. The use of the log transformation is usual with positively skewed microarray gene expression data. If skewed data is used then plots (and analyses) will tend to be dominated by the typical presence of a few very large or very small gene expression values. The log transform has the effect of stretching the lower values close to zero, and sometimes an 'offset' is use to translate the measurements away from zero.
Spectral and singular value decompositions form the basis of many data analysis methods developed for microarray analysis [18–20]. Principal Components analysis (PCA) is well known, but its usefulness for displaying gene expression data depends on using the appropriately scaled and organized gene expression matrix . Other Biplots which use SVD in their algorithms include correspondence analysis (CA)  which was developed for frequency tables and Aitchinson's compositional biplot  that is most suited to data where the 'variables' or vectors involve an implicit constraint – such as sum. They differ in the scaling applied to the matrix on which the SVD is performed.
Both plots described in this paper have close connections to other methods based on spectral decompositions and singular value decomposition algorithms when used with suitable transformed matrices. However it is only with the factorizations and transformation described in the Methods section that the interpretations in terms of variances and correlations are possible. To maintain these interpretations it is also important to keep the scales on each axis of the display the same, which is not the default in some biplot software.
In this paper we have concentrated on the use of time series plots, because the h-profile plot and covariance biplot was demonstrated on cell cycle data. However other thumbnails in other settings are possible including for example, residuals from statistical and time series models, fitted values from models, smoothed representations, survival curves etc. Another, very important, use of the methodology presented in this paper is the association of annotation information with similar gene expression profiles. Annotation can be added to the profiles, such as chromosome number and/or location, gene ontology information etc. Alternatively suitable thumbnails summarizing annotation information can be displayed at the coordinates described.
Gene selection is important for producing meaningful plots or analyses. The gene selection method demonstrated in this paper has the advantage of providing an estimate of amplitude and its standard error. In this paper we were more interested in selecting genes whose profiles matched a 2 cycle periodic function of any amplitude and to use the plots to explore the profiles further. In a comparison of computational methods for the identification of cell cycle-regulated genes, de Lichtenberg et al.  concluded that methods that did not take into account the amplitude of the response did not perform well. The placement of the thumbnails in the h-profile plot also provides a visual estimate of the standard deviation of the genes, which was used by de Lichtenberg et al.  to define regulation. They argued that heavily regulated genes would have large standard deviations, whereas genes without significant regulation would display little deviation from the mean.
If the fit value is not high then the variances and covariances (correlation) will not be well represented and in practice it is often useful to modify the selection criteria, using the visualisation tools demonstrated in this paper in an iterative process.
We have demonstrated the plots here on data generated in order to study cell cycle regulated genes, but the utility of the plots extends beyond these types of applications. In other microarray studies, correlation between gene expression profiles has often been used as a criteria for selecting co-expressed genes. Pearson's product moment correlation coefficient is arguably the most popular measure of similarity between genes used in microarray data analysis. However clustering methods using correlation are often unable to detect subsets of genes which have the same profiles in a subset of the microarrays, but are uncorrelated in other subsets. This is because genes which show highly correlated patterns of expression in one biological state, but not in another (i.e. different subsets of microarrays), may not be highly correlated across the entire data set. The GE-biplot and the h-profile plots facilitate the finding of such subsets.
We have not addressed the question of simultaneous hypothesis testing. Many of the approaches developed to address this challenge are not suited to highly correlated data, and here we are interested in finding correlated gene expression patterns.
In this paper we demonstrate the usefulness of h-profile plots to simultaneously visualize many gene profiles. By using the coordinates of the h-plot [2, 14] to position genes in k-dimensional space, genes with similar gene expression profiles are placed in the same sector, those with high and low variance profiles being separated according to the distance from the origin. Further, positively correlated profiles are placed on the opposite sides of the plot to negatively correlated profiles. Using thumbnail versions of profile plots, the profiles can be readily identified.
Visualization methods are an important component of data analysis, both supporting and being supported by bioinformatic and statistical analyses. The methods described here are easy to implement and could become part of routine analyses (R source code for the GE-Biplot and Gene-plot is available at http://dayhoff.anu.edu.au/software). h-Profile plots can be used for detecting novel or unexpected groups of genes with similar profiles, or for displaying profiles following selection of genes, or following application of a machine learning algorithm or a statistical method. They facilitate the finding of genes which have similar or reversed profiles.
The h-plot and covariance Biplot
Gabriel  introduced the biplot as a method for displaying the elements of a matrix as inner products of vectors corresponding to the rows and columns of a matrix. Gower and Hand  have extended and generalized the ideas of Gabriel. Two earlier descriptions of the biplot for microarray data are Chapman  and Pittelkow and Wilson . There is some confusion about the definition of biplots  and a variety of biplots have been proposed, some more suited than others to the analysis of microarray data.
Any matrix, Z, of rank r and size N c by N g can be factored as
Z = CG T
where C is a matrix of size N c by r, and G is N g by r. Any element of Z can therefore be written as z ij = , where c i is row of C and g j is a row of G.
The biplot, as described by Gabriel , is the plot of all the N c + N g vectors, c i , i = 1, ..., N c and g j , j = 1, ..., N g . Originally, both the rows and the columns were represented on the biplot by vectors (lines from the origin or rays), then it became common practice to use the vector representation only for the columns (variables) but, with microarray data, the vector representation is used for the rows (microarrays).
To be practical as a display, r would need to be two or three. The matrix of gene expression indices is not, in general, of sufficiently low rank to plot usefully. The usual approach, in this case, is to find a rank k approximation to Z where k is usually two or three. An exact biplot representation of the rank-k approximation matrix is known as an 'approximate biplot', but for ease of notation we drop the qualifier 'approximate'.
It is known that one may use the singular value decomposition (SVD) to approximate any rectangular matrix by a matrix of the same size but of lower rank, such that the sums of squares of the differences between the elements of the matrix and its approximation is minimized. Let the SVD of a rank-r matrix, Z, be Z = U Λ V T where U, size N c by N c , and V T , size N g by N g , are orthogonal matrices such that U T U = I and V T V = I. The notation I is used to denote a conformable identity matrix.
where U k and V k are matrices comprising the first k columns of U and V respectively and Λ k is a sub matrix of Λ formed from the first k columns and rows of Λ.
where 0 ≤ α ≤ 1, can be used to factorize . If C k and G k are defined as C k = U k , and G k = V k , then .
To ease notation in the following, assume that k is given, say 2, and write and = G k .
Different scalings affect the position of the points in the plane. If the matrix Z is gene centered, i.e. if for each gene, the average gene expression measurement is subtracted from its individual measurements, and in the biplot factorization with α = 0 and W = I or W = I, the configuration of the gene points is determined from the variance-covariance matrix of the genes. (Or, more correctly, an approximation to the variance-covariance matrix of the genes.)
Given W = I (or W = I), , and , .
The covariance Biplot uses the coordinates, (microarray points) and (gene points). The h-plot uses the coordinates .
where is the cosine of the angle subtended at the origin between the vectors, and and ||·|| denotes length. Thus if the point represented by is close to the point then it can be deduced that the gene represented by the point is relatively up regulated in the microarray represented by the point . If is on the opposite side of the plot to then the gene is relatively down regulated on the microarray.
How accurate these reductions are depends on how good the approximation is in the lower ranked space. Two measures of goodness-of-fit for a k-dimensional display, which we call I1 and I2, are available. They range from zero to one, with a high value assuring close approximation. is an absolute goodness of fit statistic, and is a goodness of fit statistic for the variances and covariances between the genes.
An alternative geometric interpretation of the inner product is as the product of the signed length of one of the vectors and the length of the projection of the other vector onto it. Since , row points can be projected onto a line joining two column points to obtain a ranking of the rows in terms of column differences. Further since one can project column points onto a row difference vector to estimate the difference in elements between two rows. Since is the interaction between gene j and j' on chips i and i', row difference vectors which lie at right angles to column difference vectors are indicative of no interaction. Further interpretations can be found in [3, 13] and .
Since log(z ij /z i'j ) ≃ () g j , gene points can be projected onto a chip difference vector to estimate the log ratio of gene expression between two microarrays. Projection of gene points onto a vector joining two microarrays points provides a visual ranking of the genes in terms of their log ratios on the two microarrays.
There is an indeterminacy in this model which is approached in different ways. Wentzell et al , for example, use multiple curve resolution methods from chemistry to resolve the indeterminacy in the solution space. Note that Wentzell et al's profiles are another example of profiles to those described here. In the h-plot and covariance-plot, the indeterminacy is resolved so that the inner-products of the gene vectors () approximates the variance-covariance matrix of the genes.
We thank referees for insightful comments on an earlier draft of the paper.
- Chambers JM, Cleveland WS, Kleiner B, Tukey PA: Graphical Methods for Data Analysis. 1983, statistics/probability, WadsworthGoogle Scholar
- Corsten L, Gabriel K: Graphical exploration in comparing variance matrices. Biometrics. 1976, 32 (4): 851-863. 10.2307/2529269.View ArticleGoogle Scholar
- Gabriel K: The biplot graphical display of matrices with applications to principal component analysis. Biometrika. 1971, 8: 453-467. 10.1093/biomet/58.3.453.View ArticleGoogle Scholar
- Wichert S, Fokianos K, Strimmer K: Identifying periodically expressed transcripts in microarray time series data. Bioinformatics. 2004, 20: 5-20. 10.1093/bioinformatics/btg364.View ArticlePubMedGoogle Scholar
- Spellman PT, Sherlock G, Zhang MQ, Iyer VR, Anders K, Eisen MB, Brown PO, Botstein D, Futcher B: Comprehensive identification of cell cycle-regulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization. Molecular Biology of the Cell. 1998, 9: 3273-3297.PubMed CentralView ArticlePubMedGoogle Scholar
- Luan Y, Li H: Model-based methods for identifying periodically expressed genes based on time course microarray gene expression data. Bioinformatics. 2004, 20 (3): 332-339. 10.1093/bioinformatics/btg413.View ArticlePubMedGoogle Scholar
- Liu H, Tarima S, Borders AS, Getchell TV, Getchell ML, Stromberg1 AJ: Quadratic regression analysis for gene discovery and pattern recognition for non-cyclic short time-course microarray experiments. BMC Bioinformatics. 2005, 6: 106-10.1186/1471-2105-6-106.PubMed CentralView ArticlePubMedGoogle Scholar
- Venezia1 TA, Merchant AA, Ramos CA, Whitehouse NL, Young AS, Shaw CA, Goodell MA: Molecular signatures of proliferation and quiescence in hematopoietic stem cells. PLoS Biol. 2004, 2 (10): e301-. 10.1371/journal.pbio.0020301.View ArticleGoogle Scholar
- Liu D, Weinberg CR, Peddada SD: A geometric approach to determine association and coherence of the activation times of cell-cycling genes under differing experimental conditions. Boininformatics. 2004, 20 (16): 2521-2528. 10.1093/bioinformatics/bth274.View ArticleGoogle Scholar
- Glynn EF, Chen J, Mushegian AR: Detecting periodic patterns in unevenly spaced gene expression time series using Lomb-Scargle periodograms. Bioinformatics. 2006, 22 (3): 310-316. 10.1093/bioinformatics/bti789.View ArticlePubMedGoogle Scholar
- de Lichtenberg U, Jensen LJ, Fausboll A, Jensen TS, Bork P, Brunak S: Comparison of computational methods for the identification of cell cycle-regulated genes. Bioinformatics. 2005, 21 (7): 1164-1171. 10.1093/bioinformatics/bti093.View ArticlePubMedGoogle Scholar
- Cho RJ, Campbell MJ, Winzeler EA, Steinmetz L, Conway A, Wodicka L, Wolfsberg TG, Gabrielian AE, Landsman D, Lockhart DJ, Davis RW: A Genome-Wide Transcriptional Analysis of the Mitotic Cell Cycle including control of mRNA transcription. Molecular Cell. 1998, 2: 65-73. 10.1016/S1097-2765(00)80114-8.View ArticlePubMedGoogle Scholar
- Aitchison J, Greenacre M: Biplots of Compositional Data. Applied Statistics. 2002, 4 (51): 375-392.Google Scholar
- Pittelkow YE, Wilson SR: Visualisation of gene expression data - the GE-biplot, the Chip-plot and the Gene-plot. Stat Appl Genet Mol Biol. 2003, 2: Article6-PubMedGoogle Scholar
- Mitotic Cell Cycle Data. [http://genomics.stanford.edu]
- Bar-Joseph Z, Farkash S, Gifford DK, Simon I, Rosenfeld R: Deconvolving cell cycle expression data with complementary information. Bioinformatics. 2004, 20 (suppl 1): i23-30. 10.1093/bioinformatics/bth915.View ArticlePubMedGoogle Scholar
- Johansson D, Lindgren P, Berglund A: A multivariate approach applied to microarray data for identification of genes with cell cycle-coupled transcription. Bioinformatics. 2003, 19: 467-473. 10.1093/bioinformatics/btg017.View ArticlePubMedGoogle Scholar
- Fellenberg K, Hauser N, Brors B, Neutzner A, Hoheisel J, Vingron M: Correspondence analysis applied to microarray data. PNAS. 2001, 98: 10781-10786. 10.1073/pnas.181597298.PubMed CentralView ArticlePubMedGoogle Scholar
- Chapman S, Schenk P, Kazan K, Manners J: Using biplots to interpret gene expression patterns in plants. Bioinformatics. 2001, 18 (1): 202-204. 10.1093/bioinformatics/18.1.202.View ArticleGoogle Scholar
- Alter O, Brown P, Botstein D: Singular value decomposition for genome-wide expression data processing and modeling. Proceedings of the National Academy of Sciences USA. 2000, 97 (18): 10101-6. 10.1073/pnas.97.18.10101.View ArticleGoogle Scholar
- Pittelkow Y, Wilson S: Use of Principal Component Analysis and of the GE-biplot for the Graphical Exploration of Gene Expression Data. Biometrics. 2005, 61 (2): 630-632. 10.1111/j.1541-0420.2005.00366.x.View ArticlePubMedGoogle Scholar
- Greenacre J: Theory and Application of Correspondence Analysis. 1984, Academic Press, New YorkGoogle Scholar
- Gower J, Hand D: Biplots. 1996, Chapman and Hall, LondonGoogle Scholar
- Wentzell PD, Karakachm TK, Roy S, Martinez MJ, Allen C, Werner-Washburne M: Multivariate curve resolution of time course microarray data. BMC Bioinformatics. 2006, 7 (343):Google Scholar
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.