Skip to content

Commit

Permalink
20230602_readme
Browse files Browse the repository at this point in the history
  • Loading branch information
LawrenceLiu023 committed Jun 2, 2023
1 parent 9912a41 commit e7bd2ec
Showing 1 changed file with 0 additions and 40 deletions.
40 changes: 0 additions & 40 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -308,43 +308,3 @@ if api_key:
if s:
s.close()
```

R API代码示例:

```r
#安装disgenet2r包
library(devtools)
install_bitbucket("ibi_group/disgenet2r")
library(disgenet2r)

#获取API密钥
library(disgenet2r)
disgenet_api_key <- get_disgenet_api_key(
email = "user@gmail.com",
password = "myspwd" )
Sys.setenv(DISGENET_API_KEY= disgenet_api_key)
#检索疾病
library(disgenet2r)
dq <- disease2gene(disease = "C0028754",
database = "ALL",
score = c(0.3,1))

#检索与基因列表相关的疾病
results <- gene2disease( gene = c( "KCNE1", "KCNE2", "KCNH1", "KCNH2", "KCNG1"), verbose = TRUE)

#检索与疾病列表相关的基因
results2 <- disease2gene( disease = c("C0036341", "C0002395", "C0030567","C0005586"), database = "CURATED", verbose = TRUE )

#检索与变异列表相关的疾病
results3 <- variant2disease( variant= "rs121913279", database = "CURATED")

#检索与疾病列表相关的变体
results4 <- disease2variant(disease = c("C3150943", "C1859062", "C2678485", "C4015695"), database = "CURATED", score = c(0.75, 1))

#疾病富集功能接收一系列基因,并在DisGeNET中对疾病进行富集分析。输入的基因列表应由HGNC符号或Entrez基因标识符识别。词汇表应该使用参数词汇表来指定。默认情况下,vocabulary = "HGNC"。该函数还有另一个可选参数:源数据库(默认情况下,database = CURATED)
#由多次Fisher检验得到的p值,使用benjamin - hochberg方法对错误发现率进行了校正

res_enrich <-disease_enrichment( entities ="KCNE1,KCNE2", vocabulary = "HGNC", database = "CURATED")

#更多用法参见相关链接
```

0 comments on commit e7bd2ec

Please sign in to comment.