TilePlot (version 1.3.1)

tileplot.single: TilePlot Single Microarray Processing Function

Description

This function is for processing a single functional gene tiling DNA microarray. As input it takes array data, and its output is the hybridization plot for each gene on the array, the mean probe intensity, the median probe intensity, the bright probe fraction, the bright segment length dependent score, the mean bright probe intensity, and the median bright probe intensity. The output is both graphical (as a set of tex and encapsulated-postscript format files that can be assembled using latex) and a tab-delimited text file with rows corresponding to genes and columns corresponding to output parameters.

Usage

tileplot.single(genesonchip, array1data, annotationslist, cutoff,
cutoff_multiplier = 3,outputfile, graphdirectory, outputtable,
array1name, smoothing_factor = 6)

Arguments

genesonchip

Path to a text file containing a list of all gene identifiers (I try to use IMG/M identifiers) on the microarray.

array1data

Path to a tab-delimited text file containing all probe identifiers in the first (left-hand) column and probe intensities in the second (right-hand) column. File must be ordered in a specific way (see below).

annotationslist

Path to a text file containing a list of all gene annotations. Each annotation must contain the gene identifier mentioned in genesonchip

cutoff

An integer cutoff value for distinguishing dark probes from bright probes based on probe intensity value. This is optional - the default value is three times the median value for all probes on the array.

cutoff_multiplier

Alternatively the cutoff can be expressed as a multiple of the median value for all probes on the array. The default is 3, but in case you wish to trade off specificity for sensitivity (or vice versa) in complex samples then you can set the multiplier lower (or higher) using this parameter.

outputfile

Path leading to desired output filename (a .tex file for input into latex).

graphdirectory

A directory for output of plots in encapsulated postscript format.

outputtable

Path leading to desired output filename (a .tdt file for further analysis of statistics for each individual gene).

array1name

Identifier for Array 1 (Default is "Array 1").

smoothing_factor

Smoothing factor for BPF recommendation. The default of 6 is usually OK.

Details

All probe names must have the following format: geneidentifier-probenumber. array1data must be sorted so that the probes for each gene are in the correct order (ie from probe 1 to the final probe). This is easily achieved with the unix sort command (sort -n -k2 -t- filename > sorted\_filename). The output file should be processed with LaTeX (http://www.latex-project.org/) to obtain a PDF.

References

http://www.stanford.edu/~ianpgm/tileplot

Examples

#This example code will deposit tex, tdt, and eps output into your R session directory
tileplot.single(genesonchip=system.file("allgenesonchip.ID", package="TilePlot"),
array1data=system.file("array1", package="TilePlot"),
annotationslist=system.file("all_annotations", package="TilePlot"),
outputfile=file.path(tempdir(), "example_output.tex"),
graphdirectory=tempdir(), outputtable=file.path(tempdir(), "example_output.tdt"))