91超碰碰碰碰久久久久久综合_超碰av人澡人澡人澡人澡人掠_国产黄大片在线观看画质优化_txt小说免费全本

溫馨提示×

溫馨提示×

您好,登錄后才能下訂單哦!

密碼登錄×
登錄注冊×
其他方式登錄
點擊 登錄注冊 即表示同意《億速云用戶服務條款》

Python數據可視化:冪律分布實例詳解

發布時間:2020-09-22 18:55:37 來源:腳本之家 閱讀:545 作者:墨竹 | kevinelstri 欄目:開發技術

1、公式推導

對冪律分布公式:

Python數據可視化:冪律分布實例詳解

對公式兩邊同時取以10為底的對數:

Python數據可視化:冪律分布實例詳解

Python數據可視化:冪律分布實例詳解

所以對于冪律公式,對X,Y取對數后,在坐標軸上為線性方程。

2、可視化

從圖形上來說,冪律分布及其擬合效果:

Python數據可視化:冪律分布實例詳解

對X軸與Y軸取以10為底的對數。效果上就是X軸上1與10,與10與100的距離是一樣的。

Python數據可視化:冪律分布實例詳解

對XY取雙對數后,坐標軸上點可以很好用直線擬合。所以,判定數據是否符合冪律分布,只需要對XY取雙對數,判斷能否用一個直線很好擬合就行。常見的直線擬合效果評估標準有擬合誤差平方和、R平方。

3、代碼實現

#!/usr/bin/env python
# -*-coding:utf-8 -*-

import matplotlib.pyplot as plt
import numpy as np
from sklearn import linear_model
from scipy.stats import norm

def DataGenerate():
 X = np.arange(10, 1010, 10) # 0-1,每隔著0.02一個數據 0處取對數,會時負無窮 生成100個數據點
 noise=norm.rvs(0, size=100, scale=0.2) # 生成50個正態分布 scale=0.1控制噪聲強度
 Y=[]
 for i in range(len(X)):
  Y.append(10.8*pow(X[i],-0.3)+noise[i]) # 得到Y=10.8*x^-0.3+noise

 # plot raw data
 Y=np.array(Y)
 plt.title("Raw data")
 plt.scatter(X, Y, color='black')
 plt.show()

 X=np.log10(X) # 對X,Y取雙對數
 Y=np.log10(Y)
 return X,Y

def DataFitAndVisualization(X,Y):
 # 模型數據準備
 X_parameter=[]
 Y_parameter=[]
 for single_square_feet ,single_price_value in zip(X,Y):
  X_parameter.append([float(single_square_feet)])
  Y_parameter.append(float(single_price_value))

 # 模型擬合
 regr = linear_model.LinearRegression()
 regr.fit(X_parameter, Y_parameter)
 # 模型結果與得分
 print('Coefficients: \n', regr.coef_,)
 print("Intercept:\n",regr.intercept_)
 # The mean square error
 print("Residual sum of squares: %.8f"
  % np.mean((regr.predict(X_parameter) - Y_parameter) ** 2)) # 殘差平方和

 # 可視化
 plt.title("Log Data")
 plt.scatter(X_parameter, Y_parameter, color='black')
 plt.plot(X_parameter, regr.predict(X_parameter), color='blue',linewidth=3)

 # plt.xticks(())
 # plt.yticks(())
 plt.show()

if __name__=="__main__":
 X,Y=DataGenerate()
 DataFitAndVisualization(X,Y)

Python數據可視化:冪律分布實例詳解

Python數據可視化:冪律分布實例詳解

以上這篇Python數據可視化:冪律分布實例詳解就是小編分享給大家的全部內容了,希望能給大家一個參考,也希望大家多多支持億速云。

向AI問一下細節

免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。

AI

瓦房店市| 德清县| 阳曲县| 寻甸| 元氏县| 香格里拉县| 旺苍县| 奉化市| 新龙县| 新晃| 武威市| 兴安盟| 乐昌市| 鸡西市| 县级市| 大悟县| 天祝| 五大连池市| 成安县| 南汇区| 栾城县| 连州市| 南靖县| 临清市| 庄河市| 喀喇| 宜良县| 深圳市| 大城县| 贵德县| 聊城市| 博乐市| 色达县| 红桥区| 富民县| 阿城市| 临高县| 安图县| 河源市| 玉门市| 焉耆|