similar to: Custom Plot - means, SD & 5th-95th% (Plotmeans or Boxplot)?

Dear R users, I've had no joy finding a solution to this online or in any of my R books. Many thanks in advance for any help you can give. I'm seeking to output a data frame (or matrix - it doesn't matter which for my purposes) to a .txt file, but omit any row or column names. The data frame that I'm using doesn't actually have column or row names to start with as it has been

Hi all, I'm just getting started with R and I would appreciate some help. I'm having trouble creating a boxplot with whiskers at the 95th and 5th percentiles instead of at 1.5 * IQR. I have read the relevant documentation, and checked existing mails on this topic. I found a small modification that should work : https://stat.ethz.ch/pipermail/r-help/2001-November/016817.html and tried to

Hi, I want to plot means with standard deviations of Total Nitrogen (TN) across 4 stations (S1-S4) and over 3 years (2007-2009). I want this plot in one panel. I do not want medians (bwplot, boxplot). I have tried a few different packages and it seems that ggplots with plotmeans was the fastest (I am not extremely skilled in writing my own scripts). Unfortunately, there is no grouping

I have a 20x3 matrix as follows: > m <- replicate(3, matrix(rnorm(20),20,1)) I need to compute, say, 95th and 99th percentiles of each column such that the resulting matrix becomes 2x3 with each row representing the respective percentile. My "best effort" is to compute one column at a time as follows: > quantile(m[,1], c(0.95, 0.99)) To do the same for columns 2 and 3, I

R users, I've been searching around for a solution to this for a few days, any help would be appreciated. I have working R code (a custom function) that transforms 3D coordinate data into angles (between any three sets of coordinates). The code outputs a single scalar (an angle). For my research I will have several hundred .csv files containing coordinate data - each file is a set of

Hello everyone, I have tried to look for this everywhere and so far have no luck. I have a plotMeans(DV, IV1, IV2, error.bars="se") graph that plots my data (DV-continuous, IVs are factors, IV1 - two levels, IV2-four levels). I am trying to increase a scale of my y-axis (to be consistent with my other graphs), but unfortunately nothing works with "plotMeans" function, which

I know that simple line plots showing a point representing the mean + an error bar are difficult in R, but I am clearly missing something. I simply want to show how a the mean and 95% CI of a chemical concentration (y) vary over space (x) under four temperature regimes (T; a categorical factor). plotmeans, in the gplots library, should allow for such a figure, right? After all, it invokes lm()

Hello I'm having trouble being able to plot more than one set of data using the plotmeans function. Plotmeans is where you can plot means and standard errors easily. Here is part of my code plotmeans(Mmeans ~ Month, data=BW.hydat.05HG001, n.label=FALSE)#one set of data par(new=TRUE) plotmeans(Mmeans ~ Month, data=Mouth.hydat.05HH001, col="red")#the second set of data Basically

Hello, I am using library(gplots) to do something like data(state) x1 <- state.area/10000 x2 <- x1+round((rnorm(length(state.area),3,3))) plotmeans(x1 ~ state.region) Is it possible to plot x2 to x1 in the same graph, something like: linesmeans(x2 ~ state.region) Best wishes, Cheba [[alternative HTML version deleted]]

I can not seem to get a ylim argument into plotMeans() is there anyway to do this. thanks in advance Stephen -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis [[alternative

Dear R fans, I am preparing a lecture discussing the paper of Wilkinson and APA Task Force on Statistical Inference (1999, American Psychologist, 54, 594-604.) I'd like to demo their Figure 3B (p. 602) with R. That is to add plotMeans(...) with dots of BHH2::dotPlot(...) . Is there any elegant script to do the job? Thanks --------------- LI, Xiaoxu

Hi, anybody know why this not works for several plots ? When i set onefile=T the plots are stacked one about another , onefile=F only the first plot is shown in AllbusMeansPlots.svg. [h2 is a data.frame] ......hist and sapply works for several plots nice with RSvgDevice ! Maybe setting the title after apply is a problem, but until yet i didn't found a better solution ? library(RSvgDevice)

Hello, I want to create a plot whit plotMeans. How can I change the font? Thanks in advance. Teresa Boca __________________________________________________ [[alternative HTML version deleted]]

Hi, I often use plotmeans() from the gplots package to quickly visualize a pattern of change. I would like to be able to plot separate lines for different groups, but the function gives an error when a grouping variable is included in the formula argument. For instance, > require(gplots) > x <- data.frame(Score=rnorm(100), Time=rep(1:10, 10),

Dear All, I am trying to graph a simple scatter plot where the x axis is year and the y axis is a percentage (percentage of infant death). Instead of plotting the raw data, I want to plot summary statistics such as mean and median. Here is the problem: the value range of y is between 0 and 1, but since infant death is a rare event, the mean and median is very low (something like 5%), which shows

Hi everyone, i'm having a problem extracting objects out of functions i've created, so i can use them for further analysis. Here's a small example: # --------------------------- test <- function(i, j){ x <- i:j y <- i*j z <- i/j return(x,y,z) } # --------------------------- This returns the 3 objects as $x, $y and $z. I cannot, however, access these objects

When the binary search returns 0 (exact match), the target key will necessarily be at slot 0 of all nodes below the current one, so in this case the binary search is not needed because it will always return 0, and we waste time doing it, holding node locks for longer than necessary, etc. Below follow histograms with the times spent on the current approach of doing a binary search when the

Dear list, I am calculating the 95th percentile of a set of values with R and with SPSS In R: > normal200<-rnorm(200,0,1) > qnorm(0.95,mean=mean(normal200),sd=sd(normal200),lower.tail =TRUE) [1] 1.84191 In SPSS, if I use the same 200 values and select Analyze -> Descriptive Statistics -> Frequencies and under "Statistics", I type in '95' under Percentiles,

How do I add to a trellis plot the best fit line from a robust fit? I can use panel.lm to add a least squares fit, but there is no panel.rlm function. -- Alan S Barnett <asb at mail.nih.gov> NIMH/CBDB

Is there a way of getting a list of indices or something indicating which set of points lqs when doing a regression? The $bestone, seems to return only 2 pts. no matter how good or bad the fit. I checked the archive for the past few months, but didn't find anything. thanks in advance, Mark Hall -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help