diff --git a/R/SRCGrob.R b/R/SRCGrob.R
index 832f127..38d8653 100644
--- a/R/SRCGrob.R
+++ b/R/SRCGrob.R
@@ -26,7 +26,11 @@ SRCGrob <- function(
     label.nodes = TRUE,
     disable.polygons = FALSE,
     length.from.node.edge = TRUE,
-    size.units = 'npc'
+    size.units = 'npc',
+    scale.bar = FALSE,
+    scale.bar.coords = c(0.5, 1),
+    scale.size.1 = NA,
+    scale.size.2 = NA
     ) {
 
     add.node.text <- !is.null(node.text);
@@ -107,7 +111,11 @@ SRCGrob <- function(
         main = main,
         main.cex = main.cex,
         main.y = main.y,
-        size.units = size.units
+        size.units = size.units,
+        scale.bar = scale.bar,
+        scale.bar.coords = scale.bar.coords,
+        scale.size.1 = scale.size.1,
+        scale.size.2 = scale.size.2
         );
 
     out.tree <- gTree(
diff --git a/R/calculate.clone.polygons.R b/R/calculate.clone.polygons.R
index 1caed88..7ef9061 100644
--- a/R/calculate.clone.polygons.R
+++ b/R/calculate.clone.polygons.R
@@ -331,7 +331,6 @@ compute.clones <- function(
 	v <- v[is.na(v$parent) | v$parent != -1, ];
 	v <- rbind(root, v);
 	v <- count.leaves.per.node(v);
-
 	if (no.ccf) {
 	    tree$angle <- if ((is.null(fixed.angle) && nrow(v) > 6) || any(table(v$parent) > 2)) {
     		tau <- -(pi / 2.5);
diff --git a/R/make.clone.tree.grobs.R b/R/make.clone.tree.grobs.R
index 0ccee07..219225c 100644
--- a/R/make.clone.tree.grobs.R
+++ b/R/make.clone.tree.grobs.R
@@ -37,42 +37,46 @@ make.clone.tree.grobs <- function(
     main.cex,
     main.y,
     size.units,
+    scale.bar,
+    scale.bar.coords,
+    scale.size.1,
+    scale.size.2,
     ...
     ) {
 
-	#initializing dataframe for subclones
-	if ('excluded' %in% colnames(ccf.df)) {
-		v <- ccf.df[!ccf.df$excluded,];
-	} else {
-		v <- ccf.df;
-		v$excluded <- FALSE;
-	    }
-
-	v <- v[order(v$id), ];
-	no.ccf <- FALSE;
-
-	if (!('ccf' %in% colnames(ccf.df)) || all(is.na(ccf.df$ccf)) || add.polygons == FALSE) {
-	    v$vaf <- NULL;
-	    v$vaf[v$parent == -1] <- 1;
-	    no.ccf <- TRUE;
-	} else {
-	    v <- v[order(v$id),]
-	    v$vaf[!v$excluded] <- v$ccf[!v$excluded] / max(v$ccf[!v$excluded]);
-	    }
-
-	if (all(is.null(ccf.df$colour))) {
-	    v$colour <- node.col
-	    }
+    #initializing dataframe for subclones
+    if ('excluded' %in% colnames(ccf.df)) {
+        v <- ccf.df[!ccf.df$excluded,];
+    } else {
+        v <- ccf.df;
+        v$excluded <- FALSE;
+        }
+
+    v <- v[order(v$id), ];
+    no.ccf <- FALSE;
+
+    if (!('ccf' %in% colnames(ccf.df)) || all(is.na(ccf.df$ccf)) || add.polygons == FALSE) {
+        v$vaf <- NULL;
+        v$vaf[v$parent == -1] <- 1;
+        no.ccf <- TRUE;
+    } else {
+        v <- v[order(v$id),]
+        v$vaf[!v$excluded] <- v$ccf[!v$excluded] / max(v$ccf[!v$excluded]);
+        }
+
+    if (all(is.null(ccf.df$colour))) {
+        v$colour <- node.col
+        }
 
     extra.len <- if (no.ccf) node.radius else node.radius * 4;
 
-	v$x <- v$y <- v$len <- v$x.mid <- numeric(length(nrow(v)));
-	v <- v[order(v$tier, v$parent), ];
+    v$x <- v$y <- v$len <- v$x.mid <- numeric(length(nrow(v)));
+    v <- v[order(v$tier, v$parent), ];
 
-	#initializing line segment dataframe and adjusting lengths to accomodate the node circles
-	tree$angle <- numeric(length = nrow(tree));
-	tree$angle[tree$parent == -1] <- 0;
-	if ('length2' %in% colnames(tree)) {
+    #initializing line segment dataframe and adjusting lengths to accomodate the node circles
+    tree$angle <- numeric(length = nrow(tree));
+    tree$angle[tree$parent == -1] <- 0;
+    if ('length2' %in% colnames(tree)) {
         tree$length2.c <- tree$length2 / scale1 * scale2;
 
         tree$length <- apply(
@@ -82,103 +86,123 @@ make.clone.tree.grobs <- function(
                 max(x[c(3, 6)]);
                 }
             );
-	} else {
-	    tree$length <- tree$length1;
-	    }
-
-	if (length.from.node.edge == TRUE) {
-	    tree <- adjust.branch.lengths(v, tree, node.radius, scale1);
-    	}
-
-	extra.len <- extra.len * (1 / scale1);
-
-	clone.out <- make.clone.polygons(
-	    v,
-	    tree,
-	    wid,
-	    scale1,
-	    scale2,
-	    extra.len,
-	    node.col,
-	    spread = spread,
-	    sig.shape = sig.shape,
-	    fixed.angle = fixed.angle,
-	    add.polygons = add.polygons,
-	    no.ccf = no.ccf
-	    );
-
-	clone.out$no.ccf <- no.ccf;
-	plot.size <- calculate.main.plot.size(
-	    clone.out,
-	    scale1,
-	    wid,
-	    min.width,
-	    node.radius
-	    );
-
-	if (!no.ccf) {
-		get.CP.polygons(clone.out);
-    	}
-
-	add.tree.segs(clone.out, node.radius, default.branch.width, scale1, seg1.col, seg2.col);
-
-	if (!is.null(cluster.list)) {
-	    message(paste(
-	        'Clustered pie nodes will be supported in a future version.',
-	        'Plain nodes will be used.'
-	        ));
-	    # TODO Implement pie nodes
-		# add.pie.nodes(clone.out, node.radius, cluster.list);
-    	}
+    } else {
+        tree$length <- tree$length1;
+        }
+
+    if (length.from.node.edge == TRUE) {
+        tree <- adjust.branch.lengths(v, tree, node.radius, scale1);
+        }
+
+    extra.len <- extra.len * (1 / scale1);
+
+    clone.out <- make.clone.polygons(
+        v,
+        tree,
+        wid,
+        scale1,
+        scale2,
+        extra.len,
+        node.col,
+        spread = spread,
+        sig.shape = sig.shape,
+        fixed.angle = fixed.angle,
+        add.polygons = add.polygons,
+        no.ccf = no.ccf
+        );
+
+    clone.out$no.ccf <- no.ccf;
+    plot.size <- calculate.main.plot.size(
+        clone.out,
+        scale1,
+        wid,
+        min.width,
+        node.radius
+        );
+
+    if (!no.ccf) {
+        get.CP.polygons(clone.out);
+        }
+
+    add.tree.segs(clone.out, node.radius, default.branch.width, scale1, seg1.col, seg2.col);
+
+    if (!is.null(cluster.list)) {
+        message(paste(
+            'Clustered pie nodes will be supported in a future version.',
+            'Plain nodes will be used.'
+            ));
+        # TODO Implement pie nodes
+        # add.pie.nodes(clone.out, node.radius, cluster.list);
+        }
 
     add.node.ellipse(clone.out,node.radius, label.nodes, label.cex, scale1);
 
-	if (add.normal == TRUE) {
-		add.normal(clone.out,node.radius,label.cex, normal.cex)
-	    }
-
-	if (yaxis.position != 'none') {
-		add.axes(
-		    clone.out,
-		    yaxis.position,
-		    scale1 = scale1,
-		    scale2 = scale2,
-		    yat = yat,
-		    axis.label.cex = axis.label.cex,
-		    axis.cex = axis.cex,
-		    no.ccf = no.ccf,
-		    xaxis.label = xaxis.label,
-		    yaxis1.label = yaxis1.label,
-		    yaxis2.label = yaxis2.label
-		    );
-	    }
-
-	if (add.node.text == TRUE & !is.null(text.df)) {
-	    node.text.grobs <- add.text2(
-	        clone.out$tree,
-	        text.df,
-	        label.nodes = text.on.nodes,
-	        line.dist = node.text.line.dist,
-	        main.y = clone.out$height,
-	        panel.height = clone.out$height,
-	        panel.width = clone.out$width,
-	        xlims = clone.out$xlims,
-	        ymax = clone.out$ymax,
-	        cex = node.text.cex,
-	        v = clone.out$v,
-	        axis.type = yaxis.position,
-	        node.radius = node.radius,
-	        scale = scale1,
-	        clone.out = clone.out,
-	        alternating = FALSE
-	        );
-
-	    clone.out$grobs <- c(clone.out$grobs, list(node.text.grobs));
-	    }
-
-	if (!is.null(main)) {
-		add.main(clone.out, main, main.cex, main.y, size.units);
-    	}
-
-	return(clone.out);
+    if (add.normal == TRUE) {
+        add.normal(clone.out,node.radius,label.cex, normal.cex)
+        }
+
+    if (yaxis.position != 'none' & scale.bar == FALSE) {
+        add.axes(
+            clone.out,
+            yaxis.position,
+            scale1 = scale1,
+            scale2 = scale2,
+            yat = yat,
+            axis.label.cex = axis.label.cex,
+            axis.cex = axis.cex,
+            no.ccf = no.ccf,
+            xaxis.label = xaxis.label,
+            yaxis1.label = yaxis1.label,
+            yaxis2.label = yaxis2.label
+            );
+        }
+
+    if (scale.bar) {
+        scale.lengths <- prep.scale.length(
+            tree,
+            scale.size.1,
+            scale.size.2
+            );
+
+        add.scale.bar(
+            clone.out,
+            scale1,
+            scale2,
+            yaxis1.label = yaxis1.label,
+            yaxis2.label = yaxis2.label,
+            scale.length = scale.lengths,
+            main.cex = axis.label.cex$y,
+            label.cex = axis.cex$y,
+            pos = scale.bar.coords
+            );
+        }
+
+    if (add.node.text == TRUE & !is.null(text.df)) {
+        node.text.grobs <- add.text2(
+            clone.out$tree,
+            text.df,
+            label.nodes = text.on.nodes,
+            line.dist = node.text.line.dist,
+            main.y = clone.out$height,
+            panel.height = clone.out$height,
+            panel.width = clone.out$width,
+            xlims = clone.out$xlims,
+            ymax = clone.out$ymax,
+            cex = node.text.cex,
+            v = clone.out$v,
+            axis.type = yaxis.position,
+            node.radius = node.radius,
+            scale = scale1,
+            clone.out = clone.out,
+            alternating = FALSE
+            );
+
+        clone.out$grobs <- c(clone.out$grobs, list(node.text.grobs));
+        }
+
+    if (!is.null(main)) {
+        add.main(clone.out, main, main.cex, main.y, size.units);
+        }
+
+    return(clone.out);
     }
diff --git a/R/scale.bar.R b/R/scale.bar.R
new file mode 100644
index 0000000..071239b
--- /dev/null
+++ b/R/scale.bar.R
@@ -0,0 +1,150 @@
+prep.scale.length <- function(
+    tree,
+    scale.size.1,
+    scale.size.2
+    ) {
+
+    scale.lengths <- c(scale.size.1, scale.size.2);
+    tree.lengths <- c(
+        auto.scale.length(tree$length1),
+        if ('length2' %in% names(tree)) auto.scale.length(tree$length2) else NA
+        );
+
+    # if scale.length is NA, replace with tree.lengths
+    scale.lengths[is.na(scale.lengths)] <- tree.lengths[is.na(scale.lengths)];
+
+    return(scale.lengths);
+    }
+
+auto.scale.length <- function(edge.lengths) {
+    scale.length <- median(edge.lengths[edge.lengths > 0], na.rm = TRUE);
+    adjusted.length <- 10 ** floor(log10(as.numeric(scale.length)));
+    return(adjusted.length);
+    }
+
+add.scale.bar <- function(
+    clone.out,
+    scale.length,
+    scale1,
+    scale2,
+    yaxis1.label,
+    yaxis2.label,
+    pos,
+    ...
+    ) {
+
+    # Generate the first scale bar
+    scale.bar1.glist <- create.scale.bar(
+        main = yaxis1.label,
+        scale.length = list(
+            label = scale.length[1],
+            length = scale.length[1]
+            ),
+        edge.col = most.common.value(clone.out$v$edge.colour.1),
+        edge.width = most.common.value(clone.out$v$edge.width.1),
+        edge.type = most.common.value(clone.out$v$edge.type.1),
+        left.x = pos[1],
+        top.y = pos[2],
+        ...
+        );
+    clone.out$grobs <- c(clone.out$grobs, scale.bar1.glist);
+
+    # Create second scalebar if specified
+    if (!is.null(yaxis2.label)) {
+        scale.bar2.glist <- create.scale.bar(
+            main = yaxis2.label,
+            scale.length = list(
+                label = scale.length[2],
+                length = scale.length[2] * scale1 / scale2
+                ),
+            edge.col = most.common.value(clone.out$v$edge.colour.2),
+            edge.width = most.common.value(clone.out$v$edge.width.2),
+            edge.type = most.common.value(clone.out$v$edge.type.2),
+            left.x = pos[1],
+            top.y = pos[2] + 0.1,
+            ...
+            );
+        clone.out$grobs <- c(clone.out$grobs, scale.bar2.glist);
+        }
+    }
+
+most.common.value <- function(x) {
+    if (is.null(x)) {
+        return(NULL);
+        }
+    n.table <- table(x);
+    return(names(n.table)[which.max(n.table)]);
+    }
+
+create.scale.bar <- function(
+        main,
+        scale.length,
+        left.x,
+        top.y,
+        edge.col,
+        edge.width,
+        edge.type,
+        main.cex,
+        label.cex
+        ) {
+
+    edge.width <- unit(edge.width, "points");
+    left.x <- unit(left.x - 1, "npc");
+    top.y <- unit(top.y, "npc");
+    xat <- left.x + unit(c(0, scale.length$length), 'native');
+    main.size <- unit(main.cex * 12, 'points');
+    label.size <- unit(label.cex * 12, 'points');
+
+    title <- textGrob(
+        label = main,
+        x = left.x + unit(scale.length$length / 2, 'native'),
+        y = top.y,
+        hjust = 0.5,
+        vjust = 1,
+        gp = gpar(
+            cex = main.cex
+            )
+        );
+
+    scale.bar.y <- top.y - (main.size * 2)
+    scale.line <- segmentsGrob(
+        x0 = xat[1],
+        x1 = xat[2],
+        y0 = scale.bar.y,
+        y1 = scale.bar.y,
+        gp = gpar(
+            col = edge.col,
+            lwd = edge.width,
+            lty = edge.type,
+            lineend = "butt"
+            )
+        );
+
+    tick.length <- edge.width + (label.size / 4);
+    ticks <- segmentsGrob(
+        x0 = xat,
+        x1 = xat,
+        y0 = scale.bar.y + (edge.width / 2.5),
+        y1 = scale.bar.y - tick.length,
+        default.units = 'native',
+        gp = gpar(
+            lineend = "butt"
+            )
+        );
+
+    tick.labels <- textGrob(
+        label = c(0, scale.length$label),
+        x = xat,
+        y = scale.bar.y - tick.length * 2,
+        gp = gpar(
+            cex = label.cex
+            )
+        );
+
+    return(gList(
+        title,
+        scale.line,
+        ticks,
+        tick.labels
+        ));
+    }
\ No newline at end of file