Think that we need to plot many pairs, e.g. comparing A vs. G. What we can do in R:
#generate two correlated variables
a=rnorm(1000000)
b=a+rnorm(a)
plot is rather horrible
![](https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhVeWNa8GJ7GxldB0Mfqysr7vlGa100FZ4piFXRxT4pBJS9d6NPPXpvrGUjNAlmpPxTzzhSfO1Sh_YYZKzUiMH9mviOjNCvJG9bGytn9j6PPpjhMb6fv1B2IYy7W8NvM3WCk5zbTrElV6U/s320/Untitled3.png)
plot(a,b)
We can draw a map with level curves
require(MASS)
contour(kde2d(a,b))
![](https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEg6hwg33EMtTS_pBL3o-VB3nNlVl4suPJi_qaBgY4Xx1fGur_jIyYA3VZrN_EA8rCn0X3g0aDNws8a2k5wPgYjLHJAkGSS60lL9vvruILUvx4MRNuLgq2pdF7x1lrafqFBEB4SzHHojevU/s320/Untitled2.png)
But I find that hexagonal binning is even nicer (and very fast !!)
require(hexbin)
tt=hexbin(a,b)
plot(tt)
![](https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEg9Od6LC6cXduhe0Mq7rASSgr7yYvHDr8Kp7dD8H28sv4LO5QRUr6z2TtU-vtmHAJsOaFZSASixEbXJjGicRxXo7w_Xezj4woBGgkEUkuyd7cWRZwzSSj0ct_mLqs3_NCN5vjRQ0dn2Jro/s320/Untitled1.png)
Finally, we have persp:
persp(kde2d(a,b))
No comments:
Post a Comment