有一个软件asreml的R包：
asreml(protein~snp10001, random=~ped(ID),data=SNPs,ginverse=list(ID=ainv))
这样可以出结果
ASReml: Thu Oct 22 17:52:33 2015

     LogLik         S2      DF      wall     cpu
  -1648.7841500142039.1348   156  17:52:33     0.0
  -1648.5315487038930.9418   156  17:52:33     0.0
因为我有好几百个snp，从snp00001-snp10001，每一个都要跑一边，我想写个循环，却发现无法进行，比如：

我把a <- snp10001 #把a赋值
然后运行程序：
asreml(protein~a, random=~ped(ID),data=SNPs,ginverse=list(ID=ainv))
报错：
Error in asreml.mf.default(formula = protein ~ a + ID, data = list(ID = 1:157,  :
  变数的长度不一样('a')


在R中怎么调用，这个功能在perl中可以实现，在R中怎么调用a呢？
ps，我把文件里面的snp名称换成1234，也不能调用，这不是引号的问题。
请大神指导！