R中的循环分析

Question

I am trying to do some analysis on several variables. I have a dataset called "data". Thanks to someone's suggestion here I have managed to make histograms in a loop like this:

my_vars <- c("Var1","Var2","Var3","Var4") # I have changed the names of the variables for "privacy"

par(mfcol=c(1,1))
for(v in my_vars) {
  hist(data[[v]], probability = TRUE, main=v, xlab="")
  lines(density(data[[v]]), col=2)
}

Now, I am trying to do the same thing with some analysis like this:

for(v in my_vars) {
  shapRes= shapiro.test(data[[v]])
  shapRes
}

However, when I run this latest part this is all I get in the console:

> for(v in my_vars) {
+   shapRes= shapiro.test(data[[v]])
+   shapRes
+ }
> 

If I open "shapRes" this is what I see:

Shapiro-Wilk normality test

data: data[[v]]
W = 0.9561, p-value = 0.09105

I think that this is just the result of the last one. Does anybody know how can I fix this?

Thank you in advance for your time.

Answer 1

You need to store each shapRes. I would recommend to use a list to store them:

mylist <- list()

for(v in my_vars) {
  shapRes= shapiro.test(data[[v]])
  mylist[[v]] <- shapRes
}