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
如何使用treeWAS
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
