[go: up one dir, main page]
Scribd
Upload
92 views5 pages

R Optimization and Plotting Guide

The document describes using R code to solve an optimization problem graphically and numerically using linear programming. It: 1) Defines constraints using functions and plots the constraint lines along with axes breaks and labels to visualize the feasible region. 2) Adds the feasible region polygon and optimal solution point to the plot. 3) Defines the objective function and overlays it on the plot. 4) Uses the lpSolve package to solve the linear program numerically, outputting the optimal solution values.

Uploaded by

Trivikraman Srinivasan
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
92 views5 pages

R Optimization and Plotting Guide

The document describes using R code to solve an optimization problem graphically and numerically using linear programming. It: 1) Defines constraints using functions and plots the constraint lines along with axes breaks and labels to visualize the feasible region. 2) Adds the feasible region polygon and optimal solution point to the plot. 3) Defines the objective function and overlays it on the plot. 4) Uses the lpSolve package to solve the linear program numerically, outputting the optimal solution values.

Uploaded by

Trivikraman Srinivasan
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 5

# Define constraints

cons.1 <- function(x) 6 - x

cons.2 <- function(x) 5 - 0.5*x

# cons.3 x1 = 4 (plotted using geom_vline)

# cons.4 x1 = 0 (defined by setting x1 axis limits)

# cons.5 x2 = 0 (defined by setting x2 axis limits)

# Import ggplot2 package

library(ggplot2)

# Build plot

p = ggplot(data = data.frame(x = 0), aes(x = x)) +

# Add axes

geom_vline(xintercept = 0) +

geom_hline(yintercept = 0) +

# Add constraints lines

stat_function(colour = "Red", fun = cons.1) +

stat_function(colour = "Blue", fun = cons.2) +

geom_vline(xintercept = 4, colour = "Green") +

# Specify axes breaks and limits

scale_x_continuous(breaks = seq(0, 10, 1), lim = c(0, 10)) +

scale_y_continuous(breaks = seq(0, 10, 1), lim = c(0, 10)) +

# Define labels

labs(title = "Optimization Problem",

subtitle = "Graphical Method",

x = "x1",

y = "x2") +
# Add black and white theme

theme_bw()

# Print plot

print (p)

# Define feasible region polygon

feasible_region = data.frame(x = c(0, 4, 4, 2, 0), y = c(0, 0, 2, 4, 5))

# Add feasible region to current plot

p <- p + geom_polygon(data = feasible_region, mapping = aes(x = x, y = y), fill = "#dddddd") +

geom_point(data = feasible_region, aes(x = x, y = y), color = "Black")

# Print plot

print(p)

# Define objective function

z1 <- function(x) -0.8*x + 2 # (the + 2 constant is just for displacement visualization purposes)

# Visualize objective function displacement

p <- p + stat_function(colour = "Black", fun = z1, lty=2)

# Print plot

print(p)

library(lpSolve)

# Set coefficients of the objective function

f.obj <- c(4, 5)


# Set matrix corresponding to coefficients of constraints by rows

f.con <- matrix(c(1, 2,

6, 6,

1, 0), nrow = 3, byrow = TRUE)

# Set unequality/equality signs

f.dir <- c("<=", "<=", "<=")

# Set right hand side coefficients

f.rhs <- c(10, 36, 4)

# Final value (z)

lp("max", f.obj, f.con, f.dir, f.rhs)

# Variables final values

solution <- lp("max", f.obj, f.con, f.dir, f.rhs)$solution

# Highlight optimum solution in plot

p <- p + geom_point(aes(x = solution[1], y = solution[2]), color = "red", size = 4)

# Print plot

print(p)
Code for north west corner rule in R

First install TransP packages


install.packages(TransP)
Rcode
library(TransP)
#Input matrix where last row is the Demand and last
column is the Supply
ex_matrix=data.frame(M1=c(13,10,25,17,210),M2=c(25,
19,10,24,240),

M3=c(8,18,15,18,110),M4=c(13,5,14,13,80),M5=c(20,12
,18,19,170),
Supply=c(430,150,100,130,810),
row.names =
c("W1","W2","W3","W4","Demand"))

>nwc(ex_matrix)
2) One more problem in North West corner
rule

You might also like