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
[email protected]
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.
_______________________________________________
R-br mailing list
[email protected]
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.
