similar to: lowess + turnpoints = doubling integers?

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, ...) >

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 :

I am working on a plot and would be like to click on a few points and then have a line connect them. Could anyone help me with this or advise me in a direction that would suit this. I know I would be using ievent.wait in iplot but not sure about this. thank you. -- View this message in context: http://www.nabble.com/ievent.wait-tf3891095.html#a11030568 Sent from the R help mailing list archive

Hi, How can I get the '2+' into superscript in the following title: '[Ca2+]i onsets' I tried the command below, but it doesn't work. What am I missing? hist(X, main=expression(([Ca*]i^2+) 'onsets'), xlab = 'sec') -- View this message in context: http://r.789695.n4.nabble.com/superscript-characters-in-title-with-tp3006981p3006981.html Sent from the R

Hello, To me as a beginner, every problem looks big. Below is what I was asked to do as part of a code that will solve my problem. I have used read.table to read my data into R and assigned the column names with colnames(dat)<- ... . But to go from txt<-" to the last " at the bottom of the table below is what I am unable to do. Could anybody please tell me how to type or arrive at

Good morning once more. My problem of yesterday has been addressed. Having learned a few tricks from that, I wish to ask another question in connection with that. My data is a cosmic ray data consisting of dates and 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.

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

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

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

Good evening to 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

Dear S-PLUS/R users: Do you know any default function or a user contributed function that can draw an ellipse with given axes and origin? Thanks for any help. Paul. -- Romance, like alcohol, should be enjoyed, but should not be allowed to become necessary. -- Edgar Friedenberg -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list --

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 |

Dear all, I need to maximize the v: v= D' W D D is a column vector ( n , 1) W is a given matrix (n, n) subject to: sum D= 1 (BTW, n is less than 300) I´ve tried to use maxBFGS, as follows: ##################################### objectiveFunction<-function(x) { return(t(D)%*%W%*%D) } Amat<-diag(nrow(D)) Amat<-rbind((rep(-1, nrow(D))), Amat) bvec<-matrix( c(0), nrow(D)+1,

Hello again, First I want to thank all the people who answered my question about line width in graphs. I promise I will learn the 'par' help page by heart for the end of the month ! I now want to trace some ellipses to emphasize groups of data. I found how to trace circles with 'symbols()', but no ellipse. I'm planning on writing my own function based on

Jean-Marc Valin wrote: > The noise suppressor will only attempt to remove stationary noise, > such as thermal noise, fans, ... The AGC can indeed do strange > things in these cases, but it's been improved in svn (compared to > 1.2beta1). OK, then the problem is that I misunderstood the feature. I assumed that dynamic squelch was part of it, but it's really something more

Hi list!, this is my first mail, nice to meet you ;-) The first, sorry for my english, this is not good, i am from Spain :-( The second THANKS in advance ;-) I am writing because I have a problem with the virtualization of M$ Windows on Xen. All system works OK when I virtulize systems based on GNU/Linux. When I start my Windows Machine for make installation with the CD (iso). The Windows

How can I use R to record the time of a mouse click? -- Seth Roberts Associate Professor Psychology Department University of California Berkeley [[alternative HTML version deleted]]

Full_Name: Erik Lukac Version: R version 2.8.1 (2008-12-22) OS: Mac OS 10.5 Submission from: (NULL) (141.84.28.167) I tried to install JGR, but somehow it doesnt work. > install.packages("JGR") Warning in install.packages("JGR") : argument 'lib' is missing: using '/Users/erikl/R/i386-apple-darwin9.6.0-library/2.8' trying URL

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 30 turnpoints? Second question: the extract function

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