Repository 'shm_csr'
hg clone https://toolshed.g2.bx.psu.edu/repos/davidvanzessen/shm_csr

Changeset 25:8d1c4c75f81b (2016-12-05)
Previous changeset 24:9955b23db68c (2016-12-02) Next changeset 26:ce25fb581ca3 (2016-12-05)
Commit message:
Uploaded
modified:
shm_csr.r
b
diff -r 9955b23db68c -r 8d1c4c75f81b shm_csr.r
--- a/shm_csr.r Fri Dec 02 05:51:13 2016 -0500
+++ b/shm_csr.r Mon Dec 05 08:26:01 2016 -0500
[
@@ -462,6 +462,12 @@
 
 frequency_bins_data_by_class = frequency_bins_data
 
+frequency_bins_data_by_class = frequency_bins_data_by_class[order(frequency_bins_data_by_class$best_match_class, frequency_bins_data_by_class$frequency_bins),]
+
+frequency_bins_data_by_class$frequency_bins = gsub("-", " to ", frequency_bins_data_by_class$frequency_bins)
+frequency_bins_data_by_class[frequency_bins_data_by_class$frequency_bins == "20", c("frequency_bins")] = "20 or higher"
+frequency_bins_data_by_class[frequency_bins_data_by_class$frequency_bins == "0", c("frequency_bins")] = "0 or lower"
+
 write.table(frequency_bins_data_by_class, "frequency_ranges_classes.txt", sep="\t",quote=F,row.names=F,col.names=T)
 
 frequency_bins_data = data.frame(data.table(dat)[, list(frequency_count=.N), by=c("best_match", "best_match_class", "frequency_bins")])
@@ -472,6 +478,11 @@
 
 frequency_bins_data$frequency = round(frequency_bins_data$frequency_count / frequency_bins_data$class_sum * 100, 2)
 
+frequency_bins_data = frequency_bins_data[order(frequency_bins_data$best_match, frequency_bins_data$frequency_bins),]
+frequency_bins_data$frequency_bins = gsub("-", " to ", frequency_bins_data$frequency_bins)
+frequency_bins_data[frequency_bins_data$frequency_bins == "20", c("frequency_bins")] = "20 or higher"
+frequency_bins_data[frequency_bins_data$frequency_bins == "0", c("frequency_bins")] = "0 or lower"
+
 write.table(frequency_bins_data, "frequency_ranges_subclasses.txt", sep="\t",quote=F,row.names=F,col.names=T)