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


17

我的问题的实质是:

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



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

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

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

f(y1,…,yn)=g(t1−1(y),…,tn−1(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为负时,如何获取日志Y?loglog⁡YY
— 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.