的R - 顏色地塊根據p值

Question

我有數百個箱圖的個從我的巨大的數據集自動打印和我已製造的代碼來着色地塊「珊瑚」如果變量具有跨越選定的組的統計顯着性和'海藍寶石'如果沒有檢測到重要性。這是我用col=ifelse(...< 0.05,'coral','aquamarine')做的。

我以爲我的輸出會更容易滑過，如果箱線圖是根據p值着色。因此，我想告訴col=...將p值從0.05 - 0.01變爲黃色，0.01 - 0.005變爲桔黃色和< 0.005變爲紅色的變量着色。

我試圖與

boxplot(a~b, data=df, col=if(pv1<0.05)"yellow" else if (pv1<0.01) "orange" else if (pv1<0.005) "red" else "green") 

正如你可以看到，這不會做，因爲值重疊。我需要做的是寫在if句子裏面的範圍，但我只是不知道如何去做。

這裏有一些數據可以使用。請注意，組的實際p值與此處給出的值不同。我只顯示pv1,2和3，以便測試代碼。

a <- c(23,24,64,12,4,75,12,65,86,76) 
b <- c(1,2,2,1,2,1,3,3,1,3) 
df <- data.frame(a,b) 
pv1 <- c(0.05) 
pv2 <- c(0.01) 
pv3 <- c(0.005) 

Answer 1

要修改if...else說法，你應該認爲它周圍（在這裏，你問的其他方式，「如果它不是< 0.05，那麼如果是< 0.01 ......」這是不可能的），所以你可以寫：

col=if(pv1<0.005)"red" else if (pv1<0.01) "orange" else if (pv1<0.05) "yellow" else "green" 

例

pv1 <- 0.04 ; if(pv1<0.005)"red" else if (pv1<0.01) "orange" else if (pv1<0.05) "yellow" else "green" # "yellow" 
pv1 <- 0.004 ; if(pv1<0.005)"red" else if (pv1<0.01) "orange" else if (pv1<0.05) "yellow" else "green" # "red" 
pv1 <- 0.06 ; if(pv1<0.005) "red" else if (pv1<0.01) "orange" else if (pv1<0.05) "yellow" else "green" # "green" 
pv1 <- 0.006 ; if(pv1<0.005) "red" else if (pv1<0.01) "orange" else if (pv1<0.05) "yellow" else "green" # "orange" 

Anothe -r選項，如果你的p值是一個矢量，你可以嘗試cut：

mycolours <- as.character(cut(pv1, ,c(1,0.05,0.01,0.005,0), right=F, labels=c("red","orange","yellow","aquamarine"), include.lowest=T)) 

然後

boxplot(a~b, data=df, col=mycolours)