沉闷科学的掘墓人分享 http://blog.sciencenet.cn/u/Bearjazz

博文

将fastr格式序列转化为整体字体串(使用R语言)

已有 3182 次阅读 2012-3-31 14:22 |个人分类:我的研究|系统分类:科研笔记| R语言, fastr格式, 序列转化, 字体串

fastr格式序列转化为整体字体串(使用R语言)

 

 

熊荣川

六盘水师范学院生物信息学实验室

xiongrongchuan@126.com

http://blog.sciencenet.cn/u/Bearjazz

 

之前我们写的很多R语言程序在处理序列输入时,推荐较为简单的fastr格式

http://bbs.sciencenet.cn/home.php?mod=space&uid=508298&do=blog&id=553655

但是这种序列导入R语言平台是,由于对每一行字符数的限制,整体序列被分割成很多单元(每单元最多70个字符),如果直接进行字符串操作则会由于这些分割而不能顺利运算。这里我们通过一段代码将fastr格式的序列成功连接成一个整体字符串。

 

data <- readLines("D:\ziliao\zhuanye\R bear\isk2.fastr")

导入fastr格式序列

x=data

data赋值个x,避免由于运算而变化data内容,再次需要的时候还有输入那段冗长命令

y=x[1]

初始化我们的y值,使得y为一个合适的字符变量

for (i in 2:length(x))

{

y=paste(y, x[i], sep = "")

}

x的第二个单元开始逐一合并到y上,最终成为一个整体

length(y)

[1] 1

检验y的整体性

nchar(y)

[1] 349

查看y的完成性

 

 



https://blog.sciencenet.cn/blog-508298-553858.html

上一篇:统计碱基组成的R语言小程序Rjuan 1.20(序列输入格式.fastr)
下一篇:在基因组中搜索目标片段(使用R语言)
收藏 IP: 210.75.236.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-6-30 06:07

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部