lapply(rg2, function(x) x$par.ests) there is no slot residuals! The function gev() does return a S3-object with class attribute 'gev', see ?gev.
> > Dr. Pfaff: > > After using str; can you give an example on data extration > (e.g. for $par.ests and @residuals) > > > > ----- Original Message ---- > From: "Pfaff, Bernhard Dr." <bernhard_pf...@fra.invesco.com> > To: Peter Maclean <pmaclean2...@yahoo.com>; Dr. Bernhard Pfaff > <bernh...@pfaffikus.de> > Cc: "r-help@r-project.org" <r-help@r-project.org> > Sent: Wed, July 6, 2011 8:17:12 AM > Subject: AW: [R] BY GROUP in evir R package > > Hello Peter, > > str(rg2) > > us quite revealing for this; by() returns a list and hence > lapply() can be > employed, e.g.: > lapply(rg2, rlevel.gev, k.blocks = 5) > > By the same token, you can extract the relevant bits and > pieces and put them > together in a data.frame. > > Best, > Bernhard > > > -----Ursprüngliche Nachricht----- > > Von: r-help-boun...@r-project.org > > [mailto:r-help-boun...@r-project.org] Im Auftrag von Peter Maclean > > Gesendet: Mittwoch, 6. Juli 2011 09:25 > > An: Dr. Bernhard Pfaff > > Cc: r-help@r-project.org > > Betreff: Re: [R] BY GROUP in evir R package > > > > Dr. Pfaff: > > How do we pass the "by" results to "rlevel.gev" function to > > get the return level and also save the results (both > > rg2(par.ests and $par.ses) and rl) as.data.frame? > > > > #Grouped vector > > Gdata <- data.frame(n = rep(c(1,2,3), each = 100), y = rnorm(300)) > > library(evir) > > require(plyr) > > > > #Model for Grouped > > rg2<- by(Gdata,Gdata[,"n"], function(x) gev(x$y, 5, method = > > "BFGS", control =list(maxit = 500))) # rl <- rlevel.gev(rg2, > > k.blocks = 5, add = TRUE) > > > > > > > > > > ----- Original Message ---- > > From: Dr. Bernhard Pfaff <bernh...@pfaffikus.de> > > To: Peter Maclean <pmaclean2...@yahoo.com> > > Sent: Fri, June 3, 2011 2:45:28 PM > > Subject: Re: BY GROUP in evir R package > > > > Hello Peter, > > > > many thanks for your email. Well, as you might have guessed, > > there is also a > > function by() in R that does the same job. See help("by") for > > more information. > > > > Best, > > Bernhard > > > > Peter Maclean schrieb: > > > Hi, > > > I am new in R and I want to use your package for data > > analysis. I usually use > > >SAS. I have rainfall data for different points. Each point > > has 120 observations. > > >The rainfall data is in the first column (RAIN) and the > > categorical variable > > >that group the data is in the second column (GROUP). The > > data frame is > > >rain.data. How can I use the gev function to estimate all > > three parameters by > > >GROUP variable group? In SAS there is a by() function that > > estimate the model by > > >group. However, I would like to move to R. > > > With thanks, > > > Peter Maclean > > > Department of Economics > > > University of Dar -es- Salaam, Tanzania > > > > > > > > > > ______________________________________________ > > R-help@r-project.org mailing list > > https://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting guide > > http://www.R-project.org/posting-guide.html > > and provide commented, minimal, self-contained, reproducible code. > > > ***************************************************************** > Confidentiality Note: The information contained in this message, > and any attachments, may contain confidential and/or privileged > material. It is intended solely for the person(s) or entity to > which it is addressed. Any review, retransmission, dissemination, > or taking of any action in reliance upon this information by > persons or entities other than the intended recipient(s) is > prohibited. If you received this in error, please contact the > sender and delete the material from any computer. > ***************************************************************** > ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
