Python數據分析--Numpy常用函數介紹(2)

来源:https://www.cnblogs.com/codingchen/archive/2022/05/21/16290862.html
-Advertisement-
Play Games

摘要：本篇我們將以分析歷史股價為例，介紹怎樣從文件中載入數據，以及怎樣使用NumPy的基本數學和統計分析函數、學習讀寫文件的方法，並嘗試函數式編程和NumPy線性代數運算，來學習NumPy的常用函數。 一、文件讀入 ：讀寫文件是數據分析的一項基本技能 CSV（Comma-Separated Valu ...


摘要：本篇我們將以分析歷史股價為例，介紹怎樣從文件中載入數據，以及怎樣使用NumPy的基本數學和統計分析函數、學習讀寫文件的方法，並嘗試函數式編程和NumPy線性代數運算，來學習NumPy的常用函數。

一、文件讀入 ：讀寫文件是數據分析的一項基本技能

CSV（Comma-Separated Value，逗號分隔值）格式是一種常見的文件格式。通常，資料庫的轉存文件就是CSV格式的，文件中的各個欄位對應於資料庫表中的列。

NumPy中的 loadtxt 函數可以方便地讀取CSV文件，自動切分欄位，並將數據載入NumPy數組。

1、保存或創建新文件 

import numpy as np

i = np.eye(3) #eye(n)函數創建n維單位矩陣
print(i)
np.savetxt('test.txt', i) #savetxt()創建並保存test.txt文件
savetxt()函數，如果有已經文件則更新，如目錄中沒有，則創建並保存test.txt文件

運行結果如下：

[[1. 0. 0.]
 [0. 1. 0.]
 [0. 0. 1.]]

2、讀取csv文件的函數loadtxt

1）先在保存程式的目錄下創建一個名稱為data.csv的文件，並設置數據如下圖：

2）讀取文件，如下：

c,v=np.loadtxt('data.csv', delimiter=',', usecols=(6,7), unpack=True)

usecols 的參數是一個元組，以獲取第7欄位至第8欄位的數據，也就是上述文件中 股票的收盤價和成交量數據。 unpack 參數設置為 True ，是分拆存儲不同列的數據，即分別將收盤價和成交量的數組賦值給變數c和v。

3、常見的函數

成交量加權平均、時間加權、算術平均值、中位數、方差等

import numpy as np

i = np.eye(3) #eye(n)函數創建n維單位矩陣
print(i)
np.savetxt('test.txt', i) #savetxt創建並保存test.txt文件

#讀取csv文件
c,v=np.loadtxt('data.csv', delimiter=',', usecols=(6,7), unpack=True)
"""usecols 的參數為一個元組，以獲取第7欄位至第8欄位的數據，也就是股票的收盤價和成交量數據。 unpack 參數設置為 True ，是分拆存儲不同列的數據，即分別將收
盤價和成交量的數組賦值給變數c和v"""
vwap = np.average(c, weights=v)  #調用了average函數,將v作為權重參數使用，
print(vwap)
print('\n')
print( np.mean(c)) #算術平均值
print('\n')
t = np.arange(len(c))
print( t )
print('\n')
twap =np.average(c, weights=t) #按時間權重
print( twap )
print('\n')
h,l=np.loadtxt('data.csv',delimiter=',', usecols=(4,5), unpack=True)
# 獲取第4欄位至第5欄位的數據，即股票的最高價和最低價

print ( np.max(h)) #獲取最大值max()
print ( np.min(l)) #獲取最小值min()
print('\n')
print( np.ptp(h) ) # 用ptp()函數計算了極差，即最大值和最小值之間的差值
print( np.ptp(l) )
print('\n')
print( np.median(c)) # 中位數median()函數，即多個數據中，處於中間的數
print( np.msort(c))#msort(( ))函數對價格數組進行排序,可以驗證上述中位數
#方差的計算
variance = np.var(c) #方差函數var()
print(variance)

用代碼、excel進行相關計算，運行結果如下：

為後面計算，將data.csv中的數據多增加幾行，修改如下並保存(為後面日期讀寫與修改，日期形式修改成如下)：

603112,2022-4-1,,13.56,13.97,13.55,13.87,3750000
603112,2022-4-2,,13.75,14.25,13.69,14.03,4003500
603112,2022-4-3,,13.69,14.11,13.61,13.95,3956500
603112,2022-4-4,,14.3,14.3,13.73,13.89,4250000
603112,2022-4-5,,14.1,14.5,13.93,14,4013500
603112,2022-4-6,,14.5,15.4,14.35,15.4,9056500
603112,2022-4-7,,16,16.94,15.85,16.94,3750000

4、股票的收益率等

股市中最常見的就是漲幅，也就是今日收盤價相對昨日漲跌的比例，即  （今日收盤價-昨天收盤價）/昨日收盤價*100，numpy中的 diff() 函數可以返回一個由相鄰數組元素的差值構成的數組，由於相鄰數據相減，因此diff()數組數據較原數組少一個。

如上述修改後，有7天的收盤價，diff()計算出的結果就只有6位，

import numpy as np

#讀取csv文件
c,v=np.loadtxt('data.csv', delimiter=',', usecols=(6,7), unpack=True)

#股票的簡單收益率
# diff 函數可以返回一個由相鄰數組元素的差值構成的數組
results = np.diff(c)
print(results)
print('\n')
results1 = np.diff(c)/c[:-1]*100  #相對前一天的漲幅
print(results1)
print('\n')
Standard_deviation =np.std(results) # 計算出標準差
print(Standard_deviation)

運行結果，代碼、excel進行相比較：

5、對數收益與波動率

1）對數收益：log 函數得到每一個收盤價的對數，再對結果使用 diff 函數即可，

logreturns = np.diff( np.log(c) )
print(logreturns)

運行結果：

[ 0.01146966 -0.00571839 -0.00431035  0.00788817  0.09531018  0.09531018]

2) where的作用

 where 函數可以根據指定的條件返回所有滿足條件的序列索引值，比如上述logreturns中有兩個小於0的數據。

posretindices = np.where(results1 > 0) 
print('Indices with positive returns1',posretindices)

運行結果：

Indices with positive returns1 (array([0, 3, 4, 5], dtype=int64),)

3）波動率：波動率=對數收益率的標準差除以其均值，再除以交易周期倒數的平方根。下麵代碼分別為以年、月進行統計的波動率.

annual_volatility =(np.std(logreturns)/np.mean(logreturns))/np.sqrt(1./252.)#使用浮點數才能得到正確的結果
print ( annual_volatility )
#月波動率
month_volatility =(np.std(logreturns)/np.mean(logreturns))/np.sqrt(1./12.)
print ( month_volatility )

 6、日期分析

處理日期總是很煩瑣。NumPy是面向浮點數運算的，因此需要對日期做一些專門的處理。

通過上述代碼，我們知道，修改函數np.loadtxt('data.csv', delimiter=',', usecols=(6,7), unpack=True)中的參數 usecols=(6,7)就可以讀取不同的列，日期是在第2列，即下標應該為1（數列下標是從0開始的），可以重新定義新日期數列並獲取後存入。

代碼如下：

dates, c=np.loadtxt('data.csv', delimiter=',', usecols=(1,6), unpack=True) #讀取下標為1、6的數據，分別存入到dates和c數列中。

但實際運行過程中會報錯，

代碼需要作如下修改：

import numpy as np
from datetime import datetime

def datestr2num(s): #定義一個函數
    return datetime.strptime(s.decode('ascii'),"%Y-%m-%d").date().weekday()  
#decode('ascii') 將字元串s轉化為ascii碼

#讀取csv文件
dates,close=np.loadtxt('data.csv',delimiter=',', usecols=(1,6),converters={1:datestr2num},unpack=True)
print(dates)

運行結果：[4. 5. 6. 0. 1. 2. 3.]，也是從0開始，到6結束。為了更好地說明數據，可以採用真實的數據，即從通信達軟體直接下載真實的交易數據，如下圖所示：

（註意：較原來少了一列空格列）

修改代碼如下：

import numpy as np
from datetime import datetime

def datestr2num(s): #定義一個函數
    return datetime.strptime(s.decode('ascii'),"%Y-%m-%d").date().weekday()  
#decode('ascii') 將字元串s轉化為ascii碼

#讀取csv文件
dates,c=np.loadtxt('data.csv',delimiter=',', usecols=(1,5),
                       converters={1:datestr2num},unpack=True)
print(dates)

print(len(dates)) #統計導出的天數

運行結果：

如上圖，導出有420天數據。

按周一到周五，統計相關數據：

averages = np.zeros(5) #創建包含5個元素的數組,保存交易日收盤價，0-4分別代表周一到周五五個交易日
for i in range(5):  #遍歷0到4的日期標識
    indices =np.where(dates==i)   #where函數得到各工作日的索引值並存儲在 indices 數組
    prices=np.take(c,indices)   #take函數獲取各個工作日的收盤價。
    avg= np.mean(prices) #每個工作日計算出平均值存放在 averages 數組
    averages[i] = avg  #每個工作日計算出平均值存放在 averages 數組
    print('day', i)
    #print('prices', prices)
    print("Average", avg)

print(averages)

當然除了上述外，還可以求得420天里的最大值、最小值以及交易日平均值中最大值、最小值等，對代碼進行如下修：

import numpy as np
from datetime import datetime

def datestr2num(s): #定義一個函數
    return datetime.strptime(s.decode('ascii'),"%Y-%m-%d").date().weekday()  
#decode('ascii') 將字元串s轉化為ascii碼

#讀取csv文件
dates,c=np.loadtxt('data.csv',delimiter=',', usecols=(1,5),
                       converters={1:datestr2num},unpack=True)

averages = np.zeros(5) #創建包含5個元素的數組,保存交易日收盤價，0-4分別代表周一到周五五個交易日
for i in range(5):  #遍歷0到4的日期標識
    indices =np.where(dates==i)   #where函數得到各工作日的索引值並存儲在 indices 數組
    prices=np.take(c,indices)   #take函數獲取各個工作日的收盤價。
    avg= np.mean(prices) 
    averages[i] = avg  #每個工作日計算出平均值存放在 averages 數組,共有5個數值
    print('day', i)
    #print('prices', prices)
    print("Average", avg)

print(averages)
print('\n')

print('the top close price:',np.max(c)) #最高收盤價
print('the low close price:',np.min(c)) #最低收盤價
print('\n')

top = np.max(averages)  #找出averages數列中的最大值
print ("Highest average", top)
print ("Top day of the week", np.argmax(averages)) #argmax函數返回的是averages數組中最大元素的索引值
print('\n')

bottom = np.min(averages) #找出averages數列中的最小值
print ("Lowest average", bottom)
print ( "Bottom day of the week", np.argmin(averages))#argmin函數返回的是averages數組中最小元素的索引值

運行結果如下：

 本篇初步導入了真實的股票交易信息，並利用numpy常見函數對進行了初步的計算，列舉了下列常用函數：

loadtxt() 函數可以方便地讀取CSV文件，自動切分欄位，並將數據載入NumPy數組。
savetxt()創建並保存test.txt文件
np.loadtxt('data.csv', delimiter=',', usecols=(6,7),)usecols參數用來選擇讀取的數列
np.average(c, weights=v) 加權平均,將v作為權重參數使用，
np.mean(c)) #算術平均值
np.max(h)) #獲取最大值max()
np.min(l)) #獲取最小值min()
np.ptp(h) ) 用ptp()函數計算了極值差，
np.median(c)) 中位數median()函數，即多個數據中，處於中間的數
np.msort(c))函數對價格數組進行排序,
np.var(c) 方差函數var()
np.diff(c) 函數可以返回一個由相鄰數組元素的差值構成的數組
np.std(results) # 標準差
np.diff( np.log(c) )
np.where(results1 > 0) 選擇
np.sqrt()#平方根sqrt(),浮點數
s.decode('ascii') 將字元串s轉化為ascii碼
np.take(c,indices) #take函數獲取各個工作日的收盤價。
np.argmax(averages)) #argmax函數返回數組中最大元素的索引值
np.argmin(averages))#argmin函數返回數組中最小元素的索引值

下一篇將介紹股市中常見的幾個指標線，如移動均線、布林線等。


您的分享是我們最大的動力!

-Advertisement-
Play Games
更多相關文章
  • MySQL使用步驟
    出現mysqld: Can't create directory 'D:\Environment\mysql-5.7.37 \data' (Errcode: 2 - No such file or directory)極大原因是my.ini的配置有問題 my.ini配置 [mysqld] based ...
  • 阿裡媽媽牟娜：定向廣告新一代點擊率預估主模型——深度興趣演化網路
    分享嘉賓：牟娜 阿裡巴巴 高級演算法工程師 編輯整理：孫鍇 內容來源：DataFun AI Talk《定向廣告新一代點擊率預估主模型——深度興趣演化網路》 出品社區：DataFun 導讀： 本次帶給大家分享是阿裡媽媽在2018年做的模型上的創新——深度興趣演化網路（Deep Interest Evol ...
  • 大數據Hadoop之——任務調度器Oozie（Oozie環境部署）
    一、概述 Oozie是一個基於工作流引擎的開源框架，依賴於MapReduce來實現，是一個管理 Apache Hadoop 作業的工作流調度系統。是由Cloudera公司貢獻給Apache的，它能夠提供對Hadoop MapReduce和Pig Jobs的任務調度與協調。Oozie需要部署到Java ...
  • MySQL深入學習--day3
    B+樹索引的正確使用 索引並不是越多越好，索引創建越多，MySQL維護的代價越高，如果SQL未能完全使用到索引，創建索引的意義是不大的。 適用條件 表x，創建索引a,b,c。主鍵y。 全值匹配 select * from x where a = '' and b = '' and c = '' 當我 ...
  • 徐童：視頻人物社交關係圖生成與應用
    導讀： 線上社交媒體平臺的發展，帶來了細粒度檢索、視頻語義摘要等媒體智能服務的巨大需求。現有的視頻理解技術缺乏深入的語義線索，結合視頻中人物的社交關係才能更完整、準確地理解劇情，從而提升用戶體驗，支撐智能應用。這裡主要介紹我們將動態分析和圖機器學習相結合，圍繞視頻中的人物社交關係網路所開展的兩個最新 ...
  • 處理mysql主從不同步問題
    問題描述：發現主庫操作數據從庫沒有變動問題，可能原因是從庫重啟導致的無法同步問題。 排查思路： 1、查看主從複製狀態 發現從庫的IO和SQL進程都是no(正常狀態應該是yes) 註意：mysql replication中slave機器上有兩個關鍵進程，死一個都不行，一個是slave_sql_runn ...
  • Material Desgin - Toolbar 和 滑動菜單的簡單使用
    Toolbar 新建項目會預設顯示ActionBar，ActionBar根據項目指定的主題來顯示，android:theam指定主題，在res/values/theams.xml中定義主題，因為用toolbar代替ActionBar，所以將parent主題設置為不帶ActionBar的主題，如下： ...
  • 在瀏覽器測試JavaScript的方法
    在瀏覽器測試JavaScript的方法 JSFiddle Cross Browser Testing Tool Karam + Jasmine + Google Chrome CodePen JSBin Liveweave 測試JavaScript是一件很痛苦的事情。很多工具、技術和框架已經被開發出 ...
一周排行
    -Advertisement-
    Play Games
  • Dapr Outbox 執行流程
    Dapr Outbox 是1.12中的功能。 本文只介紹Dapr Outbox 執行流程，Dapr Outbox基本用法請閱讀官方文檔 。本文中appID=order-processor，topic=orders 本文前提知識：熟悉Dapr狀態管理、Dapr發佈訂閱和Outbox 模式。 Outbo ...
  • Advanced .Net Debugging 7：托管堆與垃圾收集
    引言 在前幾章我們深度講解了單元測試和集成測試的基礎知識，這一章我們來講解一下代碼覆蓋率,代碼覆蓋率是單元測試運行的度量值,覆蓋率通常以百分比表示，用於衡量代碼被測試覆蓋的程度，幫助開發人員評估測試用例的質量和代碼的健壯性。常見的覆蓋率包括語句覆蓋率（Line Coverage）、分支覆蓋率（Bra ...
  • C#S7.NET實現西門子PLCDB塊數據採集的完整步驟
    前言 本文介紹瞭如何使用S7.NET庫實現對西門子PLC DB塊數據的讀寫，記錄了使用電腦模擬，模擬PLC，自至完成測試的詳細流程，並重點介紹了在這個過程中的易錯點，供參考。 用到的軟體： 1.Windows環境下鏈路層網路訪問的行業標準工具（WinPcap_4_1_3.exe）下載鏈接：http ...
  • WPF控制項：密碼框綁定MVVM
    從依賴倒置原則（Dependency Inversion Principle, DIP）到控制反轉（Inversion of Control, IoC）再到依賴註入（Dependency Injection, DI）的演進過程，我們可以理解為一種逐步抽象和解耦的設計思想。這種思想在C#等面向對象的編 ...
  • Python中的私有屬性與私有方法
    關於Python中的私有屬性和私有方法 Python對於類的成員沒有嚴格的訪問控制限制，這與其他面相對對象語言有區別。關於私有屬性和私有方法，有如下要點： 1、通常我們約定，兩個下劃線開頭的屬性是私有的（private）。其他為公共的（public）； 2、類內部可以訪問私有屬性（方法）； 3、類外 ...
  • C++ 訪問說明符詳解：封裝數據，控制訪問，提升安全性
    C++ 訪問說明符 訪問說明符是 C++ 中控制類成員（屬性和方法）可訪問性的關鍵字。它們用於封裝類數據並保護其免受意外修改或濫用。 三種訪問說明符： public：允許從類外部的任何地方訪問成員。 private：僅允許在類內部訪問成員。 protected：允許在類內部及其派生類中訪問成員。 示 ...
  • static_cast與dynamic_cast到底是什麼？
    寫這個隨筆說一下C++的static_cast和dynamic_cast用在子類與父類的指針轉換時的一些事宜。首先，【static_cast，dynamic_cast】【父類指針，子類指針】，兩兩一組，共有4種組合：用 static_cast 父類轉子類、用 static_cast 子類轉父類、使用 ...
  • C語言數據結構：雙向鏈表的增刪操作
    /******************************************************************************************************** * * * 設計雙向鏈表的介面 * * * * Copyright (c) 2023-2 ...
  • @ComponentScan註解的實現，Spring掃描包的過程
    相信接觸過spring做開發的小伙伴們一定使用過@ComponentScan註解 @ComponentScan("com.wangm.lifecycle") public class AppConfig { } @ComponentScan指定basePackage，將包下的類按照一定規則註冊成Be ...
  • opensips開啟python支持
    操作系統 ：CentOS 7.6_x64 opensips版本： 2.4.9 python版本：2.7.5 python作為腳本語言，使用起來很方便，查了下opensips的文檔，支持使用python腳本寫邏輯代碼。今天整理下CentOS7環境下opensips2.4.9的python模塊筆記及使用 ...
所有分類