rplot_combo.R

library(dplyr, warn.conflicts=FALSE)
library(ggplot2); theme_set(theme_bw())
library(forcats)
library(shellpipes)
rpcall("rplot_combo.Rout .pipestar rplot_combo.R exp.egf_single.rds logistic.egf_single.rds series.tsv")

startGraphics()

dat <- (bind_rows(rdsReadList())
	|> left_join(tsvRead(),by=c("loc"="varname"))
	|> mutate(loc_final = ifelse(loc == "Tokyo1", "Tokyo",loc)
		, loc_final = ifelse(loc_final == "Tokyo2", "Tokyo",loc_final)
   	, loc_final = paste0(loc_final, "\n", year)
	)

#	rbind(readRDS("exp.egf_single.rds"),readRDS("logistic.egf_single.rds"))
	## |> mutate(loc = factor(loc, levels=rev(levels(loc))))
)


print(dat)

print(
       ggplot(dat)
	+ aes(
	    # to simply reverse:
	    # x=forcats::fct_rev(loc)
	    # This is reordering by the first est for each loc:
	    # x = forcats::fct_reorder(loc, est)
	    # This is reordering by the mean est:
	    y = forcats::fct_reorder(loc_final, est, .fun = mean)
	    , color=phase)
        + geom_pointrange(aes(x=est,xmin=lwr,xmax=upr,lty=method,
                          shape = method
                          )
		, position = position_dodge(width=-0.4)
	)
	+ scale_linetype_manual(values=c("dotted","solid"))
	+ scale_color_manual(values=c("black","red", "orange"))
	+ labs(x = "r (1/month)", y = NULL)
)