最近我在B站的視頻頁面下發(fā)現(xiàn)了這樣的一條評論，面對質(zhì)譜數(shù)據(jù)分析領(lǐng)域內(nèi)的初學(xué)者的求教，其實自己也是非常的誠惶誠恐的。因為在視頻中所使用的腳本語言是自己開發(fā)的一門新語言，所以可能給@徐冬冬阿哈哈哈哈同學(xué)造成了一部分的困擾哈哈。首先先對這位粉絲說一聲抱歉哈。抱歉抱歉抱歉抱歉，，，
針對上述的提問，我的回答大概是有以下的幾點：
(資料圖)
在B站的教程視頻《【BioNovoGene Mzkit教程】代謝組學(xué)原始數(shù)據(jù)處理基礎(chǔ)》中，所提到的腳本語言其實并不是真正的R語言，而是基于R語言的向量化編程特性而仿照著開發(fā)出來的一門國產(chǎn)編程語言R#。所以有一些視頻中所提到的軟件包在真正的R語言環(huán)境之中是無法找到的。
所以netCDF.utils程序包模塊其實是R#編程語言的環(huán)境庫之中自帶的一個用于讀寫CDF文件的模塊，在真正的R語言環(huán)境之中是找不到的哈哈
因為R#語言的開發(fā)結(jié)合了多種語言的特性，所以看起來可能會有一些四不像。但是腳本中所展示的imports之類的語句，其實并不是來自于python，而是來自于javascript語言。代碼中using相關(guān)的代碼語法其實是來自于C#語言或者VB.NET語言之中的using。
那說完了上面的針對評論回復(fù)的幾點，下面我們開始進(jìn)行一些拓展性的知識講解吧：關(guān)于如何使用上R#語言進(jìn)行質(zhì)譜數(shù)據(jù)分析。
首先，先向大伙簡單介紹一下MZKit桌面工作站軟件的一些背景信息。在視頻中所介紹的MZKit軟件，其實是帕諾米克從2017年啟動的一個質(zhì)譜原始數(shù)據(jù)開源項目，相關(guān)的項目源代碼大家可以從Github上下載得到：https://github.com/xieguigang/mzkit
目前呢，經(jīng)過多年的發(fā)展，MZKit開源項目目前在質(zhì)譜原始數(shù)據(jù)處理相關(guān)領(lǐng)域內(nèi)的應(yīng)用其實已經(jīng)比較廣泛了：就帕諾米克而言，MZKit軟件以及R#腳本編程語言已經(jīng)非常廣泛的應(yīng)用于各種商業(yè)項目中的質(zhì)譜原始數(shù)據(jù)處理的自動化流程搭建之中。那現(xiàn)在廢話不多說，回到我們的正題：我們該如何通過R#編程語言腳本調(diào)用MZKit進(jìn)行質(zhì)譜數(shù)據(jù)的一些處理操作呢？假若各位同學(xué)已經(jīng)安裝有MZKit桌面工作站軟件的話，實際上你已經(jīng)獲取得到了一個R#編程語言的運行環(huán)境：因為MZKit的一些數(shù)據(jù)分析功能和數(shù)據(jù)可視化功能是基于R#腳本語言編寫的，R#語言編程環(huán)境已經(jīng)集成在了MZKit桌面工作站軟件之中了。
既然R#編程語言環(huán)境已經(jīng)集成在了視頻中所提到的MZKit桌面工作站軟件之中，那么我們怎么打開它呢？其實這個非常的簡單，我們只需要按照下面的操作步驟就可以進(jìn)入R#編程環(huán)境之中：
首先，通過軟件左上角的【File】菜單入口打開軟件的主菜單；接著可以在主菜單之中看見有【MZKit Data Toolkit】菜單欄；就是在這個菜單欄之中，展開后可以發(fā)現(xiàn)有一個命令行圖標(biāo)的東西【Developer CommandLine Tools】；對，就是它了，點擊這個東西我們就可以進(jìn)入到R#編程語言環(huán)境之中了。假設(shè)點擊這個菜單之后，跳出來了下面所示的這樣的一個命令行界面，那么就說明你已經(jīng)成功打算開了質(zhì)譜數(shù)據(jù)分析的開發(fā)工具環(huán)境：
那現(xiàn)在我們來編寫一個進(jìn)行CDF數(shù)據(jù)可視化的腳本來初步講解MZKit軟件自帶的質(zhì)譜數(shù)據(jù)處理函數(shù)庫的使用吧。
像傳統(tǒng)的R語言編程一樣，一般編程語言環(huán)境之中是不包含有針對具體的實際問題解決工具的。所以在進(jìn)行相關(guān)的實際問題處理之前，我們需要導(dǎo)入面向具體實際問題進(jìn)行數(shù)據(jù)處理的相關(guān)的軟件包。在這里我們先進(jìn)行接下來所需要解決的具體問題的描述：通過讀取CDF文件之中的色譜圖數(shù)據(jù)，然后分別通過ggplot軟件包或者mzkit軟件包繪制出來。
那根據(jù)上面的問題描述，實際上我們就可以分解出三個腳本操作相關(guān)的模塊：CDF文件數(shù)據(jù)讀取模塊，ggplot繪圖模塊或者mzkit質(zhì)譜數(shù)據(jù)可視化模塊。那么首先，我們向腳本執(zhí)行環(huán)境之中導(dǎo)入CDF文件數(shù)據(jù)讀取模塊netCDF.utils進(jìn)行數(shù)據(jù)讀取操作：
CDF文件數(shù)據(jù)讀取
具體的來講，上面的語句其實就是從base.dll這個程序模塊文件之中導(dǎo)入了CDF文件讀取相關(guān)的函數(shù)內(nèi)容，這樣子我們就可以在R#腳本之中打開并且讀取CDF文件之中的數(shù)據(jù)啦。接下來我們就可以通過下面所列舉的幾個函數(shù)來從我們的目標(biāo)CDF文件之中讀取數(shù)據(jù)了：
open.netCDF: 按照一個所設(shè)定的文件路徑，打開CDF文件流
var: 從CDF文件流之中獲取得到指定名稱的變量信息
getValue: 基于特定的變量信息讀取CDF文件之中所存儲的向量數(shù)據(jù)
對于氣質(zhì)色譜數(shù)據(jù)分析而言，CDF文件是一種比較通用的標(biāo)準(zhǔn)格式，CDF文件里面主要是存儲很多個向量數(shù)據(jù)，而氣質(zhì)色譜的總離子流相關(guān)的向量變量信息在CDF文件中分別為：
scan_acquisition_time 向量即為總離子流圖的保留時間向量信息
total_intensity 顧名思義就是總離子流圖的每一個保留時間點上所對應(yīng)的總離子信號響應(yīng)強度信息
那，現(xiàn)在了解到了上面的主要相關(guān)信息之后，我們就可以寫出如下腳本進(jìn)行數(shù)據(jù)的讀取了：
現(xiàn)在保存一下我們的腳本文件，在命令行窗口中輸入R#；接著呢，復(fù)制文件路徑，然后粘貼我們的文件路徑在命令行輸出窗口中，例如像下面這樣子：
最起始的R#就是代表我們所使用的R#腳本的解釋器程序，后面的路徑就是我們上面所保存的R腳本了?，F(xiàn)在在命令行上按下回車執(zhí)行我們的腳本，經(jīng)過一會兒之后，在我們的命令行窗口上應(yīng)該會出現(xiàn)和下圖類似的dataframe打印結(jié)果，出現(xiàn)這個就說明我們通過腳本從CDF文件之中讀取總離子流色譜圖相關(guān)的向量數(shù)據(jù)成功了！
ggplot繪制總離子流色譜圖
現(xiàn)在，我們已經(jīng)擁有了繪圖相關(guān)的數(shù)據(jù)，那現(xiàn)在我們可以做繪圖操作了。首先我們祭出R語言的繪圖利器，首先展示一下在R#編程語言環(huán)境之中通過ggplot進(jìn)行總離子流色譜圖的繪制：
通過執(zhí)行上面的一段ggplot的繪圖腳本代碼，我們得到了如下所示的TIC圖數(shù)據(jù)可視化輸出結(jié)果：
MZKit程序包繪制總離子流色譜圖
上面我們展示的是基于ggplot程序包按照一般的數(shù)據(jù)繪制流程做的一個數(shù)據(jù)可視化，那現(xiàn)在假若我們想要基于MZKit程序包進(jìn)行上面所展示的類似的TIC色譜圖的繪制操作呢？當(dāng)然也是沒問題的了，因為MZKit軟件包就是專門進(jìn)行質(zhì)譜數(shù)據(jù)處理操作的工作，繪制這種基礎(chǔ)原始數(shù)據(jù)圖當(dāng)然不在話下了，首先來看繪圖代碼：
從上面的繪圖代碼可以看見，相關(guān)的繪圖操作和ggplot相比較，自定義方面的內(nèi)容（例如坐標(biāo)軸，標(biāo)題之類的信息）會較少一些，但是相關(guān)的繪圖操作卻簡化了一些：我們只需要通過chromatogram對象相關(guān)的操作函數(shù)，將讀取出來的原始數(shù)據(jù)轉(zhuǎn)換為chromatogram對象，接下來就可以直接通過plot函數(shù)繪制出來了！除了操作簡單這個特點，默認(rèn)輸出的圖形樣式也要漂亮了許多。
下面截圖的是在這篇文章中所展示的基礎(chǔ)數(shù)據(jù)作圖功能的全部代碼了：
上面的所有代碼都是基于MZKit桌面工作站軟件最新發(fā)布的版本所完成的，如果各位同學(xué)需要更新或者安裝最新版本的MZKit桌面工作站軟件，可以訪問BioDeep網(wǎng)站進(jìn)行下載：http://www.biodeep.cn/downloads?lang=zh-CN
希望這篇評論回復(fù)文章可以解開粉絲的一部分疑惑哈