mstep: M-step to maximize theta
In sonicLength: Estimating Abundance of Clones from DNA Fragmentation Data
Description
Usage
Arguments
Details
Value
Author(s)
See Also
Examples
Description
This is a utility function for maxEM
Usage
1
mstep(x, theta, phi)
Arguments
x
a zero-one indicator matrix whose rows correspond to unique lengths with rownames
indicating those lengths
theta
a vector of the abundance estimates. length(theta)==ncol(x)
phi
a vector of the probabilities of sonicant lengths. length(phi)==nrow(x)
Details
The M-step for theta is computed. Probably, there is no need to use
this function directly, but just in case it is here.
Value
a vector like the input theta
Author(s)
Charles C. Berry
See Also
Examples
1
2
3
4
5
mat <- diag(10)
mat[ ,10 ] <- 1.0
phi1 <- prop.table( rep(1,10))
theta1 <- 1:10
sonicLength:::mstep( mat, theta1,phi1)
sonicLength documentation built on Sept. 20, 2021, 9:06 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Description Usage Arguments Details Value Author(s) See Also Examples
Description
This is a utility function for maxEM
Usage
1 | mstep(x, theta, phi)
|
Arguments
x |
a zero-one indicator matrix whose rows correspond to unique lengths with rownames indicating those lengths |
theta |
a vector of the abundance estimates. |
phi |
a vector of the probabilities of sonicant lengths. |
Details
The M-step for theta is computed. Probably, there is no need to use this function directly, but just in case it is here.
Value
a vector like the input theta
Author(s)
Charles C. Berry
See Also
Examples
1 2 3 4 5 | mat <- diag(10)
mat[ ,10 ] <- 1.0
phi1 <- prop.table( rep(1,10))
theta1 <- 1:10
sonicLength:::mstep( mat, theta1,phi1)
|
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.