にrssnp変換する私は、これは明らかに問題がある場合ので、私を許してください、私は生物学で何のバックグラウンドを持っていませんR.Biomartは、遺伝子名
library(biomaRt)
snp_mart = useMart("ENSEMBL_MART_SNP", dataset="hsapiens_snp")
snp_attributes = c("refsnp_id", "chr_name", "chrom_start",
"associated_gene", "ensembl_gene_stable_id", "minor_allele_freq")
getENSG <- function(rs, mart = snp_mart) {
results <- getBM(attributes = snp_attributes,
filters = "snp_filter", values = rs, mart = mart)
return(results)
}
getENSG("rs144864312")
refsnp_id chr_name chrom_start associated_gene ensembl_gene_stable_id
1 rs144864312 8 20254959 NA ENSG00000061337
minor_allele_freq
1 0.000399361
に次のコードを持っています。私はrs144864312が遺伝子名 "LZTS1"と一致するはずだと言われました。 私が上に書いたコードは主にインターネットからのものです。私の質問はどこから遺伝子名を抽出するのですか?私はlistAttributes(snp_mart)がすべての可能な出力のリストを提供するが、上記の "遺伝子名"を私に与えるものは見当たらない。バイオマーカーを使用することからこの遺伝子名を抽出することができます(そしてrs番号が与えられています)?前もって感謝します。
PS:これは、500件のエントリ(1件だけではない)のために行う必要があります。それでなぜ遺伝子名を抽出するために上記のような単純な関数を作りましたか?言っ