请教各位大佬：
1）
length(gene$PHE)  可以
length(gene$AX.108768500) 可以
length(gene$marker[1]) 不可以
我将AX.108768500，输到一个名为marker的字符串中，然后用marker[1]的方式调用，不行。

2）
for(i in c(1:5)){
       a = ddply(.data = gene,
              .variables = marker,
              .fun = function(x){
                mean(x$PHE)
        })
       write.csv(a,"a.csv")
}

只能输出最后一次的结果。
3）
其实和第一个很像，因为我准备做一个for循环：
for(i in c(1:23)){
  a = ddply(.data = gene,
            .variables = marker,
            .fun = function(x){
              kruskal.test(x$PHE~marker,data=gene,method = "holm")
            })
  print(a)
}

请教大佬有什么好的方法吗？

数据：

rs#	AX.108768500	AX.110999783	AX.108750670	AX.111672328	PHE
alleles	T/C	C/G	A/C	C/G	NA
S43	TT	GG	AA	GG	29.3417
S151	TC	CG	AC	GG	32.0381
S189	CC	CC	AC	GG	37.3938
S172	CC	CC	CC	GG	35.4436
S292	CC	CC	AC	GG	38.068
S246	CC	CC	CC	GG	31.0997