R语言如何计算GC/AT含量
发表于:2025-02-03 作者:千家信息网编辑
千家信息网最后更新 2025年02月03日,本篇文章给大家分享的是有关R语言如何计算GC/AT含量,小编觉得挺实用的,因此分享给大家学习,希望大家阅读完这篇文章后可以有所收获,话不多说,跟着小编一起来看看吧。今天给大家介绍一个可以处理FASTA
千家信息网最后更新 2025年02月03日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安全错误
数据库的锁怎样保障安全
数据库扩容量计算
规范软件开发方法
宜宾市计算机网络技术学院
如何清除服务器账号和密码
数据库使用场景 ios
手机闹钟软件开发
科技互联网简介
高校买知网数据库账务处理
网络技术型公司
软件开发可在哪些部门工作
网络安全论文博客
台州哲商网络技术有限公司
非关系型数据库简答题
世界网络安全产业的发展趋势
江西联想服务器维修调试哪家好
网络安全工程师出路友凉
虹口区一站式软件开发大概费用
科技互联网独角兽企业
用dhcp服务器管理
当兵后学习软件开发
珠海恒任互联网科技
王牌战争服务器怎么关
从前端获取数据储存到数据库
管家婆电脑数据库安装教程
有什么服务器安全防护的建议
数据库安全管理特点
无限流量用什么服务器
qq炫舞广东专区服务器
网络安全性类型怎么选
鹤壁app软件开发费用