dmultinomialACGT.Rd 1.41 KB
Newer Older
Eric Coissac committed
1 2
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/proba.R
3 4 5
\name{dmultinomialACGT}
\alias{dmultinomialACGT}
\title{Density probability function of a multinomial distribution for nucleotides}
Eric Coissac committed
6
\usage{
7
dmultinomialACGT(A, C, G, T, pA, pC, pG, pT, log = FALSE)
Eric Coissac committed
8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
}
\arguments{
\item{A}{numeric; number of observed *A* at a given site.}

\item{C}{numeric; number of observed *C* at a given site.}

\item{G}{numeric; number of observed *G* at a given site.}

\item{T}{numeric; number of observed *T* at a given site.}

\item{pA}{numeric; probability to observed a *A* at a given site.}

\item{pC}{numeric; probability to observed a *C* at a given site.}

\item{pG}{numeric; probability to observed a *G* at a given site.}

\item{pT}{numeric; probability to observed a *T* at a given site.}

\item{log}{logical; if TRUE, probabilities p are given as log(p).}
}
\value{
numeric; the density of probability of the log that probability
        of the observations given the probabilities of observation
}
\description{
Computes the probalility to observe a given number of `A`, `C`, `G`, `T`
when probalilities of each nucleotide are `pA`, `pC`, `pG`, `pT`.
}
\examples{
37 38 39 40 41
  dmultinomialACGT(c(3,4),c(3,6),c(3,2),c(3,5),
                   0.25,0.25,0.25,0.25)
  dmultinomialACGT(c(3,4),c(3,6),c(3,2),c(3,5),
                   pA = 0.25, pC = 0.25, pG = 0.25, pT = 0.25,
                   log=TRUE)
Eric Coissac committed
42 43 44 45 46

}
\author{
Eric Coissac
}