QtoJackson: Q to Jackson: transform a Q-residual into a confidence...
In MEET: MEET: Motif Elements Estimation Toolkit

Description Usage Arguments Value Author(s) See Also Examples

Given Q-residual, this function uses the parameters calculated in JAcksonParameters to convert a Q-residual into a confidence interval

1	QtoJackson(Q,h0, x1, x2, x3)

`Q`	Q-residual
`h0`	h0 parameter
`x1`	x1 parameter
`x2`	x2 parameter
`x3`	x3 parameter

output:confidence interval of the given residual.

Erola Pairo <epairo@ibec.pcb.ub.es>

kfold.PCA, Predict.PCA, JacksonParameters

data(iicc)
    TFBS<-iicc$Transcriptionfactor
     Prob<-iicc$background
     missing<-iicc$missing
     NumericalMatrix<-numericalDNA(Prob)
    suma<-apply(TFBS,2,function(y){sum(y=="-")})
    threshold<-floor(nrow(TFBS)*missing/100)
    TFBS<-TFBS[, suma<=threshold]
    ncolTFBS<-ncol(TFBS)
    TFBSnum<-apply(TFBS,1,function(x){as.vector(t(NumericalMatrix[x,]))})
    TFBSnum<-t(TFBSnum)
    nPCs<-1
    model<-pca(TFBSnum, nPcs=nPCs, method="svd", center=TRUE)
    JacksonPars<-JacksonParameters(nPCs,TFBSnum)
    Qres<-3.45
    h0<-JacksonPars$h0
    x1<-JacksonPars$x1
    x2<-JacksonPars$x2
    x3<-JacksonPars$x3
    confidence<-QtoJackson(Qres,h0,x1,x2,x3)