Форумы GIS-Lab.info

Есть два грида (кригинг), из них вырезается участок (число пикселей No data одинаково) и считается регрессия (Multiple regression analyses (Grid/Grid). Все замечательно, коэффициенты хорошие, но... на карте регрессии число пикселей No data уменьшилось, а на карте остатков в этих пикселях очень большие отрицательные числа. Как с этим бороться?

SONIA писал(а):Есть два грида (кригинг), из них вырезается участок (число пикселей No data одинаково) и считается регрессия (Multiple regression analyses (Grid/Grid). Все замечательно, коэффициенты хорошие, но... на карте регрессии число пикселей No data уменьшилось, а на карте остатков в этих пикселях очень большие отрицательные числа. Как с этим бороться?

судя по всему, NA пропадают в силу интерполяции, там по умолчанию сплайны, если NA на границе, могут дать "свечку" (для точного ответа смотрите исходники). Поставьте метод ближайшего соседа, должно помочь.

P.S. Да, и если это реальная работа, то для регрессии лучше использовать боле подходящие средства (тот же R), и соответствующие протоколы анализа, поскольку критерий нравится/не нравится ненаучный

Спасибо, буду пробовать!

В продолжение

R2 для регрессии равен 0.57, что совсем неплохо. остатки хотелось бы посмотреть на предмет нормальности.

.57 -- ето ниочём.

Ну, это как посмотреть. К сожалению, зависимости с R2=0.99 часто оказываются артефактами. Здесь R2 - в долях, а не в процентах, так что это очень о чем.

SONIA писал(а):R2 для регрессии равен 0.57, что совсем неплохо. остатки хотелось бы посмотреть на предмет нормальности.

у нас народ традиционно путает fit и model, вот R2 - это про fit, а нормальность остатков, отсутствие зависимости дисперсии остатков от среднего, расстояние Кука, H-в-шляпе, и прочие мелочи - это про model. Проблема только в том, что для выводов нужен именно model

P.S. а R2 получить близкий к вашему не проблема, хватит одной точки где-нибудь вдали от основной группы.

Вот поэтому и хочется посмотреть на остатки. Как-то мне подозрительна эта хорошая зависимость, скорее всего она действительно обусловлена отскоком. Если после обрезки остаются пикселы с no DATA (а это -99999), то вполне возможно, что это и так. Наверное, я как-то неправильно вырезаю кусок, хотя размер пиксела, и число точек с no data в обоих гридах совпадают.