similar to: Peak finding algorithm

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

I'm just getting started with R, so feel free to point me to the appropriate documentation if this is already answered somewhere (though I've been unable to find it myself). This does seem like a rather basic question. I want to fold a table into a matrix. The table is formatted like so: Column_Index Value 1 486 2 688 3 447 4 555 5

Hi! My Platform: Windows 2000 R version >1.5 Active Perl 5.6 Rtools at the start of the PATH. I am trying to build a package. The functions are written only in R (no C code). I am trying to proceed as described in R-exts.pdf But instead to get a build package i am getting (see below). At first glance it seems to me that i understand the output. At something similar i am looking most of my

Hi, I was wondering if anyone knows of a package that contains the ability for peak/valley detection. Here is an example of what I'm looking for, only problem is that it's written in Matlab. http://www.billauer.co.il/peakdet.html Thanks for any help in advance. -Joe

Dear R-gurus, We are planning an experiment to test if plants produced by selfing are less fit than those produced by outcrossing. We have plants from three different alpine valleys, picked randomly among all the possible valleys. In each valley, we have a number of individuals, also picked at random. seeds from this individuals were brought back to the green house and sawned. when they

I have a data frame with three columns, x, y, and a. I want to create a matrix from these values such that for matrix m: m[x,y] == a Obviously, I can go row by row through the data frame and insert the value a at the correct x,y location in the matrix. I can make that slightly more efficient (perhaps), by doing something like this: > for (each.x in unique(df$x)) m[each.x, df$y[df$x ==

I'm currently using the below function from some library (MASS?) for writing my data out to file. I'm using it instead of plain old "write" because it does buffering. The problem that I'm having is that the numbers are space padded, but I need true tab-delineated files. It looks like the spaces are coming from 'format', but I don't see an option for format 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 all, I'm looking for smooth monotone estimation packages, preferably using splines. I downloaded the 'cobs' package and intend to use it, but since it offers only quadratic splines based on L1 minimization, I'd like to compare its performance to that of a more 'mainstream' cubic-spline, L2-norm minimizing spline. Preferably a smoothing spline. Does anyone know of such

Hi I went through some extensive search to find suitable method (package, function) to fit multiple peaks. The best I found is ALS package but it requires rather complicated input structure probably resulting from GC-MS experimental data and seems to be an overkill to my problem. I have basically simple two column data frame with columns time and sig dput(temp) structure(list(time = c(33, 34,

Hi all, I think this question could be quite trivial, but I can?t find out the solution... How can you compute the statistic "mode" of a sample, in case it exists (as mode() returns the mode of an object)? I tried help.search("mode") but I couldn't find a clue... Any help would be much appreciated. Regards, Aurora

R : Copyright 2002, The R Development Core Team Version 1.5.1 (2002-06-17) > install.packages("O:/Rpack/mscalib.zip", .libPaths()[1], CRAN = NULL) updating HTML package descriptions Warning message: error -1 in extracting from zip file > install.packages("O:/Rpack/zzz.zip", .libPaths()[1], CRAN = NULL) updating HTML package descriptions Warning message: error -1 in

Hi, I'm trying to match peaks between chromatographic runs. I'm able to match peaks when they are chromatographed with the same method, but not when there are different methods are used and spectra comes in to play. While searching I found the ALS package which should be usefull for my application, but I couldn't figure it out. I made some dummy chroms with R, which mimic my actual

Hello! I have a dataset with the following two vectors: year<-c(1967,1968,1969,1970,1971,1972,1973,1974,1975,1976,1977,1978,1979,1980,1981,1982,1983,1984,1985,1986,1987,1988,1989,1990,1991,1992,1993,1994,1995,1996,1997,1998,1999,2000,2001,2002,2003,2004,2005,2006,2007,2008,2009)

Dear all, Another newbie just got attracted to this mailing list. I am a biologist currently working my way through R, had sort play around with python earlier this year. I have some data exhibiting periodicity ** my data consists of peaks and valleys, with peaks arising due to the presence of a repetitive structural unit,** with x being a reference grid (position along a chromosome) and y

Hi all Is there a function that can find the start and end position of peaks in a set of numbers. eg. x=c(14,15,12,11,12,13,14,15,16,15,14,13,12,11,14,12) y=somefunction(x) y 4 14 Thanks John [[alternative HTML version deleted]]

Grrr ... new trial with code here: http://pastebin.com/RjHNNG9J Maybe the amount of inline-code prevented posting? Hello, I am writing a simple peak detector and it works quite well ... however there's one special case below, that I can't get my head wrapped around ... the problem is in the "Deal with not fully qualified peaks at the sequence extremes" section, but I cannot

Dear All, I have a binomial response with one continuous predictor (d) and one factor (g) (8 levels dummy-coded). glm(resp~d*g, data, family=binomial) Y=b0+b1*X1+b2*X2 ... b7*X7 how can I get the inflection point per group, e.g., P(d)=.5 I would be grateful for any help. Thanks in advance, Ren?

Sorry to all that are angry about the form of my previous mail. I didn't realise what would happen :((. Here it is in (hopefully) plain text (if my mailer doesn't spoil it again): ############## Dear developers, I have a problem with some discrepancy between R 1.2.1 for Windows and R 1.2.2 (and less) for Linux. While trying to correct the wilcox.test (see my previous bug report) I

Hi Is there a function which determines the location, i.e., index of the all minimums or maximums of a numeric vector. Which.min(x) only finds the (first) of such. > x <- c(1:4,0:5, 4, 11) > x [1] 1 2 3 4 0 1 2 3 4 5 4 11 > which.min(x) [1] 5 > which.max(x) [1] 11 > but I need which.min(x) to be 5 11 which.max(x) to be 4 10