plot.single.slide: M vs. A Plot

Description Usage Arguments Details Value Author(s) See Also Examples

Description

For a single slide, this function produces a scatter plot of log intensity ratios M = log_2(R/G) versus average log intensities A = log_2(R*G)/2, where R and G represent the fluorescence intensities in the red and green channels respectively, with the addition of coloured lines indicating the cutoffs given by each of the three single slide methods.

Usage

1
plot.single.slide(x,layout,norm="p",image.id=1,...)

Arguments

x

A list with at least 4 elements. Each element of the list being a matrix with p rows for p genes and n columns for n slides. The first element 'R' contains the raw red intensities, the second element 'G' contains the raw green intensities, the third element 'Rb' contains the background red intensities and the 4th element 'Gb' contains the background green intensities. This data structure can be generated by an interactive function init.data.

layout

a list specifying the dimensions of the spot matrix and the grid matrix. This can be generated by calling init.grid.

norm

character string, one of "n", "m", "l", "p" or "s". This argument specifies the type of normalization method to be performed: "n" no normalization between the 2 channels; "m" median normalization, which sets the median of log intensity ratios to zero; "l" global lowess normalization; "p" print-tip group lowess normalization and "s" scaled print-tip group lowess normalization.

image.id

integer value; the index of the slide which is considered.

...

additional graphical parameters

Details

Plots an M vs A plot with results from each of the three single slide methods, the methods of Newton, Sapir-Churchill and Chen, indicated by colored lines. Solid black lines correspond to Newtons method at 1:1, 1:10 and 1:100 log odds, note that the points are colored red and green depending on where they fall according to this method. Dashed lines are Chens method at the 95% and 99% confidence levels. Dotted lines are Sapir-Churchil method at 95% and 99% posterior probabilites.

Value

A plot is created on the current graphics device. The plot is an M vs A plot with lines indicating cuttoff for the three single slide methods

Author(s)

Ben Bolstad, bolstad@stat.berkeley.edu

See Also

stat.Chen, stat.ChurSap, stat.Chen

Examples

1
2
3
4
data(MouseArray)
# mouse.setup <- init.grid()
# mouse.data <- init.data() ## see \emph{init.data}
plot.single.slide(mouse.data, mouse.setup, norm="p", image.id=3)

gnyamundanda/sma documentation built on May 3, 2019, 5:17 p.m.