Muszę narysować złożoną grafikę do wizualnej analizy danych. Mam 2 zmienne i dużą liczbę przypadków (> 1000). Na przykład (liczba wynosi 100, jeśli dyspersja jest mniej „normalna”):
x <- rnorm(100,mean=95,sd=50)
y <- rnorm(100,mean=35,sd=20)
d <- data.frame(x=x,y=y)
1) Muszę wykreślić surowe dane z rozmiarem punktu, odpowiadającym względnej częstotliwości zbieżności, więc plot(x,y)
nie ma opcji - potrzebuję rozmiarów punktów. Co należy zrobić, aby to osiągnąć?
2) Na tym samym wykresie muszę wykreślić elipsę 95% przedziału ufności i linię reprezentującą zmianę korelacji (nie wiem jak poprawnie ją nazwać) - coś takiego:
library(corrgram)
corrgram(d, order=TRUE, lower.panel=panel.ellipse, upper.panel=panel.pts)
ale z dwoma wykresami na jednym wykresie.
3) Na koniec muszę narysować wynikowy model regresji liniowej nad tym wszystkim:
r<-lm(y~x, data=d)
abline(r,col=2,lwd=2)
ale z zakresem błędów ... coś jak na wykresie QQ:
ale w przypadku błędów dopasowania, jeśli jest to możliwe.
Pytanie brzmi:
Jak to wszystko osiągnąć na jednym wykresie?
źródło
df.new <- data.frame(x = seq(min(x), max(x), 0.1))
s size is also strange (too small). Also tryed
library(car) cr.plots(m0)
car::dataEllipse
ellipse
corrgram
pakiecie: pokazuje ona 95% obszar ufności parami, zakładając dwuwymiarowy rozkład normalny wyśrodkowany na średniej i skalowany przez SD (x) i SD (y). Jednak nie jestem wielkim fanem tego, gdy używa się go na wykresie rozrzutu. Ale patrz Murdoch i Chow, Graficzny wyświetlacz dużych macierzy korelacji , Am Stat (1996) 50: 178 lub Friendly, Corrgrams: Wyświetlacze eksploracyjne dla macierzy korelacji , Am Stat (2002) 56: 316.Dla punktu 1 wystarczy użyć
cex
parametru na wykresie, aby ustawić rozmiar punktu.Na przykład
Aby mieć wiele wykresów na jednym wykresie, użyj
par(mfrow=c(numrows, numcols))
równomiernie rozmieszczonego układu lublayout
wykonaj bardziej złożone.źródło
cex
, ale myślę, że OP chce wszystkich rzeczy na tym samym obszarze kreślenia, a nie na osobnych.curve
lubpoints
przesłonić trzy wykresy;)