> 

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

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