許 龍，葉 樺，楊建偉，溫瑞睿，朱文靜

（1. 國家藥品監(jiān)督管理局執(zhí)業(yè)藥師資格認證中心，北京 100061； 2. 復旦大學藥學院，上海 201203）

執(zhí)業(yè)藥師職業(yè)資格考試報考人數(shù)多、社會影響大，考試成績和合格率均受到公眾的廣泛關(guān)注。一方面，國家規(guī)定零售藥店需配備執(zhí)業(yè)藥師［1］，為滿足零售藥店的發(fā)展需求，各地藥學專業(yè)人員報名參考踴躍；另一方面，參加考試的準入條件較寬泛，藥學專業(yè)（或中藥學專業(yè)）及醫(yī)學類、化學類、生物學類等相關(guān)專業(yè)，且具備一定專業(yè)工作年限的高等職業(yè)學校（大專）、大學本科及以上學歷考生均可報考。大學、高職高專和中專藥學類專業(yè)的課程設置差異較大，但《藥事管理與法規(guī)》（簡稱《藥事法規(guī)》）作為各層次藥學專業(yè)的必修課程，受專業(yè)設置影響較小，且有代表性。執(zhí)業(yè)藥師應了解和遵守相關(guān)的藥事管理法規(guī)，為公眾提供安全、有效和可靠的藥學服務，促進醫(yī)藥產(chǎn)業(yè)的快速健康發(fā)展［2］。本研究中通過提取全國2021年度執(zhí)業(yè)藥師職業(yè)資格考試藥學類考試科目的相關(guān)成績，采用隨機森林回歸模型，分析并揭示了與《藥事法規(guī)》考試成績密切相關(guān)的影響因素，以幫助藥學技術(shù)人員根據(jù)自身特點進行復習備考，并為推動培養(yǎng)藥學專業(yè)人才的高等學校開展專業(yè)教育改革、促進考試制度與完善相關(guān)政策提供建議?，F(xiàn)報道如下。

1 資料與方法

1.1 研究對象與數(shù)據(jù)提取

根據(jù)執(zhí)業(yè)藥師職業(yè)資格藥學類考試相關(guān)規(guī)定［3］，考試分為藥學、中藥學2 個專業(yè)類別；藥學類考試科目包括《藥事法規(guī)》（2 個專業(yè)均有）及《藥學專業(yè)知識（一）》《藥學專業(yè)知識（二）》《藥學綜合知識與技能》（后3科以下簡稱《藥學（一）》《藥學（二）》《藥學綜合》）；考試成績有效期為4年，參考人員可在期間不同年份參加不同科目考試，也可多次參考同一科目。本研究中選取全國2021 年度同時參加藥學類4 個科目考試的人員為研究對象。通過國家執(zhí)業(yè)藥師職業(yè)資格考試成績數(shù)據(jù)庫提取4個科目的考試成績，以及參考人員的性別、學歷、培養(yǎng)方式、年齡、專業(yè)、職稱、工作年限、單位類型等資料。

1.2 數(shù)據(jù)清洗

由于基于數(shù)學函數(shù)構(gòu)建的機器學習算法內(nèi)部無法識別性別、學歷、培養(yǎng)方式、專業(yè)、單位性質(zhì)等文本數(shù)據(jù)，故需對其進行預處理。數(shù)據(jù)處理環(huán)境為Jupyter Notebook（Python 3.9），分別采用Scikit-learn庫中的Ordinal Encoder和Label Encoder對有序和無序分類變量進行編碼，字符串轉(zhuǎn)換為整數(shù)或浮點值，矢量化文本［4］；再利用Standard Scaler 工具對數(shù)據(jù)進行標準化處理，將數(shù)據(jù)特征縮放到均值為0、方差為1的標準正態(tài)分布中，以消除特征之間的量級差異，提高模型的穩(wěn)定性和準確性［5］。

1.3 分析方法

隨機森林回歸分析［6］：構(gòu)建執(zhí)業(yè)藥師《藥事法規(guī)》考試成績影響因素隨機森林回歸模型（以下簡稱模型）。采用模型中的feature_importances_函數(shù)計算特征的重要性，獲取每個特征的重要性評分，評分越高表明該特征貢獻度越高，對模型的預測性能影響越大。預試驗中剔除重要性評分極低的特征后，最終確定了執(zhí)業(yè)藥師其他3 個科目考試成績、年齡、工作年限、專業(yè)、單位類型、職稱、培養(yǎng)方式、學歷8種特征作為《藥事法規(guī)》考試成績的影響因素。按各個特征進行分組，采用方差分析進行多組間比較，確定各因素對《藥事法規(guī)》科目成績的影響；進一步采用Tukey 真實顯著差異檢驗進行成對比較［7］。以F值在隨機分布下的概率（PR>F）值接近0為差異有統(tǒng)計學意義。

Pearson 相關(guān)性分析：采用Pearson 相關(guān)性分析不同科目考試成績間的相關(guān)性。相關(guān)系數(shù)范圍為［- 1，0），0，（0，1］分別為負相關(guān)、不相關(guān)、正相關(guān)。

2 結(jié)果

2.1 參考人員基本情況

共納入?yún)⒖既藛T82 922 人。其中，女性考生遠多于男性；年齡主要為26～45歲；學歷主要為大學本科和大學?？啤Ｔ斠姳?。

表1 2021年參加執(zhí)業(yè)藥師職業(yè)資格藥學類考試全部科目的人員基本情況［人（%），n=82 922］Tab.1 Basic information of examinees participated in all pharmacy subjects of the professional qualification examination for licensed pharmacists in 2021［examinee（%），n = 82 922］

2.2 各科目考試成績

2021 年度《藥事法規(guī)》《藥學（一）》《藥學（二）》《藥學綜合》平均分分別為59.97 分、55.66 分、61.80 分、63.51 分；合格率分別為17.96%、19.82%、34.52%、34.74%，以《藥事法規(guī)》最低，《藥學（二）》與《藥學綜合》相對較高。

2.3 模型效果分析

通過均方根誤差（RMSE）、R2分數(shù)及袋外估計準確率全面評估模型的性能。其中，RMSE反映了模型預測值與實際值間的平均差異，R2分數(shù)衡量了模型對總方差的解釋程度，袋外估計準確率通過模型在未分析過的數(shù)據(jù)上的性能表現(xiàn)來評估其泛化能力［8］。上述指標的綜合考量有助于全面了解模型的效果和可靠性。本研究中，訓練集（80%數(shù)據(jù)集）的RMSE（2.957）較小，提示模型在訓練數(shù)據(jù)上的預測效果較好；R2分數(shù)為0.939，表明該模型能解釋因變量變異性的93.9%，提示模型在訓練數(shù)據(jù)上的擬合效果好（擬合度為0.861 7）。測試集（20%數(shù)據(jù)集）的RMSE（8.030）、R2分數(shù)（0.551）及袋外估計準確率評分（0.552）結(jié)果顯示，模型的預測精度較高，但存在一定誤差。綜上所述，所選模型的預測結(jié)果與實際值間的誤差處于可接受范圍，即該模型性能較好。

2.4 各因素對《藥事法規(guī)》考試成績的影響

另3 科考試成績對《藥事法規(guī)》考試成績的影響較大（隨機森林回歸影響系數(shù)相對較大），提示個人的學習能力是影響成績的關(guān)鍵因素，其次為年齡和工作年限，影響較小的因素為單位類型、專業(yè)、職稱、培養(yǎng)方式、學歷。詳見表2。

表2 《藥事法規(guī)》成績影響因素Tab.2 Factors affecting the score of Pharmacy Administration and Regulations

各科目間的相關(guān)性：Pearson 相關(guān)性分析結(jié)果顯示，《藥事法規(guī)》與《藥學（一）》考試成績的相關(guān)性強于其與他科目考試成績的相關(guān)性（見圖1）。從各科目知識結(jié)構(gòu)來看，可能由于《藥事法規(guī)》與《藥學（一）》均屬基礎知識較多且較難掌握的科目，故更具挑戰(zhàn)性，學生的表現(xiàn)更加一致。整體來看，上述結(jié)果表明了藥學各傳統(tǒng)課程間的關(guān)聯(lián)度較高，相輔相成，有助于學生在多個方面的綜合發(fā)展。

圖1 執(zhí)業(yè)藥師職業(yè)資格藥學類考試科目成績Pearson相關(guān)系數(shù)熱力圖Fig.1 Heatmap of Pearson correlation coefficient for scores of pharmacy subjects in the professional qualification examination for licensed pharmacists

年齡：參考人員主要集中于26～45歲。各年齡段參考人員《藥事法規(guī)》平均分數(shù)在50 歲之前及50 歲以后隨年齡增加分別呈升高和下降趨勢（見圖2）。

圖2 各年齡段參考人員《藥事法規(guī)》成績分布Fig.2 Distribution of scores in Pharmacy Administration and Regulations among examinees of different ages

工作年限：參考人員工作年限主要集中于0～15年。參考人員工作年限為0～5 年時《藥事法規(guī)》分數(shù)較高；6～10 年時出現(xiàn)較明顯下降，但此后直至40 年，總體呈升高趨勢；超過40年后該科分數(shù)出現(xiàn)顯著下降。當工作年限達到35 年之后，參考人員群體分數(shù)分布存在正偏斜，可能是由于較高分數(shù)的參考人員數(shù)量較多，較低分數(shù)的數(shù)量較少，且由于樣本量相對較少，數(shù)據(jù)代表性不夠，引入了一定的偏差，導致該年齡段的分數(shù)表現(xiàn)與其他年齡段不同（見圖3）。

圖3 各工作年限參考人員《藥事法規(guī)》成績分布Fig.3 Distribution of scores in Pharmacy Administration and Regulations among examinees of different years of work experience

專業(yè)：參考人員的專業(yè)主要為藥學類和醫(yī)學類。中藥學類、中醫(yī)學類專業(yè)的參考人員的平均分數(shù)均高于醫(yī)學類專業(yè)（見圖4）。

圖4 各專業(yè)參考人員《藥事法規(guī)》成績分布Fig.4 Distribution of scores in Pharmacy Administration and Regulations among examinees of different majors

單位類型：參考人員主要分布于藥品零售、批發(fā)企業(yè)。藥品零售、藥品批發(fā)、藥品生產(chǎn)企業(yè)參考人員的平均分數(shù)較醫(yī)療衛(wèi)生機構(gòu)和其他領(lǐng)域高（見圖5）。

圖5 各單位類型參考人員《藥事》平均成績Fig.5 Distribution of scores in Pharmacy Administration and Regulations among examinees from different types of employing organizations

職稱：多數(shù)參考人員無職稱，與多數(shù)執(zhí)業(yè)藥師在藥品零售企業(yè)不參與職稱評定有關(guān)?！端幨路ㄒ?guī)》平均分數(shù)與參考人員的職稱呈正相關(guān)；但無職稱參考人員群體平均分數(shù)也較高（見圖6）。

圖6 不同職稱參考人員《藥事法規(guī)》成績分布Fig.6 Distribution of scores in Pharmacy Administration and Regulations among examinees with different professional titles

培養(yǎng)方式：全日制教育和非全日制教育參考人數(shù)相當，但接受全日制教育參考人員的平均分數(shù)相對更高（見圖7）。

圖7 不同培養(yǎng)方式參考人員《藥事法規(guī)》成績分布Fig.7 Distribution of scores in Pharmacy Administration and Regulations among examinees with different training methods

學歷：參考人員以大學?？茷橹?，大學本科次之。該科目平均分數(shù)與學歷基本呈正相關(guān)（見圖8）。

圖8 不同學歷參考人員《藥事法規(guī)》成績分布Fig.8 Average scores of Pharmacy Administration and Regulations in examinees with different educational backgrounds

方差分析：《藥學（一）》《藥學（二）》《藥學綜合》考試成績、年齡、工作年限、專業(yè)、單位類型、職稱、培養(yǎng)方式、學歷分組數(shù)據(jù)F值較大且均大于F臨界，PR（>F）顯著性水平接近于0，提示上述因素對《藥事法規(guī)》分數(shù)的影響均顯著（見表3）。

表3 方差分析結(jié)果Tab.3 Results of analysis of variance

3 討論

《藥事法規(guī)》是一門綜合性學科，覆蓋藥品的生產(chǎn)、經(jīng)營和使用環(huán)節(jié)，內(nèi)容包括藥品監(jiān)管法律法規(guī)和政策，以考試作為導向，要求參考人員掌握必備的藥事管理與法規(guī)知識，合法合規(guī)地從事執(zhí)業(yè)活動。學好藥事管理與法規(guī)不但需要參考人員認真學習和記憶，還需結(jié)合工作經(jīng)驗加以理解掌握，故系統(tǒng)的藥學教育及適當?shù)墓ぷ髂晗抻兄趯λ幨鹿芾砼c法規(guī)知識融會貫通。

4 個考試科目之間關(guān)系密切，各科目知識點相輔相成，故參考人員備考時需系統(tǒng)學習，尤其是未接受過藥學全日制教育的參考人員更需全面綜合掌握各科目的內(nèi)容?！端幨路ㄒ?guī)》分數(shù)與參考人員的工作年限和年齡相關(guān)。隨著工作年限和年齡的增加，工作經(jīng)驗不斷豐富，考生對藥事管理與法規(guī)的理解程度逐步加深。但工作年限滿40 年或年齡滿50 歲以后，多數(shù)人員即將退休，備考積極性減弱，其學習能力和記憶力均有所下降。此外，執(zhí)業(yè)藥師職業(yè)資格藥學類考試對不同學歷、專業(yè)有0～5年的工作年限要求。由于藥學本專業(yè)、學歷較高的參考人員比相關(guān)專業(yè)、學歷較低的人員工作年限要求相應縮短，故工作年限為0～5 年的人員組成專業(yè)度和學歷均更高，其《藥事法規(guī)》分數(shù)也更高。

中藥學類和中醫(yī)學類參考人員通常是在已具備中藥類執(zhí)業(yè)藥師資格的情況下報名藥學類執(zhí)業(yè)藥師職業(yè)資格考試，其中取得藥學類執(zhí)業(yè)藥師職業(yè)資格的，可注冊為藥學與中藥學類執(zhí)業(yè)藥師（俗稱“雙證執(zhí)業(yè)藥師”）［9］。這類參考人員已參加過《藥事法規(guī)》科目考試并取得合格成績，故平均分數(shù)比藥學類和醫(yī)學類參考人員略高。

醫(yī)學類和在醫(yī)療機構(gòu)工作的參考人員由于工作崗位需求、專業(yè)背景等方面原因，對《藥事法規(guī)》知識的掌握較薄弱，故在備考執(zhí)業(yè)藥師職業(yè)資格藥學類考試時要重點加強該方面學習，從而進一步提高藥學服務水平［10］。

根據(jù)本研究中所得結(jié)果，建議開展考前輔導的施教機構(gòu)根據(jù)受教育人群的特點因材施教，針對不同學歷、專業(yè)等人員進行課程設計，從而提升課程的針對性和有效性［11］。各藥學專業(yè)高等學校要按實際需要開展藥學專業(yè)教育，進一步加強藥學服務型人才培養(yǎng)，從而為公眾提供更優(yōu)質(zhì)的藥學服務。