similar to: How to identify the two largest peaks in a trimodal distribution

Hi Elisa, I hope this is what you wanted. dat1<-read.csv("peaks.csv",sep=",") #Subset dat2<-dat1[1:5,] res1<-do.call(cbind,lapply(seq_len(nrow(dat2)),function(i) do.call(rbind,lapply(split(rbind(dat2[i,],dat2[-i,]),1:nrow(rbind(dat2[i,],dat2[-i,]))), function(x) {x1<-rbind(dat2[i,],x);

Hello, I need to calculate the correlation for all pairwise combinations in a very large matrix. I have 25,000 elements and need to calculate the pairwise correlation with a different set of 5,000 elements. I have written code that works, but it is extremely slow. At the current rate, it will take a few weeks to finish running. I'm looking for suggestions on performing the calculations

Hi in matlab I defined a function (double gamma, parameters at the end of this mail) as h(i)=((t/d1)^a1)*exp(-(t-d1)/b1)-c*((t/d2)^a2)*exp(-(t-d2)/b2); h=h/norm(h); I do know that norm() in matlab is equal to: sqrt(sum(x^2)) in R so in R I do it like: #function (double gamama) h <- expression((t/d1)^a1*exp(-(t-d1)/b1)-c*(t/d2)^a2*exp(-(t-d2)/b2)) # plot it t <- seq(0, 20000,

Hello everyone Not strictly an R question but close... hopefully someone will be able to help. I wish to compare the switchpoints in two switchpoint regressions. The switchpoints were estimated using the segmented library running in R, and I have standard errors for the estimates. I initially thought I could just bootstrap confidence intervals for the difference between the switchpoints,

Hi, is there any package/function that can tell if a numeric vector (continuous data) has a bimodal or trimodal distribution and caluclate the location of the corresponding modes? Thanks

Hello people I am doing some thinking about how to analyse data on dimorphic animals - where different individuals of the same species have rather different morphology. An example of this is that some male beetles have large horns and small wings, and rely on beating the other guys up to get access to mates, whereas others have smaller horns and larger wings, and rely on mobility to

Hi folks I wonder if anyone can help me. I want to run some simulations to see how big a sample size might be necessary to distinguish a curved bivariate relationship (e.g. something that might be best described by a quadratic model) from a relationship that is two straight lines with a sudden change in slope (e.g. something best described by a breakpoint regression). I am using

People I apologise for asking a general stats question, but I'm at a bit of a loss as to what to do following some hostile referees' comments. If I have a fully randomised blocked design, with only three blocks, should I treat block as a random or fixed effect? I have read comments about not treating block as a random effect if the number of blocks is less than 6 or 7: is this

Hi there The default option for saving graphics from R (1.9.1) on my Mac is as a pdf file. If I open the file in Acrobat reader it looks really good and crisp, and is obviously saved as vector graphics, since I can zoom in as much as I like and it continues to look really nice. If I import it into MS Word (from office 2000), or Textedit, however, it imports it as a bitmap and unless I save

> From: Christoph Lehmann > > Hi > > in matlab I defined a function (double gamma, parameters at > the end of > this mail) as > h(i)=((t/d1)^a1)*exp(-(t-d1)/b1)-c*((t/d2)^a2)*exp(-(t-d2)/b2); > h=h/norm(h); > > I do know that norm() in matlab is equal to: > > sqrt(sum(x^2)) > in R > so in R I do it like: > > #function (double

Hello Can anyone let me know how I might change the size of the numbers on the axes in a pairs plot? The normal cex.axis call doesn't do anything, cex.labels only seems to change the font size in the diagonal labels. Using R 1.8.0 on a Mac with OS X.3. Thanks Rob Knell

Dear R-helpers Can anyone direct me to a method for plotting what you might call a partial wireframe plot? I have two explanatory variables in a dataset which give me a significant interaction term when I fit a model. The two variables are correlated with each other to a moderate degree, and if I plot the predicted values from the model as a surface in a 3D wireframe plot there are some

what is the largest file size limitation in ext3 ? Thank you

Hi, If I use the eigen() function to find the eigenvalues of a matrix, how can I find the eigenvector corresponding to the largest eigen value? Thanks! [[alternative HTML version deleted]]

Hello all, I was wondering if there is a function in R that only computes the eigenvector corresponding to the largest/smallest eigenvalue of an arbitrary real matrix. Thanks Minh -- Living on Earth may be expensive, but it includes an annual free trip around the Sun.

We were chatting at work and a couple of my co-workers thought we (Substantial) might be the largest rails shop in Seattle. We''ve got 30 full-time Ruby on Rails developers. Anybody else nearby have more? - Mike -- You received this message because you are subscribed to the Google Groups "Ruby on Rails: Talk" group. To view this discussion on the web visit

Hi Jason, > Are there any Linux users out there using the likes of > RAID'ed-NBD, CODA or > Intermezzo for a similar effect? > The NBD (network block device) looks interesting, it allows > you to mount a remote raw partition - so you can effectively > RAID over the network. I'd recomend drbd over a nbd + raid solution: * drbd knows to read from local device only so

Hi One of my packages needs a look-up table of pre-calculated numbers in the data directory. I would like to have the matrix as large as possible. What is the largest size matrix that would be an acceptable datafile in an R package? [ The table is a square, upper triangular matrix consisting of logs of Stirling numbers calculated by Maple. As discussed on the List a few days ago (thanks

I am attempting to use seek() to navigate through a large binary file (3GB). I would like to move to a file position which would require my "where" value to be greater than 2*10^9. This doesn't appear to be working, however. Is there any way around this problem on a 32-bit OS? If not I plan to modify my code to move through the binary file in smaller steps. I was curious as to

For me, I've found that I can easily work with 1 GB datasets. This includes linear models and aggregations. Working with 5 GB becomes cumbersome. Anything over that, and R croaks. I'm using a dual quad core Dell with 48 GB of RAM. I'm wondering if there is anyone out there running jobs in the 100 GB range. If so, what does your hardware look like? --Nathan [[alternative HTML