similar to: Cumulative Points and Confidence Interval Manipulation in barplot2

Hi all, I'm building a plot of the values in tmeant (below) against positions 1 to 5, using matplot. tmeant looks like this: case1 case2 pos1 861.8466 818.5909 pos2 961.2841 976.3466 pos3 878.6080 1262.8523 pos4 950.8011 1129.6080 pos5 968.1080 1063.3920 I also have lower (object tl) and upper (object tu) bounds on the confidence intervals as follows: tl: pos1

I'm using barplot2 to plot some data. Is there any way to determine the width of the bars in the generated plot? I know that barplot2 returns a list of the coordinates of the center of each bar, but since there is some white space between each bar, I don't know how to get the width of each bar. Jamie

Hi list, I'm using barplot2 form the gplots package to plot a few numbers (I want to add SD bars later). However, I would like the y-axis not to start from 0 but 500. When I add the parameters YLIM, something goes wrong. The graph is not 'cut off' at 500. Instead the bars seems to sink trough the bottom of the graph. Because its a little hard to explain, here is a self-containing

Hello, I would like to know if it's possible to modify the name of groups of bar because on my barplot2, I have 5 groups of bars and one of them is called "User Contributes" and when I save the plot "User contributes" is to big so I don't have it on my plot! Is it pssible to put the name vertically! Thanks! Sabine --------------------------------- ils, photos

Hello, I get a nice looking barplot using the barplot2 function in the gregmisc package: body2 <- barplot2(hh3, beside = TRUE, col = c("mistyrose", "lightcyan"), .... cex.names = 1.0, plot.ci = TRUE, ci.l = cil, ci.u = ciu, plot.grid = TRUE) box() However, obviously I lose the collors when converting from ps to a pdf (outside of R)

Dear R users, I am trying to build a barplot2 graph however I can't find a way of defining the scale for the x-axis. I would like to show in my x-axis only the numbers 0, 25, 50, 75 etc. (so far R is giving me a random scale hard to interpret and it doens't look nice...). Could anyone advise me on how to do this please, it would be a great help! Thank you. Below I show the code I have

Hi R-users, I inted to make multiple plots using for loop. The question is how can I adjust the left hand side margin of the plot according to the names.arg argument in barplot2. In every plot I have different annotations in the y axis and they vary in length. Now when I have fixed margins opar <- par(mar=c(3,15,0,2)... I get the same margins in all of the plots. That leaves lots of white

Dear All, I have the following data coming out from s <- with(final, summarize(norm, llist(gtt,fdiab), function(norm) { n <- sum(!is.na(norm)) s <- sum(norm, na.rm=T) binconf(s, n) }, type='matrix') ) ie gtt fdiab norm.norm norm.norm2 norm.norm3 18

Hello, I try to handle a simple bar-plot, but it turns out to be not as simple as I thought. 1) I have created a .dat-File, e.g. test.dat: DATA DATA-SEM 2.2 0.32 6.2 1.30 12.7 1.61 48.6 3.08 4.1 0.86 4.5 0.32 1.5 1.13 1.2 1.08 The first row is the data represented by bars. The second row deals with the Standard Error of Mean. The lines correspond to time-intervals of experiments. 2) I now

A few random observations: 1. Clang could do better with large but boring switches like this: http://embed.cs.utah.edu/embarrassing/jan_10/harvest/source/E8/E88C5111.shtml Performance of clang's output will be fine but this is a major code size lose. 2. Destruction of stupid loops is incomplete, sometimes due to phase ordering problems:

First, I love R and am grateful to be using this free and extremely high quality software. Recently I have been comparing two algorithms and naturally I programmed in R first. It is so slow that I can almost feel its pain. So I decided to do a comparison with Java. To draw 500,0000 truncated normal, Java program takes 2 second and R takes 72 seconds. Not a computer science major, I find it hard

> Umm, can you find one that isn't a popcount implementation? Ok. MMX psadbw instruction: http://embed.cs.utah.edu/embarrassing/jan_10/harvest/source/CE/CE3DA132.shtml Position of first set bit: http://embed.cs.utah.edu/embarrassing/jan_10/harvest/source/1F/1F4003C7.shtml Log2 floor: http://embed.cs.utah.edu/embarrassing/jan_10/harvest/source/83/837A80E9.shtml Pixel format

On Tue, Jan 26, 2010 at 12:36 PM, John Regehr <regehr at cs.utah.edu> wrote: > 2. > Sometimes not: > > http://embed.cs.utah.edu/embarrassing/jan_10/harvest/source/EC/ECC74C0C.shtml The primary issue here is that scalar evolution doesn't know how to deal with loops using "sle" for the exit condition. Shouldn't be too hard to fix now that we have overflow flags

>> Repetitive code with lots of bitwise operations is compiled by LLVM into >> much larger code than the other compilers: >> >> http://embed.cs.utah.edu/embarrassing/jan_10/harvest/source/ED/ED37DAF5.shtml >> http://embed.cs.utah.edu/embarrassing/jan_10/harvest/source/1F/1F4003C7.shtml >> >> Note that this is straight-line code, so LLVM's output will

On Tue, Jan 26, 2010 at 5:55 PM, John Regehr <regehr at cs.utah.edu> wrote: >>> Repetitive code with lots of bitwise operations is compiled by LLVM into >>> much larger code than the other compilers: >>> >>> >>> http://embed.cs.utah.edu/embarrassing/jan_10/harvest/source/ED/ED37DAF5.shtml >>> >>>

Hi there I'm trying to find the mle fo a multinomial model ->*L(N,h,S?x)*. There is only *N* I want to estimate, which is used in the number of successes for the last cell probability. These successes are given by: p^(N-x1-x2-...xi) All the other parameters (i.e. h and S) I know from somewhere else. Here is what I've tried to do so far for a imaginary data set:

Hi, I have written a function "harvest" and I would like to run the function for each value in a vector c(1:1000). The function returns 4 list objects (obj_1, obj_3, obj_3, obj_4) using the following code at the end of the function: return(list(obj_1 = obj_1, obj_2 = obj_2, obj_3 = obj_3, obj_4 = obj_4)). Since I am connecting with the web in the function and the connection sometimes

Hello everybody, In am trying to analyse a BACI experiment and I really want to do it with R (which I find really exciting). So, before moving on I though it would be a good idea to repeat some known experiments which are quite similar to my own. I tried to reproduce 2 published examples but without much success. The first one in particular is a published dataset analysed with SAS by

There's some interest in my "auto-simplifier", which is nice :), so let me explain a bit about it. On 19/01/11 19:35, Sandeep Patel wrote: > You've mentioned your auto-simplifier a few times now and curiosity is > getting the better of me. Can you explain it a bit more? On 20/01/11 00:32, Nuno Lopes wrote: > Just out of curiosity, what's this auto-simplifier?

Dear list, Given a linear mixed model (from lme4) I want to 1) first change the input dataset and then 2) change the model formula. I want this to happen in a function call; Please see below. Options 1) and 2) below work whereas 3) fails with the message > foo() Error in is.data.frame(data) : object 'beets2' not found Question: What is it one must to in case 3) to have R look