prot2D.R
In prot2D: Statistical Tools for volume data from 2D Gel Electrophoresis

### R code from vignette source 'prot2D.Rnw'

###################################################
### code chunk number 1: prot2D.Rnw:67-68
###################################################
library(prot2D)


###################################################
### code chunk number 2: prot2D.Rnw:71-73
###################################################
data(pecten)
data(pecten.fac)


###################################################
### code chunk number 3: prot2D.Rnw:90-91
###################################################
RIplot(pecten, n1=6, n2=6)


###################################################
### code chunk number 4: prot2D.Rnw:97-98
###################################################
pecten.norm <- Norm.qt(pecten, n1=6, n2=6, plot=TRUE)


###################################################
### code chunk number 5: prot2D.Rnw:106-107
###################################################



###################################################
### code chunk number 6: prot2D.Rnw:112-113
###################################################
ES.p <- ES.prot(pecten.norm, n1=6, n2=6, f=pecten.fac)


###################################################
### code chunk number 7: prot2D.Rnw:137-139
###################################################
ES.diff <- modT.Prot(ES.p, plot=TRUE, fdr.thr=0.1, 
                     method.fdr="BH" )


###################################################
### code chunk number 8: prot2D.Rnw:142-143
###################################################
featureNames(ES.diff)


###################################################
### code chunk number 9: prot2D.Rnw:146-147
###################################################
head(fData(ES.diff))


###################################################
### code chunk number 10: prot2D.Rnw:150-151
###################################################
head(exprs(ES.diff))


###################################################
### code chunk number 11: prot2D.Rnw:166-168
###################################################
Sim.data <- Sim.Prot.2D(data=pecten, nsp=700,
                        nr=10, p0=0.1, s2_0=0.2, d0=3)


###################################################
### code chunk number 12: prot2D.Rnw:171-173
###################################################
res.stud <- ttest.Prot(Sim.data, fdr.thr=0.1, plot=FALSE)
res.mo <- modT.Prot(Sim.data, fdr.thr=0.1, plot=FALSE)


###################################################
### code chunk number 13: prot2D.Rnw:176-177
###################################################
featureNames(res.stud) 


###################################################
### code chunk number 14: prot2D.Rnw:180-181
###################################################
featureNames(res.mo) 


###################################################
### code chunk number 15: prot2D.Rnw:184-185
###################################################
notes(Sim.data)$SpotSig 


###################################################
### code chunk number 16: prot2D.Rnw:190-191
###################################################
sessionInfo()