Plotting single points ?

[Thomas Petit-Jean]

Hello,

I have an issue when trying to plot a single value with mf_map. For instance:

mtq<-mf_get_mtq()
mf_map(mtq)
mf_map(mtq[1,], var = c("POP", "STATUS"), type = "prop_typo")

Returns

Error in xy.coords(x, y, xlab = deparse1(substitute(x)), ylab = deparse1(substitute(y))) : 
  'x' is a list, but does not have components 'x' and 'y'

The problem actually remains when only a single row has a non-zero POP:

tmp<- mtq[0:2,]
tmp$POP<-c(10,0)
mf_map(mtq)
mf_map(tmp, var = c("POP", "STATUS"), type = "prop_typo")

returns the same error.

Is there something I'm doing wrong ? Thanks for the help.

[Timothée Giraud]

I'm AFK right now, I'll look into it next week.

[Denis de Oliveira Rodrigues]

Hello,

I'm not sure if classify intervals works with just one record or two records. Is it really necessary or is it a just instance?

[Timothée Giraud]

Sure, it may not be very useful to plot just one proportional symbol, it's an edge case. But at least it should fail gracefully (or work :-))

[Thomas Petit-Jean]

Well in my case this situation occurred quite naturally when iterating over multiple territories: some of them had only one point to plot, and it would be quite natural to be able to plot them as the other ;)

[Timothée Giraud]

Well in my case this situation occurred quite naturally when iterating over multiple territories: some of them had only one point to plot, and it would be quite natural to be able to plot them as the other ;)

That's a perfectly legitimate case. I'll be able to work on this bug on Thursday.

[Timothée Giraud]

Commit f702f26 should have fixed the problem.
reprex:

library(mapsf)
mtq <- mf_get_mtq()
opar <- par(mfrow = c(3,4))
mf_theme(x = "dark")
for( i in 9:18){
  mf_map(mtq)
  mf_map(x = mtq[i, ], var = c('POP', 'STATUS'), type = "prop_typo", 
         symbol = "circle",inches = .2,        
         val_order = c('Prefecture', "Sub-prefecture", "Simple municipality"),
         pal = c("red", "blue", "green"),
         val_max = max(mtq$POP), leg_pos = c(NA,NA))
  mf_title(mtq[i,"LIBGEO", drop = TRUE])
}
# legends
mf_theme(bg = NA)
mf_init(mtq)
mf_legend(type = "prop", pos = "top", title = "Population",
          val = c(min(mtq$POP), 10000, 40000, max(mtq$POP)),
          inches = .2, col = "grey80",
          cex = 2, title_cex = 1.5, val_cex = 1)
mf_init(mtq)
mf_legend(type = "typo", pos = "top", title = "Status",
          val =  c('Prefecture', "Sub-prefecture", "Simple municipality"), 
          pal = c("red", "blue", "green"),
          cex = 1.5, title_cex = 1.5, val_cex = 1)
par(opar)

^{Created on 2022-11-03 with reprex v2.0.2}