Partek Flow Documentation

Page tree

Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

Principal components (PC) analysis (PCA) is an exploratory technique that is used to describe the structure of high dimensional data by reducing its dimensionality. It is a linear transformation that converts n original variables (typically: genes or transcripts) into n new variables, which are called PCs, they have three important properties:

  • PCs are ordered by the amount of variance explained
  • PCs are uncorrelated
  • PCs explain all variation in the data

PCA is a principal axis rotation of the original variables that preserves the variation in the data. Therefore, the total variance of the original variables is equal to the total variance of the PCs.

If read quantification (i.e. mapping to a transcript model) was performed by Partek EPartek E/M algorithm, PCA can be invoked on a quantification output data node (Gene counts or Transcript counts) or, after normalisationnormalization, on a Normalized counts data node. Select a node on the canvas and then PCA in the Visualization Exploratory analysis section of the toolbox. context sensitive menu.

There are two options for  features contribute (Figure 1):

equally: all the features are standardized to mean of 0 and standard deviation of 1 .  This option will give all the features equal weight in the analysis, this is the default option for e.g bulk RNA-seq data.

by variance: the analysis will give more emphasis to the features with higher variances. This is the default option for e.g. single cell RNA-seq data

If the input data node is in linear scale, you can perform log transformation on PCA calculation. 


Numbered figure captions
SubtitleTextPCA setup dialog
AnchorNamepac_config
Image Added


The PCA task creates a new task node, and to open it and see the result, do one of the following: select the PCA task node, proceed to the toolbox context sensitive menu and go to the Task result; or double-click on the PCA task node.

Numbered figure captions
SubtitleTextRNA-Seq workflow with principal components analysis (PCA) task invoked on the Gene counts data node
AnchorNamernaseq-pca

Image Removed

On the other hand, the Cufflinks quantification node is not associated with a detailed report like the Partek E/M nodes, so opening the node by double clicking (or going via the toolbox) brings up the PCA plot directly.

Also note that PCA is available for ERCC assessment; open the QA/QC on ERCC controls task node and push View PCA in the lower left corner.

PCA in Partek Flow is using correlation matrix to find the components, which means all the  features are standardized to mean of 0 and standard deviation of 1 (the standardization is performed during the computation and does not modify the values in the data node). The results are presented as scatterplot  (Figure 2), with each dot on the plot being a sample, while the axes represent the PCs and the axes values correspond to the respective PC values. By default, The report containing eigenvalues, PC projections, component loadings, and mapping error information for the first three PCs. 

When the PCA node is opened in Data viewer, by default, it is the 3D scatterplot, Scree plot with Eigenvalues, and Component loadings table  (Figure 2). Each dot on the 3D scatter plot represents an observation, while the first three PCs are shown on the X-, Y-, and Z-axis respectively, with the information content of an individual PC is in the parenthesis. 

As an exploratory tool, the the PCA scatterplot is applied to view any groupings in the data set and generate hypotheses based on the outcome, or to spot possible outliers. 

Numbered figure captions
SubtitleText

...

Double-click on the PCA result node to open the Principal components analysis result
AnchorNamepca

...

Image Added

 To rotate the 3D scatter plot left click & drag. To zoom in or out, use the mouse wheel. Click and drag the legend can move the legend to different location on the viewer.

The Detailed configuration on PCA plot can be customised by using the controls on the left (Figure 3). Color by shows the sample attributes as listed in the Data tab or you can set it to Fixed to have all the dots of the same color. Size by option works in the same way, but affects the dot sizes.found by clicking Help>How-to videos>Data viewer section.

In the Data viewer, when a PCA data node is selected from Get Data under Setup (left panel), the node can be dragged and dropped to the screen (Figure 3), then you will have the option to select a scree plot and tables. 

Numbered figure captions
SubtitleText

...

Drag PCA data node to plot
AnchorName

...

dragPCA

...

Image Removed

...

Connect by option is particularly useful for dependent study designs, where you can highlight the samples based on the same biological source by the connecting lines. Example on Figure 4 depicts results of a study where each RNA sample was processed by both RNA-seq and  gene expression array; the lines connect the same samples.

Numbered figure captions
SubtitleTextPrincipal components analysis plot in 3D. Each dot is a sample and samples originating from the same biological source (dependent study design) are connected by lines. The axes show the first three principal components, with the fraction of explained variance in the parenthesis
AnchorNamepca-3d

Image Removed

...

Image Added

When choose Scree plot icon Image Added, it will plot a 2D viewer, X-axis represents PCs, Y-axis represents eigenvalues (Figure 4)

...

Image Removed

...

Next, Show legend turns the legend on or off. Select all selects all the dots, while Show axis turns the coordinate axis on or off. To change the size of the dots, use the Dot size slider. Grid cells increases or decreases the size of the cells in the plot grid and the Major grid interval specifies the frequency of major grid lines (fat lines). E.g. setting the Major grid interval to 4 highlights every 4th grid line.

You can also reduce the number of dimensions of the plot by switching the Plot style to 2D (Figure 6 is based on the same data as Figure 2, but is plotted in 2D).

 

...

Numbered figure captions
SubtitleText

...

Scree plot
AnchorNamescree

Image Added

When mouse over on a point on the line, it will display detailed information of the PC. The scree plot shows how much variation each PC represents,  so it is often used to determine the number of principal components to keep for downstream analysis (e.g. tSNE, UMAP, graph-base clustering). The "elbow" point of the graph where the eigenvalues seem to level off should be considered as a cutoff point for downstream analysis.

PCA data node can also be draw as tables, when choose Table icon (Image Added), it will display the component loadings matrix in the viewer (Figure 5). The Content can be modified using the Content configuration option; the table can be paged through here or from the lower right corner. 

...

Image Removed

...

Although first three PCs are shown by default, you can plot any of the first nine PCs, by using the X, Y, and Z drop-down lists.

Once you are pleased with the appearance of the dot plot, push Save image button to save it to the local machine. The resulting dialog (Figure 7) controls the resolution of the image file. The image will be saved in png format, and the default filename is PCA plot.png. Or, if you are not happy with your edits, you can always revert to the initial view by pushing the Reset button.

Numbered figure captions
SubtitleText

...

Component loadings are the correlation coefficients between the features and PCs.
AnchorNamecomponet
Image Added

In the table, each row is a feature, the column represent PCs, the value is the correlation coefficient. Under Content, there is a PCA projections option,  change to this option to display the projection table (Figure 6). In this table, each row is an observation, each column is a PC, the values are the PC scores. 

Numbered figure captions
SubtitleText

...

PCA project table
AnchorName

...

projection

...

Image Added




Additional assistance


 

Rate Macro
allowUsersfalse

...