Mercurial > repos > greg > ks_distribution

changeset 65:ad406d4702b2 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author greg
date Tue, 27 Jun 2017 11:15:45 -0400
parents b7468c75c4a5
children 02de4eb54f27
files ks_distribution.R
diffstat 1 files changed, 5 insertions(+), 8 deletions(-) [+]
line wrap: on
line diff
--- a/ks_distribution.R	Tue Jun 27 11:00:36 2017 -0400
+++ b/ks_distribution.R	Tue Jun 27 11:15:45 2017 -0400
@@ -14,8 +14,7 @@
 args <- parse_args(parser, positional_arguments=TRUE)
 opt <- args$options
 
-set_component_colors = function(colors, number_comp)
-{   
+set_component_colors = function(colors, number_comp) {
     # Handle colors for components.
     if (is.null(colors)) {
         # Randomly specify colors for components.
@@ -55,8 +54,7 @@
     return(component_colors)
 }
 
-get_pi_mu_var = function(components_data, number_comp)
-{
+get_pi_mu_var = function(components_data, number_comp) {
     if (number_comp == 1) {
         pi <- c(components_data[1, 9])
         mu <- c(components_data[1, 7])
@@ -86,8 +84,7 @@
     return(results)
 }
 
-plot_ks<-function(kaks_input, component_colors, output, pi, mu, var)
-{
+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)
@@ -116,8 +113,7 @@
     }
 }
 
-calculate_fitted_density <- function(pi, mu, var, max_ks)
-{
+calculate_fitted_density <- function(pi, mu, var, max_ks) {
     comp <- length(pi)
     var <- var/mu^2
     mu <- log(mu)
@@ -172,6 +168,7 @@
     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)