######ESTIMATE POLYNOMIAL MODEL FOR RSA
RSA_NSfit <- RSAmodel(formula= engagement ~ needs*supplies,
data= sim_NSfit, model= c("FM26_PARALLELASYMWEAK"))
######PLOT EXTREMA OVER LOC AND LOIC
EXTsim <- plotting.ext(RSA_NSfit,model="FM26_PARALLELASYMWEAK",
xlim=c(-3,3),zlim=c(-3,3),acceleration=c(0,-0.3),
text_size=0.7,elabel_size=3,e_size=2,breaks_x=seq(-10, 10, 0.5),breaks_z=seq(-10, 10, 0.5))
ggpubr::ggarrange(EXTsim[["LOC"]], EXTsim[["LOIC"]],
labels = c("A. Response over LOC", "B. Response over LOIC"),
nrow=1,ncol=2,font.label = list(size = 11))
Run the code above in your browser using DataLab