[R-br] Não consigo fazer uma ANOVA com arquivo com 1e+06 registros.

[Luiz Roberto Martins Pinto]

Ok.
Vou enviar por meio do seu e-mail particular.
Luiz Roberto

Luiz Roberto Martins Pinto
Prof. Pleno/DCET/UESC
Laboratório de Estatística Computacional
Universidade Estadual de Santa Cruz
Ilhéus-Bahia

luizroberto.uesc em gmail.com
skype: lrmpinto
http://lattes.cnpq.br/2732314327604831




Em 24 de outubro de 2013 00:13, Mauro Sznelwar <sznelwar em uol.com.br>escreveu:

> **
> Eu não consegui abirir o arquivo de dados que colocou no datafilehost.
> Poderia me enviar anexado? Ou colocar no Mega Upload para baixar?
>
> O código abaixo toma partes de  scripts do PDF dos pacotes ff e ffbase.
>
> Eu testei com seus dados  e deu certo. Mas na hora de copiar e colar para
> o e-mail, houve erros... Acho que agora está certo.
> Mas ressalto, isso é só uma regressão linear. Não sei como calcular o
> two-way anova a partir disso.. somente uma anova comum. Mas estou certo de
> que tem jeito...
>
>
>  #####################################
>
> require(ff)
> require(ffbase)
> require(biglm)
>
> load("RCBD_Data.Rdata")
>
> dados.ff = as.ffdf(Data)
> rm(Data)
> gc()
>
> form = y ~ factor(Treat) + factor(block)
>
> for (i in chunk(dados.ff, by=25000)){
>   if (i[1]==1){
>     message("first chunk is: ", i[[1]],":",i[[2]])
>     biglmfit <- biglm(form, data=dados.ff[i,,drop=FALSE])
>   }else{
>     message("next chunk is: ", i[[1]],":",i[[2]])
>     biglmfit <- update(biglmfit, dados.ff[i,,drop=FALSE])
>   }
> }
>
> summary(biglmfit)
>
>
> produzível.
>
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131024/da75a20c/attachment.html>