EBSeq/man/DenNHistTable.Rd

   1 \name{DenNHistTable}
   2 \alias{DenNHistTable}
   3 %- Also NEED an '\alias' for EACH other topic documented here.
   4 \title{
   5 9 Density plots to compare the empirical q's and the simulated q's from the fitted beta distribution.
   6 }
   7 \description{
   8 Check the beta fit of 9 different groups on isoform level data.
   9 }
  10 \usage{
  11 DenNHistTable(QList, Alpha, Beta,  AList = "F")
  12 }
  13 %- maybe also 'usage' for other objects documented here.
  14 \arguments{
  15   \item{QList}{
  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.
  18
  19 }
  20
  21   \item{Alpha}{
  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.
  24 }
  25   \item{Beta}{
  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.
  28
  29 }
  30
  31   \item{AList}{
  32 Whether a list of alpha's are used
  33
  34 }
  35
  36   }
  37 \details{
  38
  39 }
  40 \value{
  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"
  45 }
  46 \references{
  47 DenNHist, beta.mom, QQP, NBBetaEB.bias.uniqueP_PoolVarSpeedUp_MDFPoi_NoNormVar
  48 }
  49 \author{
  50 Ning Leng
  51 }
  52 \note{
  53 %%  ~~further notes~~
  54 }
  55
  56 %% ~Make other sections like Warning with \section{Warning }{....} ~
  57
  58 \seealso{
  59 %% ~~objects to See Also as \code{\link{help}}, ~~~
  60 }
  61 \examples{
  62 # Generate Data
  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
  65
  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])))
  73
  74 #Run EBSeq
  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)
  77
  78 DenNHistTable(QList=IsoEBres$QList1, Alpha=IsoEBres[[1]][5,], Beta=IsoEBres[[2]][5,], AList="F")
  79
  80
  81 }
  82 % Add one or more standard keywords, see file 'KEYWORDS' in the
  83 % R documentation directory.
  84 \keyword{ ~kwd1 }
  85 \keyword{ ~kwd2 }% __ONLY ONE__ keyword per line