找不到对象错误

Question

首先我需要解释一下，我在 R 方面接受过最低限度的培训，并且对编码语言或 R 等程序的了解为 0，所以如果我提出愚蠢的问题或不了解基本知识，请见谅。

另外，我试图查看过去的主题/答案，但我很难将答案与我的数据联系起来，所以如果这个问题已经得到解答，我深表歉意。

基本上我有一个数据集，我试图在这个集合中找到两个变量的平均值（在寒冷中行走之前的峰值流量，以及在寒冷中行走之后的峰值流量）。这是我到目前为止使用的全部代码：

 drugs <- read.table(file = "C:\\Users\\Becky\\My Documents\\Asthmadata.txt", header = TRUE)

 drugs
 str(drugs)

 mean.Asthmadata <- tapply (Asthmadata$trial1, list(Asthmadata$PEFR1), mean)
 mean.Asthmadata

在 mean.Asthmadata 之前它工作正常。数据在 R 中与其他代码一起出现得很好，但是当我达到平均值并执行 mean.Asthmadata [...] 代码时，我不断收到相同的错误：“找不到对象'mean.Asthmadata'”

我的朋友使用了与我相同的代码，并且对他有用，所以我很困惑。难道我做错了什么？

谢谢

编辑： @BenBolker

这是我的数据集

trial1  PEFR1   trial2  PEFR2
Before  310 After   299
Before  242 After   201
Before  340 After   232
Before  388 After   312
Before  294 After   221
Before  251 After   256
Before  391 After   327
Before  401 After   331
Before  287 After   231

这是我用过的所有代码：

drugs <- read.table(file = "C:\\Users\\Becky\\My Documents\\Asthmadata.txt", header = TRUE)

drugs
str(drugs)

mean.drugs <- tapply (drugs$trial1, list(drugs$PEFR1), mean)
mean.drugs

我拥有的 R 版本有两个版本：i386 3.1.3 和 x64 3.1.3 – 我都尝试过，但似乎都没有达到我想要的效果。我也在使用 Windows 7 Home Premium 64 位。希望我已经包含了您需要的所有内容，如果我的格式关闭，我深表歉意——我还不能完全弄清楚如何在此处正确格式化。

当运行 Roland 提供的代码时，我现在得到的错误是：“split.default(X, group) 中的错误：第一个参数必须是一个向量”。所以我每次尝试都会遇到不同的错误——这一定是我做错了。

希望我已正确格式化所有内容并包含您需要的所有内容。谢谢:)

Answer 1

drugs <- read.table(header=TRUE,text="
trial1  PEFR1   trial2  PEFR2
Before  310 After   299
Before  242 After   201
Before  340 After   232
Before  388 After   312
Before  294 After   221
Before  251 After   256
Before  391 After   327
Before  401 After   331
Before  287 After   231")

在当前格式中，您可以通过以下方式计算之前和之后的平均值

mean(drugs$PEFR1)

和

mean(drugs$PEFR2)

你可能想到的是这个形状：

drugs2 <- with(drugs,
                data.frame(trial=c(as.character(trial1),
                                   as.character(trial2)),
                           PEFR=c(PEFR1,PEFR2)))

• 为方便起见，我使用了with()——这是一种临时附加数据框的方法，因此您可以直接引用其中的变量。）
• 在组合trial1 和trial2 时会遇到一些问题，因为它们会被强制转换为数字代码，这两种情况下都是1，除非您使用as.character() ...
• 您有要聚合的变量和要向后分组的变量的顺序（您希望通过 trial 聚合 PEFR，而不是相反）

mean.drugs <- with(drugs2,
                   tapply (PEFR, list(trial), mean))
##    After   Before 
## 267.7778 322.6667