[R] A question from a newbee...

Petr Pikal petr.pikal at precheza.cz
Fri Aug 2 10:59:40 CEST 2002



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.00294

My results

> summ_summary(anov)
> summ
            Df 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 Lamoral

anova() 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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list