[valse.git] / R / initSmallEM.R

vec_bin = function(X,r)
{
	Z = c()
	indice = c()
	j = 1
	for (i in 1:length(X))
	{
		if(X[i] == r)
		{
			Z[i] = 1
			indice[j] = i
			j=j+1
		} else
			Z[i] = 0
	}
	return (list(Z=Z,indice=indice))
}

initSmallEM = function(k,X,Y,tau)
{
	n = nrow(Y)
	m = ncol(Y)
	p = ncol(X)

	betaInit1 = array(0, dim=c(p,m,k,20))
	sigmaInit1 = array(0, dim = c(m,m,k,20))
	phiInit1 = array(0, dim = c(p,m,k,20))
	rhoInit1 = array(0, dim = c(m,m,k,20))
	piInit1 = matrix(0,20,k)
	gamInit1 = array(0, dim=c(n,k,20))
	LLFinit1 = list()

	require(MASS) #Moore-Penrose generalized inverse of matrix
	for(repet in 1:20)
	{
		clusters = hclust(dist(y)) #default distance : euclidean
Commit	Line	Data
39046da6 BA	1	vec_bin = function(X,r)
	2	{
	3	Z = c()
	4	indice = c()
	5	j = 1
	6	for (i in 1:length(X))
	7	{
	8	if(X[i] == r)
	9	{
	10	Z[i] = 1
	11	indice[j] = i
	12	j=j+1
	13	} else
	14	Z[i] = 0
	15	}
	16	return (list(Z=Z,indice=indice))
35b42a4b B	17	}
35b42a4b B	18
39046da6 BA	19	initSmallEM = function(k,X,Y,tau)
	20	{
	21	n = nrow(Y)
	22	m = ncol(Y)
	23	p = ncol(X)
35b42a4b	24
39046da6 BA	25	betaInit1 = array(0, dim=c(p,m,k,20))
	26	sigmaInit1 = array(0, dim = c(m,m,k,20))
	27	phiInit1 = array(0, dim = c(p,m,k,20))
	28	rhoInit1 = array(0, dim = c(m,m,k,20))
	29	piInit1 = matrix(0,20,k)
	30	gamInit1 = array(0, dim=c(n,k,20))
	31	LLFinit1 = list()
	32
	33	require(MASS) #Moore-Penrose generalized inverse of matrix
	34	for(repet in 1:20)
	35	{
	36	clusters = hclust(dist(y)) #default distance : euclidean
	37