Анализ нескольких временных рядов

Вопросы по статистическому пакету R. Не обязательно гео.
_taras_
Активный участник
Сообщения: 186
Зарегистрирован: 28 июл 2018, 08:40
Репутация: 12
Откуда: Киев

Анализ нескольких временных рядов

Сообщение _taras_ » 17 янв 2021, 20:14

Приветствую всех!
Помогите разобраться с вопросом.
Имеются многолетние данные по продуктивности вида в разных условиях местообитания. Я сделал базовую статистику, сравнение по годам и анализ временных изменений в каждому местообитанию.
Однако не пойму как сделать сравнение между несколькими временными рядами. Т.е. сказать, что в таких-то условиях виду хорошо на протяжении длительного времени несмотря на изменение условий, а в других ему "хуже"

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Анализ нескольких временных рядов

Сообщение Константин Силкин » 17 янв 2021, 20:52

Добрый вечер!
Постройте сначала графики зависимости состояния разных объектов от факторов, влияющих на него. Там видно будет что дальше делать. А так, не видя материала, насоветовать можно много чего, да не в тему. Ведь нет же гарантии, что вы верно описали свою ситуацию и стоящие перед вами задачи

_taras_
Активный участник
Сообщения: 186
Зарегистрирован: 28 июл 2018, 08:40
Репутация: 12
Откуда: Киев

Re: Анализ нескольких временных рядов

Сообщение _taras_ » 17 янв 2021, 23:40

Данные собирались в течение 3-х лет 3 раза за вегетацию (2-4 квартал т.е. Июнь, Сентябрь, Ноябрь).
Графики
1- изменение продуктивности от глубины воды
2- сводный ts.plot
Вложения
depth_massa.png
depth_massa.png (16.73 КБ) 6525 просмотров
massa_ts.png
massa_ts.png (14.82 КБ) 6525 просмотров

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Анализ нескольких временных рядов

Сообщение Константин Силкин » 18 янв 2021, 00:06

ph — это водородный показатель воды? Чего-то он у вас в ноль уходит…

_taras_
Активный участник
Сообщения: 186
Зарегистрирован: 28 июл 2018, 08:40
Репутация: 12
Откуда: Киев

Re: Анализ нескольких временных рядов

Сообщение _taras_ » 18 янв 2021, 01:42

Константин Силкин писал(а):
18 янв 2021, 00:06
ph — это водородный показатель воды? Чего-то он у вас в ноль уходит…
Нет это сокращенно Phragmites australis -тростник

gamm
Гуру
Сообщения: 4046
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1050
Ваше звание: программист
Откуда: Казань

Re: Анализ нескольких временных рядов

Сообщение gamm » 18 янв 2021, 06:13

С таким количеством данный у вас нет временных рядов. Максимум использовать mixed model для корреляции внутри года. Можно конечно и корреляцию AR1 привинтить, но данных очень мало
В общем, см. R, gls, lme, lmer. Если данные выложите, попробую.

_taras_
Активный участник
Сообщения: 186
Зарегистрирован: 28 июл 2018, 08:40
Репутация: 12
Откуда: Киев

Re: Анализ нескольких временных рядов

Сообщение _taras_ » 18 янв 2021, 13:50

gamm писал(а):
18 янв 2021, 06:13
Если данные выложите, попробую.
Ссылка на данные https://drive.google.com/file/d/1ZeoFnB ... sp=sharing

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Анализ нескольких временных рядов

Сообщение Константин Силкин » 18 янв 2021, 13:55

Ваши данные не соответствуют вашему вопросу. Вы спрашивали "Почему в таких-то условиях виду хорошо на протяжении длительного времени несмотря на изменение условий". А где условия-то? Только глубина водоёма и всё?

gamm
Гуру
Сообщения: 4046
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1050
Ваше звание: программист
Откуда: Казань

Re: Анализ нескольких временных рядов

Сообщение gamm » 18 янв 2021, 14:39

_taras_ писал(а):
18 янв 2021, 13:50
Ссылка на данные
есть вопросы (чтобы из этого данные соорудитить):
1) Кто у нас зависимая переменная? Тростник (про него больше всего), осока, каждый из них? Что моделируется - покрытие, число растений, вес, все сразу? Кто из них был исходным (все?), а кто получен расчетом?
2) Обычно проективное покрытие, число растений на 1м2, вес, и т.д. измеряют на нескольких плотах (на земле - бросают рамку), потом усредняют. Или это глазомерная оценка всего участки? Нужна первичка, если есть.
3) Мы участки считаем связанными ( в смысле реакции тростника на внешние условия), или анализируем каждый участок отдельно?
4) Откуда глубина, тем более целым числом. Сколько реально было замеров? Нужна первичка, если есть.
Константин Силкин писал(а):
18 янв 2021, 13:55
А где условия-то?
условия - это год и глубина

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Анализ нескольких временных рядов

Сообщение Константин Силкин » 18 янв 2021, 14:52

gamm писал(а):
18 янв 2021, 14:39
условия - это год и глубина
Таки не понятно. Ведь вопрос звучит "Почему в таких-то условиях виду хорошо на протяжении длительного времени несмотря на изменение условий". В вопросе упоминаются некие условия в разное время. Значит время (сезон) — это влияющий фактор, но оно не то условие, которое интересует. Остаётся только глубина. Возможно ещё пространственные координаты, ландшафтное положение. Но в исходных данных об этом ни гу-гу.

_taras_
Активный участник
Сообщения: 186
Зарегистрирован: 28 июл 2018, 08:40
Репутация: 12
Откуда: Киев

Re: Анализ нескольких временных рядов

Сообщение _taras_ » 18 янв 2021, 15:05

gamm писал(а):
18 янв 2021, 14:39
Нужна первичка, если есть.
Дневники на работе сам на карантине. А сохранил только усредненные данные.
gamm писал(а):
18 янв 2021, 14:39
Обычно проективное покрытие, число растений на 1м2, вес, и т.д. измеряют на нескольких плотах (на земле - бросают рамку), потом усредняют
Да именно так и делалось. Как правило 2-3 повтора рамкой 0,25 м2. Проективное покрытие видов пишем с площадки 100м2. Глазомерное здесь только проективное покрытие.
gamm писал(а):
18 янв 2021, 14:39
Откуда глубина, тем более целым числом. Сколько реально было замеров?
Глубина воды измерялась линейкой, данные усреднялось по нескольким (2-3) повторам. Значения округляли.
gamm писал(а):
18 янв 2021, 14:39
Кто у нас зависимая переменная? Тростник (про него больше всего), осока, каждый из них? Что моделируется - покрытие, число растений, вес, все сразу? Кто из них был исходным (все?), а кто получен расчетом?
Все параметры зависимы от глубины воды. Моделируется влияние снижения уровня воды на тростник (как пример). Остальные по аналогии. Так же хотелось бы ответить на вопрос как ведет себя тростник в разных ценозах при снижении уровня обводненности.
Так же не спорю, что данных может быть мало или не те, которые нужны для ответа на вопрос. Буду благодарен если поможете подкорректировать схему опыта. Единственно тонкие измерения - минерализация и т.п. недоступна (в ин-те нет рабочих приборов).

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Анализ нескольких временных рядов

Сообщение Константин Силкин » 18 янв 2021, 15:10

Вы сами указали на направление, по которому, в том числе, нужно двигаться: "как ведет себя тростник в разных ценозах при снижении уровня обводненности". Так дайте эти ценозы!

gamm
Гуру
Сообщения: 4046
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1050
Ваше звание: программист
Откуда: Казань

Re: Анализ нескольких временных рядов

Сообщение gamm » 18 янв 2021, 15:26

_taras_ писал(а):
18 янв 2021, 15:05
Моделируется влияние снижения уровня воды на тростник (как пример).
нет у нас тростника, есть набор переменных. Какие из них моделируем?

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

P.S. А что за площадки, озера эвтрофицированные или прибрежные отмели рек?

_taras_
Активный участник
Сообщения: 186
Зарегистрирован: 28 июл 2018, 08:40
Репутация: 12
Откуда: Киев

Re: Анализ нескольких временных рядов

Сообщение _taras_ » 18 янв 2021, 15:51

gamm писал(а):
18 янв 2021, 15:26
А что за площадки, озера эвтрофицированные или прибрежные отмели рек?
Это т.н. постоянно заливаемые или длительно заливаемые участки болотной растительности в Жебриянских плавнях дельты Дуная. Соленость водоемов скачет как из-за погодных условий (нагон воды с Черного моря) и накопление дождевой воды (летом бывают очень сильные ливни). За последние годы многие участки болотной растительности заменились луговой из-за падения уровня воды.

gamm
Гуру
Сообщения: 4046
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1050
Ваше звание: программист
Откуда: Казань

Re: Анализ нескольких временных рядов

Сообщение gamm » 18 янв 2021, 17:39

_taras_ писал(а):
18 янв 2021, 15:51
За последние годы многие участки болотной растительности заменились луговой из-за падения уровня воды.
это мы за 3 года не поймаем, у нас (в смысле у Вас) и тростник, и осока многолетние, просто видим угнетение, и предполагаем, что лимитирующий фактор - уровень воды. Поскольку питание скорее всего дождевое, то надо бы осадки еще посмотреть. Координаты скажите примерные, я станцию WMO попробую найти поблизости. А соленость, кстати, не измеряли? она влияет?

Ответить

Вернуться в «R»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость