R語言小作業(yè)-中級

首先需要完成R語言練習題-初級,在[link]http://www.bio-info-trainee.com/3793.html

作業(yè) 1
請根據R包org.Hs.eg.db找到下面ensembl 基因ID 對應的基因名(symbol)

ENSG00000000003.13
ENSG00000000005.5
ENSG00000000419.11
ENSG00000000457.12
ENSG00000000460.15
ENSG00000000938.11

提示:

library(org.Hs.eg.db)
g2s=toTable(org.Hs.egSYMBOL)
g2e=toTable(org.Hs.egENSEMBL)

作業(yè) 2

根據R包hgu133a.db找到下面探針對應的基因名(symbol)

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

提示:

library(hgu133a.db)
ids=toTable(hgu133aSYMBOL)
head(ids)

作業(yè) 3

找到R包CLL內置的數據集的表達矩陣里面的TP53基因的表達量矗愧,并且繪制在progres.-stable分組的boxplot圖
提示:

suppressPackageStartupMessages(library(CLL))
data(sCLLex)
sCLLex
exprSet=exprs(sCLLex) 
library(hgu95av2.db)

想想如何通過 ggpubr進行美化御吞。

作業(yè) 4

找到BRCA1基因在TCGA數據庫的乳腺癌數據集(Breast Invasive Carcinoma (TCGA, PanCancer Atlas))的表達情況

提示:使用[link]http://www.cbioportal.org/index.do 定位數據集:[link]http://www.cbioportal.org/datasets

作業(yè) 5

找到TP53基因在TCGA數據庫的乳腺癌數據集的表達量分組看其是否影響生存

提示使用:[link]http://www.oncolnc.org/

作業(yè)6

下載數據集GSE17215的表達矩陣并且提取下面的基因畫熱圖
ACTR3B ANLN BAG1 BCL2 BIRC5 BLVRA CCNB1 CCNE1 CDC20 CDC6 CDCA1 CDH3 CENPF CEP55 CXXC5 EGFR ERBB2 ESR1 EXO1 FGFR4 FOXA1 FOXC1 GPR160 GRB7 KIF2C KNTC2 KRT14 KRT17 KRT5 MAPT MDM2 MELK MIA MKI67 MLPH MMP11 MYBL2 MYC NAT1 ORC6L PGR PHGDH PTTG1 RRM2 SFRP1 SLC39A6 TMEM45B TYMS UBE2C UBE2T
提示:根據基因名拿到探針I(yè)D,縮小表達矩陣繪制熱圖瓦哎,沒有檢查到的基因直接忽略即可砸喻。

作業(yè)7

下載數據集GSE24673的表達矩陣計算樣本的相關性并且繪制熱圖,需要標記上樣本分組信息

作業(yè)8

找到 GPL6244 platform of Affymetrix Human Gene 1.0 ST Array 對應的R的bioconductor注釋包蒋譬,并且安裝它割岛!

options()$repos
options()$BioC_mirror 
options(BioC_mirror="https://mirrors.ustc.edu.cn/bioc/")
options("repos" = c(CRAN="https://mirrors.tuna.tsinghua.edu.cn/CRAN/"))
BiocManager::install("請輸入自己找到的R包",ask = F,update = F)
options()$repos
options()$BioC_mirror

作業(yè)9

下載數據集GSE42872的表達矩陣,并且分別挑選出 所有樣本的(平均表達量/sd/mad/)最大的探針犯助,并且找到它們對應的基因癣漆。

作業(yè)10

下載數據集GSE42872的表達矩陣,并且根據分組使用limma做差異分析剂买,得到差異結果矩陣

本文作者:生信技能樹團隊

最后編輯于
?著作權歸作者所有,轉載或內容合作請聯系作者
  • 序言:七十年代末惠爽,一起剝皮案震驚了整個濱河市,隨后出現的幾起案子瞬哼,更是在濱河造成了極大的恐慌婚肆,老刑警劉巖,帶你破解...
    沈念sama閱讀 219,539評論 6 508
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件坐慰,死亡現場離奇詭異较性,居然都是意外死亡,警方通過查閱死者的電腦和手機结胀,發(fā)現死者居然都...
    沈念sama閱讀 93,594評論 3 396
  • 文/潘曉璐 我一進店門赞咙,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人糟港,你說我怎么就攤上這事攀操。” “怎么了着逐?”我有些...
    開封第一講書人閱讀 165,871評論 0 356
  • 文/不壞的土叔 我叫張陵崔赌,是天一觀的道長意蛀。 經常有香客問我,道長健芭,這世上最難降的妖魔是什么县钥? 我笑而不...
    開封第一講書人閱讀 58,963評論 1 295
  • 正文 為了忘掉前任,我火速辦了婚禮慈迈,結果婚禮上若贮,老公的妹妹穿的比我還像新娘。我一直安慰自己痒留,他們只是感情好谴麦,可當我...
    茶點故事閱讀 67,984評論 6 393
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著伸头,像睡著了一般匾效。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上恤磷,一...
    開封第一講書人閱讀 51,763評論 1 307
  • 那天面哼,我揣著相機與錄音,去河邊找鬼扫步。 笑死魔策,一個胖子當著我的面吹牛,可吹牛的內容都是我干的河胎。 我是一名探鬼主播闯袒,決...
    沈念sama閱讀 40,468評論 3 420
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼游岳!你這毒婦竟也來了政敢?” 一聲冷哼從身側響起,我...
    開封第一講書人閱讀 39,357評論 0 276
  • 序言:老撾萬榮一對情侶失蹤胚迫,失蹤者是張志新(化名)和其女友劉穎堕仔,沒想到半個月后,有當地人在樹林里發(fā)現了一具尸體晌区,經...
    沈念sama閱讀 45,850評論 1 317
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...
    茶點故事閱讀 38,002評論 3 338
  • 正文 我和宋清朗相戀三年通贞,在試婚紗的時候發(fā)現自己被綠了朗若。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 40,144評論 1 351
  • 序言:一個原本活蹦亂跳的男人離奇死亡昌罩,死狀恐怖哭懈,靈堂內的尸體忽然破棺而出,到底是詐尸還是另有隱情茎用,我是刑警寧澤遣总,帶...
    沈念sama閱讀 35,823評論 5 346
  • 正文 年R本政府宣布睬罗,位于F島的核電站,受9級特大地震影響旭斥,放射性物質發(fā)生泄漏容达。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 41,483評論 3 331
  • 文/蒙蒙 一垂券、第九天 我趴在偏房一處隱蔽的房頂上張望花盐。 院中可真熱鬧,春花似錦菇爪、人聲如沸算芯。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,026評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽熙揍。三九已至,卻和暖如春氏涩,著一層夾襖步出監(jiān)牢的瞬間届囚,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 33,150評論 1 272
  • 我被黑心中介騙來泰國打工削葱, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留奖亚,地道東北人。 一個月前我還...
    沈念sama閱讀 48,415評論 3 373
  • 正文 我出身青樓析砸,卻偏偏與公主長得像昔字,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子首繁,可洞房花燭夜當晚...
    茶點故事閱讀 45,092評論 2 355