[R-sig-ME] Export several lme outputs to a single excel file
Emmanuel Curis
emm@nue|@cur|@ @end|ng |rom p@r|@de@c@rte@@|r
Wed Sep 25 22:55:48 CEST 2019
Dear Mario,
Since anova( ) results are basically data.frames, a simple way could
be to add a column to identify the model, then row-bind the results.
Something like (using examples in the lme man page)
# First model
fm1 <- lme(distance ~ age, data = Orthodont)
anova( fm1 ) -> r1
r1$Md <- 1
# Second model
fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1)
anova( fm2 ) -> r2
r2$Md <- 2
# Bind results
r <- rbind( r1, r2 )
# Export
write.xlsx( r, "your_file.xlsx" )
For more model, can be adapted with do.call( rbind, list_of_results)
and either a for or lapply to generate the additional
columns... Depends on how you obtain your several models.
Best regards,
NB: note that, to avoid duplicate row names, terms with the same name
will be incremented - see the results of the example above. You can
avoid this by explicely add the term names in the data.frame ---
r1$Terms <- rownames( r1 ) --- and removing the row names of the final
result.
On Wed, Sep 25, 2019 at 01:33:46PM +0200, Mario Garrido wrote:
« Dear users,
« it is not such an statistical question but how to export results
« I fit several *lme* model s and I know how can I export the data to excel,
« this way:
« lme1 <- lme(Mean ~ x*y, data = Data, random = ~ 1|factor(ID))
« anova(lme1)->resultslme1
« And to export to excel since keep the data in ordered in rows/columns
« write.xlsx( resultslme1 ,"C:/Users/Desktop/resultslme1.xlsx")
«
« What I want to do is to export output from several lme at the same time,
« and adding the name/reference of the model. Something like this
« lme1
« numDF denDF F.value p.value
« (Intercept) 1 78 653,6152 0
« x 2 78 0,822612 0,443057
« y 2 39 0,479357 0,622781
« x:y 4 78 2,593825 0,042787
«
« lme2
« numDF denDF F.value p.value
« (Intercept) 1 78 653,6152 0
« x 2 78 0,822612 0,443057
« y 2 39 0,479357 0,622781
« x:y 4 78 2,593825 0,042787
«
«
« I tried with merge, abind, paste,... but I did not find the solution
«
« Thanks in advance
«
« --
« Mario Garrido Escudero, PhD
« Dr. Hadas Hawlena Lab
« Mitrani Department of Desert Ecology
« Jacob Blaustein Institutes for Desert Research
« Ben-Gurion University of the Negev
« Midreshet Ben-Gurion 84990 ISRAEL
«
« gaiarrido using gmail.com; gaadio using post.bgu.ac.il
« phone: (+972) 08-659-6854
«
« [[alternative HTML version deleted]]
«
« _______________________________________________
« R-sig-mixed-models using r-project.org mailing list
« https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
--
Emmanuel CURIS
emmanuel.curis using parisdescartes.fr
Page WWW: http://emmanuel.curis.online.fr/index.html
More information about the R-sig-mixed-models
mailing list