代码之家 › 专栏 › 技术社区 › neo4k

Python中的渐近回归?

modeling curve-fitting statistics machine-learning python

neo4k · 技术社区 · 9 年前

我有 this y 值最终必须达到渐近极限。以下是前300个条目的曲线图:

我看到了 asymptotic regression y = b0 + b1*(1-exp(-exp(lrc) * x) )我认为这将是适合这些数据的最佳模型。R似乎已经有了 function 在它的图书馆里,但我在蟒蛇

(1) 是否存在以下Cleb所回答的渐近回归或相关饱和动力学模型的库函数?如果没有,如何使用scipy.optimize.curve_fit对其建模?

(2) 另外,如何预测曲线的最大整数值?示例,如果 y=2454.1234 是处的值 X=10**20 y=2545.5678 在 X=10**50 ,对于某个模型,我想获得 2454

谢谢你的帮助。

1 回复 | 直到 9 年前

Cleb 9 年前

scipy's curve_fit function .我使用了不同的模型,但你可以很容易地改变它。

这是一段非常不言自明的代码:

import numpy as np
from scipy.optimize import curve_fit
import matplotlib.pyplot as plt
import pandas as pd


def f(x, a, b, n):
    return a * x ** n  / (x ** n + b)


data = pd.read_csv('data.txt.txt', sep='\t')

y = data['y'].astype(float)
x = data['X'].astype(float)

popt, pcov = curve_fit(f, x, y, p0=[1800., 20., 1.])

plt.scatter(x, y)
plt.plot(x, f(x, *popt), 'r-')

plt.show()

推荐文章

Google User · Django管理员在`list_display中未显示`creation_date`字段`

1 年前

user29747013 · 如何创建一个新的数据框架,其中包含原始数据框架中列的聚合列?

1 年前

ÎÎÎ½Î· ÎÎ®Î¹Î½Î¿Ï · Python lxml.html语法错误:使用lxml find时XPATH的谓词无效

1 年前

user29715306 · from_users=和chats=电视节目中的差异

1 年前

Redshoe · 当执行numpy.genfromtxt()时,python是否会读取文件的所有行?

1 年前

RASEL MAHMUD · 为什么以及如何在is_even()函数内的IF条件中递归X变量在满足0后递增?[副本]

1 年前

prayner · 更新嵌套字典包含列表中的项

1 年前

Bringo Jr · 我可以在O(n)中解决这个问题吗?

1 年前

Dave · 如何在for循环中修改列表值

1 年前

Shukurullox Komiljonov · 从记录中获得相互和解。使用SQL

1 年前