Agreemat <- function(clusters)
{
    cdim = dim(clusters)
    nm = cdim[1]
    n = cdim[2]

    Amat = matrix(0,n,n)
    jmax=n-1
    for (j in 1:jmax)
    {   
        for (k in (j+1):n)
        {
            for (l in 1:nm)
            {
                if (clusters[l,j]==clusters[l,k])
                {
                    Amat[j,k] = Amat[j,k] + 1
                }
            }
        }
    }
    Agreemat = Amat
}