ggplot2 - 2D Normal PDF with ggplot and R -


i'm trying make contour/image plot using ggplot no success til now.

consider following piece of code in r creates matrix z pdf of bivariate normal:

require(mvtnorm) x1 = seq(-3, 3, length.out=200) x2 = seq(-3, 3, length.out=200) z = matrix(0, length(x1), length(x2)) (i in 1:length(x1)) {     = x1     b = x2[i]     z[,i] = dmvnorm(cbind(a,b)) } image(x1,x2,z)

2d normal pdf image plot

is possible plot matrix z using ggplot?

thanks!

# reshape data require(reshape2) dat <- melt(z)  # use geom_raster mimic image  gg <- ggplot(dat, aes(x=var2, y=var1, fill=value)) gg <- gg + labs(x="", y="") gg <- gg + geom_raster() gg <- gg + coord_equal() gg <- gg + scale_fill_gradient(low="red", high="yellow") gg <- gg + scale_x_continuous(expand = c(0, 0)) gg <- gg + scale_y_continuous(expand = c(0, 0)) gg <- gg + theme_bw() gg

enter image description here

you can change axis labels pretty if need to.


Comments

Post a Comment

Popular posts from this blog

ruby - How do I merge two hashes into a hash of arrays? -

http://codepad.org/wmyzqlis i'd merge these 2 arrays: a = { :a => 'a', :b => 'b', :d => 'd' } b = { :a => '1', :b => 'b', :c => 'c' } into hash looks this: {:a => ["a", "1"], :b => ["b", "b"], :c => [nil, "c"], :d => ["d", nil] } this doesn't work: p a.merge(b) { |k, v1, v2| [v1, v2] } # {:c=>"c", :a=>["a", "1"], :b=>["b", "b"], :d=>"d"} it's because hash#merge call block duplicate keys. a = { :a => 'a', :b => 'b', :d => 'd' } b = { :a => '1', :b => 'b', :c => 'c' } hash[(a.keys|b.keys).map { |key| [key, [a.fetch(key, nil), b.fetch(key, nil)]] }] # => {:a=>["a", "1"], :b=>["b", "b"], :d=>["d", nil], :c=>...
Read more

java - the value of local variable is not used -

this question has answer here: java error “value of local variable not used” 7 answers i wrote small program check given movie name exist in folder. getting warning on eclipse "the value of local variable not used" here main method public static void main(string[] args) { // todo auto-generated method stub boolean movie_exist=false; system.out.println("hellow world!"); try{ filewriting newfolrder = new filewriting("h:\\breakinbad2\\breaking bad season 2 complete 720p.brrip.sujaidr"); //movie_exist=newfolrder.checkfilesname("movie name"); system.out.println(newfolrder.checkmovieexist("breaking bad")); movie_exist = newfolrder.checkmovieexist("breaking bad"); }catch(nullpointerexception e){ system.out.println("exce...
Read more