sgeostat

Вопросы по статистическому пакету R. Не обязательно гео.
Ответить
Аватара пользователя
Игорь Черниенко
Активный участник
Сообщения: 137
Зарегистрирован: 28 мар 2009, 01:05
Репутация: 11
Откуда: Хабаровск, Южно-Сахалинск

sgeostat

Сообщение Игорь Черниенко » 07 янв 2010, 03:04

Рад приветствовать!
После применения fix к объектам этого модуля он перестает воспринимать их классы как свои (скажем, point после открытия как point уже не воспринимается). Это так и должно быть?

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Re: sgeostat

Сообщение Максим Дубинин » 07 янв 2010, 05:10

не покажете ли фрагмент кода?
пристегивайтесь, турбулентность прямо по курсу

Аватара пользователя
Игорь Черниенко
Активный участник
Сообщения: 137
Зарегистрирован: 28 мар 2009, 01:05
Репутация: 11
Откуда: Хабаровск, Южно-Сахалинск

Re: sgeostat

Сообщение Игорь Черниенко » 07 янв 2010, 05:18

Примерно так:
> sj.xvg_o<-est.variogram(sj.pt,sj.pr,"cmj2") ' sj.pt - точки; оценивается вариограмма, при желании - строится график
> fix(sj.pt) открываем, смотрим, ничего не меняем
> sj.xvg_o<-est.variogram(sj.pt,sj.pr,"cmj2","cmf2") то же что в первой строке
Ошибка в est.variogram(sj.pt, sj.pr, "cmj2", "cmf2") :
Point.obj must be of class, "point".

т.е sj.pt больше не поинт. если выполнить команду
sj.pt<-point(sj.pt) все встает на свои места.

такая же ситуация с вариограммами и прочим.

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Re: sgeostat

Сообщение Максим Дубинин » 07 янв 2010, 07:53

да, fix() может поменять класс объекта, так как ему приходится "засовывать" его в некоторые слишком узкие для него рамки, в итоге обкусываются некоторые атрибуты и меняется класс. А иногда может и не менять, смотрите например:

Код: Выделить всё

x <- ts(1:5)
class(x)
fix(x)
class(x)
Результат - ts, класс не изменился. А вот, если mts:

Код: Выделить всё

y <- ts(cbind(1:5,1:5)) 
class(y)
fix(y)
class(y)
Результат был mts, cтал matrix. Странновато, но некоторая логика есть. Из справки:
fix invokes edit on x
Editing an R object may change it in ways other than are obvious
пристегивайтесь, турбулентность прямо по курсу

Аватара пользователя
Игорь Черниенко
Активный участник
Сообщения: 137
Зарегистрирован: 28 мар 2009, 01:05
Репутация: 11
Откуда: Хабаровск, Южно-Сахалинск

Re: sgeostat

Сообщение Игорь Черниенко » 07 янв 2010, 22:00

Спасибо, буду иметь ввиду.

Аватара пользователя
Игорь Черниенко
Активный участник
Сообщения: 137
Зарегистрирован: 28 мар 2009, 01:05
Репутация: 11
Откуда: Хабаровск, Южно-Сахалинск

Re: sgeostat

Сообщение Игорь Черниенко » 08 янв 2010, 17:06

Кстати, то же самое происходит с объектами хтабс в модуле гмоделс.

Ответить

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

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

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