Molbiol.ru | Project | Protocols | Programs | Literature Web | Companies | Marketplace | Labor exchange Today's active topics [ Log In* | Register* ] Forum: | |
23432432432432dd | Posted 20.01.2024 15:44 |
ИНО | Posted 07.09.2022 22:41 |
Увы с картой при помощи Вашего кода ничего не вышло. Делал на новом рабочем пространстве. Прикрепляю фрагмент с матюками: CODE > plot(bb, cols=0) Warning messages: 1: In plot.window(...) : "cols" is not a graphical parameter 2: In plot.xy(xy, type, ...) : "cols" is not a graphical parameter 3: In axis(side = side, at = at, labels = labels, ...) : "cols" is not a graphical parameter 4: In axis(side = side, at = at, labels = labels, ...) : "cols" is not a graphical parameter 5: In box(...) : "cols" is not a graphical parameter 6: In title(...) : "cols" is not a graphical parameter > points(bb$x, bb$y, pch=1, col=cc[rr]) # next step: deal with legend Error in bb$x : $ operator is invalid for atomic vectors При том выдается не карта, а обычный скаттерплот с незакрашенными точками. Возможно играют роль версии ПО, в моем случае это R 3.3.1 и spatstat 1.46-1. В то время как plot(сн_2017Pgppp) работает вполне корректно у меня, никаких error. И проблемы с кириллицей видел только в заголовках столбцов таблиц. Наверное, дело в операционке. Ну, уберите из названия объекта эти две незнакомые буржуям буквы, если мешают, это ж просто пример, мне б принцип понять. С деревом можно никаких файлов не грузить, а обойтись простым синтетическим примером: CODE y<-c(rnorm(10, 0, 1), rnorm(10, 0, 100)) x<-factor(c(rep(1,10), rep(2,10))) xy<-cbind(x, y) ctree(y~x, data=xy)#ветвление дерева нет, так как проверяется разница в средних, а надо чтобы проверялась разница в дисперсиях (или любом другом параметре масштаба). xy$z<-c(abs(y[1:10]-mean(y[1:10])), abs(y[1:20]-mean(y[1:20])))#модуль отклонений от среднего в каждой группе ctree(z~x)#теперь различия уловлены правильно Надо чтоб преобразование y в z проходило автоматически в зависимости от формулы, переданной ctree(). По идее для этого есть в ctree() аргумент ytrafo, но как его правильно задавать, я так и не понял. |
|
plantago | Posted 07.09.2022 07:21 |
(ИНО @ 06.09.2022 10:49) Как заставить ctree() из {party} работать не с положением (средними), а с масштабом (дисперсией) зависимой переменной? Насколько понял из руководства, для этой цели можно использовать трансформацию при помощи аргумента ytrafo, но как именно им пользоваться так и не уразумел. Что только не перепробовал и с готовыми ytrafo из {coin}, на которые дается ссылка в том же руководстве и с самопиской функцией, считающей абсолютное отклонение от среднего, - либо матерится разными матюками, либо молча считает различия в средних, точно так же как с настройками по умолчанию, будто бы никаких манипуляций с ytrafo и не было. К сожалению, без конкретных примеров и ссылок не смог ничего понять. |
|
plantago | Posted 07.09.2022 07:20 |
(ИНО @ 06.09.2022 11:37) Прицепом еще один вопрос. Есть у меня точки на карте (файлы прилагаются): CODE library(spatstat) load("сн_2017Pgppp.Rda") load("сн_2017Pg_lL.Rda") plot(сн_2017Pgppp) Задача - залить их цветом в зависимости от значений сн_2017Pg_lL, и приделать сбоку цветовую шкалу в виде б. м. непрерывного градиента и подписями значений соответствующего показателя через равные интервалы. CODE library(spatstat) load("сн_2017Pgppp.Rda") load("сн_2017Pg_lL.Rda") # plot(сн_2017Pgppp) # error! aa <- (get(ls()[1])) # does not always work; if not, use rm(list=ls()) _before_ load() bb <- (get(ls()[2])) # and please do not name R objects with anything except [A-z0-9_.] rr <- rank(aa, ties.method="last") cc <- terrain.colors(max(rr)) plot(bb, cols=0) points(bb$x, bb$y, pch=1, col=cc[rr]) # next step: deal with legend Про шкалу можно, наверное, посмотреть сюда |
|
ИНО | Posted 06.09.2022 18:37 |
Прицепом еще один вопрос. Есть у меня точки на карте (файлы прилагаются): CODE library(spatstat) load("сн_2017Pgppp.Rda") load("сн_2017Pg_lL.Rda") plot(сн_2017Pgppp) Задача - залить их цветом в зависимости от значений сн_2017Pg_lL, и приделать сбоку цветовую шкалу в виде б. м. непрерывного градиента и подписями значений соответствующего показателя через равные интервалы. |
|
ИНО | Posted 06.09.2022 17:49 |
Как заставить ctree() из {party} работать не с положением (средними), а с масштабом (дисперсией) зависимой переменной? Насколько понял из руководства, для этой цели можно использовать трансформацию при помощи аргумента ytrafo, но как именно им пользоваться так и не уразумел. Что только не перепробовал и с готовыми ytrafo из {coin}, на которые дается ссылка в том же руководстве и с самопиской функцией, считающей абсолютное отклонение от среднего, - либо матерится разными матюками, либо молча считает различия в средних, точно так же как с настройками по умолчанию, будто бы никаких манипуляций с ytrafo и не было. | |
plantago | Posted 27.11.2021 12:47 |
Вот и замечательно, что Вы все сами нашли. Только paste0(...) -- без пробела между скобками и именем команды (это для красоты и именно paste0() а не paste(), потому что второй без дополнительных аргументов вставляет пробел, в то время как Excel -- нет. | |
Алекс3212 | Posted 26.11.2021 20:42 |
Очень долгие ответы тут, если вдруг кому то нужно, то используйте paste () , где в кавычках перечисляете столбцы которые надо присоединить |
|
Алекс3212 | Posted 26.11.2021 15:53 |
Подскажите пожалуйста аналог в R , знака & с помощью которого я могу присоединить одну яцейку к другой. Пока что приходится умножать на н значение и складывать с другой ячейкой. Мож знает кто более легкий способ. Например A B C 10 3 13 в R C <- ((A*100)+B) в excel C= A&B |
|
plantago | Posted 08.03.2021 02:54 |
Review the complete topic (launches new window) | |