根據組成員資格向循環數據添加顏色

Question

我正在嘗試根據組成員身份為循環數據中的特定點添加顏色（我分為兩組：一組具有一定的醫療狀況的個人，另一組只是健康的對照組）。 我已經將其數據從度轉換為弧度，並將其放在繪圖上，但是我還沒有能夠根據我擁有的因子變量有選擇地更改點的顏色）。

知道我已經加載了庫（循環庫），這不允許我使用ggplot。 這是我一直在使用的語法：

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

如果您已經注意到，我在顏色部分指定了factor變量（具有兩個級別），但它只是將所有內容都放在一種顏色中。 有什么建議么？

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)

Answer 2

受到Paul Regular的示例的啟發，這是一個使用相同數據的版本，其中一個條件向內堆疊，另一個條件向外堆疊。

library(circular)

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

## start with blank plot, then add group-specific points
dev.new(); par(mai = c(1, 1, 0.1,0.1))
plot(circular(subset(bcirc, condition == 1)$angles, units = 'degrees'), stack=T, bins=60, shrink= 1, col=1,sep = 0.005, tcl.text = -0.073,#text outside
     axes=T, xlab ="Basal sCORT", ylab = "Basal sAA")
par(new = T)
plot(circular(subset(bcirc, condition == 2)$angles, units = 'degrees'), stack=T, bins=60, shrink= 1.05, col=2,
     sep = -0.005, axes=F)#inner circle, no axes, stacks inwards