R/print.CA3variants.R

In CA3variants: Three-Way Correspondence Analysis Variants

print.CA3variants <-function(x, printall=FALSE, digits=3,...) {
if (printall==FALSE){   
#--------------------------------------------
#cat("The three-way contingency table considered\n")
#print(x$Data)
#--------------------------------------------
if ((x$ca3type=="CA3")||(x$ca3type=="OCA3")) {

cat("Percentage contributions of the components to the total inertia for column-tube biplots\n")
    print(round(x$inertiacoltub,digits=digits))
cat("\n")
#---------------------------------------------------
cat("Percentage contributions of the components to the total inertia for row-tube biplots\n")
    print(round(x$inertiarowtub,digits=digits))
cat("\n")
#---------------------------------------------------
cat("Percentage contributions of the components to the total inertia for row-column biplots\n")
    print(round(x$inertiarowcol,digits=digits))
cat("\n")
#---------------------------------------------------
}#end for symmetric
#-----------------------------------------------------------
if ((x$ca3type=="NSCA3")||(x$ca3type=="ONSCA3")) {
if (x$resp=="row"){
cat("Percentage contributions of the components to the total inertia for pred biplots\n")
    print(round(x$inertiacoltub,digits=digits))
cat("\n")
}
#---------------------------------------------------
if ((x$resp=="column")||(x$resp=="col")){
cat("Percentage contributions of the components to the total inertia for pred biplots\n")
    print(round(x$inertiarowtub,digits=digits))
cat("\n")
 }
#---------------------------------------------------
if (x$resp=="tube"){
cat("Percentage contributions of the components to the total inertia for pred biplots\n")
    print(round(x$inertiarowcol,digits=digits))
cat("\n")
}
#---------------------------------------------------
}#end for non-symmetric
#-----------------------------------index partition
if ((x$ca3type=="OCA3")||(x$ca3type=="ONSCA3")){
cat("\n Index partition\n\n")
print(round(x$index3res$z,digits=digits))    
cat("\n")
cat("\n Partition of the Term-IJK using polynomials \n\n")
print(round(x$index3,digits=digits))    
cat("\n")
cat("\n Partition of the Term-IJ using polynomials \n\n")
print(round(x$index3res$zij,digits=digits))
cat("\n")
cat("\n Partition of the Term-IK using polynomials \n\n")
    print(round(x$index3res$zik,digits=digits))
cat("\n")
cat("\n Partition of the Term-JK using polynomials \n\n")
    print(round(x$index3res$zjk,digits=digits))
      cat("\n\n")
}
if ((x$ca3type=="CA3")||(x$ca3type=="NSCA3")){
cat("\n Index partition\n\n")
    print(round(x$index3,digits=digits))
    cat("\n\n")}
   
}#end printall false
###############################################################################
if (printall==TRUE){   
#cat("The three-way contingency table \n\n")
#print(x$Data)
# cat("Number of iteration steps \n")
 #   print(x$iteration)
    cat("\n    RESULTS for 3-way Correspondence Analysis\n")
      #  if ((x$ca3type=="CA3")|(x$ca3type=="OCA3")){
      #  cat("Three-way Pearson standardised residuals \n")
      #  print(round(x$xs, digits = 2))
    #} 
#else {
 #       cat("Three-way (weighted) centered column profile table \n")
  #      print(round(x$xs, digits = 2))
   # }
cat("\n Row marginals\n\n")
     print(round(x$pi, digits = digits))
cat("\n Column marginals\n\n")
     print(round(x$pj, digits = digits))
cat("\n Tube marginals\n\n")
     print(round(x$pk, digits = digits))
cat("\n Row-Column marginals\n\n")
     print(round(x$pij, digits = digits))
cat("\n Column-Tube marginals\n\n")
     print(round(x$pjk, digits = digits))
cat("\n Row-Tube marginals\n\n")
     print(round(x$pik, digits = digits))
exin<-x$inertiatot/x$n
orin<-x$inertiaorig/x$n
    cat("Explained inertia (reduced dimensions)", exin, "\n\n")
    cat("Total inertia (complete dimensions)", orin, "\n\n")
#    cat("Percentage of explained inertia on total inertia  \n\n")
# cat("Percentage contribution of components to the total variation\n\n")
#   print(round(x$inertiapcsum,digits=digits))
#   print(round(x$inertiapc,digits=digits))
    #cat("Proportion of explained inertia (when reducing dimensions)\n\n")
    #print(x$prp)
    cat("\n Rows in principal coordinates\n\n")
    print(round(x$fi, digits = digits))
    cat("\n Rows in standard coordinates\n\n")
    print(round(x$fiStandard, digits = digits))
    cat("\n Column-tubes in principal coordinates\n\n")
    print(round(x$gjk, digits = digits))
 cat("\n Column-tubes in standard coordinates\n\n")
    print(round(x$gjkStandard, digits = digits))
    cat("\n Inner Product of coordinates (row x interactive)\n\n")
    print(round(x$iproductijk, digits = digits))
    cat("\n Inner Product of coordinates (col x interactive)\n\n")
    print(round(x$iproductjik, digits = digits))
    cat("\n Inner Product of coordinates (tube x interactive)\n\n")
    print(round(x$iproductkij, digits = digits))
  
#--------------------------------------------
if ((x$ca3type=="CA3")||(x$ca3type=="OCA3")) {
  cat("Core array i.e. Generalised singular values  \n\n")
    print(round(x$g/sqrt(x$cost), digits = digits))
cat("Percentage contributions of the components to the total inertia for row biplots\n\n")
    print(round(x$inertiarow,digits=digits))
cat("Percentage contributions of the components to the total inertia for column-tube biplots\n\n")
    print(round(x$inertiacoltub,digits=digits))
#---------------------------------------------------
cat("Percentage contributions of the components to the total inertia for column biplots\n\n")
    print(round(x$inertiacol,digits=digits))
cat("Percentage contributions of the components to the total inertia for row-tube biplots\n\n")
    print(round(x$inertiarowtub,digits=digits))
#---------------------------------------------------
cat("Percentage contributions of the components to the total inertia for tube biplots\n\n")
    print(round(x$inertiatube,digits=digits))
cat("Percentage contributions of the components to the total inertia for row-column biplots\n\n")
    print(round(x$inertiarowcol,digits=digits))
}#end for symmetric
#--------------------------------------
if ((x$ca3type=="NSCA3")||(x$ca3type=="ONSCA3")) {
if (x$resp=="row"){
  cat("Core array i.e. Generalised singular values  \n\n")
    print(round(x$g/sqrt(x$cost), digits = digits))
cat("Percentage contributions of the components to the total inertia for resp biplots\n\n")
    print(round(x$inertiarow,digits=digits))
cat("Percentage contributions of the components to the total inertia for pred biplots\n\n")
    print(round(x$inertiacoltub,digits=digits))
}
#---------------------------------------------------
if ((x$resp=="column")||(x$resp=="col")){
cat("Percentage contributions of the components to the total inertia for resp biplots\n\n")
    print(round(x$inertiacol,digits=digits))
cat("Percentage contributions of the components to the total inertia for pred biplots\n\n")
    print(round(x$inertiarowtub,digits=digits))
 }
#---------------------------------------------------
if (x$resp=="tube"){
cat("Percentage contributions of the components to the total inertia for resp biplots\n\n")
    print(round(x$inertiatube,digits=digits))
cat("Percentage contributions of the components to the total inertia for pred biplots\n\n")
    print(round(x$inertiarowcol,digits=digits))
}
#---------------------------------------------------
}#end for non-symmetric

#--------------------------------------
    if ((x$ca3type=="OCA3")||(x$ca3type=="ONSCA3")){
cat("\n Index partition\n\n")
print(round(x$index3res$z,digits=digits))    
cat("\n")
cat("\n Partition of the Term-IJK using polynomials \n\n")
print(round(x$index3,digits=digits))    
cat("\n")
cat("\n Partition of the Term-IJ using polynomials \n\n")
print(round(x$index3res$zij,digits=digits))
cat("\n")
cat("\n Partition of the Term-IK using polynomials \n\n")
    print(round(x$index3res$zik,digits=digits))
cat("\n")
cat("\n Partition of the Term-JK using polynomials \n\n")
    print(round(x$index3res$zjk,digits=digits))
      cat("\n\n")
}
if ((x$ca3type=="CA3")||(x$ca3type=="NSCA3")){
cat("\n Index partition\n\n")
    print(round(x$index3,digits=digits))
    cat("\n\n")}
}
}