[R密度圖：如何從x軸的實線增加了密度曲線

Question

我有一個密度圖的頂部繪製使用：

plot(density(x)) 

我感興趣的東西做的是創造一個從x軸到曲線上相應的點，從x = 5行開始。

像這樣：

Answer 1

可以通過第一存儲密度值中的對象，然後檢索從該對象的x和y元件做到這一點。在下面的例子中，我使用findInterval檢索y值對於給定的x值：

x <- rnorm(1000) # Sample data 
y <- density(x) # Calculate and store density 

x0 <- 2 # Desired position on x-axis 
y0 <- y$y[findInterval(x0, y$x)] # Corresponding y value 

plot(density(x)) 
segments(x0, 0, x0, y0) 

Answer 2

如果它真的z分數那麼就實際上畫出密度函數dnorm()。它也像你想有你的實際x軸爲0

zmax <- 4 
curve(dnorm, -zmax, zmax, xaxt = 'n', bty = 'n') 
axis(1, -zmax:zmax, pos = 0) 

要在您的線畫，你可以再次使用dnorm功能。

zscore <- 1.65 
segments(zscore, 0, zscore, dnorm(zscore)) 

你也甚至遮陽出來很好... :)

x <- seq(zscore, zmax, 0.01) 
y <- c(0, dnorm(x)[1:(length(x)-2)],0) 
polygon(x,y, density = 20) 

您還可以使用segments和text命令，將正確的標籤在你的什麼陰影圖和無陰影地區的意思。