R：如何將線與多邊形對齊

Question

我無法讓多邊形與曲線對齊。 我將代碼放在這里：

critvalmax <- qt(0.975,df=4)
critvalmin <- qt(0.025,df=4)
xvals <- seq(-5, 5, length=100) 
fx.samp.t <- dt(xvals, df=4)

plot(xvals, dnorm(xvals), col="white") 
lines(xvals, fx.samp.t, lty=1, lwd=2)
abline(v= critvalmin, lty=2)
abline(v= critvalmax, lty=2)
abline(h=0, lty=3)
polygon(cbind(c(critvalmin, xvals[xvals>=critvalmin & xvals<=critvalmax], critvalmax, critvalmax), c(0, dt(critvalmin, df=4), fx.samp.t[xvals>=critvalmin & xvals<=critvalmax], 0)), density=10, lty=3)

結果是多邊形被畫到了右邊，我自己找不到解決方案。 此外，左下角似乎沒有正確填充。

Answer 1

我刪除了polygon中的一些變量，例如第二個critvalmax和dt(critvalmin, df=4)和cbind 。 dt(critvalmin, df=4)似乎將您的多邊形向右移動了0.02558082 。 您可以使用以下代碼：

critvalmax <- qt(0.975,df=4)
critvalmin <- qt(0.025,df=4)
xvals <- seq(-5, 5, length=100) 
fx.samp.t <- dt(xvals, df=4)

plot(xvals, dnorm(xvals), col="white") 
lines(xvals, fx.samp.t, lty=1, lwd=2)
abline(v= critvalmin, lty=2)
abline(v= critvalmax, lty=2)
abline(h=0, lty=3)
polygon(c(critvalmin, xvals[xvals>=critvalmin & xvals<=critvalmax], critvalmax), c(0, fx.samp.t[xvals>=critvalmin & xvals<=critvalmax], 0), density=10, lty=3)

^{使用reprex v2.0.2創建於 2022-08-31}

polygon只需要你的 x 和 y：

包含多邊形頂點坐標的向量。

和density ：

陰影線的密度，以每英寸線數為單位。 NULL 默認值表示不繪制陰影線。 密度為零意味着沒有陰影或填充，而負值和 NA 抑制陰影（因此允許顏色填充）。