Advertisement

Python Statistical model uses Scipy's curve_fit for curve fitting

阅读量:

点赞发Nature
关注中Science

近期正在学习一些与计量相关的知识,我发现他人通常采用State或者Eviews来进行计量分析,我也尝试使用Python进行操作,以下展示的是利用Python完成的基础线性回归分析过程。

复制代码
    from scipy.optimize import curve_fit
    
    def cost_func(x_data, c1):
    return c1 * x_data   #目标函数形式
    
    def cost_curve_fit(x_data, y_data):
    
    c, cov = curve_fit(cost_func, x_data, y_data)
    
    n = len(x_data)
    fit_y = np.empty(n)
    for i in range(n):
        fit_y[i] = cost_func(x_data[i], c[0])
    return c, cov, fit_y 
    # c, cov, fit_y分别表示你和参数,拟合方差,拟合曲线
    
    
      
      
      
      
      
      
      
      
      
      

全部评论 (0)

还没有任何评论哟~