Teria um set de dados pra rodar o CMR? De antemão parece haver problema na execução da rnorm() cuja saída está sendo alocada em um objeto indexado que ainda não existe ('z'). A cada ciclo do for será gerada uma linha com 'npart' elementos e imagino que você queira acumular as saídas, mas essa não é a forma correta. Pode funcionar definindo z como uma lista e adicionando a saída em cada ciclo do for.

Trabalhando só no primeiro trecho pensei em algo assim...