I think part of your problem with your chosen method of attack arises
from your misspelling of the logLik function. The other part arises
because you have not yet digested FAQ 7.20:
http://cran.r-project.org/doc/FAQ/R-FAQ.html#How-can-I-turn-a-string-into-a-variable_003f
I haven't have much success assembling things in eval(substitute(...))
constructions, so if I were using a loop, I would have gone this way:
> xlist <- list(x1,x2,x3)
> for (i in 1:3 )
+ {
+ aic.table$Log.lik[i] <- logLik(lm(y~xlist[[i]]))[1]
+ }
> aic.table
model.number Log.lik
1 1 -43.39382
2 2 -41.24109
3 3 -34.20999
--
David Winsemius
On Feb 25, 2009, at 6:58 PM, Mark Drever wrote:
hi folks,
I'm trying to build a table that contains information about a series
of General Linear Models in order to calculate Akaike weights and
other measures to compare all models in the series.
i have an issue with indexing models and extracting the information
(loglikehood, AIC's, etc.) that I need to compile them into the
table. Below is some sample code that illustrates my approach so far
and my problem. I realize that somehow i need to provide actual
object for the AIC or LogLik call, and I've tried a couple of
different ideas (e.g., as.name, expression), but have come up empty.
Thanks so much. Any help much appreciated.
Mark.
Sample code begins here....
x1 <- rnorm(10,1,1)
x2 <- rnorm(10,2,2)
x3 <- rnorm(10,3,3)
y <- x1*2 + x2*3 + x3*4 + rnorm(10,0,0.4)
model1 <- lm(y ~ x1)
model2 <- lm(y ~ x2)
model3 <- lm(y ~ x3)
aic.table <- data.frame(model.number = 1:3, Log.lik = NA)
for (i in 1:3)
{
aic.table$Log.lik[i] <- Loglik(paste("model",i,sep=""))[1]
}
______________________________________________
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.
______________________________________________
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.
