Supponiamo di eseguire una regressione lineare semplice bayesiana. Vorrei visualizzare i risultati tracciando più linee di regressione in base alle distribuzioni posteriori di (intercetta) e b (pendenza). Mi chiedo come visualizzare i risultati in stile heatmap o in alternativa utilizzare la trasparenza per evitare sovrapposizioni. Ecco un semplice approccio ggplot.Heatmap delle linee di regressione
library(ggplot2)
set.seed(123)
N = 1000
x = 1:80
a = rnorm(N,10,3)
b = rnorm(N,5,2)
y = vector("list",length=N)
for(i in 1:N) {y[[i]] = a[i]+b[i]*x}
df = data.frame(x=rep(x,N),y=unlist(y))
df$f = rep(1:N,each=80)
(plt <- ggplot(df, aes(x, y,group=f)) +
geom_jitter(alpha=1/30,width=5,col="blue") + theme_classic())
Ci sono modi migliori per farlo? Sarebbe bello se il colore cambierebbe a seconda della quantità di sovrapposizione (come avviene in heatmaps).
io pensi di avere un errore di battitura: 'df $ f = rep (1: N, ciascuna = 80)' – csgillespie
Grazie! Corretto – beginneR