Learn R Programming

FeedbackTS (version 1.1)

map.statistic: Mapping of a spatial index

Description

Mapping of a spatial index with circles whose sizes and colors vary with the values of the index

Usage

map.statistic(coordinates, statistic, region, cex.circles = c(3, 0.2), legend, 
	main = NULL, add = FALSE)

Arguments

coordinates
a 2-column matrix with latitudes and longitudes of observation sites.
statistic
a numeric vector specifying the values, at observation sites, of the index to be mapped.
region
a list of arguments defining the region over which the index is mapped: [object Object],[object Object],[object Object]

Value

  • a plot.

item

  • cex.circles
  • legend
  • main
  • add

describe

  • x:a numeric vector specifying the longitudes where the legend circles are drawn.
  • y:a numeric vector specifying the latitudes where the legend circles are drawn.
  • xtext:a numeric vector specifying the longitudes where the values corresponding to the sizes of the legend circles are drawn.
  • ytext:a numeric vector specifying the latitudes where the values corresponding to the sizes of the legend circles are drawn.
  • digits:integer indicating the number of decimal places to be used.

code

add = FALSE

References

Soubeyrand, S., Morris, C. E. and Bigg, E. K. (in press). Analysis of fragmented time directionality in time series to elucidate feedbacks in climate data. Environmental Modelling and Software.

See Also

map.statistic rain.feedback.stats.

Examples

Run this code
#### load data of feedback and change-in-feedback indices in 88 sites across Australia
data(rain.feedback.stats)

#### spatial coordinates of the 88 sites and corresponding feedback index 
#### computed from the whole data series
coord=rain.feedback.stats[,3:4]
stat1=rain.feedback.stats[["Feedback.whole.period"]]

#### map of feedback index 
map.statistic(coord,stat1,cex.circles=c(3,0.2),
   region=list(border='Australia',xlim=c(110,155)),
   legend=list(x=c(rep(114,3),rep(123,2)),y=-c(37,39.5,42,37,39.5),
      xtext=c(rep(114,3),rep(123,2))+1,ytext=-c(37,39.5,42,37,39.5),digits=2),
   main="Feedback")

Run the code above in your browser using DataLab