PSSM矩阵（位置权重矩阵）

定义

-参考维基百科
-同时参考这篇优秀的博客： http://www.nohup.cc/article/112/

可以反映出每个位置上不同碱基出现的频率，矩阵每一列表示相应位置上碱基出现的频率。
构造PSSM的第一步：通过计算每个位置上每个碱基出现的次数来创建一个基本频率矩阵（PFM)
第二步：标准化，用每个位置的原核苷酸计数除以序列数。构建位置频率矩阵
给定一个长度为l的序列集合X (N),
第三步：构建位置比重矩阵

获得PSSM矩阵
参考了以下博客：

-http://www.nohup.cc/article/115/
-https://blog.csdn.net/weixin_43405130/article/details/103716040
-http://www.nohup.cc/article/110/（CKSAAP）
-PSSM特征-从生成到处理
-如何生成PSSM矩阵
-https://blog.csdn.net/cpc784221489/article/details/88879650
-https://www.cnblogs.com/cong3Z/p/12775414.html

第一步：下载blast，使用psiblast
第二步：下载数据库，下载swissprot，网址为：

https://ftp.ncbi.nlm.nih.gov/blast/db/或者ftp://ftp.ncbi.nlm.nih.gov/blast/db/
进去后选择：swissprot，
下载链接为：https://ftp.ncbi.nlm.nih.gov/blast/db/swissprot.tar.gz或者
ftp://ftp.ncbi.nlm.nih.gov/blast/db/swissprot.tar.gz
注意下载到blast/bin路径下
解压缩：tar zxvf filename.tar

第三步：批量处理序列：

i = 0
fw = open('/blast/bin/0.txt', 'w')
for line in open('filename.fasta', 'r'):fw.write(line)i += 1if i % 2 == 0:fw.close()fw = open(str(i) + '.txt', 'w')
fw.close()

第四步：格式化数据库（可以不用做这一步）：

makeblastdb -in swissprot -parse_seqids -hash_index -dbtype prot(nulc)

第五步：批量对序列进行处理，获得其pssm矩阵：

import os
os.chdir('/blast/bin/sequence/')
for i in range(0,num(sequence),2):os.system("/blast/bin/psiblast -query /blast/bin/sequence/"+ str(i) + ".txt"+" -db /blast/swissprot -evalue 0.001 -num_iterations 3 -num_threads 38 -out /blast/bin/pssm_file/"+str(i)+".pssm")

注意-db前的空格，否则会报错“too many positional arguments”
第六步：处理得到的PSSM矩阵

pssm = []
with open('out.pssm') as f:lines = f.readlines()[3:-6]pssm = np.array([line.split()[2:22] for line in lines], dtype=int)```

PSSM矩阵（位置权重矩阵）相关推荐

stata 导出相关系数表_Stata: 空间权重矩阵的构建
作者:潘星宇 (清华大学) Stata 连享会:知乎 | 简书 | 码云 | CSDN Stata连享会计量专题 || 公众号合集连享会-内生性专题现场班-2019.11.14-17 特别说明 ...
白话空间统计之二十五：空间权重矩阵（四）R语言中的空间权重矩阵(5)完结篇：自然临近关系
自然临近是R语言中spdep中内置的最后一种临近关系. 所谓的自然临近,指的是不进行任何的预设关系,通过其空间位置来判断是否属于相互临近,那么这个空间位置指的是什么呢?众所周知,在几何图形中,三角形是 ...
空间计量 python_一文读懂空间计量入门之空间数据、权重矩阵、空间统计、空间模型等简介...
原标题:一文读懂空间计量入门之空间数据.权重矩阵.空间统计.空间模型等简介来源:计量经济学服务中心编辑整理,转载请联系 1 前言空间分析起源于地理学.空间经济学和相关学科.到目前为止,它被认为是专 ...
Python计算经济距离权重矩阵
前情回顾:空间权重矩阵一年前的博客介绍了基本权重矩阵的定义,今天在其基础上分享Python构建的方式目录一.经济权重矩阵定义二.数据准备三.全部代码一.经济权重矩阵定义当然,这里YiY_ ...
[MATLAB] 主成分分析法求矩阵的权重
主成分分析法求矩阵的权重 clc clear close all R = ...[ 0.12 0.41 0.33 0.12 0.03; ...0.14 0.52 0.21 0.10 0.03; .. ...
空间权重矩阵与相关性检验(Stata)
空间权重矩阵与相关性检验(Stata) 文章目录空间权重矩阵与相关性检验(Stata) @[toc] 1 空间相关性检验 1.1 全局空间相关性检验 1.2 局部空间自相关检验 1.3 散点图 2 ...
计算机视觉：2.3.1、梯度下降法优化权重矩阵
二.优化方法和正则方法 "Nearly all of deep learning is powered by one import algorithm :Stochastic Gradien ...
二进制空间权重矩阵_白话空间统计之二十五：空间权重矩阵（三）解构空间权重矩阵...
原标题:白话空间统计之二十五:空间权重矩阵(三)解构空间权重矩阵中国古代就有"天圆地方"一说,所谓的"方"就是所谓的矩阵-- 在军事上面,最为讲究就是团体的力 ...
二进制空间权重矩阵_空间权重矩阵（SWM）
基本原理如何利用数学(如用面积.距离等)建立空间数据和非空间数据之间的关系?对于空间统计而言,空间权重矩阵是一种有效的表达空间关系的方式.因此,它是用量化的方法表示了数据之间的"空间结构& ...

PSSM矩阵（位置权重矩阵）

PSSM矩阵（位置权重矩阵）

PSSM矩阵（位置权重矩阵）相关推荐

最新文章

热门文章