我有关于某些鱼类捕获率的数据。
fish 1 fish 2 fish 3
0.000 3.265 9.872
2.147 1.013 0.000
我使用以下方法计算了每条鱼的平均捕获率:
a <- colMeans(df)
我还计算了标准误差:
stdError <- (sapply(df,sd))/sqrt(length(df))
我使用以下方法创建了一个dotplot:
dotplot(a, xlab="mean catch", ylab = "species",las =2,)
如何在此图中添加误差线?如果可能的话,我宁愿不使用ggplot。我目前正在使用R中的内置函数,但可以访问莱迪思。
对不起可能是一个基本问题,我对R中的情节完全是新的。
答案 0 :(得分:7)
dotplot
是一个点阵函数,大多数默认的格函数都不具备对置信区间的强大支持。 Hmisc
包扩展了大多数晶格函数,以更好地结合置信区间。
以下是您如何使用它的示例。请注意,我们将您想要绘制的数据组合到data.frame中,以便我们可以使用正确的公式symtax
mm<-data.frame(a,stdError, fish=names(a))
library(lattice)
library(Hmisc)
Dotplot(fish~Cbind(a, a-stdError, a+stdError), mm,
xlab="mean catch", ylab = "species",las =2)
这会产生
请注意,Hmisc
版本的函数名为Dotplot
,lattice
版本称为dotplot
;资本化很重要。
在这里,我只是从平均值中添加/减去了一个标准误差。您可以随意计算置信区间。