Как сравнить статистически 2 набора данных?

[rhot]

Есть 2 разновременных растра, содержат в себе информацию о местоположении сосняков (доля сосны > 70%) на одну и ту же территорию. Разница во времени между растрами 80 лет.

Как сравнить эти наборы данных статистически?

Первое, что приходит на ум, это использовать взаимнокорреляционную функцию, например, stats::ccf из R. Но я не уверен, что это подойдёт для пространственных временных рядов.

[bolotoved]

Думаю, что вопрос "Как сравнить эти наборы данных статистически?" не очень правильный с точки зрения эколога. Нужно задать осмысленные вопросы типа: "Изменилось ли расположение сосняков в пространстве?", "Увеличилась ли фрагментированность местообитаний?", "Изменилась ли форма и размер сосняков?" и т.п.

На такие вопросы уже можно отвечать. Например можно использовать разнообразные ландшафтные метрики о которых можно почитать для начала, например в руководстве по программе FragStats. Думаю, что большинство из этих метрик уже реализовано в R.

[rhot]

Спасибо за наводку.

Моя гипотеза, что боры никуда не ушли, т.к. это статус ООПТ уже очень давно. Поэтому хочу просто оценить качество старого набора данных. Фрагментированность, форма и размер не интересуют.

[bolotoved]

Чтобы утверждать, что сосняки не изменились, вам нужно знать их возраст. За 80 лет сосна вырастает очень здорово.

Если же вы уверены, что сосняки где были, там и остались, то задача у вас действительно сводится к оценке качества старых данных. В этом случае нужно использовать стандартные методы оценки качества категоризованных данных типа матрицы ошибок, каппы и т.п. Подробнее о методах тут: Congalton, R.G., and Green, K. (2008). Assessing the accuracy of remotely sensed data: principles and practices (CRC).

[nadiopt]

что за снимок возрастом 80 лет назад?
а если это лесные карты, "в лоб" их сравнить со снимком не получится. Вам ли не знать, что доля сосны по снимку и по лесоустройству - разные вещи

[rhot]

В общем, нашёл 2 способа:
1) посчитать перекрывающиеся категории 2 растров (на обожке книги, предложенной bolotoved'ом, как раз про это )
2) вычислить proximal alignment (не знаю как по-русски будет...)

Источник

PS: nadiopt, я про снимок ни разу не упомянул. Это не ДДЗЗ.

[nadiopt]

две специальные карты чтоль? интересно, расскажете?
а есть эталон? мы же не знаем, кто больше наврал при составлении карты

[rhot]

Сравнение можно сделать с помощью недавно вышедшего пакета diffeR.