Hello . I have a little question .. I work with Anova (aov) and Manova (manova) and want to print some result in a text file. I see the results with the summary function or the summary.aov function but i do not know how to save the differents probabilities of error in a value. Sorry for my poor English Thanks a lot. -- Lamoral Julien -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
>>>>> "jlamoral" == jlamoral <jlamoral at ulb.ac.be> writes:jlamoral> Hello . I have a little question .. I work with jlamoral> Anova (aov) and Manova (manova) and want to print jlamoral> some result in a text file. I see the results with jlamoral> the summary function or the summary.aov function jlamoral> but i do not know how to save the differents jlamoral> probabilities of error in a value. Sorry for my jlamoral> poor English If you use myAov <- aov(.......) str(myAov) and even Saov <- summary(myAov) str(Saov) should `see' the STRucture (typically list components) of the resulting object. {but I have to admit I do not fully understand what you mean by ``differents probabilities of error''} jlamoral> Thanks a lot. jlamoral> -- Lamoral Julien You're welcome Martin Maechler <maechler at stat.math.ethz.ch> http://stat.ethz.ch/~maechler/ Seminar fuer Statistik, ETH-Zentrum LEO C16 Leonhardstr. 27 ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND phone: x-41-1-632-3408 fax: ...-1228 <>< -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
On 2 Aug 2002 at 9:46, jlamoral wrote:> Hello > Thanks for this reponse it is perfect for the sink function > but if i use per exemple : > i_c(1,1,2,2) > j_c(0.1,0.2,0.8,0.6) > anov_aov(j~i) > summ_summary(aov)do you have an object with name "aov" in your directory??? Check with ls() I recommend to remove it as it is not recommended to give objects built in function names.> summ > the answer is ... > > Df Sum Sq Mean Sq F value Pr(>F) > j 1 0.99412 0.99412 338 0.002946 ** > Residuals 2 0.00588 0.00294My results> summ_summary(anov) > summDf Sum Sq Mean Sq F value Pr(>F) i 1 0.3025 0.3025 24.2 0.03893 * Residuals 2 0.0250 0.0125 --- Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1> --- > Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 > > > for summ[1] > the answer is ... > [[1]] > Df Sum Sq Mean Sq F value Pr(>F) > j 1 0.99412 0.99412 338 0.002946 ** > Residuals 2 0.00588 0.00294 > --- > Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 > > and i want to put the value (Pr(>F)) 0.002946 into a variables but i > can not acces to this varible Re thanks a lot Julien Lamoralanova() gives you the result which is maybe more suitable for subsetting then summary try unlist(anova(anov)[5])[1] Pr(>F).1 0.03892554 it works just with your example, for more dimensions go through manual pages and some trial and error.> > >Hope this helps. Petr Pikal petr.pikal at precheza.cz p.pik at volny.cz -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._