[R-sig-Geo] How to show the color bar legend in R?

Tue Jan 27 17:08:56 CET 2015

These are the information I got:
R version 3.1.2 (2014-10-31)Platform: x86_64-w64-mingw32/x64 (64-bit)
locale:[1] LC_COLLATE=English_United States.1252 [2] LC_CTYPE=English_United States.1252   [3] LC_MONETARY=English_United States.1252[4] LC_NUMERIC=C                          [5] LC_TIME=English_United States.1252    
attached base packages:[1] stats     graphics  grDevices utils     datasets  methods   base     
other attached packages:[1] plotrix_3.5-11     ks_1.9.3           rgl_0.95.1201      mvtnorm_1.0-2     [5] misc3d_0.8-4       KernSmooth_2.23-13
loaded via a namespace (and not attached):[1] tools_3.1.2

     On Monday, January 26, 2015 6:33 PM, Robert J. Hijmans <r.hijmans at gmail.com> wrote:

 > It seems there is a problem with my R !

It is for a good reason that you are requested to include your

sessionInfo()

information in messages to this list.  But even more important is to
first make sure that you are using the current version of R and the
current versions of all relevant packages.

Robert

On Mon, Jan 26, 2015 at 12:28 PM, Ferra Xu via R-sig-Geo
<r-sig-geo at r-project.org> wrote:
> Thank you for your email Zack. When I run the code that you provided, it gives me this error:
> Error in if (axes) { : argument is not interpretable as logicalIn addition: Warning message:In if (axes) { :  the condition has length > 1 and only the first element will be used
> It seems there is a problem with my R !
>
>      On Monday, January 26, 2015 11:02 AM, Zack Holden <zaholdenfs at gmail.com> wrote:
>
>
>  Hi Ferra,Without reproducible code, it's difficult to help. You might have a look at image.plot() in the fields library, and image()
> I'm not sure if this is exactly what you're looking for, but here is a quick example of one way to control the image and legend in a plot, using the layout command. You can create a small panel below your main plot and then add your legend there with more control.
> I hope this is helpful.
> ############################################ example code for image plot of volcano datarequire(fields)require(raster)
> dat <- raster(volcano)
>
> zc <- as.matrix(seq(minValue(dat), maxValue(dat), length.out=100))
> mycolor1 <- two.colors(n=100, start="darkblue", end="red", middle="yellow1", alpha=.9)
> layout(matrix(c(1,2), 2, 1, byrow = TRUE), heights=c(4,1), widths=c(4) )layout.show(2)image(dat, axes=F, col=mycolor1, xlab="", ylab="")
> par(mar=c(2,4,.1,2))image(zc, col=mycolor1, axes=F, xlab="", ylab="" )
> axis(1, at=seq(0,1, length.out=5),  labels=round(seq(minValue(dat), maxValue(dat),length.out=5), digits=0), las=1,  tcl=-.15, padj=-1, cex=1.5) # every .1 corr# end code#############################################
>
>
> On Mon, Jan 26, 2015 at 10:36 AM, Ferra Xu via R-sig-Geo <r-sig-geo at r-project.org> wrote:
>
> I am still struggling with this question. I appreciate any help from you!
>
>
>      On Friday, January 23, 2015 11:00 AM, Ferra Xu <ferra.xu at yahoo.com> wrote:
>
>
>  Based on the dataset, I computed a 3D kernel (longitude, latitude, depth) and displayed the magnitude with the sphere radius and time with color as follows. I am wondering how I can show the color bar legend.
>
> Also I have problem with the scale of the graph. I need to see if the balls touch or overlap each other from the graph but since the scales are not real I cannot recognize it. How can I correct it?
>
> library(ks)
> x <- read.csv("HK1.csv")
> y <- x[,1:3]
> fhat1 <- kde(x=y)
> n <- nrow(x)
> col <- grey((n:1)/n)
> plot(fhat1, xlim=c(0,500), ylim=c(0,250), zlim=c(0, 4.4))
> spheres3d(y, radius=15*(x$magnitude), alpha=1, col= col)
>
> I tried this code (which didn't work):
>
> legend_image <- as.raster(matrix(grey((n:1)/n), ncol=1))
> plot(c(0,2),c(0,200),type = 'n', axes = F, xlab = '', ylab = '', main =  'time scale')
> text(x=1.5, y = seq(0,200,l=5), labels = seq(0,200,l=5))
> rasterImage(legend_image, 0, 0, 1, 1)
>
> And this one in "plotrix" package:
>
> testcol<- col
> col.labels<- seq(0,200, l=5)
> color.legend(600, 0, 601,200,col.labels,testcol,gradient="y")
>
> But it gives me this error:
>
> Error in rect(xleft, ybottoms, xright, ytops, col = col, lty = 0) :
> plot.new has not been called yet
>
> Also I used legend() in Graphics library:
>
> legend("right", legend=seq(0,200,l=20), fill= colfunc(20), title="Time Scale")
>
> But it gives me the legend in a separate window. Here is my data:
>
>  longitude        latitude        depth            time          magnitude
>
> 422.697323  164.19886  2.457419    8.083796636  0.83367586
> 423.008236  163.32434  0.5551326  37.58477455  0.893893903
> 204.733908  218.36365  1.9397874  37.88324312  0.912809449
> 203.963056  218.4808    0.3723791  43.21775903  0.926406005
> 100.727581  46.60876    1.4022341  49.41510519  0.782807523
> 453.335182  244.25521  1.6292517  51.73779175  0.903910803
> 134.909462  210.96333  2.2389119  53.13433521  0.896529401
> 135.300562  212.02055  0.6739541  67.55073745  0.748783521
> 258.237117  134.29735  2.1205291  76.34032587  0.735699304
> 341.305271  149.26953  3.718958    94.33975483  0.849509216
> 307.138925  59.60571    0.6311074  106.9636715  0.987923188
> 307.76875  58.91453    2.6496741  113.8515307  0.802115718
> 415.025535  217.17398  1.7155688  115.7464603  0.875580325
>
>
>
>
>        [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>
>
>
>
>        [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo

	[[alternative HTML version deleted]]