Добавление цвета в круговые данные на основе членства в группе

Question

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

Знайте, что я загрузил библиотеку (круговую), что не позволяет мне использовать ggplot. Вот синтаксис я работал с:

plot(bcirc, stack=FALSE, bins=60, shrink= 1, col=w$dx, axes=FALSE, xlab ="Basal sCORT", ylab = "Basal sAA") 

Если вы заметили, я указал переменный фактор (который имеет два уровня) в разделе цвета, но он просто продолжает положить все в одном цвете. Какие-либо предложения?

Answer 1

Похоже, plot.circular не любит назначать несколько цветов. Вот один потенциал обходные:

library(circular) 

## simulate circular data 
bcirc1 <- rvonmises(100, circular(90), 10, control.circular=list(units="degrees")) 
bcirc2 <- rvonmises(100, circular(0), 10, control.circular=list(units="degrees")) 
bcirc <- c(bcirc1, bcirc2) 
dx <- c(rep(1,100),rep(2,100)) 

## start with blank plot, then add group-specific points 
plot(bcirc, stack=FALSE, bins=60, shrink= 1, col=NA, 
    axes=FALSE, xlab ="Basal sCORT", ylab = "Basal sAA") 
points(bcirc[dx==1], col=rgb(1,0,0,0.1), cex=2) # note: a loop would be cleaner if dealing with >2 levels 
points(bcirc[dx==2], col=rgb(0,0,1,0.1), cex=2)