[R] a new-bie question about obtaining certain value from the print out
Rolf Turner
r.turner at auckland.ac.nz
Wed Aug 29 23:31:54 CEST 2007
On 30/08/2007, at 8:53 AM, kurt Zhao wrote:
> Hi everyone,
>
> I am quite new to R.
>
> my command is
> t.test(c(1:5,7:11), y=c(1:10),alternative = c("two.sided"), paired
> = TRUE)
>
> The output is
>
> Paired t-test
>
> data: c(1:5, 7:11) and c(1:10)
> t = 3, df = 9, p-value = 0.01496
> alternative hypothesis: true difference in means is not equal to 0
> 95 percent confidence interval:
> 0.1229738 0.8770262
> sample estimates:
> mean of the differences
> 0.5
> According to the reference:
>
> Value
> A list with class "htest" containing the following components:
> p.value the p-value for the test.
>
> How can I get only the p.value as the output instead of the whole
> print out?
t.test(c(1:5,7:11), y=c(1:10),alternative = c("two.sided"), paired =
TRUE)$p.value
Or:
rslt <- t.test(c(1:5,7:11), y=c(1:10),alternative = c("two.sided"),
paired = TRUE)
rslt$p.value
Look at names(rslt) and/or str(rstl) to see what else you can dig out.
A side-comment: You ***do not need*** the c() wrappers for 1:10 and
"two.sided";
the c() function is for ``combine'' or ``catenate''. You don't need
to ``combine'' a single
item! The wrapper does no harm, but it's a waste of key strokes and
indicates a certain
lack of insight which *could* do harm in more complex situations.
cheers,
Rolf Turner
######################################################################
Attention:\ This e-mail message is privileged and confidenti...{{dropped}}
More information about the R-help
mailing list