R/graphics.R

   1 #draw (France or...) map with all sites of colors 'cols'
   2 drawMapWithSites = function(M, cols=rep(1,nrow(M)))
   3 {
   4     xMin = range(M[,1])[1]
   5     xMax = range(M[,1])[2]
   6     yMin = range(M[,2])[1]
   7     yMax = range(M[,2])[2]
   8     par(mar=c(2,2,2,2))
   9     plot(0,0,xlim=c(xMin,xMax),ylim=c(yMin,yMax),col="white")
  10     #plot by color groups (limited to integers)
  11     maxColor = max(cols)
  12     for (i in 1:maxColor)
  13     {
  14         indices = (1:nrow(M))[cols==i]
  15         if (length(indices) > 0)
  16             points(M[indices,1],M[indices,2],col=i,xtitle=NULL)
  17     }
  18 }
  19
  20 #draw neighborhoods graph on top of a country map (or any other map)
  21 drawNeighborhoodGraph = function(M, NI)
  22 {
  23     for (i in 1:length(NI))
  24     {
  25         for (j in NI[[i]])
  26             lines(c(M[i,1],M[j,1]),c(M[i,2],M[j,2]))
  27     }
  28 }
  29
  30 #plot a matrix of curves (in rows)
  31 plotCurves = function(M, cols=rep(1,nrow(M)))
  32 {
  33     n = nrow(M)
  34     rg = c(min(M),max(M)) #range for plotting
  35     for (i in 1:n)
  36     {
  37         plot(M[i,],col=cols[i],ylim=rg,type="l")
  38         if (i < n) par(new=TRUE)
  39     }
  40 }