R/graphics.R

   1 #' Draw (France or...) map with all sites of colors 'cols'
   2 #'
   3 #' @param M Coordinates matrix (in columns)
   4 #' @param cols Vector of colors for each row of M [default: all black]
   5 #' @export
   6 #'
   7 drawMapWithSites = function(M, cols=rep(1,nrow(M)))
   8 {
   9     xMin = range(M[,1])[1]
  10     xMax = range(M[,1])[2]
  11     yMin = range(M[,2])[1]
  12     yMax = range(M[,2])[2]
  13     par(mar=c(2,2,2,2))
  14     plot(0,0,xlim=c(xMin,xMax),ylim=c(yMin,yMax),col="white")
  15     #plot by color groups (limited to integers)
  16     maxColor = max(cols)
  17     for (i in 1:maxColor)
  18     {
  19         indices = (1:nrow(M))[cols==i]
  20         if (length(indices) > 0)
  21             points(M[indices,1],M[indices,2],col=i,xtitle=NULL)
  22     }
  23 }
  24
  25 #' Draw neighborhoods graph on top of a country map (or any other map)
  26 #'
  27 #' @param M Coordinates matrix (in columns)
  28 #' @param NI Neighborhoods of M rows (list of integer vectors)
  29 #' @export
  30 #'
  31 drawNeighborhoodGraph = function(M, NI)
  32 {
  33     for (i in 1:length(NI))
  34     {
  35         for (j in NI[[i]])
  36             lines(c(M[i,1],M[j,1]),c(M[i,2],M[j,2]))
  37     }
  38 }
  39
  40 #' Plot a matrix of curves (in rows)
  41 #'
  42 #' @param M Coordinates matrix (in columns)
  43 #' @param cols Vector of colors for each row of M [default: all black]
  44 #' @export
  45 #'
  46 plotCurves = function(M, cols=rep(1,nrow(M)))
  47 {
  48     n = nrow(M)
  49     rg = c(min(M),max(M)) #range for plotting
  50     for (i in 1:n)
  51     {
  52         plot(M[i,],col=cols[i],ylim=rg,type="l")
  53         if (i < n) par(new=TRUE)
  54     }
  55 }