geWorkbench

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

Overview

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.

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 August 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:

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 set 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 expression values are inside (or outside) a user-defined range. For example, when "inside range" is selected, all expression values beteen the minimum and maximum values given will be marked missing.
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:

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 webmatrix2.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.2_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.

• Next, select Log2 Transformation.
• Click Normalize.

• In the Filtering component, select the Deviation filter.
• The picture below was obtained using a deviation bound of 1.0. However, to recreate the current example data set, please use a deviation bound of 1.2. 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.

• 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.

Using these settings, 2226 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.

• 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.