3 %- Also NEED an '\alias' for EACH other topic documented here.
5 9 Density plots to compare the empirical q's and the simulated q's from the fitted beta distribution.
8 Check the beta fit of 9 different groups on isoform level data.
11 DenNHistTable(QList, Alpha, Beta, AList = "F")
13 %- maybe also 'usage' for other objects documented here.
16 The estimated q's from the output of NBBetaEB.bias.uniqueP_PoolVarSpeedUp_MDFPoi_NoNormVar
17 . Input should be a list of different groups of transcripts. The number of lists here should be 9.
22 The fitted parameter alpha from the output of NBBetaEB.bias.uniqueP_PoolVarSpeed
23 Up_MDFPoi_NoNormVar. Input should be a number if AList is not defined.
26 The fitted parameter beta from the output of NBBetaEB.bias.uniqueP_PoolVarSpeedUp_MDFPoi_NoNormVar.
27 Input could be one single number or a vector of several numbers. The length of the input should be 9.
32 Whether a list of alpha's are used
41 A plot contains 9 dub-plots will be generated.
42 The empirical estimation of q's will be represented as blue histogram and the density of
43 the fitted beta distribution will be represented as the green line.
44 The main title of the plot will be "GroupName name"
47 DenNHist, beta.mom, QQP, NBBetaEB.bias.uniqueP_PoolVarSpeedUp_MDFPoi_NoNormVar
56 %% ~Make other sections like Warning with \section{Warning }{....} ~
59 %% ~~objects to See Also as \code{\link{help}}, ~~~
63 IsoGenerate=IsoSimu(DVDconstant=NULL, DVDqt1=.97, DVDqt2=.98, Conditions=as.factor(rep(c(1,2),each=5)), NumofSample=10, NumofIso=NULL, DEIsoProp=.1, Phiconstant=NULL, Phi.qt1=.25, Phi.qt2=.75, OnlyData="Y" )
64 IsoList=IsoGenerate$data
66 # Get Ng Vector, 5End Vector and 3End Vector
67 ngv=c(1,2,3,2,3,2,3,2,3)
68 b3v=c(1,0,0,1,1,0,0,1,1)
69 b5v=c(1,0,0,0,0,1,1,1,1)
70 NgV=unlist(sapply(1:9,function(i)rep(ngv[i],dim(IsoList[[i]])[1])))
71 Bias3V=unlist(sapply(1:9,function(i)rep(b3v[i],dim(IsoList[[i]])[1])))
72 Bias5V=unlist(sapply(1:9,function(i)rep(b5v[i],dim(IsoList[[i]])[1])))
75 IsoData=do.call(rbind,IsoList)
76 IsoEBres=NBBetaEB.bias.uniqueP_PoolVarSpeedUp_MDFPoi_NoNormVar(Data=IsoData, NgVector=NgV, Vect5End=Bias5V, Vect3End=Bias3V, Conditions=as.factor(rep(c(1,2),each=5)),sizeFactors=rep(1,10), maxround=5)
78 DenNHistTable(QList=IsoEBres$QList1, Alpha=IsoEBres[[1]][5,], Beta=IsoEBres[[2]][5,], AList="F")
82 % Add one or more standard keywords, see file 'KEYWORDS' in the
83 % R documentation directory.
85 \keyword{ ~kwd2 }% __ONLY ONE__ keyword per line