Learn R Programming

prot2D (version 1.10.0)

RIplot: Ratio-Intensity Plot for 2D Gel Volume data

Description

R-I plots consist in plotting the intensity log2-ratio (R) against mean log10 intensity (I) of spot's volume data from 2D Gel experiments.

Usage

RIplot(data, n1, n2, ...)

Arguments

data
either a dataframe or a matrix of 2D Gel Volume data. data should be spots intensities displayed with gel as columns with the name of columns corresponding to the names of the gels and spots as rows with the names of the rows corresponding to the name of the spots. The replicates for each condition should be ordered in following columns.
n1
an integer. Number of replicates in condition 1.
n2
an integer. Number of replicates in condition 2.
...
Further arguments to be passed to methods, such as graphical parameters (see par).

Value

Display the RI-plot of data, the values on the plot are the number of spots with a Ratio greater then 1 (upper values) or lesser then -1 (lower value). A list is invinsibly returns with components :
RI
a dataframe with ratio and intensity values calculated.
diff
a vector with number of spots with a Ratio greater then 1 or lesser then -1 (in that order).

Details

Dudoit et al (2002) proposed a method for visualization of artifacts in microarray datasets, called the MA-plot, which was transposed for proteomics data as the Ratio-Intensity plot by Meunier et al (2005). R-I plots allow to directly visualize artifacts in the original data set as well as the effects of normalization. The log2-ratio (R) is the log2 of the mean of volume data in condition 2 upon the mean of volume data in condition 1: $$R = log_2\frac{mean(V_{Cond2})}{mean(V_{Cond1})}$$ The intensity (I) is the log10 of the mean of volume data in condition 2 by the mean of volume data in condition 1 $$I = log_10(mean(V_{Cond2})\times mean(V_{Cond1}))$$

References

  • Artigaud, S., Gauthier, O. & Pichereau, V. (2013) "Identifying differentially expressed proteins in two-dimensional electrophoresis experiments: inputs from transcriptomics statistical tools." Bioinformatics, vol.29 (21): 2729-2734.
  • Dudoit, S., Yang, Y.H., Callow, M.J., & Speed, T.P. (2002) "Statistical methods for identifying differentially expressed genes in replicated cDNA microarray experiments" Statistica Sinica, vol. 12: 111-139.
  • Meunier, B., Bouley, J., Piec, I., Bernard, C., Picard, B., & Hocquette, J.-F. (2005) "Data analysis methods for detection of differential protein expression in two-dimensional gel electrophoresis" Analytical Biochemistry, vol. 340 (2): 226-230.

See Also

Norm.qt,Norm.vsn

Examples

Run this code
data(pecten)

RIplot(pecten, 6, 6, main="Without Normalization")

Run the code above in your browser using DataLab