我有一个实验，包括测试受试者爬上树枝所需的时间，我们在 2016 年至 2019 年测试了这些受试者（男性和女性），其中一些正在接受治疗 A 或 B。我需要构建一个条形图比较 2016 年女性接受每种治疗的时间，与 2017 年女性接受每种治疗的时间，依此类推。男性完全一样。
我的数据存储在这样的数据库中：
Year Treatment ID Sex Time(Mins)
2016 A K1 M 1,63
2016 A H4 F 5
2017 B J9 F 3,5
2018 ...
2019 ...
我尝试使用此代码，其中 BD 是我的数据库名称：
barplot(c(mean(BD$Time[bd2016$Sex=="F" & BD$Treatment=="A" & BD$Year=="2019"]),
(mean(BD$Time[bd2016$Sex=="F" & BD$Treatment=="B" & BD$Year=="2019"]]))),
ylim= c(0,40), ylab="xx",names=c("A","B"))