similar to: AIC and BIC in mixed effects model

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 My ultimate goal is to best fit time series by comparing AICs and BICs (as in Bayesian) from arima() and nnet(). I looked at the arima.R source code, but I am afraid I do not understand it. What I only miss really is the number of parameters p, where: AIC = n*log(S/n) + 2*p with S the squared residuals and n the number of observations. Can I get p

I'm trying to do the following out of sample regression with autoregressive terms and additional x variables: y(t+1)=const+B(L)*y(t)+C(1)*x_1(t)...+C(K)*x_K(t) where: B(L) = lag polynom. for AR terms C(1..K) = are the coeffs. on K exogenous variables that have only 1 lag Question 1: ----------- Suppose I use arima to fit the model:

Een ingesloten tekst met niet-gespecificeerde tekenset is van het bericht gescrubt ... Naam: niet beschikbaar Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070803/79b6292b/attachment.pl

On 01/21/10 20:42, Jan Sjodin wrote: >>> Imo, a loop is simply a special kind of region, so a "filter" is perhaps the way to >>> go if you are interested in loops. Regions containing loops will have to be inspected >>> using the PST. >> >> Except loops that have multiple exits. they are not necessarily (single >> entry single exit) region, if

Is AIC or BIC available when using the nnet package? Thank you Paul Green

Hi All, How can extract AIC,BIC from a fitted Garch model? -- SUMANTA BASAK. [[alternative HTML version deleted]]

Hello I'm with a doubt using BIC and AIC. I want to know if both of then are a way to steem the best model to use. How i know which of then to choose? Talita Perciano Costa Leite Graduanda em Ci??ncia da Computa????o Universidade Federal de Alagoas - UFAL Departamento de Tecnologia da Informa????o - TCI Constru????o de Conhecimento por Agrupamento de Dados - CoCADa

Hi all, I obtained negative values for AIC and BIC criteria for a particular model that I have developped... I don't remember to have negative values for these crietria for others applications, so I am a little suprised... Could anyone tell me if something is wrong or his conclusion concerning my model? Best regards, Olivier.

Dear all, I am trying to model number of samples from a given series. The series are modelled according Gamma function. In order to estimate the # samples, I use BIC/AIC with MLE (computed from dgamma function). Here is the code I have. __BEGIN__ mlogl <- function( x_func, theta_func, samp) { # computing log_likelihood return( - sum(dgamma(samp, shape = x_func, scale=theta_func, log

Hello everybody, Please help with connecting the AIC and BIC numbers printed by summary.gls to the logLik number. 1. is the logLik number the true ML or density scaling constants have been omitted? 2. what is the formula for calculating the AIC and BIC from logLik (and how can I see it)? I tried printing summary.gls but it says object not found. Thank you very much. Stephen [[alternative

Hi all, First: I have a small line of code I'm applying to a variable which will be placed in a matrix table for latex output of accuracy measures: acc.aarima <- signif(accuracy(forecast(auto.arima(tix_ts, stepwise=FALSE), h=365)), digits=3). The time series referred to is univariate (daily counts from 12-10-2010 until 5-8-2010 (so not 2 full periods of data)), and I'm working on

Dear All 1) I'm fitting spatial CAR models using R2Winbugs and although everything seems to go reasonably well (or I think so) the next message appears from WINBUGS 1.4 window: gen.inits() Command #Bugs: gen.inits cannot be executed (is greyed out) The question is if this message means that something is wrong and the results are consequently wrong, or Can I assume it as a simple warning

Hi, I'm slowly working through Tsay's "Analysis of Financial Time Series" 3rd ed. ?I'm trying to replicate Table 2.1 on p.47, which gives PACF, AIC, and BIC for the monthly simple returns of the CRSP value-weighted index. The data: http://faculty.chicagobooth.edu/ruey.tsay/teaching/fts3/m-ibm3dx2608.txt > da <-

Hi all, I am having issues comparing models with lmer. As an example, when I run the code below the model summaries (AIC, BIC, loglik) differ between the summary() and anova() commands. Can anyone clear up what's wrong? Thank you! Darren Ward library(lme4) data(sleepstudy) fm1<-lmer(Reaction ~ Days + (1|Subject), sleepstudy) summary(fm1) fm2<-lmer(Reaction ~ Days +

R 2.3.0, all packages up to date Linux, SuSE 10.0 Hi I want to calculate AIC or BIC from several results from mle calculation. I found the AIC function, but it does not seem to work with objects of class mle - If I execute the following: ml1 <- mle(...) AIC(ml1) I get the following error messale: Error in logLik(object) : no applicable method for "logLik" Therefore I am using the

Hi all, Could anybody please help me understand AIC and BIC and especially why do they make sense? Furthermore, I am trying to devise a new metric related to the model selection in the financial asset management industry. As you know the industry uses Sharpe Ratio as the main performance benchmark, which is the annualized mean of returns divided by the annualized standard deviation of returns.

TL;DR This RFC is a joint effort by Intel and Saarland University to bring the divergence analysis of the Region Vectorizer [1,2,3,4,5] (dubbed the vectorization analysis of RV) to LLVM. The implementation is available on github for feedback [0]. The existing divergence analysis infrastructure in LLVM has conceptual limitations (structured control, SCEV based). The new analysis resolves bugs

Dear R-users, I use glm() to do logistic regression and use stepAIC() to do stepwise model selection. The common AIC value comes out is about 100, a good fit is as low as around 70. But for some model, the AIC went to extreme values like 1000. When I check the P-values, All the independent variables (about 30 of them) included in the equation are very significant, which is impossible, because we

Dear all, I want to compare some different models for a dataset by QQ plots and AIC. I get the following AICs: - linear model: 19759.66 - GAMLSS model: 18702.7 - linear model with lognormal response: -7862.182 The QQ plots show that the lognormal model fits better than the linear model, but still much worse than the GAMLSS. So, in my opinion, the AIC of the lognormal model should be between the

Hi all. We are getting our first PRI line to use with Asterisk and one of the technical specifications is about framing, structured or unstructured. The main difference about them is almost clear for me: http://ckp.made-it.com/g704.html says: "G.704 is the framing specification for G.703. A carrier can 'steal' a 64kbps time slot (TS0) from a 2.048 Mbps line and use this to