Saya perlu menggambar grafik yang kompleks untuk analisis data visual. Saya memiliki 2 variabel dan sejumlah besar kasus (> 1000). Misalnya (angka 100 jika membuat dispersi kurang "normal"):
x <- rnorm(100,mean=95,sd=50)
y <- rnorm(100,mean=35,sd=20)
d <- data.frame(x=x,y=y)
1) Saya perlu memplot data mentah dengan ukuran titik, yang sesuai dengan frekuensi relatif dari kebetulan, jadi plot(x,y)
bukan pilihan - saya perlu ukuran titik. Apa yang harus dilakukan untuk mencapai ini?
2) Pada plot yang sama saya perlu memplot interval kepercayaan 95% elips dan garis yang mewakili perubahan korelasi (tidak tahu bagaimana cara menamainya dengan benar) - sesuatu seperti ini:
library(corrgram)
corrgram(d, order=TRUE, lower.panel=panel.ellipse, upper.panel=panel.pts)
tetapi dengan kedua grafik pada satu plot.
3) Akhirnya, saya perlu menggambar model regresi linar yang dihasilkan di atas semua ini:
r<-lm(y~x, data=d)
abline(r,col=2,lwd=2)
tetapi dengan rentang kesalahan ... sesuatu seperti pada QQ-plot:
tetapi untuk kesalahan pemasangan, jika memungkinkan.
Jadi pertanyaannya adalah:
Bagaimana cara mencapai semua ini dalam satu grafik?