R/covCor.R

In tera-insights/gtStats: Higher Level Statistical Functions for Grokit

Covariance <- cov.data <- function(data, inputs, outputs = result) {
  if (missing(inputs))
    inputs <- convert.schema(data$schema)
  else
    inputs <- substitute(inputs)
  inputs <- convert.exprs(inputs)

  outputs <- substitute(outputs)
  check.atts(outputs)
  outputs <- convert.atts(outputs)
  if (length(outputs) != 1)
    stop("There must be exactly one output specified.")

  agg <- Aggregate(data, GLA(statistics::CovCor_Matrix, which = "cov"), inputs, outputs)

  ## if (exists("grokit.jobid") && !force.frame) {
  ##   View(agg)
  ## } else {
  ##   result <- as.object(agg)$content[[1]][[1]]
  ##   terms <- unlist(lapply(grokit$expressions[inputs], deparse))
  ##   cov <- matrix(result$covariance$data, result$covariance$n_cols, result$covariance$n_rows)
  ##   rownames(cov) <- colnames(cov) <- terms
  ##   cov
  ## }
}

Correlation <- function(data, inputs, outputs = result) {
  if (missing(inputs))
    inputs <- convert.schema(data$schema)
  else
    inputs <- substitute(inputs)
  inputs <- convert.exprs(inputs)

  outputs <- substitute(outputs)
  check.atts(outputs)
  outputs <- convert.atts(outputs)
  if (length(outputs) != 1)
    stop("There must be exactly one output specified.")

  agg <- Aggregate(data, GLA(statistics::CovCor_Matrix, which = "cor"), inputs, outputs)

  ## if (exists("grokit.jobid") && !force.frame) {
  ##   View(agg)
  ## } else {
  ##   result <- as.object(agg)$content[[1]][[1]]
  ##   terms <- unlist(lapply(grokit$expressions[inputs], deparse))
  ##   cor <- matrix(result$correlation$data, result$correlation$n_cols, result$correlation$n_rows)
  ##   rownames(cor) <- colnames(cor) <- terms
  ##   cor
  ## }
}

CovCor <- function(data, inputs, outputs = result) {
  if (missing(inputs))
    inputs <- convert.schema(data$schema)
  else
    inputs <- substitute(inputs)
  inputs <- convert.exprs(inputs)

  outputs <- substitute(outputs)
  check.atts(outputs)
  outputs <- convert.atts(outputs)
  if (length(outputs) != 1)
    stop("There must be exactly one output specified.")

  agg <- Aggregate(data, GLA(statistics::CovCor_Matrix, which = "both"), inputs, outputs)

  ## if (exists("grokit.jobid") && !force.frame) {
  ##   View(agg)
  ## } else {
  ##   result <- as.object(agg)$content[[1]][[1]]
  ##   terms <- unlist(lapply(grokit$expressions[inputs], deparse))
  ##   cov <- matrix(result$covariance$data, result$covariance$n_cols, result$covariance$n_rows)
  ##   cor <- matrix(result$correlation$data, result$correlation$n_cols, result$correlation$n_rows)
  ##   rownames(cov) <- colnames(cov) <- rownames(cor) <- colnames(cor) <- terms
  ##   list(cov = cov, cor = cor)
  ## }
}