Войти в систему

Home
    - Создать дневник
    - Написать в дневник
       - Подробный режим

LJ.Rossia.org
    - Новости сайта
    - Общие настройки
    - Sitemap
    - Оплата
    - ljr-fif

Редактировать...
    - Настройки
    - Список друзей
    - Дневник
    - Картинки
    - Пароль
    - Вид дневника

Сообщества

Настроить S2

Помощь
    - Забыли пароль?
    - FAQ
    - Тех. поддержка



Пишет mi_b ([info]mi_b)
@ 2004-08-12 21:30:00


Previous Entry  Add to memories!  Tell a Friend!  Next Entry
смотрим на кривые
Человеческий глаз очень плохо приспособлен к анализу выпуклости кривых. Выше-ниже, растет-падает на графиках видно хорошо, а с выпуклостью, знаком второй производной, нетренированный глаз справляется куда хуже. У меня, похвастаюсь, глаз к этому тренированный, потому что я 11 часов в день пялюсь в кривые, где выпуклость и даже изменения выпуклости очень важны.

Вот тут возник спор о том, как меняется содержание CO2 в атмосфере последние полвека, эспоненциально или линейно (график, данные). Некоторые обсуждающие считали, что, если линейно, то это делает маловероятным антропогенный характер изменений, потому что всякая промышленность и население растут экспоненциально. Это утверждение мне, кстати неочевидно, но это неважно, потому что ясно, что рост там суперэкспоненциальный. Один из участников дискуссии в другом треде попросил меня объяснить, почему это мне очевидно.

Update 2: удивительным образом, после обмена многими десятками комментариев выяснилось, что наиболее горячие противники идеи об экспоненциальном росте не знакомы с определнием "экспоненциального роста", которое было бы применимо к подобным экспериментальным данным.

По определению, временной ряд имеет экспоненциальный рост если его логарифм имеет линейный рост. Это определение эквивалентно постоянному периоду удвоения.

Линейный рост обычно определяют как значимость положительного коэффициента при наклоне и незначимость коэффициента при "нелинейном" члене в соответствующей случаю регрессии. Суперэкспоненциальный рост определяют через более быстрый, чем линейный рост логарифма, то есть, значимо положительный коэффициент при нелинейном выпуклов вверх члене регрессии.




Update 1: картинки заменены линками для облегчения траффика.

1. Берем данные по ссылке. Работать будем с последним столбцом, кототрый отличается от простого среднего за год тем, что авторы использовали чуть более сложную фильтрацию. Выкидываем первый год, где слишком много наблюдений пропущено.


Fig 1


2. Рисуем график Fig 1. Видим: кривая выпукла вниз, она лежит ниже красной хорды, а нарисованное внизу отклонение от ближайшей к ней (зеленой) прямой имеет ярко выраженную чашеобразную форму. Видим, что, помимо роста, есть компонента высокой частоты, поэтому ни на какую простую кривую в точности наши данные не лягут. Но можно говорить об общем характере роста. Рост со временем ускоряется, но это еще не говорит об экспоненциальном росте.


Fig 2


3. Нарисуем теперь график логарифма концентрации Fig 2. Если изменение концентрации приближенно описывается какими-нибудь медленно растущими функциями, например, линейными, то этот график будет выпуклым вверх - лежать выше достаточно длинных отрезков, соединяющих его точки. Если оно описывается двумя склеенными кусками таких функций, мы увидим две склеенных выпуклых вверх кривых. Мы же видим кривую, лежащую заметно ниже таких отрезков (показанных красным цветом). Значит, рост быстрее, чем экспоненциальный. Быстрее ли рост, чем exp(exp(время))?

Fig 3


4. Построим для этого третью картинку, график log(log(концентрации)), Fig 3. Поскольку и он более-менее выпукл вниз, и, большей частью, лежит ниже своих красных хорд, мы имеем дело точно с суперэкпоненциальным ростом, и, возможно, даже двойным экпоненциальным exp(exp(время)).



А если кто статистические тесты любит больше графиков, ему сюда. Результат там тот же самый.




(Читать комментарии) - (Добавить комментарий)


[info]sowa@lj
2004-08-14 13:22 (ссылка)
Ну по экспоненциальному росту-то вы специалист. Дайте ваше определение. Неужели, кроме ругательств, аргументов нет?

Я свое пояснил. Для тех, кто в танке. В моих аргументах изпользуется только одно обстоятельство - функция, точно заданная полиномиальной формулой небольшой степени не может рассматриваться как экспоненциально растущая. Такая, как x^2+1.1 на интервале [0,1]. Функция, заданная полиномом очень высокой степени (скажем, 100) на коротком интервале, может, в принципе, рассматриваться как приближенно экспоненциальная (на конечном интервале, как всем известно, любую функцию можно приблизить полиномом достаточно высокой степени).

Какой сделать выбор - полином или экспонента (или тригонометрический полином, скажем) - зависит от физического смысла данных. От размера коэффициентов, как я уже писал - не численного, а в смысле значимости в рассматриваемой физической задаче. Никакой анализ данных здесь не поможет.

(Ответить) (Уровень выше) (Ветвь дискуссии)


[info]mi_b@lj
2004-08-14 13:45 (ссылка)
Это - не определение. Перечитайте условие - требуется определить "экспоненциальный рост на дискретной выборке зашумленных данных". Никаких точных полиномов степени 100 в условии не упоминается. Не мне же Вас учить, как формулируются определения. Думайте дальше.

Когда будет определение, попробуйте проверить его для какой нибудь небольшой эмпирической нетривиальной выборки, которую принято называть экспоненциально растущей. Например, численности населения какой-нибудь страны или города, или ВВП какой-нибудь страны. Таких данных полно в сети.

И, я Вас прошу в третий раз, не пишите, пожалуйста, ко мне в журнал пока у Вас не будет практически применимого определения.

(Ответить) (Уровень выше) (Ветвь дискуссии)


[info]sowa@lj
2004-08-14 14:01 (ссылка)
Послушайте, я вам экзамен не сдаю. Что это вы мне какие-то условия ставите? Вы мне тут заданий надавали, будто я к вам на упражнения по "анализу данных" хожу. Оставьте это вашим студентам, если они у вас есть.

Я вам свои представления об экспоненциальном росте объяснил. Объясните теперь вы свои, как принято в дискуссиях.

(Ответить) (Уровень выше) (Ветвь дискуссии)


[info]ex_ost922@lj
2004-08-15 19:05 (ссылка)
Два взрослых человека...

(Ответить) (Уровень выше)


(Читать комментарии) -