similar to: inconsistency in return value of peaks() {splus2R} (PR#13988)

In 2004 there was this R-Help posting from Jan 2004: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/33097.html R equivalent of Splus peaks() function? The peaks function there has worked well for me on a couple of projects, but some code using "peaks" failed today, which had worked fine in the past. I was looking for a peak in a test case that was a sine curve over one

Dear All, For function: rmvnorm{mvtnorm} in (library mvtnorm, not splus2R), if I generate 2 bivariate normal samples as follows: > rmvnorm(2,mean=rep(0,2),sigma=diag(2)) [,1] [,2] [1,] 2.0749459 1.4932752 [2,] -0.9886333 0.3832266 Where is the first sample, it is stored in the first row or the first column? Does this function store samples row-wise or column-wise? Thank

Hello, Thank you to all those great folks that have helped me in the past (especially Dennis Murphy). I have a new challenge. I often generate time-series data sets that look like the one below, with a variable ("Phase") which has a series of flat-topped peaks (sample data below with 5 "peaks"). I would like to calculate the phase value for each peak. It would be great to

http://r.789695.n4.nabble.com/file/n3312061/x_and_y_values.txt x_and_y_values.txt I have the absorbance values form HPLC Chromatogram. I need to find the peaks in datapoints and area under each peak. I am not sure how how to find the peaks, I tried couple of libraries and peak function but they return me a list of values which when computed from area turns out to be huge nnumbers. In the

Hi, I apologize for not following the posting rules? Here is the text from my previous post: "I started evaluating the 'Peaks' package a couple of months ago and found it to be quite useful. Getting back to it last week I had to set up my R environment due to hardware changes again. The Peaks package loads with no problem. After successfully reinstalling all packages (RedHat

Hello all, I have some histograms of amount of DNA in some cells (DU145 cells overexpressing Bax and Bcl-xL for those who wish to know). The histograms show not only two peaks as expected, but three, indicating that some cells have more than normal amounts of DNA. I am interested in knowing how much of the cell populations are in each peak as well as between. I am not really sure how to go

Hello list! I would like to fit a distribution to each of the peaks in a histogram, such as this: http://photos1.blogger.com/blogger/7029/2724/1600/DU145-Bax3-Bcl-xL.png . The peaks are identified using Petr Pikal peaks function ( http://finzi.psych.upenn.edu/R/Rhelp02a/archive/33097.html), but after that I am quite stuck. Any idea as to how I can: Fit a distribution to each peak Integrate the

Greetings all, I am using Fourier analysis to search for periodicities in IP network traffic by generating periodograms and then visually examining them for large, distinct peaks. However, in many cases it is not readily apparent where there are periodicities. I have no experience with discrete maths so I've come up against a block here: How do I define what the "noise floor"

I've been asked in private, (and am replying BCC to the asker), >> I saw your post on the R-help archives page about the possibility of >> porting a function from S-Plus called peaks() to R. I am looking for >> some way to locate peaks in a simple x,y data set, and thought that R >> might be the way to go. "of course" it is the way to go, don't get

Hello all I'm trying to do a simulation that involves identifying the minimum point between two peaks of a (usually) bimodal distribution. I can do this easily if there are only two peaks: CnBdens<-density(Ys/Xs) #probability density function for ratio of Ys to Xs for(p in 1:512) ifelse(CnBdens$y[p]>CnBdens$y[p-1],peak1<-p,break) #identifies first peak in probability

Hi R-fellows, I am trying to simulate a multivariate correlated sample via the Gaussian copula method. One variable is a binary variable, that should be autocorrelated. The autocorrelation should be rho = 0.2. Furthermore, the overall probability to get either outcome of the binary variable should be 0.5. Below you can see the R code (I use for simplicity a diagonal matrix in rmvnorm even if it

This is probably a trivial question, but I did not see it in the FAQs. Perhaps I simply need to be pointed to some URLs. In any case, I have some sampled data that has peaks, and I would like to extract peak height, location, and width, according to almost any statistical model. Gaussian would seem to make the most sense. There can also be "steps" increasing the amplitude over a

> ?agrep > Results in: No documentation for 'agrep ' in specified packages and libraries: you could try 'help.search("agrep ")' There is white space after agrep, that ? doesn't ignore. --please do not edit the information below-- Version: platform = i486-pc-linux-gnu arch = i486 os = linux-gnu system = i486, linux-gnu status = major = 2 minor =

I am looking for a function in R that is equivalent to the function "peaks" in Splus. This function gives the local maxima of a vector. I do not find it using the help and I 'm not sure it exits in R. Could someone help me ? Thanks, Herve Cardot -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

If there something available in R that has the functionality of the S-PLUS peaks() function? Thanks, Mirka

If save() fails because an object is not found, it should not overwrite an existing file. > a <- 1:9 > save(a, file = "a.rda") > rm(a) > load("a.rda") > a [1] 1 2 3 4 5 6 7 8 9 > rm(a) > save(a, file = "a.rda") Error in save(a, file = "a.rda") : object 'a' not found > load("a.rda") Error in

I've run into a scoping problem in R. I'm calling a function that * creates a formula * calculates a weight vector * calls lm with that formula and weights This fails. Here's a simplified reproduce example: # f works, g doesn't, h is a workaround rm(w) data <- data.frame(y=runif(20), x=runif(20), z=runif(20)) f <- function(k){ w <- data$z^k coef(lm(y~x, data

On 15 Oct 2002 at 18:09, Rieckermann Joerg wrote: > Dear Petr, > > I have been fooling around with this peaks function of yours/Ripley > and don't seem to get the main idea. > > What I intend to do is locate the peaks in a vector and later use thes > values. So the position of the peaks would be of great interest. > > My questions: > * Why does peak retunr a

Hi, all. Does anyone have an equivalent to Splus "peaks", which finds local maxima (with locality defined by a parameter "span")? I thought I'd check whether anyone has done it already before trying to put something together myself. Thanks, Matt Wiener -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hi R Users, Sorry if its out of topic. I would like to ask you about twin peaks - normal distribution. How R can handle it, any example to explain it in R. Thanks, regards, Unung