Sérgio, bom dia!
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...
### <BEGIN>
### não testado!
prop <- function(x){
z <- list(NULL)
for (i in 1:nrow(x)) {
z[i]=list(rnorm(x$npart[i], mean=x$idade[i], sd=x$sd[i]))}
return(z)
}
### <END>
Desse ponto seria só prosseguir com os cálculos utilizando cada 'linha' da
lista.
Éder Comunello <c <comunello.e...@gmail.com>omunello.e...@gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
