File: multmixinit.R

package info (click to toggle)
r-cran-mixtools 2.0.0.1-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 1,944 kB
  • sloc: ansic: 781; makefile: 6
file content (19 lines) | stat: -rwxr-xr-x 415 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
multmix.init <- function(y, lambda = NULL, theta = NULL, k = 2){

n <- nrow(y)
p <- ncol(y)

    if (is.null(theta)) {
        theta = matrix(runif(p * k), k, p)
        theta = theta/apply(theta, 1, sum)
    }
    else k = nrow(theta)
    if (is.null(lambda)) {
        lambda = runif(k)
        lambda = lambda/sum(lambda)
    } else k = length(lambda)

list(lambda=lambda, theta=theta, k=k)


}