similar to: Error bars using data from bugs()

Hello, i use ggplot to plot some measures including CIs as horizontal errorbars. I get an error when the scale limits are narrower than the boundaries of the error bar and hence the CIs are not plotted. library(ggplot2) df <- data.frame(resp=c(1,2), k=c(1,2), se=c(1,2)) ggplot(df, aes(resp,y=k)) + geom_point() + geom_errorbarh(aes(xmax = resp + se, xmin = resp - se)) +

I want to obtain a barplot with error bars and a mean plot with error bars with other scale on the left in the same plot. I need help to get the same x axis (centered) when overlay two plots (barplot2 and plotCI (or lineplot.CI) with errorbars), with diferent y axis. I use par (new=T) but the X axis names are not centered with the figures. An example of my data:

Hi everyone, I'm very new to R and I like to learn a lot... actually I have a little problem concerning errorbars with xyplot. My data look like run target hemi x 1 1 Nichts Links 0.0007743240 2 2 Nichts Links -0.0008153365 3 1 Target Links -0.0015825950 4 2 Target Links 0.0088743785 5 1 Nichts Rechts 0.0015898995 6 2 Nichts Rechts -0.0011465190 7 1

Hi! Could you tell me how I can draw a graph with error bars? Sorry, I don't use R that often and I couldn't find it easily in the documentation. TIA -- myriam -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

Dear R-list I'm working with with geom_errorbar; specifically I'm trying to reproduce the example Hadley Wickham have on http://had.co.nz/ggplot2/geom_errorbar.html (all in the button of the page) where he makes an nice plot with errorbars and then draw lines between the points. What confuses me is the 'limits' he defines for the errorbars from the se variable. First he creates

Hi there, Using the code below I generated the plot attached. The error bars at day = -3, 2 and 4 appear are larger. I was thinking if there is a way I could make all the error bars to be of the same size. I don't know if that makes sense. If not, then, is there a way I can plot only these 3 largest errorbars at day axis (day = -3, 2 and 4) and omit the rest of the error bars. I would be much

Hi, Does anyone know how i get rid of the marker point in my error bars? The bars default function (it seems) is to have a dot in the middle of the bar, however i don't want it there. I am using Hmis to draw my errorbars and code is a<-as.vector(tapply(Sporangia,list(Host,Isolate),mean)) b<-as.vector(tapply(Sporangia,list(Host,Isolate),sd)) c<-sqrt(b/36)

Dear R-help members, I have several large data sets from certain simulations I did and now I want to plot the results nicely. I don't know anything about the size of the x and y values in advance. Plotting these values is not a problem. However, I want to add errorbars (errbar in the Hmisc package). 1) For this I'm factoring the data (xdata00 varies from 0 to max(xdata00)) xfactor00 =

Hi All: Greetings, and best wishes from the festive times here at Kolkata, India -- the time of Durga Puja celebrations. I seek your advice as I try plotting lines for my data. The problem: I have created a dataframe that looks like this (name: myFrame): lowest second third fourth highest significance INAS 0.107 0.115 0.123 0.115 0.166 0.000 MMA 0.091 0.107 0.115

Dear I am trying to find a 90% credible interval. I am using the following code. fit<-bugs( model.file=BUGScode, data=data, inits = list(geninits1,geninits2), parameters.to.save=keepers, n.chains=nchains, n.iter=runs, n.burnin=burn, n.thin=nthin, DIC= TRUE, bugs.directory="C:/Program Files/WINBUGS.14", \ ) But this is only giving 95%

Hi Everyone, I would like to add errorbars to a lattice plot. I already have a function that adds error bars to a plot (originally written by George Gilchrist) which takes as arguments the coordinates of the point and then the size of the error bar. I would like to integrate it in my lattice code but I seem to be stumped. I can put the error bar function inside my code as if it were a panel

The concrete problem is that I am refereeing a paper where a confidence interval is presented for the risk ratio and I do not find it credible. I show below my attempts to do this in R. The example is slightly changed from the authors'. I can obtain a confidence interval for the odds ratio from fisher.test of course === fisher.test example === > outcome <- matrix(c(500, 0, 500, 8),

I am trying to calculate Bayesian Credible Intervals for a proportion (disease prevalence values to be more specific) and am having trouble using R to do this. I am working with ncredint() function but have not had success with it. Please help! Example: Positive samples = 3 Total sampled = 10 Prevalence = 0.3 pvec <- seq(1,10,by=1) npost = dbinom(pvec,10,prob=0.3, log=FALSE) ncredint(pvec,

Hi all, I have the data from the posterior distribution for some parameter. I want to find the 95% credible interval. I think "t.test(data)" is only for the confidence interval. I did not fine function for the Bayesian credible interval. Could some one suggest me? Thanks [[alternative HTML version deleted]]

Dear all, I would like to plot credible interval for a function estimate in R. I would like to plot the credible intervals as shaded region using polygon function. Does anyone ever used that? I tried several times but I could not obtain the right figure. xis=sort(xi,decreasing=TRUE) plot(xi,fm) polygon(c(xi,xis),c(f05m,f95m)) The above piece of code produces something else. Many thanks in

This is a question of clarification. IN 2009 Higgins, Thompson and Spiegelhalter (J R Statist Soc A 172:137-159) gave WinBUGs code to get credible intervals from random effects meta analysis for the prediction interval of a new study. It appears that the predict.rma function creates approximate credible intervals (pending a function revision by the author) for that purpose. Is my assumption

I have a function "eval.delta" which does what I want but isn't very elegant. I have consulted the R documents, MASS, and S Programming. Is there a practical way to optimize the code? Thank you very much. Peter B. Function: eval.delta <- function(delta){ cat("VALUES\n") vlr <- NULL k <- 0 for(j in 1:length(delta)) if(delta[j] <= 2){

Hi all, I'm new to the list and am hoping to get some advice. I have a set of multivariate data and would like to find the densest part of the data cloud containing 95% of the data, like a 95% HPD credible volume. Is there any R code available to compute that? Thank you very much! Your help and patience are much appreciated. G.S. [[alternative HTML version deleted]]

Hi, I'm having problems using the 'width' aesthetic attribute for the geom_errorbarh. This is the same problem reported earlier here, but I'll try to write the problem more clearly: http://www.mail-archive.com/r-help@r-project.org/msg62371.html The problem I'm having is that, the 'width' attribute is supposed to set the height of the endpoints of the whiskers, and

Hi, If I have an R object UUU, where the second element is U2, based on "g" column of my.table my.table of UUU is: mmm ggg gindex map Info aaa123 U1 1 1 1 aaa124 U1 1 2 1 bbb1378 U2 2 1 1 bbb8888 U2 2 2 0 bbb1389 U2 2 3