R语言如何计算GC/AT含量
发表于:2025-02-19 作者:千家信息网编辑
千家信息网最后更新 2025年02月19日,本篇文章给大家分享的是有关R语言如何计算GC/AT含量,小编觉得挺实用的,因此分享给大家学习,希望大家阅读完这篇文章后可以有所收获,话不多说,跟着小编一起来看看吧。今天给大家介绍一个可以处理FASTA
千家信息网最后更新 2025年02月19日R语言如何计算GC/AT含量
运行,结果如下:
本篇文章给大家分享的是有关R语言如何计算GC/AT含量,小编觉得挺实用的,因此分享给大家学习,希望大家阅读完这篇文章后可以有所收获,话不多说,跟着小编一起来看看吧。
今天给大家介绍一个可以处理FASTA文件的包-Biostrings。这个包主要是处理基因组的一些序列信息,包括:序列翻译、DNA/RNA互转、统计各个碱基的含量、三连字母的含量.....这些都是一行命令可以解决的。今天就先来教大家怎样计算GC/AT含量。
首先是安装,代码如下:
source("http://bioconductor.org/biocLite.R")biocLite("Biostrings")
输入代码后需要耐心地等待几分钟。

安装完毕,只需敲几行代码,就可以实现GC/AT含量可视化。
#序列文件储存路径filepath<-"C:/Users/dell/Desktop/sequence.fasta"#读取文件(FASTA格式)x<-readDNAStringSet(filepath)chrom<-x[[1]]#每100个碱基为窗口计算AT含量at<-rowSums(letterFrequencyInSlidingView(chrom,100,c("A","T")))/100#获取描述性统计量根据此设置坐标summary(at)#画图plot(at,type='l',axes=F,xlab=NA,ylab=NA,ylim=c(0.2,0.8))axis(2,at=c(0.2,0.4,0.6,0.8),labels=c("20%","40%","60%","80%"))#纵坐标设置根据summary函数计算结果axis(1,at=c(0,2000,4000,6000,8000,10000,12000,14434),labels=c("Start","2000","4000","6000","8000","10000","12000","End"))#根据基因组显示横坐标信息
以上就是R语言如何计算GC/AT含量,小编相信有部分知识点可能是我们日常工作会见到或用到的。希望你能通过这篇文章学到更多知识。更多详情敬请关注行业资讯频道。
含量
代码
序列
文件
语言
信息
基因
基因组
更多
知识
碱基
篇文章
结果
处理
统计
实用
一行
函数
只需
命令
数据库的安全要保护哪些东西
数据库安全各自的含义是什么
生产安全数据库录入
数据库的安全性及管理
数据库安全策略包含哪些
海淀数据库安全审计系统
建立农村房屋安全信息数据库
易用的数据库客户端支持安全管理
连接数据库失败ssl安全错误
数据库的锁怎样保障安全
检察院网络安全自查工作
数据库连表和join区别
nginx 缓存服务器
长春软件开发有名吉网传媒优
公安局网络技术类公务员考试
软件开发常见面试问题回答
2016 顶级服务器
山西家业互联网科技有限公司
大公司做金融软件开发
东西湖区 网络安全中心
支付成功后数据库更新失败
报刊索引数据库帐号
地图存储 图数据库
网络安全挑战有哪些
nat代理服务器
数据库字段拼接 c
如何修改数据库的用户名字
计算机网络技术 专业特色
协议软件开发部
一般一台服务器承受多少并发数
fm2019巨型数据库
亳州市孩纸荚网络技术服务部
浙江软件开发上市公司
海南电商软件开发价格
软件开发毕业设计目录
2019网络安全公益
对于企业 网络安全法
您已与服务器断开连接
中普服务器取数
软件开发 应该转什么行业