1.07. Примеры подбора линейных моделей связи между двумя факторами. Фиктивная линейная связь
В этом разделе мы рассмотрим примеры подбора линейных моделей связи для конкретных данных.
Пример 1
В следующей таблице приведены данные об изменении потребительского спроса на куриные яйца семи семейных хозяйств в зависимости от цены на этот продукт в течение 15 недель:
I |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
Спрос |
12 |
10 |
13 |
11.5 |
12 |
13 |
12 |
12 |
12 |
13 |
Цена |
0.54 |
0.51 |
0.49 |
0.49 |
0.48 |
0.48 |
0.48 |
0.47 |
0.44 |
0.44 |
I |
11 |
12 |
13 |
14 |
15 |
Спрос |
13.5 |
14 |
13.5 |
14.5 |
13 |
Цена |
0.43 |
0.42 |
0.41 |
0.40 |
0.39 |
(спрос измерялся в дюжинах, цена — в долларах). Диаграмма рассеяния для этих данных имеет следующий вид:
Предполагая, что модель наблюдений имеет вид где — спрос в I-Ю неделю, а — цена в I-Ю неделю, мы получаем следующие оценки для неизвестных параметров И модели линейной связи между ценой и спросом: Таким образом, подобранная модель линейной связи имеет вид При этом,
Так что коэффициент детерминации оказывается равным т. е. изменчивость цен объясняет 51.4% изменчивости спроса на куриные яйца. На диаграмме рассеяния изображена прямая линия, соответствующая подобранной модели линейной связи.
Пример 2
В следующей таблице приведены данные о годовом потреблении свинины на душу населения в США (в фунтах) и оптовых ценах на свинину (в долларах за фунт) за период с 1948 по 1961 год:
Год |
Потр. |
Цена |
Год |
Потр. |
Цена | |
1948 |
67.8 |
0.5370 |
1955 |
66.6 |
0.4256 | |
1949 |
67.7 |
0.4726 |
1956 |
67.4 |
0.4111 | |
1950 |
69.2 |
0.4556 |
1957 |
61.5 |
0.4523 | |
1951 |
71.9 |
0.4655 |
1958 |
60.2 |
0.4996 | |
1952 |
72.4 |
0.4735 |
1959 |
67.6 |
0.4183 | |
1953 |
63.5 |
0.5047 |
1960 |
65.2 |
0.4433 | |
1954 |
60.0 |
0.5165 |
1961 |
62.2 |
0.4448 |
Для этих данных диаграмма рассеяния имеет вид
Предполагая, что модель наблюдений имеет вид где — потребление свинины в I-Й год рассматриваемого периода, а — оптовая цена на свинину в этом году, мы получаем следующие оценки для неизвестных параметров И модели линейной связи между оптовой ценой и потреблением: Таким образом, подобранная модель линейной связи имеет вид При этом,
, ,
Так что коэффициент детерминации здесь оказывается равным . Изменчивость оптовой цены объясняет здесь лишь 5.5% изменчивости потребления свинины.
Пример 3
Рассмотрим данные о размерах совокупного располагаемого дохода и совокупных расходах на личное потребление в США в период с 1970 по 1979 год. Обе величины выражены в текущих долларах США.
Год |
Расп. доход |
Потребление |
1970 |
695.2 |
621.7 |
1971 |
751.9 |
672.4 |
1972 |
810.3 |
737.1 |
1973 |
914.0 |
811.7 |
1974 |
998.1 |
887.9 |
1975 |
1096.2 |
976.6 |
1976 |
1194.3 |
1084.0 |
1977 |
1313.5 |
1204.0 |
1978 |
1474.3 |
1346.7 |
1979 |
1650.5 |
1506.4 |
Этим данным соответствует диаграмма рассеяния
Предполагая, что модель наблюдений имеет вид где — совокупные расходы на личное потребление в I-Й год рассматриваемого периода, а — совокупный располагаемый доход в этом году, мы получаем следующие оценки для неизвестных параметров И модели линейной связи между совокупным располагаемым доходом и совокупными расходами на личное потребление: Таким образом, подобранная модель линейной связи имеет вид При этом,
, ,
Так что коэффициент детерминации здесь оказывается равным . Изменчивость совокупного располагаемого дохода объясняет здесь более 99.95% изменчивости совокупных расходов на личное потребление.
Впрочем, не следует Слишком оптимистически интерпретировать близкие к единице значения коэффициента детерминации как указание на то, что изменения значений объясняемой переменной практически полностью определяются именно изменениями значений объясняющей переменной. В этой связи, рассмотрим следующий поучительный пример.
Пример 4
Рассмотрим динамику изменений в период с 1957 по 1966 годы трех совершенно различных по природе показателей: E — Суммарного производства электроэнергии в США (в млрд. квт-час), C — совокупных потребительских расходов в Тайланде (в млрд. бат) и H — Мирового рекорда на конец года в прыжках в высоту с шестом среди мужчин (в см). Значения этих показателей приведены в таблице:
Год |
Потребление Тайланд Млрд бат |
Эл. энергия США Млрд квт-час |
Мир. рекорд (прыжки с шестом) См |
1957 |
34.9 |
716 |
478 |
1958 |
35.9 |
724 |
478 |
1959 |
37.9 |
797 |
478 |
1960 |
41.1 |
844 |
481 |
1961 |
43.5 |
881 |
483 |
1962 |
46.7 |
946 |
493 |
1963 |
48.9 |
1011 |
520 |
1964 |
52.0 |
1083 |
528 |
1965 |
56.1 |
1157 |
528 |
1966 |
62.6 |
1249 |
534 |
Динамика изменений показателей показана на графике:
По этим данным мы можем формально, используя метод наименьших квадратов, подобрать модели линейной зависимости каждого из трех показателей от каждого из остальных показателей. Это приводит, например, к моделям
(Заметим, кстати, что произведение угловых коэффициентов двух последних прямых, соответствующих моделям линейной связи, в которых объясняемая и объясняющая переменая меняются местами, равно и совпадает со значением коэффициента детерминации в этих двух подобранных моделях.)
Мы видим, что во всех подобранных моделях значения коэффициента детерминации весьма высоки, и это формально означает, что изменчивость «объясняющих» переменных в этих моделях составляет значительный процент от изменчивости «объясняемой» переменной, стоящей в левой части уравнения. Однако, вряд ли мы всерьез можем полагать, что динамика роста суммарного производства электроэнергии в США Действительно объясняется динамикой роста мирового рекорда по прыжкам в высоту с шестом, несмотря на высокое значение 0.9 коэффициента детерминации в первом из четырех уравнений.
В ситуациях, подобных последнему примеру, принято говорить о Фиктивной (Ложной, паразитной — spurious) линейной связи между соответствующими показателями. И такие ситуации часто встречаются при рассмотрении показателей, динамика изменений которых обнаруживает Заметный тренд (убывание или возрастание) — именно такой характер имеют исследуемые показатели в последнем примере.
Чтобы понять, почему это происходит, вспомним полученное в свое время равенство
Из этого равенства вытекает, что близкие к единице значения коэффициента детерминации соответствуют Близким по абсолютной величине к единице значениям коэффициента корреляции между переменными и . Но этот коэффициент корреляции равен
Где
При фиксированных значениях и , значение Будет тем ближе к , чем большим будет значение Последнее же обеспечивается Совпадением знаков разностей и для максимально возможной доли наблюдений переменных и , что как раз и имеет место, когда в процессе наблюдения обе переменные возрастают или обе переменные убывают по величине. (В этом случае превышение одной из переменных своего среднего значения сопровождается, как правило, и превышением второй переменной своего среднего значения. Напротив, если одна из переменных принимает значение, меньшее среднего значения этой переменной, то и вторая переменная, как правило, принимает значение, меньшее своего среднего.)
Аналогичным образом, значение будет тем ближе к , чем меньшим будет значение Последнее же обеспечивается Несовпадением знаков разностей и для максимально возможной доли наблюдений переменных и , что имеет место, когда в процессе наблюдения одна из переменных возрастает, а вторая убывает. (В этом случае, если одна из переменных принимает значение, меньшее среднего значения этой переменной, то вторая переменная, как правило, принимает значение, большее своего среднего.)
Из сказанного следует, что близость к единице наблюдаемого значения коэффициента детерминации Не обязательно означает наличие причинной связи между двумя рассматриваемыми переменными, а может являться лишь Следствием тренда значений обеих переменных.
Последнее обстоятельство часто наблюдается при анализе различных экономических показателей, вычисленных Без поправки на инфляцию (Недефлированные данные). Проиллюстрируем это следующим примером.
Пример 5
Обратимся к данным о совокупном располагаемом доходе и совокупных личных расходах на местный транспорт в США за период с 1970 по 1983 год. Данные представлены как в текущих долларах США, так и в долларах 1972 года — пересчет к последним выполнен с учетом динамики индекса потребительских цен в указанном периоде. (Уровень цен в 1972 г. принят за 100%.)
Год |
Распол. доход Номинал. |
Расходы Номинал. |
Распол. доход Дефлир. |
Расходы Дефлир. |
1970 |
695.2 |
3.1 |
751.6 |
3.4 |
1971 |
751.9 |
3.3 |
779.2 |
3.4 |
1972 |
810.3 |
3.4 |
810.3 |
3.4 |
1973 |
914.0 |
3.6 |
864.7 |
3.4 |
1974 |
998.1 |
4.0 |
857.5 |
3.5 |
1975 |
1096.2 |
4.4 |
874.5 |
3.5 |
1976 |
1194.3 |
4.7 |
906.4 |
3.6 |
1977 |
1313.5 |
5.0 |
942.9 |
3.6 |
1978 |
1474.3 |
5.5 |
988.8 |
3.7 |
1979 |
1650.5 |
6.2 |
1015.7 |
3.8 |
1980 |
1828.7 |
6.3 |
1021.6 |
3.5 |
1981 |
2040.9 |
6.2 |
1049.3 |
3.2 |
1982 |
2180.1 |
6.6 |
1058.3 |
3.2 |
1983 |
2333.2 |
6.6 |
1095.4 |
3.1 |
Диаграммa рассеяния для недефлированных величин имеет вид
Соответствующая модель линейной связи: Коэффициент детерминации равен . Диаграмме рассеяния дефлированных величин
Соответствует модель линейной связи Коэффициент детерминации равен на этот раз всего лишь .
В связи с последним примером, вернемся к примеру 3 и выясним, не является ли обнаруженная там сильная линейная связь между совокупным располагаемым доходом и совокупными расходами на личное потребление лишь следствием использования недефлированных величин.
Для этого рассмотрим дефлированные значения, представленные следующей таблицей, в последнем столбце которой приведены значения индекса потребительских цен (уровень цен 1972 г. принят за 100%).
Год |
Дефлир. доход |
Дефлир. потребл. |
1970 |
695.2 |
621.7 |
1971 |
751.9 |
672.4 |
1972 |
810.3 |
737.1 |
1973 |
914.0 |
811.7 |
1974 |
998.1 |
887.9 |
1975 |
1096.2 |
976.6 |
1976 |
1194.3 |
1084.0 |
1977 |
1313.5 |
1204.0 |
1978 |
1474.3 |
1346.7 |
1979 |
1650.5 |
1506.4 |
Соответствующая этой таблице диаграмма рассеяния имеет вид
Подобранная модель линейной связи Коэффициент детерминации при переходе от номинальных величин к дефлированным остается очень высоким: . Следовательно, наличие сильной линейной связи между совокупным располагаемым доходом и совокупными расходами на личное потребление Не является только лишь следствием инфляционных процессов.
< Предыдущая | Следующая > |
---|