Hi Bruce,

One way is via apply()

# some data
set.seed(123)
X <- matrix(rnorm(100), ncol = 5)
X

# tests
t(apply(X, 2, function(x){
            sw <- shapiro.test(x)
            c(sw$statistic, P = sw$p.value)
            }))

See ?apply and ?str and ?shapiro.test for more information.

HTH,
Jorge


On Tue, Apr 26, 2011 at 4:15 PM, Bruce Kindseth <> wrote:

> I have a large amount of data which I break down into a collection of
> vectors of 100-125 values each.  I would like to test the normality of the
> vectors and compare them.  In the interactive mode I can test any one
> vector
> using the Shapiro-Wilk test or the Kolmogorov-Smirnov test.  My problem is
> that when I try to write out the results to a file, the term output is a
> mixture of alpha characters along with the number value that I really want.
> How can I get it to write out just the p-values by themselves?
>
> Thank you very much for your time.
>
> Bruce Kindseth
>
>
>        [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

	[[alternative HTML version deleted]]

