Considere os scripts de 11 a 14 aqui http://leg.ufpr.br/~walmes/cursoR/conest2013/, o 11 e o 12 aqui http://leg.ufpr.br/~walmes/cursoR/cnpaf3/.

À disposição.

Walmes.