similar to: numerical accuracy, dumb question

First, many thanks to Frank Harrell for once again helping me out. This actually relates to the next point, which is my contribution to the 'why don't social scientists use R' discussion. I am a hybrid social scientist(child psychiatrist) who trained on SPSS. Many of my difficulties in coming to terms with R have been to do with trying to apply the logic underlying SPSS, with dire

Hi, may be simle question, but a do not find it anywhere. Is there same function like max() ,but giving more results. max() give 1number-maximum I need funcion what give p bigest number. many thanks -- View this message in context: http://www.nabble.com/bigest-part-of-vector-tp22188901p22188901.html Sent from the R help mailing list archive at Nabble.com.

Dear Contributors: I have the next matrix: "X" "Y" "Z" 1 2 526 2 5 723 3 10 110 4 7 1110 5 9 34 6 8 778 7 1 614 8 4 876 9 6 249 10 3 14 I want to order the matrix from bigest Z (1110) to lower Z (14). Then I want to asign a color scale vector from blue ( bigest Z) to orange

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

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, 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

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

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

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

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

Peter, I have been doing some more testing and I find that if I modify the libusb.c file to set a report size in libusb_get_report to 16 the error messages go away. This is the code as it currently stands: static int libusb_get_report(usb_dev_handle *udev, int ReportId, unsigned char *raw_buf, int ReportSize ) { int ReportSize_test = 16; unsigned char raw_buf_test[ReportSize_test +

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

I have just setup asterisk on a debian sarge box. I am running Asterisk 1.21 with AMP and chan_capi_cm 0.6.1 using a BT Speedway (AVM Fritz) ISDN card, connected to a BT ISDN2e line. Currently we have 6 extensions (SIP) configured all using CounterPath(Xten) eyebeam softphone. After many hours of Googling I have finally got it all setup and working. We can transfer calls internally and make and

I recently bought a Belkin F6H375 UPS. It connects to the computer via USB. Searching the Internet led me to believe that the megatec_usb driver in nut supports this UPS, however I've been unable to get it to work. I have this in ups.conf: [belkinusb] driver = megatec_usb port = auto desc = "Belkin USB" Here's what happens when I try to use it: ~ #

Hello all, It seems that the next improvement to the R Graph Gallery is categorization of the graphics, that way each graph will be easier to find. That step should be done *carefully* if we want to avoid the opposite side-effect : graph not reachable through the categories. That's why the wisdom of the R community is required. Graphics will be classified in : - categories -

The following single line of code shows what I am trying to do, and the problem I am having... barplot(c(101,102,103),ylim=c(100,103)) The 'xaxis' is missing, and the grey bars 'fall off' the plot area. This is generally ugly, and I would like to trim the bars (ideally they would have a ragged appearance to show that I am 'zooming in'). I can see why what I am trying to