How do I plot the maximum value point in R?

To plot the maximum value point in R, you can use the points() function along with the which.max() function to locate the maximum value and its coordinates. Here’s a step-by-step explanation and example:

1. Create a Sample Dataset:

- First, you need some data to work with. Let's create a simple vector or data frame:

# Sample data vector
data_vector <- c(5, 12, 8, 20, 15, 9)

# Sample data frame (with x and y coordinates)
data_frame <- data.frame(x = 1:6, y = c(5, 12, 8, 20, 15, 9))

2. Find the Index of the Maximum Value:

- Use which.max() to find the index (position) of the maximum value in the vector or the 'y' column in the data frame.

# For vector data
max_index_vector <- which.max(data_vector)

# For data frame
max_index_df <- which.max(data_frame$y)

3. Extract the Coordinates of the Maximum Value:

- Get the actual x and y coordinates based on the identified index.

# For vector data, x coordinate is the index, y is the value
max_x_vector <- max_index_vector
max_y_vector <- data_vector[max_index_vector]

# For data frame, x and y are directly from the data frame
max_x_df <- data_frame$x[max_index_df]
max_y_df <- data_frame$y[max_index_df]

4. Create the Plot:

- Use the plot() function to generate the base plot.

# Plotting for vector (using index as x)
plot(1:length(data_vector), data_vector, type = "l", main = "Plot with Max Value Highlight", xlab = "Index", ylab = "Value")

# Plotting for data frame (using x and y columns)
plot(data_frame$x, data_frame$y, type = "l", main = "Plot with Max Value Highlight", xlab = "X", ylab = "Y")

5. Add the Maximum Value Point:

- Use points() to highlight the maximum value's coordinate on the plot.

# For vector data
points(max_x_vector, max_y_vector, col = "red", pch = 16, cex = 1.5)

# For data frame
points(max_x_df, max_y_df, col = "red", pch = 16, cex = 1.5)

6. Complete Example:

# Example with vector
data_vector <- c(5, 12, 8, 20, 15, 9)
max_index_vector <- which.max(data_vector)
max_x_vector <- max_index_vector
max_y_vector <- data_vector[max_index_vector]
plot(1:length(data_vector), data_vector, type = "l", main = "Max Value Point in Vector", xlab = "Index", ylab = "Value")
points(max_x_vector, max_y_vector, col = "red", pch = 16, cex = 1.5)

# Example with data frame
data_frame <- data.frame(x = 1:6, y = c(5, 12, 8, 20, 15, 9))
max_index_df <- which.max(data_frame$y)
max_x_df <- data_frame$x[max_index_df]
max_y_df <- data_frame$y[max_index_df]
plot(data_frame$x, data_frame$y, type = "l", main = "Max Value Point in Data Frame", xlab = "X", ylab = "Y")
points(max_x_df, max_y_df, col = "red", pch = 16, cex = 1.5)

This code first creates a plot with the data, then it locates the maximum value and its x and y coordinates using which.max(), and then finally plots a red circle at that point on the graph using points().

Learn more about this