最新国产好看的视频,伊人天堂AV在线,国产Aaaaaa视频,蜜臀视频在线观看一区,人妻av色图,密臀久久久精品影片,青青视频免费观看毛片,久草在线观看视,国产三级精品色情在线

python計(jì)算RPKM操作示例詳解

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

操作

np.rot90(df, 1) #矩陣,逆時(shí)針,90度
np.rot90(df, -1) #矩陣,順時(shí)針,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)實(shí)現(xiàn)矩陣除數(shù)組,類似R語(yǔ)言(按行除),

即每一行列數(shù)組中每一個(gè)個(gè)數(shù),按列除array[:,None]將數(shù)組豎著排實(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ù)點(diǎn)后保留20位

計(jì)算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)矩陣法,過(guò)于復(fù)雜不推薦
#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計(jì)算RPKM操作示例詳解的詳細(xì)內(nèi)容,更多關(guān)于python計(jì)算RPKM的資料請(qǐng)關(guān)注腳本之家其它相關(guān)文章!

相關(guān)文章

  • Python?中的?Counter?模塊及使用詳解(搞定重復(fù)計(jì)數(shù))

    Python?中的?Counter?模塊及使用詳解(搞定重復(fù)計(jì)數(shù))

    Counter 是一個(gè)簡(jiǎn)單的計(jì)數(shù)器,用于統(tǒng)計(jì)某些可哈希對(duì)象的數(shù)量。它以字典的形式存儲(chǔ)元素和它們的計(jì)數(shù),這篇文章主要介紹了Python?中的?Counter?模塊及使用詳解(搞定重復(fù)計(jì)數(shù)),需要的朋友可以參考下
    2023-04-04
  • python中最小二乘法詳細(xì)講解

    python中最小二乘法詳細(xì)講解

    在本篇文章里小編給大家整理的是一篇關(guān)于python中最小二乘法詳細(xì)講解內(nèi)容,需要的朋友們可以參考下。
    2021-02-02
  • Django框架封裝外部函數(shù)示例

    Django框架封裝外部函數(shù)示例

    這篇文章主要介紹了Django框架封裝外部函數(shù),結(jié)合Django框架表單登陸功能示例分析了封裝外部函數(shù)的相關(guān)操作步驟與實(shí)現(xiàn)技巧,需要的朋友可以參考下
    2019-05-05
  • Python利用PyMobileDevice3控制iOS設(shè)備的完整教程

    Python利用PyMobileDevice3控制iOS設(shè)備的完整教程

    PyMobileDevice3是一個(gè)純Python 3實(shí)現(xiàn)的iOS設(shè)備控制工具庫(kù),為開(kāi)發(fā)者和技術(shù)愛(ài)好者提供了強(qiáng)大的iOS設(shè)備管理能力,所以本文給大家介紹了Python利用PyMobileDevice3控制iOS設(shè)備的完整教程,需要的朋友可以參考下
    2025-12-12
  • python爬蟲(chóng)入門教程之糗百圖片爬蟲(chóng)代碼分享

    python爬蟲(chóng)入門教程之糗百圖片爬蟲(chóng)代碼分享

    這篇文章主要介紹了python爬蟲(chóng)入門教程之糗百圖片爬蟲(chóng)代碼分享,本文以抓取糗事百科內(nèi)涵圖為需求寫了一個(gè)爬蟲(chóng),,需要的朋友可以參考下
    2014-09-09
  • 一篇文章搞懂python混亂的切換操作與優(yōu)雅的推導(dǎo)式

    一篇文章搞懂python混亂的切換操作與優(yōu)雅的推導(dǎo)式

    這篇文章主要給大家介紹了如何通過(guò)一篇文章搞懂python混亂的切換操作與優(yōu)雅的推導(dǎo)式的相關(guān)資料,文中通過(guò)示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下
    2021-08-08
  • Python類定義和類繼承詳解

    Python類定義和類繼承詳解

    這篇文章主要介紹了Python類定義和類繼承詳解,本文講解了類的私有屬性、類的方法、私有的類方法、類的專有方法、類的定義、類的單繼承、類的多繼承等內(nèi)容,需要的朋友可以參考下
    2015-05-05
  • 基于Python制作簡(jiǎn)單的井字棋游戲

    基于Python制作簡(jiǎn)單的井字棋游戲

    這篇文章主要為大家詳細(xì)介紹了如何基于Python實(shí)現(xiàn)簡(jiǎn)單的井字棋游戲,文中示例代碼介紹的非常詳細(xì),具有一定的參考價(jià)值,感興趣的小伙伴們可以參考一下
    2023-04-04
  • python的構(gòu)建工具setup.py的方法使用示例

    python的構(gòu)建工具setup.py的方法使用示例

    本篇文章主要介紹了python的構(gòu)建工具setup.py的方法示例,小編覺(jué)得挺不錯(cuò)的,現(xiàn)在分享給大家,也給大家做個(gè)參考。一起跟隨小編過(guò)來(lái)看看吧
    2017-10-10
  • Python學(xué)習(xí)之直方圖均衡化原理詳解

    Python學(xué)習(xí)之直方圖均衡化原理詳解

    直方圖均衡化是以累計(jì)分布函數(shù)為核心,將原始圖像灰度直方圖從比較集中的某個(gè)灰度區(qū)間,非線性地映射為在全部灰度范圍內(nèi)的較均勻分布,從而增強(qiáng)對(duì)比度。本文將為大家詳細(xì)講解直方圖均衡化的原理,需要的可以參考一下
    2022-03-03

最新評(píng)論

馆陶县| 胶州市| 山阳县| 扎赉特旗| 沿河| 万安县| 化德县| 石屏县| 常州市| 玉龙| 龙里县| 志丹县| 阳城县| 汶川县| 虞城县| 马尔康县| 丹寨县| 寿宁县| 东光县| 定边县| 光山县| 阿尔山市| 湘潭县| 大同市| 威宁| 邹城市| 沽源县| 琼海市| 靖江市| 苍溪县| 云南省| 花莲县| 揭西县| 磐石市| 桃源县| 康平县| 平武县| 库尔勒市| 绍兴县| 新田县| 两当县|