8:37 p.m.
Dear Developers,
Today I'd found a expression that I think isn't marked for translation. In a
OLS model window
choose File -> View as equation. The expression is "R-squared".
Best regards,
--
Henrique C. de Andrade
Doutorando em Economia Aplicada
Universidade Federal do Rio Grande do Sul
www.ufrgs.br/ppge
4:22 p.m.
New subject: BIC in ols
I am a little confused about the behaviour of AIC, HQC and BIC in
the standard ols command. I estimated some AR(p) models with different p
values with the 'ols' and the 'var' commands. The last with only one
variable which I know it is not designed for, but ... it works.
I see that the BIC obtained through both commands:
1- are not the same. I suspect may be some different scale factor in the
formulas but this should not be a problem.
2- leads to different specification, and yes, this is a problem. In
concrete, the BIC in ols would suggest the maximum number of lags, and
this does not seem to be correct if one have a look at the correlograms.
<script>
#AR(p) selection
function void ARF(series y, scalar maxorder)
matrix C=zeros(maxorder,3)
loop i=1..maxorder
var i y
genr C[i,1]=i
genr C[i,2]=$bic
list xlist = lags(i,y)
ols y 0 xlist
genr C[i,3]=$bic
end loop
C
end function
open data9-3.gdt
Y=diff(reskwh)
Y=sdiff(Y)
ARF(Y, 30)
</script>
[Linux Ubuntu 10.04 and gretl from CVS build date 2010-05-14]
--
Ignacio Diaz-Emparanza
DEPARTAMENTO DE ECONOMÍA APLICADA III (ECONOMETRÍA Y ESTADÍSTICA)
UPV/EHU
Avda. Lehendakari Aguirre, 83 | 48015 BILBAO
T.: +34 946013732 | F.: +34 946013754
www.ea3.ehu.es
5:14 p.m.
New subject: BIC in ols
On Thu, 27 May 2010, Ignacio Diaz-Emparanza wrote:
I am a little confused about the behaviour of AIC, HQC and BIC in
the standard ols command. I estimated some AR(p) models with different p
values with the 'ols' and the 'var' commands. The last with only one
variable which I know it is not designed for, but ... it works.
I see that the BIC obtained through both commands:
1- are not the same. I suspect may be some different scale factor in the
formulas but this should not be a problem.
Well spotted: in fact, when we compute the information criteria for VARs,
we use as the sample size the total data span (T), whereas for OLS we use
the number of datapoints actually used (T-number of lags). To see what
I'm talking about, try this little modification of your script:
<script>
#AR(p) selection
function void ARF(series y, scalar maxorder)
matrix C=zeros(maxorder,3)
loop i=1..maxorder
var i y --quiet
genr C[i,1]=i
genr C[i,2] = $bic
list xlist = lags(i,y)
ols y 0 xlist --quiet
genr C[i,3]= $bic / $T
end loop
C
end function
open data9-3.gdt
Y=diff(reskwh)
Y=sdiff(Y)
ARF(Y, 30)
</script>
2- leads to different specification, and yes, this is a problem. In
concrete, the BIC in ols would suggest the maximum number of lags, and
this does not seem to be correct if one have a look at the correlograms.
I agree it is a problem, but only to an extent: information criteria are
likelihood-based, and as such have asymptotic validity only. Allin thinks
that this could be the outcome of the policy of "trying to agree as much as
possible with other packages". If package A uses T and package B uses T-k,
well... the issue "internal consistency vs agreement with others" has been
debated to death and I really wouldn't know what to say.
Riccardo (Jack) Lucchetti
Dipartimento di Economia
Università Politecnica delle Marche
r.lucchetti(a)univpm.it
http://www.econ.univpm.it/lucchetti
5322
days inactive
5324
days old
2 comments
3 participants
participants (3)
-
Henrique Andrade -
Ignacio Diaz-Emparanza -
Riccardo (Jack) Lucchetti