比较模型及其对数转换版本的AIC


17

我的问题的实质是:

YRn与平均值的多元正态随机变量μ和协方差矩阵Σ。让Z:=log(Y),即Zi=log(Yi),i{1,,n}。如何比较适合观察到的实现模型的AIC Y与适合观察到的实现模型的AIC Z



我最初的问题和稍长的问题:

YN(μ,Σ)是一个多变量正态随机变量。如果我想比较适合于Y的模型与适合对的模型log(Y),可以看看它们的对数似然性。但是,由于这些模型不是嵌套的,因此我无法直接比较对数可能性(以及诸如AIC之类的东西),但必须对其进行转换。

我知道如果X1,,Xn是具有联合pdf 随机变量,g(x1,,xn)并且Yi=ti(X1,,Xn)进行一对一转换tii{1,,n},则PDF的Y1,,Yn其中 J是与变换关联的雅可比行列式。

f(y1,,yn)=g(t11(y),,tn1(y))det(J)
J

我是否只需要使用转换规则进行比较

l log Y = log n i = 1 ϕ log y i; μ Σ

l(Y)=log(i=1nϕ(yi;μ,Σ))
l(log(Y))=log(i=1nϕ(log(yi);μ,Σ))

还是我还能做些什么?


[edit]忘记将对数放在最后两个表达式中。


您似乎在最后一个表达式中也失去了Jacobian。
whuber

2
我不了解转换。当Y为负时,如何获取日志YloglogYY
semibruin

Answers:


6

当拟合两个不同的数据集Z时,您无法比较AIC或BIC 。仅当适合同一数据集时,才可以比较两个基于AIC或BIC的模型。看一看模型选择和多模型推理:一种实用的信息理论方法(Burnham和Anderson,2004年)。他们在第81页(第2.11.3节“响应变量的转换”)中提到了我的答案:YZ

研究人员应确保所有假设都使用相同的响应变量建模(例如,如果整个模型集均基于log(y),则不会出现问题;这是因为响应变量混合不正确)。

顺便说一句,对于使用AIC或BIC标准,您的模型不必一定是嵌套的(参考文献,第88页,第2.12.4节非嵌套模型),实际上,这就是使用BIC的优点之一。


5

log{y(n)+1}log{y(n)+1}n=1,2,...,N

Akaike,H.,1978年。“关于时间序列模型的可能性”,《皇家统计学会杂志》,D系列(统计学家),第27(3/4)页,第217–235页。


1
R中是否有一种方法可以做到这一点?
森林生态学家
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.