###################################################
### chunk number 1: 
###################################################
library(logitT)


###################################################
### chunk number 2: 
###################################################
library(SpikeInSubset)
data(spikein95)  ## get the example data
groupvector<-c("A","A","A","B","B","B")  ## specify group affiliations
logitTex<-logitTAffy(spikein95, group=groupvector) 


###################################################
### chunk number 3: 
###################################################
logitTex[1:10]


###################################################
### chunk number 4: 
###################################################
groupvector<-c("A","A","A","B","B","B")


###################################################
### chunk number 5: 
###################################################
pvals <- (1-pt(abs(logitTex),df=4))*2              ## calculate p-values
signifgenes<-names(logitTex)[pvals<0.01]           ## find probe sets with p-values smaller than 0.01
signifgenes