R（工作室）因子与水平

Question

我是一名正在学习与R（工作室）合作的学生，为此我收到了一项任务。 我应该比较一些随机生成的数据，并据此得出结论。

但是，我遇到的问题是该数据具有5个级别的因数，我想一次比较一个级别的数据...

> str(data)
'data.frame':   275 obs. of  5 variables:  
 $ leverancier        : Factor w/ 5 levels "A","B","C","D",..: 1 1 1 1 1 1 1 1 1 1 ...   
 $ las_trekproef      : num  211 375 503 195 221 ...  
 $ las_score          : num  2.6 3.1 4.3 2.6 2.7 3.3 3.9 2.2 2.7 2.7 ...  
 $ afwijking_draaiwerk: num  0.081 0.061 0.015 0.125 0.256 0.004 0.124 0.016 0.042 0.062 ...  
 $ afwijking_freeswerk: num  0.336 0.026 0.032 0.161 0.36 0.447 0.062 0.176 0.317 0.212 ...

我想做的是将每个级别的数据放在不同的变量中，这样我就可以一次将图形和方框图设为一个级别。

A <- （此处为A级的所有数据）

代替：

summary(data)  
 leverancier las_trekproef     las_score     afwijking_draaiwerk afwijking_freeswerk
 A:55        Min.   :128.6   Min.   :2.000   Min.   :0.0000      Min.   :0.0010     
 B:55        1st Qu.:270.6   1st Qu.:2.700   1st Qu.:0.0355      1st Qu.:0.1210     
 C:55        Median :361.0   Median :3.500   Median :0.0760      Median :0.2670     
 D:55        Mean   :356.1   Mean   :3.513   Mean   :0.1268      Mean   :0.3055     
 E:55        3rd Qu.:443.6   3rd Qu.:4.250   3rd Qu.:0.1340      3rd Qu.:0.4255     
             Max.   :571.4   Max.   :5.000   Max.   :1.1390      Max.   :1.3890

提前致谢，

NH

Answer 1

你的意思是

a <- subset(data, leverancier=='A')

要获得帮助，只需查看?subset

Answer 2

您可以按因子将tapply用于摘要。 如果您使用公式界面，并且x是一个因子，则Boxplot将按因子绘制变量。 您还可以将因子索引到y子集以生成单个箱线图。 但是，为了基于因子比较y，您需要在同一图上绘制级别。 这里有些例子。

# Create example data
dat <- data.frame(leverancier=rep(c("A","A","B","B","B","A","C","D","D","C"),100),
                  las_trekproef=runif(1000,100,500),
                  las_score=runif(1000,1,4))

# Use tapply to summarize y by factor                 
tapply(dat$las_score, dat$leverancier, FUN=summary)

# Using formula interface plot y by factor 
boxplot(las_score ~ leverancier, data=dat, notch=TRUE)

# You can also index y based on a factor level to create a single boxplot of y
boxplot(dat[dat$leverancier == "A" ,]$las_score, notch=TRUE)

Answer 3

考虑使用split

split(data, data$leverancier)

将为您提供data.frame的list ，每个list对应一个级别的leverancier 。 然后，您可以一次对每个元素进行操作，也可以遍历列表依次对每个部分进行操作。

我意识到这并不能直接回答您的问题（Seb的答案就是那样），但是它应该为您提供处理R中数据的更惯用的方向。