print.rc <- function(x, digits = max(3, getOption("digits") - 4), ...) {
cat("Call:\n", deparse(x$call), "\n", sep="", fill=TRUE)
ass <- x$assoc
cat("Intrinsic association coefficients:\n")
print(format(ass$phi[1,], digits=digits, ...), quote=FALSE)
cat("\nNormalized row scores (", ass$vars[1], "):\n", sep="")
print(format(ass$row[,,1], digits=digits, ...), quote=FALSE)
cat("\nNormalized column scores (", ass$vars[2], "):\n", sep="")
print(format(ass$col[,,1], digits=digits, ...), quote=FALSE)
if(length(ass$diag) > 0) {
cat("\nDiagonal coefficients:\n")
print(format(ass$diag[1:nrow(ass$diag),], digits=digits, ...), quote=FALSE)
}
cat("\nNormalization weights:", ass$weighting)
printModelStats(x, digits=digits)
}
print.rc.symm <- function(x, digits = max(3, getOption("digits") - 4), ...) {
cat("Call:\n", deparse(x$call), "\n", sep="", fill=TRUE)
ass <- x$assoc
cat("Intrinsic association coefficients:\n")
print(format(ass$phi[1,], digits=digits, ...), quote=FALSE)
cat("\nNormalized row and column scores:\n")
print(format(ass$row[,,1], digits=digits, ...), quote=FALSE)
if(length(ass$diag) > 0) {
cat("\nDiagonal coefficients:\n")
print(format(ass$diag[1:nrow(ass$diag),], digits=digits, ...), quote=FALSE)
}
cat("\nNormalization weights:", ass$weighting)
printModelStats(x, digits=digits)
}
print.hmskew <- function(x, digits = max(3, getOption("digits") - 4), ...) {
cat("Call:\n", deparse(x$call), "\n", sep="", fill=TRUE)
ass <- x[["assoc"]]
if(length(ass) > 0) {
cat("Intrinsic symmetric association coefficients:\n")
print(format(ass$phi[1,], digits=digits, ...), quote=FALSE)
cat("\nNormalized symmetric association scores:\n")
print(format(ass$row[,,1], digits=digits, ...), quote=FALSE)
cat("\n")
}
ass <- x$assoc.hmskew
cat("Intrinsic skew association coefficients:\n")
print(format(ass$phi[1,], digits=digits, ...), quote=FALSE)
cat("\nNormalized skew association scores:\n")
print(format(ass$row[,,1], digits=digits, ...), quote=FALSE)
if(length(ass$diag) > 0) {
cat("\nDiagonal coefficients:\n")
print(format(ass$diag[1:nrow(ass$diag),], digits=digits, ...), quote=FALSE)
}
cat("\nNormalization weights:", ass$weighting)
printModelStats(x, digits=digits)
}
print.yrcskew <- function(x, digits = max(3, getOption("digits") - 4), ...) {
cat("Call:\n", deparse(x$call), "\n", sep="", fill=TRUE)
ass <- x[["assoc"]]
if(length(ass) > 0) {
cat("Intrinsic symmetric association coefficients:\n")
print(format(ass$phi[1,], digits=digits, ...), quote=FALSE)
cat("\nNormalized symmetric association scores:\n")
print(format(ass$row[,,1], digits=digits, ...), quote=FALSE)
cat("\n")
}
ass <- x$assoc.yrcskew
cat("\nIntrinsic skew association coefficients:\n")
print(format(ass$phi[1,], digits=digits, ...), quote=FALSE)
cat("\nNormalized skew association scores:\n")
print(format(ass$row[,,1], digits=digits, ...), quote=FALSE)
if(length(ass$diag) > 0) {
cat("\nDiagonal coefficients:\n")
print(format(ass$diag[1:nrow(ass$diag),], digits=digits, ...), quote=FALSE)
}
cat("\nNormalization weights:", ass$weighting)
printModelStats(x, digits=digits)
}
print.rcL <- function(x, digits = max(3, getOption("digits") - 4), ...) {
cat("Call:\n", deparse(x$call), "\n", sep="", fill=TRUE)
ass <- x$assoc
cat("Intrinsic association coefficients:\n")
print(format(ass$phi, digits=3), quote=FALSE)
if(dim(ass$row)[3] == 1) {
cat("\nNormalized row scores (" , ass$vars[1], ") for all layers:\n", sep="")
print(format(ass$row[,,1], digits=digits, ...), quote=FALSE)
}
else {
cat("\nNormalized row scores (", ass$vars[1], "):\n", sep="")
print(format(ass$row, digits=digits, ...), quote=FALSE)
}
if(dim(ass$col)[3] == 1) {
cat("\nNormalized column scores (" , ass$vars[2], ") for all layers:\n", sep="")
print(format(ass$col[,,1], digits=digits, ...), quote=FALSE)
}
else {
cat("\nNormalized column scores (", ass$vars[2], "):\n", sep="")
print(format(ass$col, digits=digits, ...), quote=FALSE)
}
if(length(ass$diag) > 0) {
cat("\nDiagonal coefficients:\n")
print(format(ass$diag[1:nrow(ass$diag),], digits=digits, ...), quote=FALSE)
}
cat("\nNormalization weights:", ass$weighting)
printModelStats(x, digits=digits)
}
print.rcL.symm <- function(x, digits = max(3, getOption("digits") - 4), ...) {
cat("Call:\n", deparse(x$call), "\n", sep="", fill=TRUE)
ass <- x$assoc
cat("Intrinsic association coefficients:\n")
print(format(ass$phi, digits=digits, ...), quote=FALSE)
if(dim(ass$row)[3] == 1) {
cat("\nNormalized row and column scores for all layers:\n")
print(format(ass$row[,,1], digits=digits, ...), quote=FALSE)
}
else {
cat("\nNormalized row and column scores:\n")
print(format(ass$row, digits=digits, ...), quote=FALSE)
}
if(length(ass$diag) > 0) {
cat("\nDiagonal coefficients:\n")
print(format(ass$diag[1:nrow(ass$diag),], digits=digits, ...), quote=FALSE)
}
cat("\nNormalization weights:", ass$weighting)
printModelStats(x, digits=digits)
}
print.hmskewL <- function(x, digits = max(3, getOption("digits") - 4), ...) {
cat("Call:\n", deparse(x$call), "\n", sep="", fill=TRUE)
ass <- x[["assoc"]]
if(length(ass) > 0) {
cat("Intrinsic symmetric association coefficients:\n")
print(format(ass$phi, digits=digits, ...), quote=FALSE)
if(dim(ass$row)[3] == 1) {
cat("\nNormalized symmetric association scores for all layers:\n")
print(format(ass$row[,,1], digits=digits, ...), quote=FALSE)
cat("\n")
}
else {
cat("\nNormalized symmetric association scores:\n")
print(format(ass$row, digits=digits, ...), quote=FALSE)
cat("\n")
}
}
ass <- x$assoc.hmskew
cat("Intrinsic skew association coefficients:\n")
print(format(ass$phi, digits=3), quote=FALSE)
if(dim(ass$row)[3] == 1) {
cat("\nNormalized skew association scores for all layers:\n")
print(format(ass$row[,,1], digits=digits, ...), quote=FALSE)
}
else {
cat("\nNormalized skew association scores:\n")
print(format(ass$row, digits=digits, ...), quote=FALSE)
}
if(length(ass$diag) > 0) {
cat("\nDiagonal coefficients:\n")
print(format(ass$diag[1:nrow(ass$diag),], digits=digits, ...), quote=FALSE)
}
cat("\nNormalization weights:", ass$weighting)
printModelStats(x, digits=digits)
}
print.rcL.trans <- function(x, digits = max(3, getOption("digits") - 4), ...) {
cat("Call:\n", deparse(x$call), "\n", sep="", fill=TRUE)
ass <- x$assoc
cat("Transition coefficients:\n")
print(format(ass$transition, digits=3), quote=FALSE)
cat("\nIntrinsic association coefficients:\n")
print(format(ass$phi, digits=3), quote=FALSE)
cat("\nNormalized row scores (", ass$vars[1], ") for first layer:\n", sep="")
print(format(ass$row[,,1], digits=digits, ...), quote=FALSE)
cat("\nNormalized row scores (", ass$vars[1], ") for last layer:\n", sep="")
print(format(ass$row[,,dim(ass$row)[3]], digits=digits, ...), quote=FALSE)
cat("\nVariation of normalized row scores\nbetween first and layer layer:\n")
print(format(ass$row[,,dim(ass$row)[3]] - ass$row[,,1], digits=digits, ...), quote=FALSE)
cat("\nNormalized column scores (", ass$vars[2], ") for first layer:\n", sep="")
print(format(ass$col[,,1], digits=digits, ...), quote=FALSE)
cat("\nNormalized column scores (", ass$vars[2], ") for last layer:\n", sep="")
print(format(ass$col[,,dim(ass$col)[3]], digits=digits, ...), quote=FALSE)
cat("\nVariation of normalized column scores\nbetween first and layer layer:\n")
print(format(ass$row[,,dim(ass$col)[3]] - ass$col[,,1], digits=digits, ...), quote=FALSE)
if(length(ass$diag) > 0) {
cat("\nDiagonal coefficients:\n")
print(format(ass$diag[1:nrow(ass$diag),], digits=digits, ...), quote=FALSE)
}
cat("\nNormalization weights:", ass$weighting)
printModelStats(x, digits=digits)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.