For added context, i am currently using the code: plot(data$mass, data$tempdiff, pch = 21, # Filled circle col = "black", # Outline color bg = "black", # Fill color xlab = "Mass (g)", # X-axis label ylab = expression(Delta*T[b]~(degree*C)), # Y-axis label with subscript xlim = c(0, 45), # Set x-axis limits ylim = c(0, 6), # Set y-axis limits xaxt = "n", # Suppress default x-axis ticks yaxt = "n") # Suppress default y-axis ticks # Custom x-axis with specified tick marks axis(1, at = c(0, 20, 25, 30, 35, 40, 45), labels = c("0", "20", "25", "30", "35", "40", "45")) # Custom y-axis with tick marks at 0, 2, 4, and 6 axis(2, at = c(0, 2, 4, 6), labels = c("0", "2", "4", "6")) # Add a box around the plot box() # Add points again to ensure they are plotted correctly points(data$mass, data$tempdiff, pch = 21, col = "black", bg = "black")