R2-30第二阶段第四次作业

土貉 2018-04-17 22:04:09 阅读: 1467
library("httr")
yeardata=c()
countdata=c()
url = "https://eutils.ncbi.nlm.nih.gov/entrez/eutils/esearch.fcgi"
for (i in 1:10) {
  year=2007 +i
  term = paste('miRNA[TI] AND ',year,'[DP]',sep = "")
  r<- POST(
    url,
    body = list(
      db='pubmed',
      term=term,
      retmode='json',
      retstart=0,
      retmax=10,
      rettype='uilist'
      )
   )
  stop_for_status(r)
  data=content(r,"parsed","application/json")
  esearchresult=data$esearchresult
  count=esearchresult$count
  yeardata=c(yeardata,year)
  countdata=c(countdata,count)
} 
data=data.frame(year=yeardata,count=countdata)
print(data)

library(ggplot2)
ggplot(data,aes(x=factor(year),y=count))+
  geom_bar(aes(fill=count),stat = "identity")+
  geom_text(label=data$count,colur='red',vjust=-0.5)+
  theme(axis.ticks = element_blank(),
      panel.grid.minor.y=element_blank(),
      panel.grid.major.y = element_blank(),
      axis.text.y =element_blank())+
  labs(x='year',y='count')

001.png

 
邀请讨论

附件

{{f.title}} 大小 {{f.file_size}} 下载 {{f.count_download}} 金币 {{f.count_gold}}
{{item.nick_name}} 受邀请回答 {{item.create_time}}
{{item.refer_comment.nick_name}} {{item.refer_comment.create_time}}

附件

{{f.title}} 大小 {{f.file_size}} 下载 {{f.count_download}} 金币 {{f.count_gold}}
切换到完整回复 发送回复
赞({{item.count_zan}}) 踩({{item.count_cai}}) 删除 回复 关闭
科研狗©2015-2024 科研好助手,京ICP备20005780号-1 建议意见

服务热线

178 0020 3020

微信服务号