如何在R中选择特定值？

米兰达·文特雷拉

我需要修改这些代码以便一次考虑一个人：

plot(breaths$breathing~as.factor(breaths$event))
plot(breaths$breathing~as.factor(breaths$event),
     names=c("control","grooming"),
     xlab="event", ylab="breathing rate (b/m)", pch=19)

在我的数据集中，这些是以下列：事件、呼吸和 id。事件要么是修饰，要么是控制。每个个体 (id) 都分为两类（梳理、控制）。我只想为每个人绘制一个图表。您认为我应该添加哪些代码？提前致谢！

r2evans

样本数据：

set.seed(42)
dat <- data.frame(event = sample(c("grooming", "control"), 100, replace=TRUE), id = sample(c("Aa", "Bb", "Cc"), 100, replace=TRUE), breathing = runif(100))
head(dat)
#      event id  breathing
# 1 grooming Bb 0.35110692
# 2 grooming Aa 0.15902238
# 3 grooming Bb 0.30409800
# 4 grooming Aa 0.01754832
# 5  control Bb 0.99655268
# 6  control Cc 0.80439331

一种vis方法：

library(ggplot2)
ggplot(dat, aes(breathing)) +
  geom_density(aes(color = event)) +
  facet_wrap(id ~ .)

或者如果您更喜欢积分：

ggplot(dat, aes(event, breathing)) +
  geom_point(aes(color = event)) +
  facet_wrap(id ~ .)

这里的关键组件是facet_wrap将数据拆分为id.

侧面要点：

使用~公式方法时，一般优先使用data=而不包含breaths$在每个变量中；虽然不是必需的，但它确实会产生轻微的审美差异（轴名称）。（我也发现代码更容易阅读。）
```
par(mfrow = c(1,2))
plot(mtcars$disp ~ mtcars$mpg)
plot(disp ~ mpg, data = mtcars)
```
不清楚names=这里是什么；缺少样本数据，我收到有关"names" is not a graphical parameter. 如果您确定没有收到该警告（多次plot调用一次），那么您的数据不仅仅是一个data.frame.