###################################################
### chunk number 1:  eval=FALSE
###################################################
## #line 56 "vignettes/CALIB/inst/doc/quickstart.Rnw"
## library(CALIB)
## calibReadMe()


###################################################
### chunk number 2: 
###################################################
#line 98 "vignettes/CALIB/inst/doc/quickstart.Rnw"
library(CALIB)
path<-system.file("arraydata", package="CALIB")
dir(system.file("arraydata", package="CALIB")) 


###################################################
### chunk number 3: 
###################################################
#line 107 "vignettes/CALIB/inst/doc/quickstart.Rnw"
datapath <- system.file("arraydata", package="CALIB")
targets <- readTargets("targets.txt",path=datapath)
targets


###################################################
### chunk number 4: 
###################################################
#line 117 "vignettes/CALIB/inst/doc/quickstart.Rnw"
RG <- read.rg(targets$FileName,columns=list(Rf="CH1_NBC_INT",Gf="CH2_NBC_INT",Rb="CH1_SPOT_BKGD",Gb="CH2_SPOT_BKGD",RArea="CH1_SPOT_AREA",GArea="CH2_SPOT_AREA"),path=datapath)


###################################################
### chunk number 5: 
###################################################
#line 123 "vignettes/CALIB/inst/doc/quickstart.Rnw"
filename <- "annotation.txt"
fullname <- file.path(datapath,filename)
annotation <- read.table(file=fullname,header=T,fill=T,quote="",sep="\t")
RG$genes <- annotation


###################################################
### chunk number 6: 
###################################################
#line 132 "vignettes/CALIB/inst/doc/quickstart.Rnw"
types<-readSpotTypes(path=datapath)
types
spotstatus<-controlStatus(types,RG$genes)
RG$genes$Status<-spotstatus


###################################################
### chunk number 7: 
###################################################
#line 141 "vignettes/CALIB/inst/doc/quickstart.Rnw"
concfile<-"conc.txt"
spike<-read.spike(RG,file=concfile,path=datapath)


###################################################
### chunk number 8: spikeci
###################################################
#line 148 "vignettes/CALIB/inst/doc/quickstart.Rnw"
arraynum <- 1
plotSpikeCI(spike,array=arraynum)


###################################################
### chunk number 9: 
###################################################
#line 164 "vignettes/CALIB/inst/doc/quickstart.Rnw"
parameter<-estimateParameter(spike,RG,bc=F,area=T,errormodel="M")


###################################################
### chunk number 10: spikehi
###################################################
#line 170 "vignettes/CALIB/inst/doc/quickstart.Rnw"
plotSpikeHI(spike,parameter,array=arraynum)


###################################################
### chunk number 11: 
###################################################
#line 186 "vignettes/CALIB/inst/doc/quickstart.Rnw"
array<-c(1,1,2,2)
condition<-c(1,2,2,1)
dye<-c(1,2,1,2)
idcol<-"CLONE_ID"
## here, we normalize the first ten genes as example.
cloneid<-RG$genes[1:10,idcol]
normdata<-normalizeData(RG,parameter,array,condition,dye,idcol=idcol,cloneid=cloneid)
normdata