学习总结

时间：2019.7.4
内容：初级题与中级题

分组

> plate <- as.data.frame(plate)
> e$plate <- plate###将plate加入e中成为plate列
> e1 <- e[e$plate=='0048']
> class(e1)
[1] "character"
> dim(e1) <- c(384,3)
> class(e1)
[1] "matrix"

> e2 <- e[e$plate=='0049']
> dim(e2) <- c(384,3)

> e1 <- e1[,-2]
> e2 <- e2[,-2]
> e1 <- as.numeric(e1)
> dim(e1 ) <- c(384,2)
> class(e1[,1])
[1] "numeric"
> colnames(e1) <- c('MBases','plate')
> e2 <- as.numeric(e2)
> dim(e2) <- c(384,2)
> colnames(e2) <- c('MBases','plate')

频数图

> hist(e1[,1])

image.png

> hist(e2[,1])

image.png

使用ggplot2画图

箱图

library(ggplot2)
class(e1)
e1 <- as.data.frame(e1)
e2 <- as.data.frame(e2)
> ggplot(e1,aes(x=plate,y=MBases))+geom_boxplot()
> ggplot(e2,aes(x=plate,y=MBases))+geom_boxplot()

image.png image.png

频数图

> ggplot(e1,aes(x=MBases))+geom_histogram(bins = 40,color="blue")
> ggplot(e2,aes(x=MBases))+geom_histogram(bins = 40,color="blue")

image.png

> ggplot(e2,aes(x=MBases))+geom_histogram(bins = 40,color="blue")

image.png

密度图

> ggplot(e1,aes(x=MBases))+geom_density()

image.png

> ggplot(e2,aes(x=MBases))+geom_density()

image.png

使用ggpubr作图

library(ggpubr)

箱图

> ggboxplot(e1,x = 'plate',y = 'MBases')

image.png

> ggboxplot(e2,x = 'plate',y = 'MBases')

image.png

频数图

> gghistogram(e1,x='MBases',bins=30)

image.png

> gghistogram(e2,x='MBases',bins=30)

image.png

密度图

> ggdensity(e1,x='MBases')

image.png

> ggdensity(e2,x='MBases')

image.png

随机取384个MBases信息，跟前面的两个plate的信息组合成新的数据框，第一列是分组，第二列是MBases，总共是384*3行数据。

> a1 <- e$MBases[1:384]
> a2 <- e$Title[1:384]
> a <- data.frame(a1,a2)
> a$plate <- as.data.frame(plate[,1][1:384])
> colnames(a) <- c('MBases','Title','plate')

中级题

作业 1

根据R包org.Hs.eg.db找到下面ensembl 基因ID 对应的基因名(symbol)

> g2s <- toTable(org.Hs.egSYMBOL)
> g2e <- toTable(org.Hs.egENSEMBL)
> ensemble_id <- c('ENSG00000000003.13','ENSG00000000005.5','ENSG00000000419.11','ENSG00000000457.12','ENSG00000000460.15','ENSG00000000938.11')

> #批量取基因名
> library(stringr)

> unlist(str_split(ensemble_id,'[.]'))
 [1] "ENSG00000000003" "13"              "ENSG00000000005" "5"              
 [5] "ENSG00000000419" "11"              "ENSG00000000457" "12"             
 [9] "ENSG00000000460" "15"              "ENSG00000000938" "11"       

> tmp <- unlist(str_split(ensemble_id,'[.]',simplify = T))###simplify = T 此参数生成为矩阵

image.png

> class(unlist(str_split(ensemble_id,'[.]',simplify = T)))
[1] "matrix"
> ensemble_id <- tmp[,1]
> ensembl_id <- as.data.frame(ensemble_id)

image.png

> colnames(ensembl_id) <- 'ensembl_id'
> merge1 <- merge(x=ensembl_id,y=g2e,by='ensembl_id')

image.png

> merge2 <- merge(x=merge1,y=g2s,by='gene_id')

image.png

作业 2

根据R包hgu133a.db找到下面探针对应的基因名(symbol)

> tmp <- c('1053_at','117_at','121_at','1255_g_at','1316_at','1320_at','1405_i_at','1431_at','1438_at','1487_at','1494_f_at','1598_g_at','160020_at','1729_at','177_at')
> probe_id <- as.data.frame(tmp)
> View(probe_id)
> colnames(probe_id) <- 'probe_id'
> View(probe_id)

image.png

> a <- toTable(hgu133aSYMBOL)
> merge <- merge(x=probe_id,y=a,by='probe_id')

image.png