R/plot_loadings.R

In cstawitz/spatiotemporalBSAIgrowth: Analyze Growth in Bering Sea Groundfish

plot_loadings <- function(L_pj,years,nfactor){
  

  clr <- growth_cols(6)[c(1,2,3,6)]
  for(i in 1:nfactor){
    lim1 <-  range(L_pj)
    Ylim = c(lim1[1]-.001,lim1[2]+.001)
    
    plot(1, type="n", xlim=c(0.5,nrow(L_pj)+0.5), ylim=Ylim, xlab="", ylab="", xaxt="n", xaxs="i",las=1)
    mtext( text=paste("Factor",i), side=3, line=0.1, adj=0)
    abline(h=0)
    Labels <- c("walleye pollock","Pacific cod", "yellowfin sole", "flathead sole")
    
    axis(labels = Labels, at=c(1:4), side=1, las=2)
    # Loop through categories and plot each
    for(p in 1:nrow(L_pj)){
      lines(y=c(0,L_pj[p,i]), x=rep(p,2), lwd=5, col=clr[p])
    }
    legend( "top", legend=paste0(round(100*sum(L_pj[,i]^2)/sum(L_pj^2),1),"%"), bty="n")
  }
}