getBranchSizes <-
function(trdata, brColors, divisions){
	if(missing(trdata) || class(trdata) != "phylo"){
		stop("A valid tree is required.")
	}
	
	edgecol <- NULL
	edgewid <- NULL
	edgelen <- trdata$edge.length
	
	noZdata <- subset(edgelen, edgelen != 0)	#get min and max without 0's in data set
	quan <- quantile(noZdata, probs = seq(0, 1, 0.2))
	
	if(missing(divisions)){ #default divisions
		divisions <- c(.1, 1, 10, 100, 1000, 10000)
	}
	sort(divisions)
	divisions <- c(quan[[1]], divisions)#always want the 0 lengths to be white
	
	if(missing(brColors)){ #default colors
		brColors <- c("red", "orange", "yellow", "green" , "cyan", "blue")
	}
	
	if(length(divisions) > (length(brColors)+1)){ #need more colors, dont care if more colors than divisons
		brColors <- c(brColors, rep(brColors[length(brColors)], (length(divisions) - length(brColors)-1)))
	}
	
	palette(brColors)
	
	for(i in (1:length(edgelen))){ 	#loop through all tree edges
		if(edgelen[as.numeric(i)] < divisions[1]){ #0 value so make it white
			edgecol <- append(edgecol, 0)
			edgewid <- append(edgewid, 0)
		}else{
			for(j in (2:length(divisions))){
				if(edgelen[as.numeric(i)] <= divisions[j]){
					cval <- (j-1)
					len <- floor(4*j/length(divisions)) 
					edgecol <- append(edgecol, cval)
					edgewid <- append(edgewid, len)
					break
				}
			}
		}
	}
	
	bs <- list(edgecol = edgecol, edgewid = edgewid)
	return(bs)
}