similar to: Plotting very skewed data in barplot

I am trying to make a barplot with a broken axis using gap.barplot (in the indispensable plotrix package). This works well when the data is a vector: > twogrp<-c(rnorm(10)+4,rnorm(10)+20) > gap.barplot(twogrp,gap=c(8,16),xlab="Index",ytics=c(3,6,17,20),ylab="Group values",main="Barplot with gap") But when the data is an array (for a bar plot with multiple

Hello R'ers In previous versions of R (I now use 2.0.1) there was a package "bootstrap". I wrote some programs that depended heavily on this package but can unfortunately not find it on Cran today. Is there any way to find an older version of "bootstrap" and use it with the new version - 2.0.1? Sincerely Fredrik Lundgren Norrk??ping Sweden

Hello fellow R'ers, I have a simple calculation with a very large data set being generated (34.9 million values) on a somewhat unreliable XP box that will likely take ~ 74hrs. I wanted to know if there is a way to have my script automatically "save.image()" throughout the calculation in case of a crash. This could be on the basis of output generated or time elapsed. I checked the

Dear all: I have the following codes: Xdata<-c(2,3,8,9,10) Ydata<-1:5 gap.plot(Xdata, Ydata,gap=c(5,6),gap.axis="x",type="o") However, the type='o' seems only work on the first part of gap plot, the second half of the plot always just points, you can not add lines on that part, any help will be highly appreciated. I would like to have these two parts of

Dear R'ers, I use R version 2.10.1 (2009-12-14) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 and Mac OS X 10.5.8 and use to work in different dierectories. When doing so it appears as if the .Rhistory file from my home directory always is loaded even when I want the file from the directory where I'm working. I have tried to fix it through preferences

What I want to do is create a horizontal bar chart matrix for a set of data that have a common set of variables (e.g., % of As, % of Bs, etc.) listed on the Y-axis and groups (e.g., Classes) on the X-axis. The key here is that the bars for each individual class plot are "centered" rather than left or right-justified. In archaeology plots similar to this are called battleship curves or

Hello R'ers, I have a large set of data which has many y samples for each unit x. The data might look like: Seconds Response_time ---------- ---------------- 0 0.150 0 0.202 0 0.065 1 0.110 1 0.280 2 0.230 2 0.156 3 0.070 3 0.185 3 0.255 3 0.311 3 0.120 4 .... and so on When I do a basic plot with type=l or the default of points it obviously plots every

Dear R- list, I have some data to present with histograms. Therefore I used hist(...). I have few values with almost 80% of the frequencies (totaly 800) and some other values with low frequencies ( totaly 5 -10 ) that I want to emphasize. Therefore I want to "cut" the y-axis on 100, but I don't know how to deal with this. Thanks in advance, Michael Graber

I just found out that this message got scrambled with other threads, so I trying to re-send... Dear R'ers, I am trying to build a composite plot (with several plots in one figure). I have tried, but I cannot use facetting, as I need to customize each plot using grid. Since all the plots are the same (with different data, but same layout and categories), I would like to have only one legend,

Dear R'ers, I've scanned available documentation and the web, but I can't seem to figure out where I've gone wrong in adding numbers to scale the colorkey in levelplot (Lattice package). For example, levelplot(z, contour=T, labels=T, cuts=10, region=T,... colorkey=list(T, space="bottom", at=seq(0,2000,length=20),

Hi R'ers- I'll believe this can be done and that I'm just not clever enough. I'm trying to do this in 2D, but i think we have the same problem in 1D. #Assume you have some 1D time series d<-rnorm(100) #Say you want to get the average over intervals of different lengths #like from time 10 to time 24 and time 51 to time 86 starts <- c(10,51) ends <- c(24,86) results

Dear R'ers, I am wondering what is the smallest geographicterritorial unit available for formatting in Canada. Provinces? I know that in the US it is the county so that I can color US counties any way I want, for example: ### Example for coloring US counties ### Creating an ARTIFICIAL criterion for coloring US counties: library(maps) allcounties<-data.frame(county=map('county',

Dear R'ers, I was not able to change the font size of axis label using cex.   Tried R archive solutions, though not solved yet. My code is: print(dotplot(ss ~ Response | bb*sr*tr, pp, pch = 21,               strip = FALSE, strip.left = TRUE,layout = c(3,8),               scales = list(y = list(relation = "free")), cex.axis=.2, cex.label=.2,                       ylab =

Dear R'ers, I have trouble installing tikzDevice in Ubuntu.  When I use install.packages("tikzDevice"), it gives error message: ERROR: dependency ‘filehash’ is not available for package ‘tikzDevice’ * removing ‘/usr/local/lib/R/site-library/tikzDevice’ Then I tried filehash installation, I get the message: "package ‘filehash’ is not available (for R version 2.13.1)"

Dear R'ers, I have a very large data frame (over 4000 rows and 2,500 columns). My task is very simple - I have to replace all NAs with a zero. My code works fine on smaller data frames - but I have to deal with a huge one and there are many NAs in each column. R runs out of memory on me ("Reached total allocation of 1535Mb: see help(memory.size)"). Is there any other, more efficient

R'ers: I was hoping I could get some direction on this. I have a dataset of the form: Y1,Y2,...,YM = f(X1,X2,...,XN), where N is >>> M The response data (Y1,Y2,...,YM) is frequency data, such that the sum of all Yi = 1.0. Both Xj and Yi are continuous variables. I'm trying to figure out the best approach(es) to solving for the model f() -- any ideas? I could solve

Dear R'ers, I am fighting with a problem that is driving me crazy. I use "lm" in my user-defined function, but it seems to be looking for weights outside of my function's environment: ### Generating example data: x<-data.frame(y=rnorm(100,0,1),a=rnorm(100,1,1),b=rnorm(100,2,1)) myweights<-runif(100) data.for.regression<-x[1:3] ### Creating function

Greetings, fellow R'ers. How can I get this frame in R: ID meas ID.1 meas.1 1 1.1 3 1.2 2 2.1 4 2.2 to look like this (stacking): ID meas 1 1.1 2 2.1 3 1.2 4 2.2 It's not really the reshape function (or is it?) because we can consider the additional columns, viz., ID.1 and meas.1, as independent of ID and meas so it is basically a stacking

Hello, dear R'ers, index<-expand.grid(1:7,1:4,1:4,1:4,1:4,1:4,1:4,1:4,1:4,1:4,1:4) In this case, dim(index) is 7,340,032 (!) and 11. I realize it's huge. Then, I am trying to get rid of the undesired combinations of columns. They should not contain identical values in any 2 columns. Also if column 1 has a value of 5, there should be no 2 in any other column, if column 1 has a value

Dear R'ers, When I used S-plus i wrote a small program for a Mantel-Haenszel test for trend (I think it worked). Unfortunately I can't get it working in R. It appears as if my use of 'el' is the problem but I can't sort it out. Error in apply(array, c(, 2, 3), function(el) el * 1:s) : argument is missing, with no default Further down in the program I use 'el'