[R] How to create list of objects?
Dimitris Rizopoulos
dimitris.rizopoulos at med.kuleuven.be
Tue Jun 6 16:44:57 CEST 2006
try something like:
lapply(f, summary)
sapply(f, function(x) AIC(logLik(x)))
Best,
Dimitris
----
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
Web: http://med.kuleuven.be/biostat/
http://www.student.kuleuven.be/~m0390867/dimitris.htm
----- Original Message -----
From: "Rainer M Krug" <rkrug at sun.ac.za>
To: "R help list" <r-help at stat.math.ethz.ch>
Sent: Tuesday, June 06, 2006 4:38 PM
Subject: [R] How to create list of objects?
> Hi
>
> I am doing several mle and want to store them in a list (or whatever
> is
> the right construct) to be able to analyse them later.
>
> at the moment I am doing:
>
> f <- list()
> f$IP <- mle(...)
> f$NE <- mle(...)
>
> but when I say:
>> summary(f)
> I get:
>
> Length Class Mode
> IP 0 mle list
> NE 0 mle list
>
> I don't get the output I would have, i.e. the one from
>> summary(f$IP)
> summary(f$IP)
> Maximum likelihood estimation
>
> Call:
> mle(minuslogl = IPNeglogPoisL, method = "L-BFGS-B", fixed = list(),
> control = list(maxit = 1e+08, factr = 1e-20))
>
> Coefficients:
> Estimate Std. Error
> a 1242.0185506 44.92341097
> b 0.8802538 0.01685811
>
> -2 log L: 145.3509
>
>
> What I want to do is something like:
>
> AICs <- AIC(logLik(f))
>
> and then have all the AICs in the vector AICs.
>
> It must be possible or is this again a namespace issue?
>
> Rainer
>
> --
> Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation
> Biology (UCT)
>
> Department of Conservation Ecology and Entomology
> University of Stellenbosch
> Matieland 7602
> South Africa
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
>
Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm
More information about the R-help
mailing list