Блог

Модули Python доступны в R через reticulate пакет. К сожалению, это не работает для GEKKO, формулы не распознаются как таковые.

Я получаю такие сообщения об ошибках при попытке передать уравнения, как в m$Equation(x1 x2 == 7) : Error in x1 x2 : non-numeric argument to binary operator .

 library(reticulate)
gekko <- reticulate::import("gekko")
m <- gekko$GEKKO(remote=TRUE) 
x1 = m$Var(value=1, lb=1, ub=5)
x2 = m$Var(value=5, lb=1, ub=5)
x3 = m$Var(value=5, lb=1, ub=5)
x4 = m$Var(value=1, lb=1, ub=5)
m$Equation(x1   x2 == 7)
m$Equation(x1 * x2 * x3 * x4 >= 25)
m$Obj(x1 * x4 * (x1   x2   x3)   x3)  
m$solve(disp=TRUE)  
 

Я бы очень хотел использовать GEKKO в R, можем ли мы что-нибудь изменить, чтобы сделать это возможным?

Повторная публикация вопроса из проблемы GitHub в репозитории Gekko, чтобы узнать, можно ли найти дополнительную помощь здесь. Также возможен интерфейс Python Gekko в MATLAB, если это поможет.

 % start Matlab from Anaconda prompt
close all; clear;
% Solve linear equations
% Initialize Model
m = py.gekko.GEKKO();
% Initialize Variables
x = m.Var();            % define new variable
y = m.Var();            % default=0
% Define Equations
m.Equation(3*x 2*y==1);
m.Equation(x 2*y==0);  
% Solve
m.solve();
% Extract values from Python lists using curly brackets
disp(['x: ' num2str(x.VALUE{1})]);
disp(['y: ' num2str(y.VALUE{1})]);
 

Одним из способов было бы запустить код Python в R с помощью команды py_run_string() . Я разбил код, используя несколько строк, чтобы показать шаги.

 Sys.setenv(RETICULATE_PYTHON = "set the path for python")

library(reticulate)
#Load the package 
str1 = " 
from gekko import GEKKO
m = GEKKO()"

py_run_string(str1)

#Assign values with lower and upper bound along with equations
str2 = "
x1 = m.Var(value=1, lb=1, ub=5)
x2 = m.Var(value=5, lb=1, ub=5)
x3 = m.Var(value=5, lb=1, ub=5)
x4 = m.Var(value=1, lb=1, ub=5)

m.Equation(x1   x2 == 7)
m.Equation(x1 * x2 * x3 * x4 >= 25)
m.Obj(x1 * x4 * (x1   x2   x3)   x3)
m.solve(disp=True)"

py_run_string(str2)

#Print the results
str3="
print('')
print('Results')
print('x1: '   str(x1.value))
print('x2: '   str(x2.value))
print('x3: '   str(x3.value))
print('x4: '   str(x4.value))
"
py_run_string(str3)

#To access an object
py$x1$VALUE