Displaying 2 results from an estimated 2 matches for "mixm".
Did you mean:
mix
2009 May 22
0
EM algorithm mixture of multivariate
...#############
library(mvtnorm)
libray(scatterplot3d)
library(MASS)
n=100
m1=c(5,-5)
m2=c(-3,3)
s1=matrix(c(2,1,1,3), 2,2)
s2=matrix(c(4,1,1,6), 2,2)
alpha=0.3
c1 <- mvrnorm(round(n*alpha),m1,s1)
c2 <- mvrnorm(round(n*(1-alpha)),m2,s2)
allval <- rbind(c1,c2)
x <- allval[sample(n,n),]
mixm<-
function(x,m1,m2,s1,s2, alpha)
{
f1<-dmvnorm(x, m1, s1, log=FALSE)
f2<-dmvnorm(x, m2, s2, log=FALSE)
f=alpha*f1+(1-alpha)*f2
f
}
plot(x)
den<-mixm(x,m1,m2,s1,s2,alpha)
scatterplot3d(x[,1],x[,2],den, highlight.3d=TRUE, col.axis="blue",
col.grid="ligh...
2009 May 22
0
EM algorithm mixture of multivariate gaussian
...#############
library(mvtnorm)
libray(scatterplot3d)
library(MASS)
n=100
m1=c(5,-5)
m2=c(-3,3)
s1=matrix(c(2,1,1,3), 2,2)
s2=matrix(c(4,1,1,6), 2,2)
alpha=0.3
c1 <- mvrnorm(round(n*alpha),m1,s1)
c2 <- mvrnorm(round(n*(1-alpha)),m2,s2)
allval <- rbind(c1,c2)
x <- allval[sample(n,n),]
mixm<-
function(x,m1,m2,s1,s2, alpha)
{
f1<-dmvnorm(x, m1, s1, log=FALSE)
f2<-dmvnorm(x, m2, s2, log=FALSE)
f=alpha*f1+(1-alpha)*f2
f
}
plot(x)
den<-mixm(x,m1,m2,s1,s2,alpha)
scatterplot3d(x[,1],x[,2],den, highlight.3d=TRUE, col.axis="blue",
col.grid="ligh...