similar to: Is this an artifact of using "which"?

Dear All, I do apologise if this question is out of place for this list but I've tried searching mailing lists and read "Introductory Statistics with R" by Peter Dalgaard, but couldn't find any hints on solving my question below: I have a data frame (d) of values which I will rank in decreasing order of "val". Each value belongs to a group, either 'A',

Dear all, I have a collection of 5000 entries which represent the evolutionary rates of 3 animals. I would like to show the differences between the rates of all 3 animals and have tried using the function parallel (from the lattice package) and pairs() function. The parallel function would have been perfect save for the large number of data (5000). The pairs() function doesn't show

Dear all, I want to read in 1000 files which contain varying number of columns. For example: file[1] contains 8 columns (mixture of characters and numbers) file[2] contains 16 columns etc I'm reading everything into one big data frame and when I try rbind, R returns an error of "Error in rbind(deparse.level, ...) : numbers of columns of arguments do not match" Below is my

Dear all, I am a newbie in R and need some help please. (I do apologise if my email is not as informative as it should be, I've tried to include the relevant details without overcrowding it with the rest of the code) I would like to sample (without replacement) Y objects based on the number of objects in X in 5 different bins. I'm having trouble because the list object in which

Hi all, I have the following x1<-paste("A", 1:6, sep = "") x2<- round(rgamma(6,2,1)) x3<-paste("B", 1:6, sep = "") x4<- round(rgamma(6,2,1)) data1 <- data.frame(x1,x2,x3,x4) I would like to get data2 <- c(A1=4, A2=1, A3=0,...) Is there any standard for such a case? Thank you very much in advance, Diego

I want to run deconvolution of a time series by an impulse or point-spread function through Wiener filter, regularized filter, Lucy-Richardson method, or any other approaches. I searched the CRAN website and the mailing list archive, but could not find any package for such a deconvolution analysis. Does anybody know an existing R function for deconvolution? TIA, Gang

Hello everyone, I get some data in the following format and I would like to combine them to form a dataframe. The data is like: cbcname1 = 0.1, cbcname2= 0.2, cbcname3=0.3,... name1, name2, name2 are just some random names. I would like to achieve sth like: (cbcname1=0.1, cbcname2=0.2, cbcname3=0.3,......) I am using the following codes do.call(cbind,

I've just found some behaviour which strikes me as odd, but I'm not sure whether it's a bug or a feature. If you don't mind, I'd like to explain via a couple of examples. Let x = 1:10. Then intuitively, to me at least, the command x[-integer(0)] should leave x untouched. However the actual output under R2.4.0 is integer(0). A slightly more involved example demonstrates

Hi, Can we get the code for calculating Oja median for multivariate data Thanks and Regards Rahul Agarwal Analyst Equities Quantitative Research UBS_ISC, Hyderabad On Net: 19 533 6363 [[alternative HTML version deleted]]

hi i am new to R language. I want to use aaMI package which calculates the amino acid mutual interaction for a given protein sequence. I had installed the package but when i run the program it gives me the error could not find function "aaMI". can anyone tell me what might be the problem.. -- View this message in context: http://www.nabble.com/aaMI-tp14915744p14915744.html Sent from

Hello, I have a survivor curve that shows account cancellations during the past 3.5 months.  Fortunately for our business, but unfortunately for my analysis, the survivor curve doesn't yet pass through 50%.  Is there a safe way to extend the survivor curve and estimate at what time we'll hit 50%? We started a new program 3.5 months ago, and I believe that this set of accounts behaves

Hi, I have to match names where names can be recorded with errors or additions. Now I am searching for a string search function which returns always the "closest" match. E.g. searching for "Washington" it should return only Washington but not Washington, D.C. But it also could be that the list contains only "Hamburg" but the record I am searching for is

Hello everybody,   I wonder if there is any swap function in R that does the following: x <- seq(1,10,12) x1 <- swap(x) x1 1 8 3 4 5 6 7 2 10 Thank you very much!   Amor __________________________________________________ Schutz gegen Massenmails. [[alternative HTML version deleted]]

Hello, I have this in excel Control 543_BU 123_AT 432_CU I want to be able to import to R so that it will read like this c<-c("543_BU","123_AT","432_CU") output: [1] "543_BU" "123_AT" "432_CU" This is just a short version. I have about 200000 rows and i need a simpler way instead of typing each one. thanks -- View this message

Hi There, There seems to be a way for calling R from .net. However, is there anyway for calling .net/C# code from R? Something similar to the RJava package for .net? Thanks a mil! Fayssal -- View this message in context: http://www.nabble.com/R-and-.net-C--tp22904715p22904715.html Sent from the R help mailing list archive at Nabble.com.

Hello all, could you tell wehere I can find information related to MCMC (Monte Carlo). And some examples about this topic. Thanks, Tania [[alternative HTML version deleted]]

Hi all, I'm looking to add a "density" smoother on top of a hist when Freq=T. In order to do this I can use the relation between count and density, but I would like to know if there is a way for me to predict it upfront. Here is an example: set.seed(242) z = rnorm(30) hist_z <- hist(z) hist_z$counts / hist_z$density # the relation is 15 # why is this 15 ?? # So I can now do:

Hi, the second step in my exercice is to calculate the sum of the amout for each class et not the frequency i have this vector x           y 1          100 2          1500 3          3250 4          6250 5          2000 6           450 i want to use the function table and cat to calculate the sum of the amount in each class [0-1000], [1000-3000],[> 3000] Thank you for your help??

Hi, I?ve written a script to run several multivariate statistical analysis automatically. As one result a biplot and screeplot is produced. Now I?d like to display the name of the inputdatset as part of the title of these graphics and I do not want to enter it each time I run the script. How can I extract the name of a dataset? An (shortened) extraction from the script: >Test1 <-

Hi all: I am having trouble dropping levels, got a few hints online?without success. Please consider the dataset below: ?I was under the inpression that subset(......drop=TRUE) would work but it doesn't library(ggplot2) ??? library(hmisc) x <- structure(list(first = c(38.2086, 43.1768, 43.146, 41.8044, 42.4232, 46.3646, 38.0813, 40.0745, 40.4889, 38.6246, 40.2826, 41.6056, 34.5353,