前幾天,有一個(gè)合作課題的朋友臨時(shí)提到絮记,現(xiàn)在咱們幾百個(gè)樣品有一個(gè)樹摔踱,但是我只想看其中30個(gè)樣品的樹,怎么辦到千?尤其是昌渤,他正好有只有這 30 個(gè)樣品的一些數(shù)值矩陣云云。
當(dāng)然憔四,最簡(jiǎn)單的做法就是直接從一棵幾百個(gè)Taxon的進(jìn)化樹里面膀息,抽取出這30個(gè)樣品的演化關(guān)系。說(shuō)起來(lái)簡(jiǎn)單了赵,做起來(lái)就比較復(fù)雜潜支。我簡(jiǎn)單檢索了,似乎并無(wú)合適實(shí)現(xiàn)柿汛。其中基本也全是冗酿,你用python某軟件云云,或者用R的什么包络断,而且也寫得不太清楚裁替,用起來(lái)很不方便。當(dāng)然貌笨,我嘗試了用chatGPT或者openAI弱判,并無(wú)解決問(wèn)題。Emmm....盡管我做過(guò)一些分析锥惋,python R 也不至于完全不會(huì)昌腰,但是這么一個(gè)小事,你就要騙我去學(xué)一個(gè)包膀跌,學(xué)一個(gè)模塊的API遭商?不接受。
不接受捅伤,那就只自己寫劫流。于是 TBtools 終于有了一個(gè)新功能。
使用簡(jiǎn)單丛忆,直接輸入一個(gè)大樹祠汇,給一下你需要的子集ID列表,點(diǎn)擊 Start 就可以看到一個(gè)小樹蘸际。
是不是覺(jué)得很實(shí)用?