IMvigor210CoreBiologies中DESeq包安裝體會(huì)

最近想研究一下IMvigor210 cohort馋袜,生信技能樹(shù)已經(jīng)把指南寫(xiě)得很清楚了IMvigor210CoreBiologies包安裝指北 (qq.com)脱衙,但是今天發(fā)現(xiàn)了一個(gè)bug开缎,花了我?guī)讉€(gè)小時(shí)時(shí)間沈贝,總算找出答案~~~~我把整個(gè)過(guò)程寫(xiě)下來(lái),研究之后發(fā)現(xiàn)珊燎,R包的安裝有時(shí)候真是一把辛酸淚庐扫。
首先我們按照流程來(lái),我的R版本是4.1

version                      
#platform       x86_64-w64-mingw32          
#arch           x86_64                      
#os             mingw32                     
#system         x86_64, mingw32             
#status                                     
#major          4                           
#minor          1.0                         
#year           2021                        
#month          05                          
#day            18                          
#svn rev        80317                       
#language       R                           
#version.string R version 4.1.0 (2021-05-18)
#nickname       Camp Pontanezen      

我們直接先本地安裝塑娇,先按照生信技能樹(shù)教程的指引澈侠,去網(wǎng)站下載IMvigor210CoreBiologies_1.0.0.tar.gz和DESeq_1.38.0.tar.gz兩個(gè)文件,然后本地安裝埋酬,接著就報(bào)錯(cuò)了哨啃。

install.packages("IMvigor210CoreBiologies_1.0.0.tar.gz",repos=NULL)
#這個(gè)是可以安裝成功的
load("cds.RData")
#Warning: namespace ‘DESeq’ is not available and has been replaced
#by .GlobalEnv when processing object ‘cds’. 載入需要的程輯包:DESeq
#Failed with error:  ‘程輯包‘DESeq’是在R版本3.0.0之前建的:你得重新安裝’
#Error in .requirePackage(package) : unable to find required package ‘DESeq’.
install.packages('DESeq')
#Warning in install.packages :package ‘DESeq’ is not available for this version of R
#A version of this package for your version of R might be available elsewhere,
#see the ideas at https://cran.r-project.org/doc/manuals/r-patched/R-admin.html#Installing-packages

可以看到有2個(gè)報(bào)錯(cuò),一個(gè)是當(dāng)我load數(shù)據(jù)集cds.RData的時(shí)候写妥,提示了DESeq包得重新安裝拳球,第二個(gè)是無(wú)法安裝DESeq包,也就是說(shuō)關(guān)鍵在于DESeq包是否成功安裝珍特。
那么有沒(méi)有可能是devtool包的問(wèn)題祝峻,我接著操作了一下,發(fā)現(xiàn)devtools是沒(méi)問(wèn)題的,我之前有安裝過(guò)了扎筒。

library(devtools)
#載入需要的程輯包:usethis
#Warning messages:
#1: 程輯包‘devtools’是用R版本4.1.2 來(lái)建造的 
#2: 程輯包‘usethis’是用R版本4.1.2 來(lái)建造的 
# devtools::find_rtools()
[1] TRUE

那么我們重新找安裝DESep包的辦法莱找,采用最官方的語(yǔ)句安裝

if (!requireNamespace("BiocManager", quietly = TRUE))
    install.packages("BiocManager")
   BiocManager::install("DESeq")
#Warning message:
#package(s) not installed when version(s) same as current; use `force = TRUE` to
#re-install: 'DESeq' 
load("cds.RData")
#Warning: namespace ‘DESeq’ is not available and has been replaced
#by .GlobalEnv when processing object ‘cds’
#載入需要的程輯包:DESeq
#Failed with error:  ‘程輯包‘DESeq’是在R版本4.0.0之前建的:你得重新安裝’
#Error in .requirePackage(package) : 
#unable to find required package ‘DESeq’

#按照它的要求,加上force=TRUE
if (!requireNamespace("BiocManager", force=TRUE)
    install.packages("BiocManager")
   BiocManager::install("DESeq")
#Warning message:
#package(s) not installed when version(s) same as current; use `force = TRUE` to
#re-install: 'DESeq' 
load("cds.RData")
#Warning: namespace ‘DESeq’ is not available and has been replaced
#by .GlobalEnv when processing object ‘cds’
#載入需要的程輯包:DESeq
#Failed with error:  ‘程輯包‘DESeq’是在R版本4.0.0之前建的:你得重新安裝’
#Error in .requirePackage(package) : 
#unable to find required package ‘DESeq’

可以看到嗜桌,這兩次安裝都有warning,但是沒(méi)報(bào)錯(cuò)奥溺,不過(guò)仍然無(wú)法順利的load("cds.RData"),所以DESeq仍沒(méi)有被正確安裝。我去看了一下bioconductor對(duì)這個(gè)包的注釋?zhuān)l(fā)現(xiàn)這個(gè)包還有一個(gè)本地的安裝包可以用骨宠,就是紅色線(xiàn)上的“DESeq_1.38.0.zip”浮定,我下載下來(lái),本地安裝一下层亿,發(fā)現(xiàn)成功安裝了壶唤,但是仍然沒(méi)作用~~


image.png
install.packages("DESeq_1.38.0.zip", repos = NULL, type = "win.binary")
#package ‘DESeq’ successfully unpacked and MD5 sums checked
load("cds.RData")
#Warning: namespace ‘DESeq’ is not available and has been replaced
#by .GlobalEnv when processing object ‘cds’
#載入需要的程輯包:DESeq
#Failed with error:  ‘程輯包‘DESeq’是在R版本4.0.0之前建的:你得重新安裝’
#Error in .requirePackage(package) : 
 #unable to find required package ‘DESeq’

我回頭去看了一下原文的代碼,發(fā)現(xiàn)作者竟然是沒(méi)有安裝DESeq這個(gè)包的(TGF-b attenuates tumor response to PD-L1 blockade by contributing to exclusion of T cells (gene.com)
)棕所,跟他很像的是DESeq2,但是不是這個(gè)啊,這個(gè)我有l(wèi)ibrary之后悯辙,仍然無(wú)效琳省。


我們不要灰心迎吵,把各種報(bào)錯(cuò)信息放到網(wǎng)上找一下答案。試一下這個(gè)针贬,不行~~
也有人說(shuō)現(xiàn)在建議使用DESeq2代替DESeq,但是我已經(jīng)試過(guò)了击费,沒(méi)效果

太難了~,該怎么辦桦他?我回頭再去細(xì)看一下waring,我總感覺(jué)肯定忽略了什么東西蔫巩,可以看到warning其中比較重要的1句話(huà)是

#Failed with error:  ‘程輯包‘DESeq’是在R版本4.0.0之前建的:你得重新安裝

這提示著這個(gè)包或許還沒(méi)有及時(shí)更新,因此用R version 4.1是沒(méi)法運(yùn)作的快压,這時(shí)候我們就得回去bioconductor看一下這個(gè)包的歷史了圆仔,我們可以看到藍(lán)色字體部分,這個(gè)包對(duì)應(yīng)的R版本是3.6蔫劣,那么是不是換成3.6版本就行了嗎坪郭?還有其他原因嗎?再仔細(xì)看一下脉幢,上面說(shuō)的這句“This package is for version 3.10 of Bioconductor; for the stable, up-to-date release version, see DESeq"歪沃,讓我恍然大悟!

我們點(diǎn)開(kāi)黃色線(xiàn)上DESeq進(jìn)去看看嫌松,總算發(fā)現(xiàn)原因了沪曙,可以看到,在Bioconductor version 3.13 萎羔,DESeq就已經(jīng)被removed液走,取而代之的是DESeq2,但是IMvigor210CoreBiologies仍然需要的依賴(lài)包是DESeq!真是淚目......

那么接下來(lái)就簡(jiǎn)單了外驱,我把R version改成3.6育灸,然后將Bioconductor version改為3.10
if (!require("BiocManager", quietly = TRUE))
    install.packages("BiocManager")
BiocManager::install(version = "3.10")
install.packages("DESeq_1.38.0.tar.gz",repos=NULL)
install.packages("IMvigor210CoreBiologies_1.0.0.tar.gz",repos=NULL)
#因?yàn)閾Q版本了,所以依賴(lài)包都要重新安裝一下昵宇,這里我不展示了
library(IMvigor210CoreBiologies)
load("cds.RData")
cds
#CountDataSet (storageMode: environment)
#assayData: 31286 features, 348 samples 
 #element names: counts 
#protocolData: none
#phenoData
#sampleNames: SAMf2ce197162ce SAM698d8d76b934 ... SAM2de7cffb5f72 (348
# total)
#varLabels: Best Confirmed Overall Response binaryResponse ... TCGA
#Subtype (25 total) 
#varMetadata: labelDescription
#featureData
# featureNames: 1 10 ... 9997 (31286 total)
#fvarLabels: entrez_id symbol ... Symbol (6 total)
# fvarMetadata: labelDescription
#experimentData: use 'experimentData(object)'
#Annotation:  

漂亮磅崭!沒(méi)報(bào)錯(cuò)了,有了這個(gè)文件瓦哎,下面的表達(dá)矩陣砸喻,臨床信息就都可以提取啦〗總之包的安裝很費(fèi)勁割岛,死活安裝不成功的時(shí)候,不妨仔細(xì)看一看waring,然后到bioconductor仔細(xì)看一下犯助,有時(shí)候或許是Rversion癣漆,bioconductor version在搞事情.....這些都是實(shí)戰(zhàn)經(jīng)驗(yàn)啊~

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個(gè)濱河市剂买,隨后出現(xiàn)的幾起案子惠爽,更是在濱河造成了極大的恐慌癌蓖,老刑警劉巖,帶你破解...
    沈念sama閱讀 206,013評(píng)論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件婚肆,死亡現(xiàn)場(chǎng)離奇詭異租副,居然都是意外死亡,警方通過(guò)查閱死者的電腦和手機(jī)较性,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,205評(píng)論 2 382
  • 文/潘曉璐 我一進(jìn)店門(mén)用僧,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái),“玉大人赞咙,你說(shuō)我怎么就攤上這事责循。” “怎么了人弓?”我有些...
    開(kāi)封第一講書(shū)人閱讀 152,370評(píng)論 0 342
  • 文/不壞的土叔 我叫張陵沼死,是天一觀的道長(zhǎng)。 經(jīng)常有香客問(wèn)我崔赌,道長(zhǎng)意蛀,這世上最難降的妖魔是什么? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 55,168評(píng)論 1 278
  • 正文 為了忘掉前任健芭,我火速辦了婚禮县钥,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘慈迈。我一直安慰自己若贮,他們只是感情好,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,153評(píng)論 5 371
  • 文/花漫 我一把揭開(kāi)白布痒留。 她就那樣靜靜地躺著谴麦,像睡著了一般。 火紅的嫁衣襯著肌膚如雪伸头。 梳的紋絲不亂的頭發(fā)上匾效,一...
    開(kāi)封第一講書(shū)人閱讀 48,954評(píng)論 1 283
  • 那天,我揣著相機(jī)與錄音恤磷,去河邊找鬼面哼。 笑死,一個(gè)胖子當(dāng)著我的面吹牛扫步,可吹牛的內(nèi)容都是我干的魔策。 我是一名探鬼主播,決...
    沈念sama閱讀 38,271評(píng)論 3 399
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼河胎,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼闯袒!你這毒婦竟也來(lái)了?” 一聲冷哼從身側(cè)響起,我...
    開(kāi)封第一講書(shū)人閱讀 36,916評(píng)論 0 259
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤搁吓,失蹤者是張志新(化名)和其女友劉穎原茅,沒(méi)想到半個(gè)月后,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體堕仔,經(jīng)...
    沈念sama閱讀 43,382評(píng)論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 35,877評(píng)論 2 323
  • 正文 我和宋清朗相戀三年晌区,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了摩骨。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 37,989評(píng)論 1 333
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡朗若,死狀恐怖恼五,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情哭懈,我是刑警寧澤灾馒,帶...
    沈念sama閱讀 33,624評(píng)論 4 322
  • 正文 年R本政府宣布,位于F島的核電站遣总,受9級(jí)特大地震影響睬罗,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜旭斥,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,209評(píng)論 3 307
  • 文/蒙蒙 一容达、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧垂券,春花似錦花盐、人聲如沸。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 30,199評(píng)論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)乙嘀。三九已至,卻和暖如春熙揍,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背近速。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 31,418評(píng)論 1 260
  • 我被黑心中介騙來(lái)泰國(guó)打工诈嘿, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留,地道東北人削葱。 一個(gè)月前我還...
    沈念sama閱讀 45,401評(píng)論 2 352
  • 正文 我出身青樓奖亚,卻偏偏與公主長(zhǎng)得像,于是被迫代替她去往敵國(guó)和親析砸。 傳聞我的和親對(duì)象是個(gè)殘疾皇子昔字,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,700評(píng)論 2 345

推薦閱讀更多精彩內(nèi)容