notes/spline_compare.R

In bbolker/McMasterPandemic: Pandemic Model

library(McMasterPandemic)
library(tidyverse);theme_set(theme_bw())

source("makestuff/makeRfuns.R")
commandEnvironments()
makeGraphics()

print(plot(ff, data=dat))

aa <- coef(ff, "fitted")

print(aa)

quit()

#spline_df <- data.frame(time = 1:nrow(X)
#	, calib_spline = exp(X[,-1] %*% matrix(aa$time_beta, ncol=1))
#	, true_spline = exp(X[,-1] %*% matrix(bb[-1], ncol=1))
#)

print(spline_df)

sdf <- (spline_df
	%>% gather(key="type",value="value",-time)
)

gg <- (ggplot(sdf,aes(x=time,y=value,color=type))
	+ geom_line()
	+ scale_color_manual(values=c("red","black"))
)

print(gg)