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＆lt; - （此处为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)

将为您提供list data.frame个，每个leverancier对应一个{{1}}级别。然后，您可以一次对每个元素进行操作，或者循环遍历列表以依次对每个元素进行操作。

我意识到这并没有直接回答你的问题（Seb的回答是这样的），但它应该指出你在R中处理数据的更惯用的方向。

R（工作室）因子与水平

3 个答案: