Re: [R] Peak Over Threshold values

2010-06-02 Thread William Dunlap
> -Original Message- > From: Tonja Krueger [mailto:tonja.krue...@web.de] > Sent: Monday, May 31, 2010 2:03 AM > To: William Dunlap > Cc: r-help@r-project.org > Subject: RE: [R] Peak Over Threshold values > > Thanks a lot for your help. That's the time period

Re: [R] Peak Over Threshold values

2010-05-31 Thread Tonja Krueger
you in advance, Tonja -Ursprüngliche Nachricht- Von: William Dunlap Gesendet: 27.05.2010 22:13:21 An: "Hutchinson,David [PYR]" ,Tonja Krueger Betreff: RE: [R] Peak Over Threshold values >> -Original Message- >> From: r-help-boun...@r-project.org >&

Re: [R] Peak Over Threshold values

2010-05-27 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of William Dunlap > Sent: Thursday, May 27, 2010 12:24 PM > To: Hutchinson,David [PYR]; Tonja Krueger > Cc: r-help@r-project.org > Subject: Re: [R] Peak

Re: [R] Peak Over Threshold values

2010-05-27 Thread William Dunlap
l Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of > Hutchinson,David [PYR] > Sent: Thursday, May 27, 2010 10:41 AM > To: Tonja Krueger > Cc: r-help@r-project.org > Subject: Re: [R] Peak Over Threshold values > > Perh

Re: [R] Peak Over Threshold values

2010-05-27 Thread Hutchinson,David [PYR]
col = 'red' ) abline(h = threshold, lty = 2, col = 'red') -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Tonja Krueger Sent: Thursday, May 27, 2010 1:47 AM To: Vito Muggeo (UniPa); Clint Bowman Cc: r-help

[R] Peak Over Threshold values

2010-05-27 Thread Tonja Krueger
I’m sorry, but that‘s not exactly what I was looking for. I obviously didn’t explain properly: Within my dataframe (df) I would like to find POT values that are not linked. In my definition two maximum values are linked if walevel does not fall below a certain value (the lower th

Re: [R] Peak Over Threshold values

2010-05-26 Thread Clint Bowman
How about? hi.rle<-rle(walevel>5.79) lo.rle<-rle(walevel<5.36) plot(walevel) abline(h=5.8,col=2,lty=3) abline(h=5.35,col=3,lty=3) hi.lo.rle<-sort(c(cumsum(hi.rle$lengths),cumsum(lo.rle$lengths))) abline(v=hi.lo.rle) You can use the $values from the rle to sort things out. Probably want to igno

Re: [R] Peak Over Threshold values

2010-05-26 Thread Vito Muggeo (UniPa)
dear Tonja, By plotting your data plot(df) it seems to me that you are looking for a piecewise linear relationships. If this is the case, have a look to the package segmented. You have to specify or not the number and the starting values for the breakpoints library(segmented) olm<-lm(waleve

[R] Peak Over Threshold values

2010-05-26 Thread Tonja Krueger
Dear List I hope you can help me: I’ve got a dataframe (df) within which I am looking for Peak Over Threshold values as well as the length of the events. An event starts when walevel equals 5.8 and it should end when walevel equals the lower threshold value (5.35). I tried “clu

[R] Peak Over Threshold values/ length of the events

2010-05-25 Thread Tonja Krueger
Dear List I hope you can help me: I’ve got a dataframe (df) within which I am looking for Peak Over Threshold values as well as the length of the events. An event starts when walevel equals 5.8 and it should end when walevel equals the lower threshold value (5.35). I tried “clu