Debug log - R - RunHarmony()

起因

  • 運行以下RunHarmony()時報錯睡蟋,可是代碼以前能跑的梅屉,唯一的不同就是換了服務(wù)器
seuratObj <- RunHarmony(sce, "orig.ident")
  • 報錯內(nèi)容
Error in harmonyObj$cluster_cpp() : 
  element-wise pow(): incompatible matrix dimensions: 100x6 and 6x1

處理過程

方案一 ?:重裝

  • 該問題比較新舞竿,根據(jù)Github解決方案 [1]搀捷,首先運行以下代碼重裝harmony:
devtools::install_github("eddelbuettel/harmony",force = TRUE)
devtools::install_github("immunogenomics/harmony",force = TRUE)
  • 不行

方案二 ?:Live Support

  • 官方在上述issue中提出與官方聯(lián)系 [2]
  • 我聯(lián)系了锅很,但是問題提交的時候一直處于Submitting…

方案三 ?:Rcpp版本問題

  • 從一開始的解決方案里其馏,有人提出是Rcpp包的版本問題 [3]
  • 于是查看Rcpp包的版本 [4]
packageVersion('Rcpp')
[1] ‘1.0.9’
  • 嗯,有可能是版本問題……
  • 安裝舊版本的Rcpp
    • 因為我查看自己電腦上的R版本是1.0.7爆安,所以安裝了一樣的以防萬一
devtools::install_version("Rcpp", version = "1.0.7", repos = "[http://cran.us.r-project.org](http://cran.us.r-project.org/)")
  • 再一次查看版本還是1.0.9
  • 推測可能是已經(jīng)加載過devtools包的緣故
    • 因為在把所有R包全部卸載的情況下叛复,首先安裝devtools,然后查看Rcpp包版本顯示1.0.9
  • 此時無法加載剛安裝的Rcpp
Package ‘Rcpp’ version 1.0.9 cannot be unloaded ……
  • 清除環(huán)境加載的R包重啟Rstudio或上邊欄Session-Restart R,此時所有R包都會unload
    • 加載并查看版本
library("Rcpp")
packageVersion('Rcpp')
[1] ‘1.0.7’
  • 皆大歡喜:职隆?г拧!
  • 然后繼續(xù)安裝其他R包撬码!
  • 此方法適用于沒安裝多少包的
  • 如果安裝了儿倒,應(yīng)該可以通過unload所有包,然后重新加載實現(xiàn)呜笑。有些依賴高版本的Rcpp的R包夫否,可能需要重新下載
  • 繼續(xù)運行RunHarmony()就沒問題了

補充

查看R包版本方法

packageVersion('Seurat')
sessionInfo('Seurat')

  1. https://github.com/immunogenomics/harmony/issues/171 ?

  2. https://devharmony.support/I/-/dev.harmony.one/index.html ?

  3. https://github.com/immunogenomics/harmony/issues/171#issuecomment-1287192826 ?

  4. http://www.reibang.com/p/4b040b29e2ac ?

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市叫胁,隨后出現(xiàn)的幾起案子凰慈,更是在濱河造成了極大的恐慌,老刑警劉巖曹抬,帶你破解...
    沈念sama閱讀 217,826評論 6 506
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件溉瓶,死亡現(xiàn)場離奇詭異,居然都是意外死亡谤民,警方通過查閱死者的電腦和手機堰酿,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,968評論 3 395
  • 文/潘曉璐 我一進店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來张足,“玉大人触创,你說我怎么就攤上這事∥梗” “怎么了哼绑?”我有些...
    開封第一講書人閱讀 164,234評論 0 354
  • 文/不壞的土叔 我叫張陵,是天一觀的道長碉咆。 經(jīng)常有香客問我抖韩,道長,這世上最難降的妖魔是什么疫铜? 我笑而不...
    開封第一講書人閱讀 58,562評論 1 293
  • 正文 為了忘掉前任茂浮,我火速辦了婚禮,結(jié)果婚禮上壳咕,老公的妹妹穿的比我還像新娘席揽。我一直安慰自己,他們只是感情好谓厘,可當(dāng)我...
    茶點故事閱讀 67,611評論 6 392
  • 文/花漫 我一把揭開白布幌羞。 她就那樣靜靜地躺著,像睡著了一般竟稳。 火紅的嫁衣襯著肌膚如雪属桦。 梳的紋絲不亂的頭發(fā)上熊痴,一...
    開封第一講書人閱讀 51,482評論 1 302
  • 那天,我揣著相機與錄音聂宾,去河邊找鬼愁拭。 笑死,一個胖子當(dāng)著我的面吹牛亏吝,可吹牛的內(nèi)容都是我干的岭埠。 我是一名探鬼主播,決...
    沈念sama閱讀 40,271評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼蔚鸥,長吁一口氣:“原來是場噩夢啊……” “哼惜论!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起止喷,我...
    開封第一講書人閱讀 39,166評論 0 276
  • 序言:老撾萬榮一對情侶失蹤馆类,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后弹谁,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體乾巧,經(jīng)...
    沈念sama閱讀 45,608評論 1 314
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 37,814評論 3 336
  • 正文 我和宋清朗相戀三年预愤,在試婚紗的時候發(fā)現(xiàn)自己被綠了沟于。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 39,926評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡植康,死狀恐怖旷太,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情销睁,我是刑警寧澤供璧,帶...
    沈念sama閱讀 35,644評論 5 346
  • 正文 年R本政府宣布,位于F島的核電站冻记,受9級特大地震影響睡毒,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜冗栗,卻給世界環(huán)境...
    茶點故事閱讀 41,249評論 3 329
  • 文/蒙蒙 一演顾、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧贞瞒,春花似錦偶房、人聲如沸趁曼。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,866評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽挡闰。三九已至乒融,卻和暖如春掰盘,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背赞季。 一陣腳步聲響...
    開封第一講書人閱讀 32,991評論 1 269
  • 我被黑心中介騙來泰國打工愧捕, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人申钩。 一個月前我還...
    沈念sama閱讀 48,063評論 3 370
  • 正文 我出身青樓次绘,卻偏偏與公主長得像,于是被迫代替她去往敵國和親撒遣。 傳聞我的和親對象是個殘疾皇子邮偎,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 44,871評論 2 354

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