Skip to main content

Using the predict() methods of fitted models with gekko

Many model-fitting Python packages have a predict() method, which outputs a prediction of the fitted model given observations of the predictor(s).

Question: How would I use these predict() methods to predict a single value when the observation is a variable in a gekko model?

Below is a very simple reproducible example:

Note: The actual model that I am fitting is a cubic spline using statsmodels.gam.smooth_basis.BSplines and statsmodels.gam.generalized_additive_model.GLMGam. However, I am hoping that this simple example with sklearn.linear_model.LinearRegression will translate to more complex model classes from other packages.

from sklearn.linear_model import LinearRegression
import numpy as np
import matplotlib.pyplot as plt
from gekko import GEKKO

# create example data
x = np.arange(100)[:, np.newaxis]
y = np.arange(100) * 2 + 10

plt.plot(x, y)  # plot x vs y data
# plt.show()

model = LinearRegression()  # instantiate linear model
model.fit(x, y)  # fit model

x_predict = np.arange(100, 200)[:, np.newaxis]  # create array of predictor observations
y_predict = model.predict(x_predict)  # use model to make prediction

plt.plot(x_predict, y_predict)  # plot prediction
# plt.show()

m = GEKKO()  # instantiate gekko model

x2 = m.FV()  # instantiate free variable
x2.STATUS = 1  # make variable available for solver

y2 = 50  # true value

# place x2 variable in numpy array to adhere to `predict()`'s argument requirements
x2_arr = np.array(x2).reshape(1, -1)

# minimize squared error between the true value and the model's prediction
m.Minimize((y2 - model.predict(x2_arr)) ** 2)

m.options.IMODE = 3
m.solve(disp=True)

print(f"x2 = {x2.value[0]:3f}")

I get the following sequential errors:

TypeError: float() argument must be a string or a real number, not 'GK_FV'

ValueError: setting an array element with a sequence.

My first thought is that I would have to create a wrapper class around the gekko.gk_parameter.GK_FV class to modify the float() method, but that's where my knowledge and skills end.



source https://stackoverflow.com/questions/76408801/using-the-predict-methods-of-fitted-models-with-gekko

Comments

Popular posts from this blog

How to show number of registered users in Laravel based on usertype?

i'm trying to display data from the database in the admin dashboard i used this: <?php use Illuminate\Support\Facades\DB; $users = DB::table('users')->count(); echo $users; ?> and i have successfully get the correct data from the database but what if i want to display a specific data for example in this user table there is "usertype" that specify if the user is normal user or admin i want to user the same code above but to display a specific usertype i tried this: <?php use Illuminate\Support\Facades\DB; $users = DB::table('users')->count()->WHERE usertype =admin; echo $users; ?> but it didn't work, what am i doing wrong? source https://stackoverflow.com/questions/68199726/how-to-show-number-of-registered-users-in-laravel-based-on-usertype

Why is my reports service not connecting?

I am trying to pull some data from a Postgres database using Node.js and node-postures but I can't figure out why my service isn't connecting. my routes/index.js file: const express = require('express'); const router = express.Router(); const ordersCountController = require('../controllers/ordersCountController'); const ordersController = require('../controllers/ordersController'); const weeklyReportsController = require('../controllers/weeklyReportsController'); router.get('/orders_count', ordersCountController); router.get('/orders', ordersController); router.get('/weekly_reports', weeklyReportsController); module.exports = router; My controllers/weeklyReportsController.js file: const weeklyReportsService = require('../services/weeklyReportsService'); const weeklyReportsController = async (req, res) => { try { const data = await weeklyReportsService; res.json({data}) console...

Confusion between commands.Bot and discord.Client | Which one should I use?

Whenever you look at YouTube tutorials or code from this website there is a real variation. Some developers use client = discord.Client(intents=intents) while the others use bot = commands.Bot(command_prefix="something", intents=intents) . Now I know slightly about the difference but I get errors from different places from my code when I use either of them and its confusing. Especially since there has a few changes over the years in discord.py it is hard to find the real difference. I tried sticking to discord.Client then I found that there are more features in commands.Bot . Then I found errors when using commands.Bot . An example of this is: When I try to use commands.Bot client = commands.Bot(command_prefix=">",intents=intents) async def load(): for filename in os.listdir("./Cogs"): if filename.endswith(".py"): client.load_extension(f"Cogs.{filename[:-3]}") The above doesnt giveany response from my Cogs ...