ReactomePAのエラー

factorがduplicatedしてるよと以下のように怒られるエラー。

(Error in levels <- ... factor level is duplicated Calls: ... fortify.enrichResult -> fortify.internal -> factor)


遺伝子セットは完璧にuniqueになっているのに、どうしても消えない。
ソースコードをみると、fortify-internalでDescriptionをfactor化しているときにエラーが入っている模様。

issueにあったけど解決してないようなので(すごい対応だなと笑った)、
https://github.com/GuangchuangYu/clusterProfiler/issues/116

無理やりEnrichの結果からDescrpitionが重複しているものをresultから除く処理を入れたよ。
https://gist.github.com/carushi/cc82d623af80d633539022aabbc32d2a

コメント

このブログの人気の投稿

再現性なんてないさ(?)

特に収穫がない日

qsubに引数