如何在R中按群组应用shapiro测试？(How to apply a shapiro test by groups in R?)

如何在R中按群组应用shapiro测试？(How to apply a shapiro test by groups in R?)

我有一个数据框，其中我的所有90个变量都有整数数据，类型如下：

代码| variable1 | 变量2 | 变量3 | ...

AB | 2 | 3 | 10 | ...

AH | 4 | 6 | 8 | ...

BC | 1 | 5 | 9 | ...

...... | ...... | ...... | ...

我想通过变量将shapiro测试（shapiro.test {stats}）应用于我的数据框，并将结果写在如下的表中：

variable_name | W | p值

有人有线索吗？

I have a dataframe where all my 90 variables have integer data, of the type:

code | variable1 | variable2 | variable3 | ...

AB | 2 | 3 | 10 | ...

AH | 4 | 6 | 8 | ...

BC | 1 | 5 | 9 | ...

... | ... | ... | ...

I want to apply a shapiro test (shapiro.test {stats}) to my dataframe by variable and write the results in a table like:

variable_name | W | p-value

Does anyone have a clue?

最满意答案

使用来自R的mtcars数据

mydata<-mtcars kk<-Map(function(x)cbind(shapiro.test(x)$statistic,shapiro.test(x)$p.value),mydata) library(plyr) myout<-ldply(kk) names(myout)<-c("var","W","p.value") myout var W p.value 1 mpg 0.9475648 1.228816e-01 2 cyl 0.7533102 6.058378e-06 3 disp 0.9200127 2.080660e-02 4 hp 0.9334191 4.880736e-02 5 drat 0.9458838 1.100604e-01 6 wt 0.9432578 9.265551e-02 7 qsec 0.9732511 5.935208e-01 8 vs 0.6322636 9.737384e-08 9 am 0.6250744 7.836356e-08 10 gear 0.7727857 1.306847e-05 11 carb 0.8510972 4.382401e-04

Using mtcars data from R

mydata<-mtcars kk<-Map(function(x)cbind(shapiro.test(x)$statistic,shapiro.test(x)$p.value),mydata) library(plyr) myout<-ldply(kk) names(myout)<-c("var","W","p.value") myout var W p.value 1 mpg 0.9475648 1.228816e-01 2 cyl 0.7533102 6.058378e-06 3 disp 0.9200127 2.080660e-02 4 hp 0.9334191 4.880736e-02 5 drat 0.9458838 1.100604e-01 6 wt 0.9432578 9.265551e-02 7 qsec 0.9732511 5.935208e-01 8 vs 0.6322636 9.737384e-08 9 am 0.6250744 7.836356e-08 10 gear 0.7727857 1.306847e-05 11 carb 0.8510972 4.382401e-04