2017-05-25 25 views
0

當我繪製了數刻度的箱線圖,晶須是基於未記錄的數據。我可以在對數級盒形圖上使用日誌晶須嗎?

x <- rlnorm(n=50, meanlog=0, sdlog=1) 
library('beeswarm') 
beeswarm(x, log=TRUE) 
boxplot(x, add = TRUE, outline = FALSE) 

Of course,如果我手動登錄數據,然後再在該晶須反映這種轉變。

beeswarm(log(x)) 
boxplot(log(x), add = TRUE, outline = FALSE) 

是否有可能具有與第一圖形與第二的晶須軸線?也就是說,我可以在日誌軸上繪製未記錄的數據,但仍有「記錄」的晶須?

回答

1

您可以將boxplot(log(x))計算得出的值轉換回原始比例x。我不知道得到的情節雖然是多麼有意義:

x <- rlnorm(n=50, meanlog=0, sdlog=1) 
library('beeswarm') 
beeswarm(x, log=TRUE) 
box = boxplot(log(x), add = FALSE, plot = FALSE, outline = FALSE) 
box$stats = exp(box$stats) 
box$conf = exp(box$conf) 
bxp(box, add=TRUE) 

enter image description here

+0

尼斯(+1)。但是,您的意思是「我不確定最終的情節有多麼有意義」? – Sparhawk

+0

按照慣例,晶須代表1.5 *四分位數間距 - 當您計算對數轉換變量上的晶須,然後顯示它們與可能引起誤解的原始值。但你比我更瞭解情節應該代表什麼。 – Marius

+0

啊,是的,好點。那麼,這是生物數據,因此(可能)是對數正態分佈。因此,理想情況下應該將其繪製在日誌軸上,而IMO對於標籤具有一些有意義的單位更具生物相關性。 FWIW [這個答案](https://stats.stackexchange.com/a/47859/20013)也表明我應該計算日誌值的鬍鬚,但是不同,因爲這些軸也被標記了記錄值。所以在我的情況下,即使軸處於對數空間中,「約定」會將晶須解釋爲算術空間嗎? [這個問題](https://stats.stackexchange.com/q/112705/20013)奮鬥... – Sparhawk

相關問題