Как метод наименьших квадратов строит линию
Номера лет становятся значениями аргумента, температуры — значениями функции. Метод подбирает такую прямую, чтобы сумма квадратов отклонений наблюдений от неё была минимальной. Наклон прямой равен отношению ковариации номера года и температуры к дисперсии номеров лет.
Наклон измеряется в градусах Цельсия за год. Положительное значение означает потепление в рамках ряда, отрицательное — похолодание. Модуль показывает скорость изменения, но не говорит, насколько оно значимо. Чтобы перейти от градусов в год к изменению за десятилетие, умножьте наклон на десять: для нашего примера это около 2,4 °C за десять лет, что намного выше типичных многолетних скоростей изменения и показывает неустойчивость короткого ряда.
Разбор двух примеров
Для ряда 8,1; 8,4; 8,2; 8,9; 9,0; 9,3 получаем наклон 0,243 °C/год. Линия, продолженная на 3 года за последнее наблюдение, даёт 9,99 °C, а R² равен 0,88: около 88 % разброса объясняется линейным ростом.
Для ряда 12,0; 11,6; 11,9; 11,2; 11,0 наклон отрицательный, -0,240 °C/год, оценка через 2 года 10,58 °C, R² = 0,77. Идеальная линия 10, 11, 12, 13 даёт наклон 1 и R² = 1: все точки лежат на прямой.
Сравнивая два примера, обратите внимание на роль разброса. У второго ряда значения скачут вверх и вниз, но общее направление вниз видно. Именно поэтому полезно смотреть не только на наклон, но и на R².
Что говорит коэффициент детерминации
R² — доля вариации ряда, которую объясняет линия. Значение, близкое к единице, говорит о почти линейном изменении, значение около нуля — о том, что тренд практически не отличим от шума.
Высокий R² на коротком ряду не доказывает климатическое изменение: при пяти-шести точках случайно получить крутую линию нетрудно. Поэтому всегда смотрите на длину ряда и на разброс значений. Во всех случаях тренд лучше сопровождать указанием периода: например, «за последние восемь лет».
Почему экстраполяция требует осторожности
Линия тренда описывает прошлое. Продолжая её вперёд, вы предполагаете, что тот же процесс сохранится. Но температура зависит от циркуляции атмосферы, океанических циклов и извержений вулканов, и межгодовые колебания могут быть значительнее тренда.
Поэтому оценку берут для небольшого горизонта и воспринимают как иллюстрацию, а не как прогноз. Для реальных выводов используют длинные ряды, проверенные станции и статистические критерии значимости.
Подготовка данных и типичные ошибки
Следите за однородностью ряда: смена места установки термометра, прибора или времени наблюдений создаёт ступеньки, которые тренд принимает за потепление. Не смешивайте среднегодовые значения с разными методами осреднения.
Не вводите пропущенные годы как нули и не меняйте порядок данных. Ряд должен идти в хронологическом порядке с шагом в один год. Если часть лет отсутствует, лучше сократить ряд до непрерывного участка, чем заполнять пропуски произвольными числами.
Наконец, сравнивайте температуры одинаковой природы: средние за календарный год с такими же значениями, а не со средними за сезон или за период с мая по сентябрь.
Формула расчёта
T(t) = a + b · t, b = Σ(t − t̄)(T − T̄) / Σ(t − t̄)²
Проверенные примеры
Ряд температур 8,1; 8,4; 8,2; 8,9; 9,0; 9,3 °C. Каков тренд и оценка через 3 года после последнего наблюдения?
Наклон b ≈ 0,243 °C/год, оценка через 3 года ≈ 9,99 °C, R² ≈ 0,88.
Температуры 12,0; 11,6; 11,9; 11,2; 11,0 °C. Найдите наклон и оценку через 2 года.
Наклон ≈ -0,240 °C/год, оценка через 2 года ≈ 10,58 °C, R² ≈ 0,77.
Как пользоваться онлайн-калькулятором
- Выпишите среднегодовые температуры по порядку без пропусков.
- Введите их через пробел.
- Задайте горизонт экстраполяции в годах.
- Нажмите «Вычислить» и оцените наклон, прогноз и R², не забывая об ограничениях короткого ряда.