Difference between revisions of "Tutorial - Filtering and Normalizing"

(Normalizers)
 
(20 intermediate revisions by the same user not shown)
Line 2: Line 2:
  
  
==Outline==
 
This tutorial will cover:
 
* Uses of filtering and normalization
 
* Types of filtering and normalization in geWorkbench
 
* Example 1 - Combined normalization and filtering
 
* Example 2 - using the Affymetrix Present/Absent/Marginal call filter
 
  
==Background==
+
The material on this page has been replaced.  Please see the following pages:
  
Filtering can be used to remove low quality data or reduce the size of the dataset by removing less interesting data.  Normalization can be used to decrease the effects of systematic differences across a set of microarrays, allowing better cross-microarray comparisons.
+
For Filtering: [[Filtering | Filtering]]
  
In geWorkbench, filtering and normalization alter the loaded dataset; the original is not retained. The effect of filtering out a value is to mark it internally as "Missing".  Many types of analysis require that missing values be dealt with properly. The '''Missing Values Filter''' allows markers that have more than a specified number of missing values to be removed.  Another option is to use the '''Missing Value Computation''' found in '''Normalizers''', which will replace missing values with imputed values.  Some analysis routines have built-in methods for replacing missing values.  Normalization results in the replacement of existing data values with new values.
+
For Normalization: [[Normalization | Normalization]]
 
 
Direct support for methods such as RMA and GCRMA, which include a normalization step, is not directly available in geWorkbench, but its addtion is planned (as of June 2006).  Affymetrix CEL files can be processed externally to geWorkbench using a program such as RMAExpress (available for Windows computers) or R/Bioconductor and then imported into geWorkbench.
 
 
 
==Filters==
 
 
 
geWorkbench comes with the following filters installed:
 
 
 
 
 
{|style="border: 1px solid lightGray"
 
!Filter||Description||
 
|-
 
|-|-
 
|-
 
 
 
|'''Affy Detection Call'''||Applicable to Affymetrix data only. Sets all measurements whose detection status is any user-defined combination of P, A or M (Present, Absent, Marginal) as missing.
 
|-
 
 
 
|'''Missing values''' ||Discards all markers that have “missing” measurements in at least n microarrays, where n is defined by the user. Missing values can arise either from the original data or from the results of another filtering step.
 
|-
 
|-
 
|-
 
 
 
|'''Deviation''' ||Marks as missing all markers whose deviation is less than a given value across all microarrays.
 
|-
 
|-
 
|-
 
 
 
|'''Expression Threshold''' ||Marks as missing all markers whose measurements are inside (or outside) a user-defined range.
 
|-
 
|-
 
|-
 
 
 
|'''2 Channel'''  ||Applicable to 2-channel arrays (Genepix) data only. Defines applicable ranges for each channel, and marks as missing all expression measurements for which either channel intensity is inside (or outside) the defined range.
 
|-
 
|-
 
|-
 
 
 
|'''GenePix Flags''' ||Remove values flagged in the Genepix software.
 
|-
 
|-
 
|-
 
|}
 
 
 
==Normalizers==
 
 
 
geWorkbench comes with the following normalization routines installed:
 
 
 
 
 
{|style="border: 1px solid lightGray"
 
!Normalizer||Description||
 
|-
 
|-
 
 
 
|'''Missing value computation''' ||Replaces every missing value with either the mean value of that marker across all microarrays or with the mean measurement of all markers in the microarray where the missing value is observed
 
|-
 
|-
 
|'''Log2 Transformation'''  ||Applies a log2 transformation to all measurements in a microarray
 
|-
 
|-
 
|-
 
|'''Threshold Normalizer'''  ||All data points whose value is less than (or greater than) a user-specified minimum (maximum) value are raised (reduced) to that minimum (maximum) value
 
|-
 
|-
 
|-
 
|-
 
|'''Marker-based centering''' ||Subtracts the mean (median) measurement of a marker profile from every measurement in the profile
 
|-
 
|-
 
|-
 
|'''Array-based centering''' ||Subtracts the mean (median) measurement of a microarray from every measurement in that microarray
 
|-
 
|-
 
|-
 
|'''Mean-variance normalizer''' ||For every marker profile, the mean measurement of the entire profile is subtracted from each measurement in the profile and the resulting value is divided by the standard deviation
 
|-
 
|-
 
|'''House-keeping genes normalizer''' || Normalize all values such that the averaged expression value of specified house-keeping markers is the same on each microarray.
 
|-
 
|-
 
|'''Quantile Normalizer''' || Adjusts expression values so that the distribution of values is the same on each microarray, though which marker has which value varies.
 
|-
 
|-
 
|}
 
 
 
==Preparation==
 
 
 
These examples use the microarray dataset file '''webmatrix.exp''', available in [[Download]]. Please refer to [[Tutorial - Projects and Data Files]] tutorial for assistance in loading a file.  It contains unfiltered, unnormalized data.
 
 
 
 
 
 
 
==Example 1: Normalization followed by filtering==
 
 
 
This example depicts several typical normalization and filtering steps.  These steps will recreate the dataset "webmatrix_quantile_log2_dev1_mv0.exp", available as part of the tutorials data download.
 
 
 
 
 
* In the '''Normalization''' component, select '''Quantile Normalization'''.  The image below shows a view of the dataset prior to normalization.
 
* Click '''Normalize'''.
 
 
 
 
 
[[Image:T_Normalization_Quantile.png]]
 
 
 
 
 
* Next, select '''Log2 Normalization'''.
 
* Click '''Normalize'''.
 
 
 
 
 
[[Image:T_Normalization_Log2.png]]
 
 
 
 
 
* In the '''Filtering''' component, select the '''Deviation''' filter.
 
* Here we set the lower bound on deviation at an arbitrarily chosen value of 1.  Markers whose deviation is less than this measured across all arrays will have each value marked missing.
 
* Click on '''Filter'''.
 
* Missing values are displayed in yellow in the '''Microarray Viewer''', as in the result shown below.
 
 
 
 
 
[[Image:T_Filtering_Deviation.png]]
 
 
 
 
 
 
 
* Select the '''Missing Values''' filter.  A cutoff of zero can be used because the markers have been set to missing on each array.
 
* Click on '''Filter'''.  This will remove all the markers whose values have been marked missing. 
 
 
 
[[Image:T_Filtering_MissingValuesResult.png]]
 
 
 
 
 
In this case, 3836 out of the original 12,600 markers remain.
 
 
 
 
 
 
 
==Example 2: Filtering out data called absent in an Affymetrix file==
 
 
 
* In the Filtering Panel, select'' Affy Detection Call Filter''.
 
 
 
[[Image:Filterpanel.gif]]
 
 
 
 
 
* Select the ‘A’ (Absent) checkbox and press '''Filter.''' Values that were called '''Absent''' in the original dataset are highlighted in yellow in the Microarray Viewer, as in the example above.  Internally, the values are now marked as '''Missing'''.
 
* In the '''Filtering''' component, select '''Missing Values Filter'''.
 
* Choose the maximum number of arrays that can have missing values before marker is removed – the default is 0.
 
* Click '''Filter'''.  Markers with more than 0 missing values are removed.
 

Latest revision as of 11:20, 21 June 2010

Home | Quick Start | Basics | Menu Bar | Preferences | Component Configuration Manager | Workspace | Information Panel | Local Data Files | File Formats | caArray | Array Sets | Marker Sets | Microarray Dataset Viewers | Filtering | Normalization | Tutorial Data | geWorkbench-web Tutorials

Analysis Framework | ANOVA | ARACNe | BLAST | Cellular Networks KnowledgeBase | CeRNA/Hermes Query | Classification (KNN, WV) | Color Mosaic | Consensus Clustering | Cytoscape | Cupid | DeMAND | Expression Value Distribution | Fold-Change | Gene Ontology Term Analysis | Gene Ontology Viewer | GenomeSpace | genSpace | Grid Services | GSEA | Hierarchical Clustering | IDEA | Jmol | K-Means Clustering | LINCS Query | Marker Annotations | MarkUs | Master Regulator Analysis | (MRA-FET Method) | (MRA-MARINa Method) | MatrixREDUCE | MINDy | Pattern Discovery | PCA | Promoter Analysis | Pudge | SAM | Sequence Retriever | SkyBase | SkyLine | SOM | SVM | T-Test | Viper Analysis | Volcano Plot



The material on this page has been replaced. Please see the following pages:

For Filtering: Filtering

For Normalization: Normalization