search for: turnpoints

Happy New Year, r-helpers! I am using lowess to smooth a scatter plot, xx<-lowess(xinput,f=.04) #defaults for other args followed by turnpoints(xx$y) #defaults for other args I plot the smoothed result as well as turnpoints (using yy$tppos) on top of raw data plot. Result is exactly as expected, graphically. For another purpose, I calcuate the difference between turnpoints (representing time intervals between turnspoints in my applicaito...

Hello, I'm trying to get a few turnpoints for a financial time series. There is a function in pastecs that does that. However, I get a large number of turnpoints: library(pastecs) data(EuStockMarkets) dax <-EuStockMarkets[,1] plot(dax) turnp <-turnpoints(dax) summary(turnp) #gives 925 peaks/pits!!! How can specify to get only 3...

...on" but rather small random fluctuations. Unfortunately, it is not possible for me to run more simulations to smooth the obtained distribution. What I am doing at the moment is the following: 1) use density() using the default parameter to obtain the estimated density distribution 2) use turnpoints() (PASTECS package) and a probability limit 0.05 to determine the turnpoints and select only the maxima for further analysis of the density distribution obtained in (1). My question is: I have several independent datasets which I want to analyse - do I have to set the width in density() to a s...

...o everybody. I have a data of four columns - three of the columns represent date while the fourth is counts. Setting 4000 as my minimum point/value, I wish to find out the number of counts that are less or equal to 4000 and the dates they occur. I have installed pastecs and have been trying to use turnpoints function to do this. I have not been much lucky. I am a still learning. I have added part of my data here where y stands for year, m month, d day and lastly the count. I would be grateful to anyone who could show me the way the do this. Best regards Ogbos y m d count 93 02 07 3974.6 93 02 08 3...

Hello, I just notice this: > x <- c(1:4,0:5, 4, 11) > library(pastecs) Loading required package: boot > tp <- turnpoints(x) > extract(tp, no.tp = FALSE, peak = TRUE, pit = FALSE) [1] FALSE FALSE FALSE TRUE FALSE FALSE FALSE FALSE FALSE TRUE FALSE FALSE Warning message: arguments after the first two are ignored in: UseMethod("extract", e, n, ...) > extract(tp) [1] 0 0 0 1 -1 0 0 0 0...

Hi, Here's the problem I ran into: the gmp package has a method for apply() so it masks the base::apply function. With gmp installed, I tried to run the function turnpoints() from the pastecs package. It fails because it calls apply() internally, like this: apply(mymatrix,1,max,na.rm=TRUE) , but the code in the gmp package which sets up the operator overload for apply() strictly limits the arguments to the first three (a matrix, a dimension, and a function). I...

Tom Blackwell wrote: >... >I summarized this to myself as "computed subscripts need explicit >rounding in R, but not in S". Here's the sample code which gave >me different results with R than with Splus. I no longer have >Splus available, so I can't check it again. >look <- (10 * seq(14)) - 76 >chk.1 <- seq(1420)[ 10 * (73.1 + look) ] #

Hi I just installed the pastecs package and I am wondering: is there an english (or german) translation of the file pastecs.pdf? If not, is there an explanation somewhere of the object of type 'turnpoints' as a result of turnpoints(), especially the "info" field? Thanks, Rainer -- NEW EMAIL ADDRESS AND ADDRESS: Rainer.Krug at uct.ac.za RKrug at sun.ac.za WILL BE DISCONTINUED END OF MARCH Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation Biology (UCT) Leslie Hill Institut...

I have a vector with about 100.000 values representing a quite regular function (sinusoid like). I would like to find all local maxima of this function (should be about 4000). Is there a native routine for R? Thanks in advance Eryk. -- _|_ \|/ \|/ Eryk Witold Wolski tel :0049-(0)30-8413-1543 w w ?v? 'v? \'v'/ MPI Moleculare Genetik fax :0049-(0)30-8413-1139 |

...857.3 93 02 23 3848.3 93 02 24 3869.8 93 02 25 3898.1 93 02 26 3920.5 93 02 27 3936.7 93 02 28 3931.9 " con <- textConnection(txt) dat <- read.table(con, header = TRUE) close(con) dat$date <- as.Date(paste(dat$y, dat$m, dat$d), format = "%y %m %d") library(pastecs) tp <- turnpoints(dat$count) [[alternative HTML version deleted]]

Hi all a week ago, I posted a newbie question in data smoothing & maximum-extraction with R. I got quite a lot of response, but I'm still kinda stuck with it... I'll restate the problem : i got a datafile with 2400 measuerements (every 250msec) of a CO2 measurement device, capturing the breath of a subject. I uploaded such a sample here :

...counts. When I plot a graph of counts versus dates, the resultant signal shows a number of maximum and minimum points. These minimum points (turning points) are of interest to me. Reading these dates and counts off from the plot is difficult as I am dealing with a large data. I have been looking at turnpoints function in pastecs library but have not been able to figure out the appropriate commands that one can use to find the minima/maxima (turning points) or pits/peaks in a series. My data is of the form shown below where y stands for year, m month, d day and finally count. Is there a way I could find...

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

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

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