如何使用treeWAS

library(treeWAS)
gene_matrix<-read.table("homolog_matrix.txt")
t_gene_matrix<-t(gene_matrix)
t_gene_matrix[which(t_gene_matrix>0)]<-1

df1<-read.table("pheno.txt")
df<-t(pheno)
phen <- as.vector(unlist(df))
names(phen) <- rownames(df)

out <- treeWAS(snps = t_gene_matrix,
              phen = phen,
              tree = "NJ*",
              plot.tree = TRUE,
              seed = 1)

str(out)
correlation<-out$terminal$corr.dat
test<-data.frame(correlation)
print(out:, sort.by.p=FALSE)


x <- sample(0:4,100,replace=T)
x[which(x==0|x==1|x==2)] <-3
x