regarding EOT figure

[soni2288]

Hi,
While generating the figures, I have faced the following issues that I hope you can help me resolve:

1. Colorbar Range Customization: For comparison, I am plotting EOT Mode results for five different years and would like to keep the color scale consistent across all figures for both the predictor (from 0 to 1) and response mode (from -1 to +1). However, I’m unsure where to modify the colorbar range in the plotting function.

2. Custom Country Boundary: In other analyses, I have defined and plotted India’s boundary only. However, in the current EOT plots, the default map includes surrounding countries as well. I would like to apply the same India-only boundary used in my other figures. How can I implement this boundary constraint in the response mode plots?

3. Lat-Lon Grid and Box: I would also like to add a latitude-longitude grid or bounding box to both the predictor and response maps for better geographic reference.
   Please see the attached figure, which I want to modify:

[soni2288]

One more query I have is whether it is possible to generate and save the outputs of the analysis. If yes, I can plot as I wish.

[tim-salabim]

One more query I have is whether it is possible to generate and save the outputs of the analysis. If yes, I can plot as I wish.

Yes, the eot function returns all relevant layers:

data(vdendool)

## claculate 2 leading modes
nh_modes <- eot(x = vdendool, y = NULL, n = 2, 
                standardised = FALSE, 
                verbose = TRUE)

str(nh_modes, max.level = 5)
Formal class 'EotStack' [package "remote"] with 2 slots
  ..@ modes:List of 2
  .. ..$ mode_01:Formal class 'EotMode' [package "remote"] with 19 slots
  .. .. .. ..@ mode              : int 1
  .. .. .. ..@ name              : chr "mode_01"
  .. .. .. ..@ eot               : num [1:50] -22.5 -4.27 5.67 3.6 8.56 ...
  .. .. .. ..@ coords_bp         : num [1, 1:2] -35 77.5
  .. .. .. .. ..- attr(*, "dimnames")=List of 2
  .. .. .. ..@ cell_bp           : int 87
  .. .. .. ..@ cum_exp_var       : num 0.215
  .. .. .. ..@ r_predictor       :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ rsq_predictor     :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ rsq_sums_predictor:Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ int_predictor     :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ slp_predictor     :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ p_predictor       :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ resid_predictor   :Formal class 'RasterBrick' [package "raster"] with 13 slots
  .. .. .. ..@ r_response        :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ rsq_response      :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ int_response      :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ slp_response      :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ p_response        :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ resid_response    :Formal class 'RasterBrick' [package "raster"] with 13 slots
  .. ..$ mode_02:Formal class 'EotMode' [package "remote"] with 19 slots
  .. .. .. ..@ mode              : int 2
  .. .. .. ..@ name              : chr "mode_02"
  .. .. .. ..@ eot               : num [1:50] 76.1 75 33.8 65.5 -44.8 ...
  .. .. .. ..@ coords_bp         : num [1, 1:2] -165 42.5
  .. .. .. .. ..- attr(*, "dimnames")=List of 2
  .. .. .. ..@ cell_bp           : int 326
  .. .. .. ..@ cum_exp_var       : num 0.382
  .. .. .. ..@ r_predictor       :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ rsq_predictor     :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ rsq_sums_predictor:Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ int_predictor     :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ slp_predictor     :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ p_predictor       :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ resid_predictor   :Formal class 'RasterBrick' [package "raster"] with 13 slots
  .. .. .. ..@ r_response        :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ rsq_response      :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ int_response      :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ slp_response      :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ p_response        :Formal class 'RasterLayer' [package "raster"] with 13 slots
  .. .. .. ..@ resid_response    :Formal class 'RasterBrick' [package "raster"] with 13 slots
  ..@ names: chr [1:2] "mode_01" "mode_02"

Maybe also have a look at the plot method, though it is using lattice and sp, which should probably be ported to ggplot at some stage...