欧美bbbwbbbw肥妇,免费乱码人妻系列日韩,一级黄片

python計算RPKM操作示例詳解

 更新時間:2023年07月06日 08:42:02   作者:小白菜學生信  
這篇文章主要為大家介紹了python計算RPKM操作示例詳解,有需要的朋友可以借鑒參考下,希望能夠有所幫助,祝大家多多進步,早日升職加薪

操作

np.rot90(df, 1) #矩陣,逆時針,90度
np.rot90(df, -1) #矩陣,順時針,90度

矩陣/數(shù)組,使用循環(huán)速度特慢

#!/usr/bin python3
# -*- coding: UTF-8 -*-
import os,sys,re
import pandas as pd
df=pd.read_csv("final.head", index_col = 0, header = 0, sep = "\t")
length=pd.read_csv("final.length", index_col = 0, header = 0, sep = "\t")
for i in range(len(df.index)):
    for j in range(len(df.columns)):
        df.iloc[i,j]=df.iloc[i,j]/length.iloc[i,0]
df.to_csv('final.norm', sep='\t', index = True)

numpy不用循環(huán)實現(xiàn)矩陣除數(shù)組,類似R語言(按行除),

即每一行列數(shù)組中每一個個數(shù),按列除array[:,None]將數(shù)組豎著排實現(xiàn)按行除或用np.rot(df,)旋轉(zhuǎn)矩陣,不推薦,容易轉(zhuǎn)暈

#!/usr/bin/env python3
import os,re,sys
import numpy as np
ms, df, length, outfile = sys.argv
df = np.loadtxt(df, dtype=np.int32, delimiter='\t')
length = np.loadtxt(length, dtype=np.int32)
out=df/length[:,None]  # [:,None]將數(shù)組豎著排
np.savetxt(outfile,out,fmt='%.8f',delimiter='\t')
#np.savetxt(outfile,out,fmt='%.20f',delimiter='\t') # 小數(shù)點后保留20位

計算rpkm

#!/usr/bin/env python3
import os,re,sys
import numpy as np
ms, df, length, mapped_reads, outfile = sys.argv
# df[gene,sample]
df = np.loadtxt(df, delimiter='\t')
length = np.loadtxt(length)
mapped_reads = np.loadtxt(mapped_reads)
#np.dtype=int32 # 這是整形32位,不適用
out=1e3*1e6*df/(length[:,None]*mapped_reads)
#使用旋轉(zhuǎn)矩陣法,過于復雜不推薦
#out=1e3*1e6*np.rot90((np.rot90(df,1)/length),-1)/mapped_reads
np.savetxt(outfile,out,fmt='%.8f',delimiter='\t')

行列求和

#!/usr/bin python
# -*- coding: UTF-8 -*-
import pandas as pd
#t=[[1,2,3],[4,5,6],[7,8,9]]
#df=pd.DataFrame(t) # List轉(zhuǎn)為dataframe
df=pd.read_csv("test.df", index_col = 0, header = 0, sep = "\t")
total=df.apply(lambda x: x.sum())
total.to_csv('test.total', sep='\t', index = True)
#df['Row_sum'] = df.apply(lambda x: x.sum(),axis=1) # 按行求和,添加為新列
#df.loc['Col_sum'] = df.apply(lambda x: x.sum()) # 各列求和,添加新的行

以上就是python計算RPKM操作示例詳解的詳細內(nèi)容,更多關于python計算RPKM的資料請關注腳本之家其它相關文章!

相關文章

  • python自動化測試三部曲之request+django實現(xiàn)接口測試

    python自動化測試三部曲之request+django實現(xiàn)接口測試

    這篇文章主要介紹了python自動化測試三部曲之request+django實現(xiàn)接口測試,文中通過示例代碼介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友們下面隨著小編來一起學習學習吧
    2020-10-10
  • python之pyinstaller組件打包命令和異常解析實戰(zhàn)

    python之pyinstaller組件打包命令和異常解析實戰(zhàn)

    前段時間在制作小工具的時候,直接在命令行用pyinstaller工具打包成功后,啟動exe可執(zhí)行文件的時候各種報錯, 今天,我們就分享一下踩坑經(jīng)過,需要的朋友可以參考下
    2021-09-09
  • Python 的 with 語句詳解

    Python 的 with 語句詳解

    這篇文章主要介紹了Python 的 with 語句,本文詳細講解了with語句、with語句的歷史、with語句的使用例子等,需要的朋友可以參考下
    2014-06-06
  • 使用Python對SQLite數(shù)據(jù)庫操作

    使用Python對SQLite數(shù)據(jù)庫操作

    本文主要介紹了Python對SQLite數(shù)據(jù)庫操作的簡單教程。SQLite是一種嵌入式數(shù)據(jù)庫,它的數(shù)據(jù)庫就是一個文件。由于SQLite本身是C寫的,而且體積很小,所以,經(jīng)常被集成到各種應用程序中,甚至在IOS和Android的APP中都可以集成。
    2017-04-04
  • python使用pygame實現(xiàn)笑臉乒乓球彈珠球游戲

    python使用pygame實現(xiàn)笑臉乒乓球彈珠球游戲

    這篇文章主要為大家詳細介紹了python使用pygame實現(xiàn)笑臉乒乓球彈珠球游戲,文中示例代碼介紹的非常詳細,具有一定的參考價值,感興趣的小伙伴們可以參考一下
    2019-11-11
  • 解決python 未發(fā)現(xiàn)數(shù)據(jù)源名稱并且未指定默認驅(qū)動程序的問題

    解決python 未發(fā)現(xiàn)數(shù)據(jù)源名稱并且未指定默認驅(qū)動程序的問題

    今天小編就為大家分享一篇解決python 未發(fā)現(xiàn)數(shù)據(jù)源名稱并且未指定默認驅(qū)動程序的問題,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧
    2018-12-12
  • 使用Python讀取Excel數(shù)據(jù)并寫入到CSV、XML和文本

    使用Python讀取Excel數(shù)據(jù)并寫入到CSV、XML和文本

    Excel工作簿是常用的表格格式,許多數(shù)據(jù)呈現(xiàn)、數(shù)據(jù)分析和數(shù)據(jù)匯報都是以Excel工作表的形式進行,本文將演示如何運用Python編程語言,將Excel工作表中的豐富數(shù)據(jù)導入到CSV、XML或文本中,需要的朋友可以參考下
    2024-03-03
  • python機器學習MATLAB最小二乘法的兩種解讀

    python機器學習MATLAB最小二乘法的兩種解讀

    這篇文章主要為大家介紹了python機器學習中MATLAB最小二乘法的兩種解讀方式,有需要的朋友可以借鑒參考下希望能夠有所幫助
    2022-02-02
  • Python程序打包exe報錯的幾種解決方法

    Python程序打包exe報錯的幾種解決方法

    本文主要介紹了Python程序打包exe報錯的幾種解決方法,文中通過幾種解決方法的介紹非常詳細,具有一定的參考價值,感興趣的小伙伴們可以參考一下
    2023-08-08
  • 使用flask如何獲取post請求參數(shù)

    使用flask如何獲取post請求參數(shù)

    近日在使用flask框架獲取前端的請求時獲取參數(shù)時,遇到了幾個問題,所以下面這篇文章主要給大家介紹了關于使用flask如何獲取post請求參數(shù)的相關資料,需要的朋友可以參考下
    2022-08-08

最新評論