similar to: Question about creating error bars

Hi, I want to use hist with non-equi-spaced breaks, picked such that the fraction of the data points falling in the cells (defined by 'breaks') is roughly equal accross all cells. Is there such a function that will automatically try to determine the breaks to fullfill this requirement? Something like.. hist( x, br = magic_function_to_pick_breaks()) For example, if x looked like

Signif. codes: 0 *** 0.001 ** 0.01 * 0.05 . 0.1 1 Signif. codes: 0 <80><98>***<80><99> 0.001 <80><98>**<80><99> 0.01 <80><98>*<80> <99> 0.05 <80><98>.<80><99> 0.1 <80><98> <80><99> 1 Signif. codes: 0 *** 0.001 ** 0.01 * 0.05 . 0.1 1 hmm... they go away when I

This has been asked before, but all the answers are hidiously complex. The legend.text=TRUE option of barplot is almost exactly what I need, except I need a legend.placement='tl' (top left) option. This option would be in contrast to the default placement which we could call 'tr' (top right). Anyone know how to edit the barplot code to make this change? Could someone like me

I found that the edit command kills my linux/R/emacs environment sometimes. How should I report this bug?

Hi, I am plotting a distribution of (ordered) values as a barplot. I would like to label groups of bars together to highlight aspects of the distribution. The label for the group should be the range of values in those bars. As this is hard to describe, here is an example; x <- rlnorm(50)*2 barplot(sort(x,decreasing=T)) y <- quantile(x, seq(0, 1, 0.2)) y plot(diff(y)) That last

I read answers to questions including the words "tapply" and "weighted.mean", but I didn't understand either the problem (data) or the solution provided. Here is my question ... > dat[1:10,] GROUP VALUE FREQUENCY 1 2 2 78 2 2 3 40 3 2 4 16 4 2 5 3 5 2 6 1 6 2 8 1 7

Is there an R wiki? Looking at the huge amount of traffic on this list, I think wiki could be an exelet outlet for all the constructive enthusiasm here. I don't think it would be too hard to port the existing R documentation (the stuff you get with the ?) onto a wiki system, then users could add their own examples and comments. Having distinct web page style organization would be easier to

Hello, supposing that I have two or three clear categories for my data, lets say pet preferece across fish, cat, dog. Lets say most people rate their preference as being mostly one of the categories. I want to do pca on the data to see three 'groups' of people, one group for fish, one for cat and one for dog. I would like to see the odd person who likes both or all three in the

I want to add values onto the end of a vector, for example... x <- vector for (i in 1:5){ add_to_end_of_vector(i,x) } I just cant find the answer to this question! Sorry for such a basic question, I tried... x <- c() for (i in 1:5) x[length(x)] <- i but it didn't work.

Hi, Since I started to make some 'final' plots of my data I found that I have tons of questions related to 'the little things'. Rather than bother the list with all the questions (ahem), or search the archives for similar questions and translate the context, I would like to find a FAQ for plotting in particular (and R programming in general). I know for sure (searching the

Full_Name: Dan B Version: na OS: na Submission from: (NULL) (80.6.127.185) I can't log into the bug tracker (I can't find where to register / login). In this way I can't add the following context diff (hopefully in the right order) for my changes to the matrix.Rd... Hmm... I guess this should be a separate report anyway... The first diff explains how the dimnames list should

I store an id as a big number, could this be a problem? Should I convert to at string when I use read.table(... example id's 1001001001001 1001001001002 ... 1002001002005 Bigest is probably 1011001001001 Ta, Dan.

I use the following code to scan a (limited) parameter space of clustering strategies ... data <- read.table(... dataTranspose <- t(data) distMeth <- c("euclidean", "maximum", "manhattan", "canberra", "binary" ) clustMeth <- c("ward",

I am trying the following code... for(i in 1:bins){ print(cat(c(i, length(var1[var1==i])))) } Where var1 is a vector with zero or more values = {1, 2, ..., bins} (or something like that). The result I get is... 1 33NULL 2 28NULL 3 39NULL 4 27NULL 5 32NULL 6 30NULL 7 23NULL 8 16NULL 9 10NULL 10 15NULL 11 13NULL 12 7NULL 13 3NULL 14 7NULL 15 2NULL 16 0NULL 17 2NULL 18 0NULL 19 0NULL 20 3NULL

This may sound crazy but... I have data like this... > results.matrix [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 949 93 2 11 26 20 7 6 10 5 0 3 [2,] 1233 124 24 35 58 57 17 21 31 19 11 21 Which is the result of binning (summing) the response variables of an underlying (nearly) continious range of predictor

First, thanks to everyone who helped me get to grips with R in (x)emacs (I get confused easily). Special thanks to Stephen Eglen for continued support. My question is about non-linear binning, or density functions over distributions governed by a power law ... y ~ mu*x**lambda # In one of its forms # (can't find Pareto in the online help) Looking at the following

Martin Maechler said... > >>>>>> "lanceh" == lanceh <lanceh at ibm.net> writes: > > lanceh> Sorry if this is a redundent question but I searched the > lanceh> archives and although I found a couple of similar questions I > lanceh> saw no responses. I am running in Linux and when I try to > >You should have looked in the

Gives strange results. I get 'weird' dendrograms with canberra / binary distance metric and median / centroid cluster methods. Is this just my data? Dan

I found this great search source("http://www.statslab.cam.ac.uk/~djw1005/Stats/Interests/search.R") helpHTML() Has (or will) this become the standard search method? Is R 'Free Software'? The dependence on Java seems a bit of a pain for 'freeness'. Did the above make it into CRAN? Cheers, Dan.

Meta: This question is somewhat long and has two parts, I would be very happy for someone just to nudge me in the right direction with the manual / tutorial, as I am somewhat lost... 1) How do I fit a curve of the form "y = k1 * x^k2" ? I want to estimate values of k1 and k2 given the x/y data I have, and I can't work out how to get R to calculate and return their estimates. 2)