Likelihood Function of the NB-Beta Model In Multiple...

'LogN' specifies the function to run (one round of) the EM algorithm for the NB-beta model.

1 2 | ```
LogN(Input, InputSP, EmpiricalR, EmpiricalRSP, NumOfEachGroup,
AlphaIn, BetaIn, PIn, NoneZeroLength)
``` |

`Input, InputSP` |
The expressions among all the samples. |

`NumOfEachGroup` |
Number of genes in each Ng group. |

`AlphaIn, PIn, BetaIn, EmpiricalR, EmpiricalRSP` |
The parameters from the last EM step. |

`NoneZeroLength` |
Number of Ng groups. |

Ning Leng

Ning Leng, John A. Dawson, James A. Thomson, Victor Ruotti, Anna I. Rissman, Bart M.G. Smits, Jill D. Haag, Michael N. Gould, Ron M. Stewart, and Christina Kendziorski. EBSeq: An empirical Bayes hierarchical model for inference in RNA-seq experiments. Bioinformatics (2013)

1 2 3 4 5 6 | ```
#Input = matrix(rnorm(100,100,1), ncol=10)
#rownames(Input) = paste("g",1:10)
#RIn = matrix(rnorm(100,200,1), ncol=10)
#res = LogN(Input, list(Input[,1:5], Input[,6:10]),
# RIn, list(RIn[,1:5], RIn[,6:10]),
# 10, .6, .7, .3, 1)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.