similar to: clusters in evd package

Hello, and please advise regarding any errors/omissions on my part. However, the documentation for R's tcplot function in the evd package appears to contain an error. I am using evd version 2.2-4 in R 2.14.0 with Windows XP. > data(portpirie) > mrlplot(portpirie) ## No Error > tlim <- c(3.6, 4.2) > tcplot(portpirie, tlim) ## Line from documentation Error in fpot(data, u[1],

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 ?clusters (?)? from ?evd package?, and varied r (see example) but it did not work for all events (again

Hi, I am using the clusters function in the evd package in order to determine storm events from a wave time series. So far I have the code working as I want it for wave height on its own but I would now like to include the period as well. The input data is in the form of: H t H t H t so every height measurement has a corresponding period. The storms are defined when the wave height exceeds a

hi, rgev function gives me random deviates and I have a data set which I am fitting to an EVD,IS there a way I can plot both observed and ideal evd on the same plot thankyou Rangesh

As I tried to draw random sample from Hüsler-Reiss density with the following code, library(evd) margins <- cbind(0, 1, seq(-0.5,0.5,0.1)) x <- rbvevd(101, dep = 1.7, model = "hr", mar1 = margins) [,1] [,2] [1,] -0.56662298 0.8448505 [2,] 0.55824918 1.0217279 [3,] 0.24741124 0.6684668 [4,] -0.31547985 0.1781680 [5,] 0.69466134 0.2299363 [6,]

Hello, I would like to sample from a Gumbell (minimum) distribution. I have installed package {evd} but the Gumbell functions there appear to refer to the maximum case. Unfortunately, setting the scale parameter negative does not appear to work. Is there a separate package for the Gumbell minimum? -- _____________________________ Rich Gwozdz Fire and Mountain Ecology Lab College of

Hi, I use the function fpot of packages evd. If I call the fit that I obtain "fit", I want to modify the name of the axis and the main title that is produced by plot(fit1). To do this, I want to create a new function where I would have the names modified. The problem is that I can't find the function that produce these plots... I tried to see in plot.uvevd but it doesn't seems

Hi, I am trying the following: library(ismev) library(evd) fit <- gev.fit(x,show=FALSE) ks.test(x,pgev,fit$mle[1],fit$mle[2],fit$mle[3]) but I am getting: Warning message: cannot compute correct p-values with ties in: ks.test(x, pgev, fit$mle[1], fit$mle[2], fit$mle[3]) where x is: [1] 239 38 1 43 22 1 5 9 15 6 1 9 156 25 3 100 6 [18] 5 100

I want to use the plot(model) function to generate Tukey-anscomb and Q-Q plots of a lm(). I manage to change all labels but the main one which apparently is neither main or sub. So far I have tried as par setting: cex (changes symbol size within the plot), cex.main (no effect), cex.sub (no effect) cex.lab (changes label size), cex.axis (changes axis label size). What I would like to change is the

hello anyone faced this error before? dtrace: processing aborted: Abort due to systemic unresponsiveness I used iosnoop -evD > iosnoop.log thanks. __________________________________________________________ Sent from Yahoo! Mail - a smarter inbox http://uk.mail.yahoo.com

I am finding real trouble in installing evd library in R for linux -- View this message in context: http://www.nabble.com/library-instal-tp19453453p19453453.html Sent from the R help mailing list archive at Nabble.com.

Hello everybody 1. If there is/ever will be a function to perform "duncan multiple range test" please inform me at once. 2. Is there a way to create a grouping as duncan does from TukeyHSD output? My experimental design contained 62 genotypes, so the pairwise comparison is not that usefull and clear to look at. How would I do that? Is there an other test that does give grouping as

Dear Rxperts, Running the following code: ======================================================= twlo=10; twhi=20; wt=154; vd=0.5; cl=0.046; tau=6; t=3; F=1; wtkg <- wt/2.2 # convert lbs to kg vd.pt <- wtkg * vd # compute weight-based vd (L) cl.pt <- wtkg * cl # compute CL (L/hr) k <- cl.pt/vd.pt # compute k (hr^-1) cave <-

Dear all: I m a newer on R.? I have some problem when I use?anova function.? I use anova function to get Type I SS results, but I also need to get Type III SS results.? However, in my code, there is some different between the result of Type I SS and Type III SS.? I don?t know why the ?seqe? factor disappeared in the result of Type III SS.? How can I do?? Here is my example and result.

Hello, I have a problem with ggplot and colors I used this function to draw somes lines and I want them to be all black (just to test) however, I dont get any black line in the figure. Do you have any experience with ggplot and manual colors?? ggplot(cmax, aes(cmax[,3], cmax[,6],colour=interaction(cmax[,1],cmax[,2]))) + geom_line() + geom_point() +

We have just published the 2002/2 issue of R News on http://cran.R-project.org/doc/Rnews where you can download the newsletter as PDF or Postscript file. It will propagate to the CRAN mirrors within a day or two. This issue focuses on the R 1.5.0 release, with key innovations in the standard add-on package 'ts' and the two new recommended packages 'grid' and 'lattice'.

We have just published the 2002/2 issue of R News on http://cran.R-project.org/doc/Rnews where you can download the newsletter as PDF or Postscript file. It will propagate to the CRAN mirrors within a day or two. This issue focuses on the R 1.5.0 release, with key innovations in the standard add-on package 'ts' and the two new recommended packages 'grid' and 'lattice'.

Dear useRs, With the following set of command i managed to create the copula density ================================================ >library(MASS) >library(evd) >X<-matrix(sample(1:3000), ncol=2) >U=cbind(rank(X[,1])/(nrow(X)+1),rank(X[,2])/(nrow(X)+1)) > mat1=kde2d(U[,1],U[,2],n=150) image(mat1$x,mat1$y,mat1$z,col= rev(heat.colors(1500)),xlab="",ylab="")

I'm trying to convert an S-Plus program to R.  Since I'm a SAS programmer I'm not facile is either S-Plus or R, so I need some help.  All I did was convert the underscores in S-Plus to the assignment operator <-.  Here are the first few lines of the S-Plus file:   sshc _ function(rc, nc, d, method, alpha=0.05, power=0.8,              tol=0.01, tol1=.0001, tol2=.005, cc=c(.1,2),

Hi, I am new to bioconductor, trying to install KEGGSOAP package, but got warnings() when installing and error message when trying to load the package, can anyone suggest what went wrong? many thanks John > source("http://bioconductor.org/biocLite.R") Bioconductor version 2.11 (BiocInstaller 1.8.3), ?biocLite for help > biocLite("KEGGSOAP") BioC_mirror: