Tutorial - Filtering and Normalizing
Contents
Outline
- 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 filter
Background
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, thus making the data more comparable.
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.
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) and then imported into geWorkbench.
Filters
Available geWorkbench filters are as follows:
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. |
Normalizers
Normalizer | Description | |
---|---|---|
Missing value calculation | 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 |
Preparation
In these examples we will use the microarray dataset file webmatrix.exp, available in Downloads. Please refer to Tutorial - Projects and Data Files tutorial for assistance in loading a file. It contains unfiltered, unnormalized data.
Apply Quantile Normalizer
1. In the Normalization Panel, select Quantile Normalizer.
2. Leave the default averaging method of Mean Profile Marker to indicate handling of missing values..
3. Click Normalize. The View Area is updated to reflect normalization (after the screen has been refreshed). Note: The first value in the second row was update from 41,394.6 to 55,779.26.
PRENORMALIZATION | NORMALIZED | |
---|---|---|
. |
Example 2: Filtering out data called absent in an Affymetrix file
The first microarray dataset, seen in the Microarray Viewer using the Absolute display setting, looks like this:
- In the Filtering Panel, select Affy Detection Call Filter.
- Select ‘A’ (Absent) checkbox and press Filter. Values that were called Absent in the original dataset are highlighted in yellow. Internally, the values are now marked as Missing.
- In the Filtering Panel, select Missing Values Filter.
- Choose the maximum number of arrays that can have missing values before marker is removed – default is 0.
- Click Filter. Markers with more than 0 missing values are removed. Notice in the below picture that the yellow values are now gone.