博文
统计碱基组成的R语言小程序Rjuan 1.04（增加比率统计）
|||
统计碱基组成的R语言小程序Rjuan 1.04（增加比率统计）
熊荣川
六盘水师范学院生物信息学实验室
http://blog.sciencenet.cn/u/Bearjazz
之前我们所开发的一段R语言代码——Rjuan 1.0可用于一段序列的各个碱基数量的统计。（详见“统计碱基组成的R语言小程序Rjuan 1.0” http://bbs.sciencenet.cn/home.php?mod=space&uid=508298&do=blog&id=553110）
我们的第一次升级解决了输入序列长度限制的问题
http://bbs.sciencenet.cn/home.php?mod=space&uid=508298&do=blog&id=553147
第二次升级在程序的后面加上一段具有打印功能的代码，并将信息保存为矩阵，并与输出为表格文件
http://bbs.sciencenet.cn/home.php?mod=space&uid=508298&do=blog&id=553163
第三次升级对核心代码部分进行函数化封装，这样的在处理新的输入序列时，只需要使用调用函数的方式就可以输出结果了。同时也避免了每次都要进行计数器清零操作。
http://bbs.sciencenet.cn/home.php?mod=space&uid=508298&do=blog&id=553282
现在，我们对程序进行升级，增加各碱基比率统计并输出。
|
|
nucount=function(x){ na=0 nt=0 ng=0 nc=0 for (i in 1:length(x)) ( if(x[i]=="A") na=na+1 else if(x[i]=="T") nt=nt+1 else if(x[i]=="G") ng=ng+1 else if(x[i]=="C") nc=nc+1 ) sum=sum(c(na,nt,ng,nc)) RateA=na/sum RateT=nt/sum RateG=ng/sum RateC=nc/sum sumR=RateA+RateT+RateG+RateC namen=c("na","nt","ng","nc","sum") numbern=c(na,nt,ng,nc,sum) Rate=c(RateA,RateT,RateG,RateC,sumR) shown=data.frame(namen,numbern,Rate) print(shown) } |
函数后端增加各个碱基比率的统计功能并输出 |
|
> |
data<-read.csv("D:\ziliao\zhuanye\R bear\test3.csv") |
导入序列 |
|
> |
length(data[,1]) [1] 17962 |
序列长度17962bp |
|
|
z=data[,1] |
序列向量化 |
|
> |
> nucount(z) namen numbern Rate 1 na 5111 0.2845452 2 nt 5272 0.2935085 3 ng 2667 0.1484801 4 nc 4912 0.2734662 5 sum 17962 1.0000000 > |
调用函数nucount（） 统计碱基数及所占比率 |
https://blog.sciencenet.cn/blog-508298-553293.html
上一篇：统计碱基组成的R语言小程序Rjuan 1.03 （核心代码函数化封装）
下一篇：统计碱基组成的R语言小程序Rjuan 1.20（序列输入格式.fastr）
