r<-read.csv('drug_info.csv', header=T)
head(r[ ,c('MFDS_CODE', 'ETC_OTC')])
drug_cnt<-table(r$ETC_OTC)
결측치 '기타'로 대체
r$ETC_OTC<- ifelse(r$ETC_OTC=='','기타',r$ETC_OTC)
table(r$ETC_OTC)
기타 원료의약품 의약외품 일반의약품
9489 469 24 13276
전문의약품 전문의약품(희귀) 한약재
33933 462 22306
fig4<-plot_ly(x=names(drug_cnt),
y=as.numeric(drug_cnt),
type='bar',
marker=list(color=hcl.colors((length(drug_cnt)),'TealGrn')))
fig4
pie(drug_cnt, main='허가의약품 분류별 원형 그래프',
col=hcl.colors((length(drug_cnt)),'TealGrn'),
labels=paste(names(drug_cnt),
'\n',drug_cnt))
r2<-aggregate(MFDS_CODE ~ ETC_OTC, r, length)
names(r2) <- c('전문일반', '의약품수')
pct<-round(r2$의약품수/length(r$MFDS_CODE)*100,1)
drug_label<-paste(r2$전문일반, ':',pct,'%' )
drug_label
[1] "기타 : 11.9 %" "원료의약품 : 0.6 %"
[3] "의약외품 : 0 %" "일반의약품 : 16.6 %"
[5] "전문의약품 : 42.4 %" "전문의약품(희귀) : 0.6 %"
[7] "한약재 : 27.9 %"
pie(drug_cnt,main='허가의약품 분류별 원형 그래프',
labels=drug_label,
col=hcl.colors((length(drug_cnt)),'TealGrn'))