R/statvc4levelt2.R

In mvctm: Multivariate Variance Components Tests for Multilevel Data

statvc4levelt2 <-
function(y,cluster,m1,n1,m2,n2,m3,n3,p,weight){

# y = Matrix Nxp (N obs, p=dimension) of (maybe standardized; signs or ranks) data 
# cluster = List of indices (each element gives the indices of the obs in a level 3 cluster nested within a level 2 cluster nested within a level 1 cluster). From the "preparecluster" function.
# N = Total number of observations
# p = Dimension
# m1 = Number of level 1 clusters
# n1 = Number of observations in each level 1 cluster (vector m1X1)
# m2 = Number of level 2 clusters within each level 1 cluster (vector m1X1)
# n2 = Number of observations in each level 2 cluster within each level 1 cluster (list of size m1)
# m3 = Number of level 3 clusters nested within each combination of level 1 and level 2 clusters (list of size m1)
# n3 = Number of observations in each level 3 cluster nested within each combination of level 1 and level 2 clusters (list of size m1)


mat=matrix(0,p,p)
w=0

for(i in 1:m1)
	{
	clusteri=cluster[[i]]
	m1i=m2[i]
	m2i=m3[[i]]
	n2i=n3[[i]]			
	mati=statvc3levelt1(y,clusteri,m1i,m2i,n2i,p,weight)
	mat=mat+mati[[1]]
	w=w+mati[[2]]
	}

list(mat,w)
}