首页 > 解决方案 > GEKKO多元非线性回归

问题描述

df = pd.read_csv("data.csv")

xm1 = np.array(df["T"]) #Dep Var 1
xm2 = np.array(df["t"])  #Dep Var 2
xm3 = np.array(df["L"]) #Dep Var 3
ym = np.array(df["S"])  #Indep Var

# GEKKO model
m = GEKKO()
a = m.FV(lb=-100.0,ub=100.0)
b = m.FV(lb=-100.0,ub=100.0)
c = m.FV(lb=-100.0,ub=100.0)
d = m.FV(lb=-100.0,ub=100.0)
e = m.FV(lb=-100.0,ub=100.0)
f = m.FV(lb=-100.0,ub=100.0)
g = m.FV(lb=-100.0,ub=100.0)
x1 = m.Param(value=xm1)
x2 = m.Param(value=xm2)
x3 = m.Param(value=xm3)
z = m.Param(value=ym)
y = m.Var()
m.Equation(y == a+x1*b+x2*c+x3*d+e*(x1**2)+f*(x2**2)+g*(x3**2)
m.Obj(((y-z)/z)**2)

我开始 SyntaxError: invalid syntaxm.Obj(((y-z)/z)**2)。为此,我遵循了 APMonitor.com 的代码。该代码非常适合那里的示例。但是当我修改它以适应更多 FV 的回归问题时,它会显示此语法错误。附上语法错误图片

不太确定这有什么问题。任何帮助,将不胜感激。

标签: pythonnon-linear-regressiongekko

解决方案


由于此行而发生错误:

m.Equation(y == a+x1*b+x2*c+x3*d+e*(x1**2)+f*(x2**2)+g*(x3**2)

y声明前有多余的括号


推荐阅读