Description Usage Arguments Details Value Author(s) References See Also Examples
A Class mixnormerr
is declared in cubfits, and this is the function
to print and summary objects.
1 2 |
x |
an object with the class attributes. |
digits |
for printing out numbers. |
... |
other possible options. |
This is an useful function for summarizing and debugging.
The results will cat or print on the STDOUT by default.
Wei-Chen Chen wccsnow@gmail.com.
https://github.com/snoweye/cubfits/
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | ## Not run:
suppressMessages(library(cubfits, quietly = TRUE))
### Get individual of phi.Obs.
GM <- apply(yassour[, -1], 1, function(x) exp(mean(log(x[x != 0]))))
phi.Obs.all <- yassour[, -1] / sum(GM) * 15000
phi.Obs.all[phi.Obs.all == 0] <- NA
### Run optimization.
X <- log(as.matrix(phi.Obs.all))
param.init <- list(K = 2, prop = c(0.95, 0.05), mu = c(-0.59, 3.11),
sigma2 = c(1.40, 0.59), sigma2.e = 0.03)
ret <- mixnormerr.optim(X, K = 2, param = param.init)
print(ret)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.