plotting standard errors in R

Question

I have data on the catch rate of certain species of fish.

fish 1  fish 2  fish 3
0.000   3.265   9.872
2.147   1.013   0.000

I have calculated the mean catch rate for each fish using:

a <- colMeans(df)

I have also calculated the standard error:

stdError <- (sapply(df,sd))/sqrt(length(df))

I have created a dotplot using:

dotplot(a, xlab="mean catch", ylab = "species",las =2,)

How do I add error bars to this plot? I would prefer not to use ggplot if possible. I am currently using the inbuilt functions in R but have access to Lattice.

Sorry for what is probably a basic question, I am entirely new to plots in R.

MrFlick · Accepted Answer · 2014-08-11 16:05:47Z

7

dotplot is a lattice function and most default lattice functions don't have great support for confidence intervals. The Hmisc package extends most lattice functions to better incorporate confidence intervals.

Here's an example of how you would use it. Note that we combine the data you want to plot here into a data.frame so we can use a proper formula 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)

and this produces

enter image description here

Note that the Hmisc version of the function is called Dotplot while the lattice version is called dotplot; capitalization matters.

Here I've just added/subtracted one standard error from the mean. You can calculate the confidence intervals however you like.

answered Aug 11, 2014 at 16:05

MrFlick

209k19 gold badges300 silver badges324 bronze badges

Sign up to request clarification or add additional context in comments.

Collectives™ on Stack Overflow

plotting standard errors in R

1 Answer 1

Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

Comments

Your Answer

Sign up or log in

Post as a guest

Related