plotting object with season dimension = 12

[marchtaylor]

Hi there,

I think there may need to be an adjustment to the way the season dimension is handled in ggplotFL. In my example, I have 12 months in the season dimension. I believe season gets originally treated as a factor, and then converted to numeric. This seems to mess up the month order in my case (e.g. "1", "10", "11", "12", "2", "3", "4", "5", "6", "7", "8", "9").
Below is a fix that re-orders the levels based on their sequence in the dimnames of the FLStock object.

direct plot with ggplotFL:

    plot(stk@stock.n)

adjusted version with ggplot2:

    DIMNAMES <- dimnames(stk@stock.n)
    df <- as.data.frame(stk@stock.n)
    df$season <- factor(df$season, levels = DIMNAMES$season)
    df$date <- as.Date(paste(df$year, 
      as.numeric(df$season), "01", sep = "-"), 
      format = "%Y-%m-%d")
    p <- ggplot2::ggplot(data = df, 
      aes(x = date, y = data, group = age)) +
      facet_wrap(~ age, ncol = 1, scales = "free_y", 
        drop=TRUE, strip.position = "right") +
      scale_y_continuous(limits = c(0,NA)) +
      geom_line()
    print(p)

[iagomosqueira]

The conversion takes place in FLCore::as.data.frame(FLQuant), with date=TRUE. The date column is of class ISOdate, and season is passed as numeric. I have added a small check for seasons being months to use the first day as season start. See here

https://github.com/flr/FLCore/blob/master/R/FLQuant.R#L1037

Does this help?

[marchtaylor]

Works perfect. I assume that since I didn't need to add any additional arguments to plot(stk@stock.n), that the default is date=TRUE? Many thanks. Feel free to close the issue. -Marc

[iagomosqueira]

The call to as.data.frame inside ggplotFL::plot sets date=TRUE, yes, although only uses it in the x axis if length(season) > 1.

[marchtaylor]

Hi there - it looks to me as if this functionality has been lost. The season dimension is no longer correctly handled.

[iagomosqueira]

Have you reinstalled from github over the last few days? I have changed a bit the engine, now using our own geom_ function. Will check this is working as it should with the new code

[marchtaylor]

Yes, I don't know exactly when it happened, but installing one of the other packages (FLBRP?) automatically reinstalled several other dependencies as well. I tried reinstalling ggplotFL today and still get the plotting problem - Could the date conversion be the issue again (i.e. ISOdate)?

[marchtaylor]

It would seem to me that season needs to be changed to a character, then numeric. As is, it looks like a direct conversion to numeric. Thus, this error is possibly not seen in an object with only 4 season levels, whereas with 12 the dimension order gets jumbled.

now: as.numeric(res$season)
should be: as.numeric(ac(res$season))

[iagomosqueira]

OK, will fix now, and will release new ggplotFL.

[iagomosqueira]

Could you possibly install the package in github and see if it has been fixed?

[marchtaylor]

Looks good now - thanks so much. One caveat: the plotting function removes empty ages / years. I think I am fine with this, but just wanted to make sure that you want this as the default.

[iagomosqueira]

Are you plotting a single FLQuant? It probably shouldn't, specially since the function offers you to plot dots at data=0 for NAs. Let me check

[marchtaylor]

Yes, a single FLQuant. The behavior you describe used to be the case, but no longer is. NAs, are converted to zeros, and dims containing all NAs are dropped completely. I can send you an example.