2013-01-15 25 views
19

Yoğunluğa orantılı regresyon arsa rengimi geliştirmek istiyorum. Örneğin güven aralığı genişse, gölge aralığı genişse, dolgu aralığı genişse dolgu rengi açık olur. Sonuç grafiği aşağıdaki gibi görünebilir:gradyan gölgeli güven aralığı

set.seed(1234) 
md <- c(seq(0.01, 1, 0.01), rev(seq(0.01, 1, 0.01))) 
cv <- c(rev(seq(0.01, 1, 0.01)), seq(0.01, 1, 0.01)) 
rv <- rnorm (length(md), 0.1, 0.05) 

df <- data.frame(x =1:length(md), F = md*2.5 + rv, L =md*2.5 -rv-cv, U =md*2.5+ rv+ cv) 
plot(df$x, df$F, ylim = c(0,4), type = "l") 

polygon(c(df$x,rev(df$x)),c(df$L,rev(df$U)),col = "cadetblue", border = FALSE) 
lines(df$x, df$F, lwd = 2) 
#add red lines on borders of polygon 
lines(df$x, df$U, col="red",lty=2) 
lines(df$x, df$L, col="red",lty=2) 
+0

sayesinde kaynak kodu vwReg fonksiyonu için, ben bir çözüm (son örnek) olan kullanılarak ... İşte – jon

+1

bulamadım temel grafik sistemi [www.alisonsinclair.ca] (http://www.alisonsinclair.ca/2011/03/shading-between-curves-in-r/) Verilerinize uyarlanabilir –

cevap

7

denstrip pakette densregion() komut istediğini yapmak görünüyor:

İşte

enter image description here

bir çalışma örneğidir. onun yardım sayfasında örnekten küçük bir adaptasyon:

require(denstrip) 
x <- 1:10 
nx <- length(x) 
est <- seq(0, 1, length=nx)^3 
se <- seq(.7,1.3,length.out=nx)/qnorm(0.975) 
y <- seq(-3, 3, length=100) 
z <- matrix(nrow=nx, ncol=length(y)) 
for(i in 1:nx) z[i,] <- dnorm(y, est[i], se[i]) 
plot(x, type="n", ylim=c(-3, 3),xlab="") 
densregion(x, y, z) 
lines(x,est,col="white") 

enter image description here

+0

Güzel yanıtlar için teşekkürler iki farklı yönlendirme açısı verir, yani uyum çizgisine doğru, güven aralığının kendisinin genişliği daha fazladır. – jon

+0

Evet, haklısınız. Bunun nedeni standart sapmanın $ x $ artırılması için artmasıdır; bu nedenle normal yoğunluğun tepe noktası daha düşüktür (ör., dnorm (0,0,1) '$> $' dnorm (0,0,2) '). Gri tonlamanın regresyon çizgisi etrafında aynı olmasını istiyorsanız, z z matrisini yeniden ölçeklendirebilirsiniz (ancak daha sonra gri tonlama, "gerçek" yoğunluğa karşılık gelmeyecektir, sadece yoğunluğun ölçeklenmiş bir katı). –