根據組成員身份向循環數據添加顏色

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)