Почему мы не можем использовать

10

Представьте, что у нас есть модель линейной регрессии с зависимой переменной . Мы находим его . Теперь мы делаем другую регрессию, но на этот раз для , и аналогично находим ее . Мне сказали, что я не могу сравнить оба чтобы увидеть, какая модель лучше подходит. Это почему? Причиной для меня было то, что мы будем сравнивать изменчивость разных величин (разных зависимых переменных). Я не уверен, что это должно быть достаточной причиной для этого.y log ( y ) R 2 log ( y ) R 2Ry2log(y)Rlog(y)2R2

Также есть ли способ формализовать это?

Любая помощь будет оценена.

Старик в море.
источник
1
Я подозреваю, что это могло обсуждаться ранее на Cross Validated. Тщательно ли вы прошли через подобные темы? Кроме того, вас интересуют разные зависимые переменные (например, ВВП относительно цены на нефть) или трансформации одной и той же переменной (ВВП против роста ВВП) или и то, и другое?
Ричард Харди
@RichardHardy Я нашел несколько, но я думаю, что они касались моего вопроса. Как этот: stats.stackexchange.com/questions/235117/… Ответ просто говорит да, не объясняя, почему.
Старик в море.
@RichardHardy Я заинтересован в преобразованиях зависимой переменной.
Старик в море.
1
R2 сравнения имеет смысл только между вложенными моделями.
LVRao
@LVRao Спасибо за ваш комментарий. Почему это так?
Старик в море.

Ответы:

8

Это хороший вопрос, потому что «разные величины» не являются чем-то вроде объяснения.

Есть две важные причины, по которым следует с осторожностью использовать для сравнения этих моделей: он слишком грубый (он на самом деле не оценивает качество соответствия ) и будет неподходящим по крайней мере для одной из моделей. Этот ответ касается второй проблемы.R2


Теоретическое лечение

сравнивает дисперсию остатков модели с дисперсией ответов. Дисперсия - среднеквадратичное аддитивное отклонение от подгонки. Таким образом, мы можем понимать R 2 как сравнение двух моделей ответа y . R2R2y

Модель «база» является

(1)yi=μ+δi

где - параметр (теоретический средний отклик), а δ i - независимые случайные «ошибки», каждая из которых имеет нулевое среднее значение и общую дисперсию τ 2 .μδiτ2

Модель линейной регрессии вводит векторы качестве объясняющих переменных:xi

(2)yi=β0+xiβ+εi.

Число и вектор β являются параметрами (точка пересечения и «уклоны»). Ε я снова независимые случайные ошибки, каждый с нулевым средним и дисперсией общей сг 2 .β0βεiσ2

оценивает уменьшение дисперсии, т 2 - σ 2 ,сравнению с исходной дисперсией т 2 .R2τ2σ2τ2

Когда вы берете логарифмы и используете наименьшие квадраты, чтобы соответствовать модели , вы неявно сравниваете отношения формы

(1a)log(yi)=ν+ζi

к одной из форм

(2a)log(yi)=γ0+xiγ+ηi.

(1)(2)(2a)

yi=exp(log(yi))=exp(γ0+xiγ)exp(ηi).

Термины ошибки теперь умножают базовые отношенияexp(ηi)yi=exp(γ0+xiγ)

Var(yi)=exp(γ0+xiγ)2Var(eηi).

xi(2)σ2

Обычно только один из этих наборов моделей может быть разумным описанием данных. Применение второго набора и ( 2 а )(1a)(2a)(1)(2)R2R2

Анализ

R2xyεiηi

Такая модель (которая обычно встречается) представляет собой метод наименьших квадратов, соответствующий экспоненциальному отношению,

(3)yi=exp(α0+xiα)+θi.

yx(2a)θi(2)R2xy

(3)300xi1.01.6(x,y)(x,log(y))

диаграммы рассеяния

R20.700.56R2R20.70


log(y)(3)

Whuber
источник
Критика на R ^ 2 несправедлива. Как и любой инструмент, его использование должно быть хорошо понято. В приведенных выше примерах R ^ 2 дает правильное сообщение. R ^ 2 выбирает лучшее соотношение сигнал / шум. Конечно, это не очевидно, когда вы кладете два графика с совершенно разными масштабами рядом. На самом деле сигнал слева очень сильный по сравнению с шумовыми отклонениями.
Кагдас Озгенц
@Cagdas Похоже, вы предлагаете противоречивое сообщение. Поскольку эти два графика неизбежно находятся в двух разных масштабах - один показывает исходные ответы, а другой - их логарифмы, - тогда умоляя о том, что что-то «не очевидно» из-за этого неизбежного факта, похоже, не поддерживается ваш случай. Жаловаться на то, что этот ответ «несправедливый», на самом деле не выдерживает явный анализ предложенных мною моделей.
whuber
В том, что я говорю, нет никаких противоречий. R ^ 2 выбирает более высокое отношение сигнал / шум. Это то, что он делает. Попытка превратить это во что-то другое и заявить, что это не работает, совершенно неправильно. Вся критика в отношении R ^ 2 также применяется к другим показателям соответствия, когда они применяются к другой переменной отклика, но по какой-то причине R ^ 2 выбран в качестве козла отпущения.
Кагдас Озгенц
R2R2R2
2
Спасибо за вашу помощь. Извините за поздний прием, у меня не было много свободного времени в последнее время. ;)
Старик в море.