similar to: error bars around a point

posting the question in r-help@ to get more feedback :-) thanks, Yun-Fang ----- Original Message ----- From: "Yun-Fang Juan" <yunfang at yahoo-inc.com> To: <s-news at wubios.wustl.edu> Sent: Tuesday, December 16, 2003 2:04 PM Subject: [S] plot stacked bar chart in R > Hi, > I am trying to plot a stacked bar chart in R but am not able to find the > documentation.

Hi, I try to use lm to fit a linear model with 600k rows and 70 attributes. But I can't even load the data into the R environment. The error message says the vector memory is used up. Is there anyone having experience with large datasets in R? (I bet) Please advise. thanks, Yun-Fang [[alternative HTML version deleted]]

I need some guidance on what is the best way to automate R I am aware of Rterm. First question: Is it true that R is not currently OLE accessible under Windows 2000? Second question: Is there an R command that echoes all R Console results to a run log system file? My first attempt to write a run log file is the example from the R Data Inport/Export documentation. > zz<-

To all I have some simple questions pertaining to weights used in regression. If the variability of the dependent variable (y) is a function of the magnitude of predictor variable (x), can the use of weights give an appropriate answer to the regression parameters and the std errors? Assume that y at x=1 and 6 has a standard deviation of 0.1 and at x=11 it is 0.4 Then according to a web page on

Any idea why I can no longer plot two graphs on the same graphics device using the par(fig=) parameter? A simpler par(mfrow=c(1,2)) does work, showing the two plots side-by-side, but I would like the first to be larger. This simple example fails: x<-c(1,1,NA,2,2,NA,3,3) y<-c(2,4,NA,3,5,NA,1,4) par(fig=c(0,2/3,0,1)) plot(x,y) par(fig=c(2/3,1,0,1)) qqnorm(x) When plotted, the last

Is there an eloquent solution to re-assign vector element values? I have a vector which contains chemical data, some of them are "flagged" as non-detected values by their negative values. I can find the statistics on the positive values in vector "v" simply by typing: >v<- c(5,5,-3,-3,7,8,10) > v[(v>0)] [1] 5 5 7 8 10 I can also convert to positive values by

The last two times I have originated message threads on R or Bioconductor I have received the message included below from someone named Patrick Connolly. Both times I was the originator of the message thread and used what I thought was a unique subject line that explained as best I could what my question was. Patrick seems to be implying that I am abusing the R and BioC help newsgroups in this

Hi, I have a two-way anova with unequal cell numbers that I want to analyze. The problem is I don't have individual observations of the data. I only have the sufficient statistics (mean, variance, # of observations) for each cell. Is there any existing function in S-plus that would allow me to do aov() with the sufficient statistics? The table is like G1 G2 G3 G5 G6 T1 T2 T3

How do I post text-string dates along the x-axis instead of the excel_date position values? I would like to plot a time-series of chemical values changing with time The first column is the internal date from excel. This should be the x-axis position of a graph of the plotted data. The second column is a text-string of date. This should be what is posted (instead of the excel_date). The third

I'm just reading Modern Applied Statistics with S, 4th Ed., Venables and Ribley I'm typing in their examples in both R and S+. I need insight in the difference in the class() statement shown in Chap. 2. Example from book: > names(powers.of.pi) <- -2:2 > powers.of.pi -2 -1 0 1 2 0.1013212 0.3183099 1.0000000 3.1415927 9.8696044 >

This is not a criticism. I'm just curious. Is there an effort to keep R comparable to S+? Or are the two languages diverging? I am doing what probably legions have done before me, and legions will after me...using R on examples from text books written with S+ code. Most of the time everything appears to be equivalent. And then there are amazing divergences in commands. For instance: S:

To All How does one get a list of functions from a loaded package so that one can then get the appropriate help for each of the functions. Currently my method is based on a lot of trial-and-error. Here's an example of what I mean... >From this forum I learn that an interesting package called "multtest" exists on Bioconductor. I then use R Console's "Packages" --

Dear R Users, I am using R on windows. how to use 'points' function to plot two curves with errbar I am doing like: x.val <- as.integer(names(co2mean)) errbar(x.val, co2mean, co2mean + co2sd, co2mean - co2sd, xaxt='n', col=1, xlab=NA,ylab=NA)# obs error bar lines(x.val, co2mean, col=1, lwd=2) errbar(x.val, co2tm3.month.mean, co2tm3.month.mean + co2sd.tm3, co2tm3.month.mean -

Hi, I'm struggling with errbar graphics. I'm trying to plot an x-y graph with correct labelling, however can't seem to get main and sub to show on my graph. They do work when I use title(main="," etc...., but this will make it look at lot messier,I'll have to blank out ylab=" " , and I need to try and get the titles to update automatically according to my

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 =

Hello R-help, I am using R version 2.15.1. I upgraded from R version 2.13 a few months back. Previously, I was able to plot error bars on an xy scatter plot using the errbar function: errbar(RAEthylene$TIME,RAEthylene$AVE,RAEthylene$AVE+RAEthylene$STD,RAEthylene$AVE-RAEthylene$STD,add = TRUE,lty=2,pch=17); Today, I went to update my plot. However, in R version 2.15.1 I get error code

Is there an easy way to make the error bars the same color as the points and lines they are plotted with. My example # fake data x=sample(1:10, 100, replace =T) y = rnorm(100) + runif(100) df=data.frame(x,y) # summarize data m = aggregate(df,list(x),mean) se = aggregate(df,list(x),sd)/sqrt(10) library(Hmisc) plot(x,y) errbar(m$x, m$y, m$y+1.96*se$y,

Dear R-users How can I get the eigenvalues out of an lda analysis? thanks a lot christoph -- Christoph Lehmann <christoph.lehmann at gmx.ch>

Hello! If I use the errbar-function and have a logarithmic scale on the x-axis, then the little horizontal bars at the end of the errbars (cap) disappear. What can I do? Thanks for helping! Ute

I''m trying to run mongrel on port 80 and switch to a different user/group from root. As root, this works: mongrel_rails start -p 3000 -u rails --group rails This: mongrel_rails start -p 80 -u rails --group rails dies with: --------------------------- ** Changing group to rails. ** Changing user to rails. ** Starting Mongrel listening at 0.0.0.0:80