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

[Elias T Krainski]

Mauro, apos fazer download do arquivo vc carrega ele no R com
    load("RCBD_Data.Rdata")
mas lembre-se que para isso funcionar o resultado de
    getwd()
deve ser o diretorio onde o arquivo se encontra.

On 10/24/2013 04:13 AM, Mauro Sznelwar wrote:
> 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.
>