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 <bernhard at pfaffikus.de>
To: Peter Maclean <pmaclean2011 at 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
>
>?
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-bounces at r-project.org > [mailto:r-help-bounces at r-project.org] Im Auftrag von Peter Maclean > Gesendet: Mittwoch, 6. Juli 2011 09:25 > An: Dr. Bernhard Pfaff > Cc: r-help at 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 <bernhard at pfaffikus.de> > To: Peter Maclean <pmaclean2011 at 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 at 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 ...{{dropped:10}}
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_Pfaff at fra.invesco.com> To: Peter Maclean <pmaclean2011 at yahoo.com>; Dr. Bernhard Pfaff <bernhard at pfaffikus.de> Cc: "r-help at r-project.org" <r-help at 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-bounces at r-project.org > [mailto:r-help-bounces at r-project.org] Im Auftrag von Peter Maclean > Gesendet: Mittwoch, 6. Juli 2011 09:25 > An: Dr. Bernhard Pfaff > Cc: r-help at 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 <bernhard at pfaffikus.de> > To: Peter Maclean <pmaclean2011 at 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 at 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. ********************************************