R-help Plot Error Bars
Contents |
here for a quick overview of the site Help Center Detailed answers to any questions you might have Meta Discuss the workings and policies of this site About Us error bars in r barplot Learn more about Stack Overflow the company Business Learn more about hiring developers error.bar function r or posting ads with us Stack Overflow Questions Jobs Documentation Tags Users Badges Ask Question x Dismiss Join the Stack scatter plot with error bars in r Overflow Community Stack Overflow is a community of 6.2 million programmers, just like you, helping each other. Join them; it only takes a minute: Sign up Scatter plot with error bars up vote
Error Bars In R Plot
21 down vote favorite 11 How can I generate the following plot in R? Points, shown in the plot are the averages, and their ranges correspond to minimal and maximal values. I have data in two files (below is an example). x y 1 0.8773 1 0.8722 1 0.8816 1 0.8834 1 0.8759 1 0.8890 1 0.8727 2 0.9047 2 0.9062 2 0.8998 2 0.9044 2 0.8960 errbar r .. ... r plot share|improve this question edited Oct 23 '12 at 15:10 Roland 74.2k463103 asked Oct 23 '12 at 14:29 sherlock85 1521313 Since you clearly don't want a boxplot, I changed the title of your question in order to reflect what you really want. –Roland Oct 23 '12 at 15:11 1 also plotrix::plotCI, gplots::plotCI, library("sos"); findFn("{error bar}") –Ben Bolker Oct 23 '12 at 17:29 add a comment| 5 Answers 5 active oldest votes up vote 52 down vote accepted First of all: it is very unfortunate and surprising that R cannot draw error bars "out of the box". Here is my favourite workaround, the advantage is that you do not need any extra packages. The trick is to draw arrows (!) but with little horizontal bars instead of arrowheads (!!!). This not-so-straightforward idea comes from the R Wiki Tips and is reproduced here as a worked-out example. Let's assume you have a vector of "average values" avg and another vector of "standard deviations" sdev, they are of the same length n. Let's make the abscissa just the number of these "measurements", so x <- 1:n. Using these, here come the plotting commands: plot(x, avg, ylim=rang
by over 573 bloggers. There are many ways to follow us - By e-mail: On Facebook: If you are an R blogger yourself you are invited to add r ggplot error bars your own R content feed to this site (Non-English R bloggers should add themselves-
R Summaryse
here) Jobs for R-usersStatistical Analyst @ Rostock, Mecklenburg-Vorpommern, GermanyData EngineerData Scientist – Post-Graduate Programme @ Nottingham, EnglandDirector, Real World Informatics &
R Arrows
Analytics Data Science @ Northbrook, Illinois, U.S.Junior statistician/demographer for UNICEF Popular Searches web scraping heatmap twitter maps time series animation boxplot shiny hadoop ggplot2 how to import image file to R trading finance latex eclipse http://stackoverflow.com/questions/13032777/scatter-plot-with-error-bars rstudio excel SQL ggplot quantmod knitr googlevis PCA market research rattle regression map tutorial coplot rcmdr Recent Posts Election 2016: Tracking Emotions with R and Python Data science for executives and managers The Worlds Economic Data, Shiny Apps and all you want to know about Propensity Score Matching! August Package Picks Slack all the things! Warsaw R-Ladies Notes from the Kölner R meeting, 14 October 2016 anytime 0.0.4: New features and https://www.r-bloggers.com/building-barplots-with-error-bars/ fixes 2016-13 ‘DOM’ Version 0.3 Building a package automatically The new R Graph Gallery Network Analysis Part 3 Exercises Annotated Facets with ggplot2 Paper published: mlr - Machine Learning in R a grim knight [cont’d] Other sites Jobs for R-users SAS blogs Building Barplots with Error Bars August 17, 2015By Chris Wetherill (This article was first published on DataScience+, and kindly contributed to R-bloggers) Bar charts are a pretty common way to represent data visually, but constructing them isn't always the most intuitive thing in the world. One way that we can construct these graphs is using R's default packages. Barplots using base R Let's start by viewing our dataframe: here we will be finding the mean miles per gallon by number of cylinders and number of gears. View(mtcars) We begin by aggregating our data by cylinders and gears and specify that we want to return the mean, standard deviation, and number of observations for each group: myData <- aggregate(mtcars$mpg, by = list(cyl = mtcars$cyl, gears = mtcars$gear), FUN = function(x) c(mean = mean(x), sd = sd(x), n = length(x))) After this, we'll need to do a little manipulation since the previous function returned matrices instead of vectors myData <- do.call(data.frame, myData) And now let's compute the standard error fo
Diet & Nutrition (28) Education (1) Evolution (35) Human Ecology (75) Infectious Disease (66) LaTeX (5) Primates (9) R (12) science (17) Social Network Analysis (17) Statistics (16) Teaching (10) Uncategorized (28) Meta Log in Entries RSS Comments http://monkeysuncle.stanford.edu/?p=485 RSS WordPress.org ← Latest Swine Flu Epidemic Curve for the United States Stanford Workshop in Biodemography → Plotting Error Bars in R August 24th, 2009 · 52 Comments · R One common frustration that I http://www.personality-project.org/r/html/error.bars.html have heard expressed about R is that there is no automatic way to plot error bars (whiskers really) on bar plots. I just encountered this issue revising a paper for submission and figured I'd share error bars my code. The following simple function will plot reasonable error bars on a bar plot. PLAIN TEXT R: error.bar <- function(x, y, upper, lower=upper, length=0.1,...){ if(length(x) != length(y) | length(y) !=length(lower) | length(lower) != length(upper)) stop("vectors must be same length") arrows(x,y+upper, x, y-lower, angle=90, code=3, length=length, ...) } Now let's use it. First, I'll create 5 means drawn from a Gaussian random variable with unit mean and variance. I error bars in want to point out another mild annoyance with the way that R handles bar plots, and how to fix it. By default, barplot() suppresses the X-axis. Not sure why. If you want the axis to show up with the same line style as the Y-axis, include the argument axis.lty=1, as below. By creating an object to hold your bar plot, you capture the midpoints of the bars along the abscissa that can later be used to plot the error bars. PLAIN TEXT R: y <- rnorm(500, mean=1) y <- matrix(y,100,5) y.means <- apply(y,2,mean) y.sd <- apply(y,2,sd) barx <- barplot(y.means, names.arg=1:5,ylim=c(0,1.5), col="blue", axis.lty=1, xlab="Replicates", ylab="Value (arbitrary units)") error.bar(barx,y.means, 1.96*y.sd/10) Now let's say we want to create the very common plot in reporting the results of scientific experiments: adjacent bars representing the treatment and the control with 95% confidence intervals on the estimates of the means. The trick here is to create a 2 x n matrix of your bar values, where each row holds the values to be compared (e.g., treatment vs. control, male vs. female, etc.). Let's look at our same Gaussian means but now compare them to a Gaussian r.v. with mean 1.1 and unit variance. PLAIN TEXT R: y1 <- rnorm(500, mean=1.1) y1 <- matrix(