Here's one way. I have used the data.table
package below, so if you do not have it, please install it first. If you really prefer that the plot be created using the regular data frame, let me know and I will adjust it accordingly.
library(data.table) # Install via install.packages(data.table)
# x is the vector of numbers for which you want to create the stem and leaf plot
# leftDigits gives the position of the '|' relative to the decimal point,
# e.g. leftDigits = 0 will result in 43.2 being represented as 43 | 0.2
# e.g. leftDigits = 1 will result in 43.2 being represented as 4 | 3.2
# I have included a rounding parameter, due to floating point problems seen when taking the difference of two numbers
# rounding = 2 means that your numbers will show to 2 decimal places. Change as desired.
myStem <- function(x, leftDigits, rounding = 2){
data = data.table("x" = x) # Create a data.table
data[, left := floor(x/10^leftDigits)] # Get the number that will be to the left of '|'
data[, right := round(x - left*10^leftDigits, rounding)] # The remainder will be on the right
data = data[, paste(sort(right), collapse = ", "), by = left] # For each 'left', Place all the 'right' values in a comma separated string.
data[, out := paste(left, " | ", V1), by = left] # Add the "|"
cat(data$out, sep = "\n") # Output
}
# Example
myStem(x = c(43.2, 45.3, 48.1, 54.2), 1)
# 4 | 3.2, 5.3, 8.1
# 5 | 4.2
Edit: If you would like a back to back plot, then the following code should work. It works in a similar manner as the previous code
myStem <- function(leftVals, rightVals, mainDigits, rounding = 10){
data = data.table("x" = leftVals, "ind" = "x")
data = rbind(data, data.table("x" = rightVals, "ind" = "y"))
data[, main := floor(x/10^mainDigits)]
data[, right := round(x - main*10^mainDigits, rounding)]
data = data[, ifelse(ind == "x", paste0(-sort(-right), collapse = ", "), paste0(sort(right), collapse = ", ")), by = c("ind", "main")] # For each 'main', Place all the 'right' values in a comma separated string.
data = dcast(data, main ~ ind, value.var = "V1")
data[, "left|" := ifelse(is.na(x), "", "|")]
data[, "right|" := ifelse(is.na(y), "", "|")]
data[, x := ifelse(is.na(x), "", x)]
data[, y := ifelse(is.na(y), "", y)]
data = data[, c("x", "left|", "main", "right|", "y"), with = F]
maxLengthY = max(nchar(data$y))
data[, y := unlist(lapply(y, function(z) paste0(z, paste0(replicate(maxLengthY - nchar(z), " "), collapse = ""))))]
colnames(data) = rep(" ", ncol(data))
data
}
# Example
myStem(leftVals = c(43.2, 45.3, 48.1, 54.2), rightVals = c(30.2, 34.5, 44.3), 1)
1: 3 | 0.2, 4.5
2: 8.1, 5.3, 3.2 | 4 | 4.3
3: 4.2 | 5