R/Plot_DiDiSTATIS_F.b...R

In michaelkriegsman/PlotDiDiSTATIS: Plot results of the package, DiDiSTATIS

# Plot_DiDiSTATIS_F.b..
#
#' Plot F.b.. (Barycentric Grand Factor Scores)
#'
#' @param res_DiDiSTATIS The output of DiDiSTATIS
#' @param axes Axes to plot
#' @param Flip_axis1 Flag to flip Component 1
#' @param Flip_axis2 Flag to flip Component 2
#' @param dev.new Flag to appease ReporteRs (set FALSE to print results to pptx)
#' @export

Plot_DiDiSTATIS_F.B.. <- function(res_DiDiSTATIS, axes = c(1,2), Flip_axis1 = FALSE, Flip_axis2 = FALSE, dev.new = TRUE){

  Flip_mat <- makeFlip_mat(Flip_axis1, Flip_axis2)

  alpha <- .8
  prettyPlot(res_DiDiSTATIS$res_BaryGrand$eig$Fb..Cond[,axes] %*% Flip_mat,
             col = add.alpha(res_DiDiSTATIS$input$DESIGN_rows$colors_B, alpha),
             cex = 3, pch=15, dev.new = dev.new,
             xlab = paste0("Component ", axes[1]," = ", round(res_DiDiSTATIS$res_BaryGrand$eig$tb..[axes][1],1), "% of Barycentric variance"),
             ylab = paste0("Component ", axes[2]," = ", round(res_DiDiSTATIS$res_BaryGrand$eig$tb..[axes][2],1), "% of Barycentric variance"),
             constraints = minmaxHelper(rbind_array_2_matrix(res_DiDiSTATIS$res_BaryGrand$Proj_B.D$F_B.D_Cond[,axes,]) %*% Flip_mat),
             main = "Barycentric Grand Compromise, F.b..")

}