根据组成员身份向循环数据添加颜色

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)