Thanks a lot, David ! Thats what I wanted. it is greatly helpful.
GP University of Guelph Guelph, ON On Tue, Nov 11, 2014 at 7:02 PM, David Winsemius <[email protected]> wrote: > > On Nov 11, 2014, at 2:53 PM, Gyanendra Pokharel wrote: > > > Thanks David, what do you mean by organized data in regular manner? > > This is what I meant by a regular manner: > > > matrix( c( rep( seq(0,1, by=.1), 11), rep( seq(0,1, by=.1),each=11) , > runif(121) ), 121,3) > [,1] [,2] [,3] > [1,] 0.0 0.0 0.048946906 > [2,] 0.1 0.0 0.332529489 > [3,] 0.2 0.0 0.967099462 > [4,] 0.3 0.0 0.565349269 > [5,] 0.4 0.0 0.024230243 > [6,] 0.5 0.0 0.421633329 > [7,] 0.6 0.0 0.965847357 > [8,] 0.7 0.0 0.719618276 > [9,] 0.8 0.0 0.948675911 > [10,] 0.9 0.0 0.180241643 > [11,] 1.0 0.0 0.804828468 > [12,] 0.0 0.1 0.713698501 > [13,] 0.1 0.1 0.991003966 > [14,] 0.2 0.1 0.936413540 > [15,] 0.3 0.1 0.941731063 > [16,] 0.4 0.1 0.373998953 > [17,] 0.5 0.1 0.988915380 > [18,] 0.6 0.1 0.500791201 > [19,] 0.7 0.1 0.070137099 > [20,] 0.8 0.1 0.968422057 > [21,] 0.9 0.1 0.827396746 > snipped the remain 100 lines > > with( df, contour( x=unique(V1), y=unique(V2), > z= matrix( V3, length(unique(V1)), length(unique(V2) ) > ) )) > > > df <- as.data.frame( matrix( c( rep( seq(0,1, by=.1), 11), rep( seq(0,1, > by=.1),each=11) , runif(121) ), 121,3)) > > str(df) > 'data.frame': 121 obs. of 3 variables: > $ V1: num 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 ... > $ V2: num 0 0 0 0 0 0 0 0 0 0 ... > $ V3: num 0.628 0.661 0.163 0.57 0.527 ... > > My original untested suggestion had some missing parentheses and a missing > comma: > > with( df, contour( x=unique(V1), y=unique(V2), > z= matrix( V3, length(unique(V1)), length(unique(V2) ) > ) )) > > > > Gyanendra Pokharel > > University of Guelph > > Guelph, ON > > > > On Tue, Nov 11, 2014 at 3:53 PM, David Winsemius <[email protected]> > wrote: > > > > On Nov 11, 2014, at 11:17 AM, Gyanendra Pokharel wrote: > > > > > Hi R users, > > > I am trying to plot two dimensional posterior likelihood surface. I > have a > > > data like > > > > > > para1 para2 likehood > > > ....... ........ ........... > > > ....... ........ ........... > > > > > > > > > > > > I looked at contour plot but it needs a shorted values of parameters > and a > > > matrix of likelihood values. Is there any way to get the plot? or how > can I > > > change my likelihood values to a matrix for the function "contour"? > > > > If the data are organized in a regular manner, then this might succeed: > > > > with( df, contour( x=unique(para1), y=unique(para2) > > z= matrix( likehood, length(unique(para1), > length(unique(para2) ) > > ) ) > > > > > > > > Any suggestions are appreciated. > > > > > > > > > GP > > > University of Guelph > > > Guelph, ON > > > > > > [[alternative HTML version deleted]] > > > > > > ______________________________________________ > > > [email protected] mailing list > > > https://stat.ethz.ch/mailman/listinfo/r-help > > > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > > > and provide commented, minimal, self-contained, reproducible code. > > > > David Winsemius > > Alameda, CA, USA > > > > > > David Winsemius > Alameda, CA, USA > > [[alternative HTML version deleted]] ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
