> 

data<-read.csv("D:\ziliao\zhuanye\R bear\test2.csv")

导入test2.csv数据存在矩阵data中

> 

length(data[,1])

[1] 602

查看行数，理论上应该和我们序列中碱基数一致

> na=0

> nt=0

> ng=0

> nc=0

将各个碱基的计数器清零

z=data[,1]

为了增加我们基础代码的通用性（较少改动），我们首先将序列数据格式又矩阵转为向量格式

for (i in 1:length(z)) {

if(z[i]=="A")

na=na+1

else if(z[i]=="T")

nt=nt+1

else if(z[i]=="G")

ng=ng+1

else if(z[i]=="C")

nc=nc+1

}

这是Rjuan 1.0核心代码

> na

[1] 139

> nc

[1] 172

> nt

[1] 182

> ng

[1] 109

查看结果

> 

sum(c(na,nt,ng,nc))

[1] 602

加总之后等于碱基总数，间接检验运算的可靠性

> 

show=function(x){

namen=c("na","nt","ng","nc")

numbern=c(na,nt,ng,nc)

shown=data.frame(namen,numbern)

print(shown)

}

加入打印函数，将数据保存为矩阵数据shown并在屏幕上输出

> 

show(1)

调用输出函数，括号中参数为任意值

  namen numbern

1    na     139

2    nt     182

3    ng     109

4    nc     172

输出结果