comparison ks_distribution.R @ 64:b7468c75c4a5 draft

Uploaded

63:99ba4a5320e0

opt <- args$options

set_component_colors = function(colors, number_comp)
{   
    # Handle colors for components.
    if (is.na(colors)) {
        # Randomly specify colors for components.
        component_colors <- c("red", "yellow", "green", "black", "blue", "darkorange")
    } else {
        # Handle selected colors for components.
        component_colors <- c()

    }
    results <- c(pi, mu, var)
    return(results)
}

plot_ks<-function(kaks_input, number_comp, component_colors, output, pi, mu, var)
{
    # Start PDF device driver to save charts to output.
    pdf(file=output, bg="white")
    kaks <- read.table(file=kaks_input, header=T)
    max_ks <- max(kaks$Ks, na.rm=TRUE)

}

# Read in the components data.
components_data <- read.delim(opt$components_input, header=TRUE)
number_comp <- as.integer(opt$number_comp)

# Set component colors.
component_colors <- set_component_colors(opt$specified_colors, number_comp)

# Set pi, mu, var.

    var <- items[11:15]
} else if (number_comp == 6) {
    pi <- items[1:6]
    mu <- items[7:12]
    var <- items[13:18]
}

# Plot the output.
plot_ks(opt$kaks_input, number_comp, component_colors, opt$output, pi, mu, var)

64:b7468c75c4a5

opt <- args$options

set_component_colors = function(colors, number_comp)
{   
    # Handle colors for components.
    if (is.null(colors)) {
        # Randomly specify colors for components.
        component_colors <- c("red", "yellow", "green", "black", "blue", "darkorange")
    } else {
        # Handle selected colors for components.
        component_colors <- c()

    }
    results <- c(pi, mu, var)
    return(results)
}

plot_ks<-function(kaks_input, component_colors, output, pi, mu, var)
{
    # Start PDF device driver to save charts to output.
    pdf(file=output, bg="white")
    kaks <- read.table(file=kaks_input, header=T)
    max_ks <- max(kaks$Ks, na.rm=TRUE)

}

# Read in the components data.
components_data <- read.delim(opt$components_input, header=TRUE)
number_comp <- as.integer(opt$number_comp)
if (number_comp == 0) {
    # Default to 1 component to allow functional testing.
    number_comp <- 1
}

# Set component colors.
component_colors <- set_component_colors(opt$specified_colors, number_comp)

# Set pi, mu, var.

    var <- items[11:15]
} else if (number_comp == 6) {
    pi <- items[1:6]
    mu <- items[7:12]
    var <- items[13:18]

# Plot the output.
plot_ks(opt$kaks_input, component_colors, opt$output, pi, mu, var)