許 龍,葉 樺,楊建偉,溫瑞睿,朱文靜
(1. 國家藥品監督管理局執業藥師資格認證中心,北京 100061; 2. 復旦大學藥學院,上海 201203)
執業藥師職業資格考試報考人數多、社會影響大,考試成績和合格率均受到公眾的廣泛關注。一方面,國家規定零售藥店需配備執業藥師[1],為滿足零售藥店的發展需求,各地藥學專業人員報名參考踴躍;另一方面,參加考試的準入條件較寬泛,藥學專業(或中藥學專業)及醫學類、化學類、生物學類等相關專業,且具備一定專業工作年限的高等職業學校(大專)、大學本科及以上學歷考生均可報考。大學、高職高專和中專藥學類專業的課程設置差異較大,但《藥事管理與法規》(簡稱《藥事法規》)作為各層次藥學專業的必修課程,受專業設置影響較小,且有代表性。執業藥師應了解和遵守相關的藥事管理法規,為公眾提供安全、有效和可靠的藥學服務,促進醫藥產業的快速健康發展[2]。本研究中通過提取全國2021年度執業藥師職業資格考試藥學類考試科目的相關成績,采用隨機森林回歸模型,分析并揭示了與《藥事法規》考試成績密切相關的影響因素,以幫助藥學技術人員根據自身特點進行復習備考,并為推動培養藥學專業人才的高等學校開展專業教育改革、促進考試制度與完善相關政策提供建議。現報道如下。
根據執業藥師職業資格藥學類考試相關規定[3],考試分為藥學、中藥學2 個專業類別;藥學類考試科目包括《藥事法規》(2 個專業均有)及《藥學專業知識(一)》《藥學專業知識(二)》《藥學綜合知識與技能》(后3科以下簡稱《藥學(一)》《藥學(二)》《藥學綜合》);考試成績有效期為4年,參考人員可在期間不同年份參加不同科目考試,也可多次參考同一科目。本研究中選取全國2021 年度同時參加藥學類4 個科目考試的人員為研究對象。通過國家執業藥師職業資格考試成績數據庫提取4個科目的考試成績,以及參考人員的性別、學歷、培養方式、年齡、專業、職稱、工作年限、單位類型等資料。
由于基于數學函數構建的機器學習算法內部無法識別性別、學歷、培養方式、專業、單位性質等文本數據,故需對其進行預處理。數據處理環境為Jupyter Notebook(Python 3.9),分別采用Scikit-learn庫中的Ordinal Encoder和Label Encoder對有序和無序分類變量進行編碼,字符串轉換為整數或浮點值,矢量化文本[4];再利用Standard Scaler 工具對數據進行標準化處理,將數據特征縮放到均值為0、方差為1的標準正態分布中,以消除特征之間的量級差異,提高模型的穩定性和準確性[5]。
隨機森林回歸分析[6]:構建執業藥師《藥事法規》考試成績影響因素隨機森林回歸模型(以下簡稱模型)。采用模型中的feature_importances_函數計算特征的重要性,獲取每個特征的重要性評分,評分越高表明該特征貢獻度越高,對模型的預測性能影響越大。預試驗中剔除重要性評分極低的特征后,最終確定了執業藥師其他3 個科目考試成績、年齡、工作年限、專業、單位類型、職稱、培養方式、學歷8種特征作為《藥事法規》考試成績的影響因素。按各個特征進行分組,采用方差分析進行多組間比較,確定各因素對《藥事法規》科目成績的影響;進一步采用Tukey 真實顯著差異檢驗進行成對比較[7]。以F值在隨機分布下的概率(PR>F)值接近0為差異有統計學意義。
Pearson 相關性分析:采用Pearson 相關性分析不同科目考試成績間的相關性。相關系數范圍為[- 1,0),0,(0,1]分別為負相關、不相關、正相關。
共納入參考人員82 922 人。其中,女性考生遠多于男性;年齡主要為26~45歲;學歷主要為大學本科和大學專科。詳見表1。

表1 2021年參加執業藥師職業資格藥學類考試全部科目的人員基本情況[人(%),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]
2021 年度《藥事法規》《藥學(一)》《藥學(二)》《藥學綜合》平均分分別為59.97 分、55.66 分、61.80 分、63.51 分;合格率分別為17.96%、19.82%、34.52%、34.74%,以《藥事法規》最低,《藥學(二)》與《藥學綜合》相對較高。
通過均方根誤差(RMSE)、R2分數及袋外估計準確率全面評估模型的性能。其中,RMSE反映了模型預測值與實際值間的平均差異,R2分數衡量了模型對總方差的解釋程度,袋外估計準確率通過模型在未分析過的數據上的性能表現來評估其泛化能力[8]。上述指標的綜合考量有助于全面了解模型的效果和可靠性。本研究中,訓練集(80%數據集)的RMSE(2.957)較小,提示模型在訓練數據上的預測效果較好;R2分數為0.939,表明該模型能解釋因變量變異性的93.9%,提示模型在訓練數據上的擬合效果好(擬合度為0.861 7)。測試集(20%數據集)的RMSE(8.030)、R2分數(0.551)及袋外估計準確率評分(0.552)結果顯示,模型的預測精度較高,但存在一定誤差。綜上所述,所選模型的預測結果與實際值間的誤差處于可接受范圍,即該模型性能較好。
另3 科考試成績對《藥事法規》考試成績的影響較大(隨機森林回歸影響系數相對較大),提示個人的學習能力是影響成績的關鍵因素,其次為年齡和工作年限,影響較小的因素為單位類型、專業、職稱、培養方式、學歷。詳見表2。

表2 《藥事法規》成績影響因素Tab.2 Factors affecting the score of Pharmacy Administration and Regulations
各科目間的相關性:Pearson 相關性分析結果顯示,《藥事法規》與《藥學(一)》考試成績的相關性強于其與他科目考試成績的相關性(見圖1)。從各科目知識結構來看,可能由于《藥事法規》與《藥學(一)》均屬基礎知識較多且較難掌握的科目,故更具挑戰性,學生的表現更加一致。整體來看,上述結果表明了藥學各傳統課程間的關聯度較高,相輔相成,有助于學生在多個方面的綜合發展。

圖1 執業藥師職業資格藥學類考試科目成績Pearson相關系數熱力圖Fig.1 Heatmap of Pearson correlation coefficient for scores of pharmacy subjects in the professional qualification examination for licensed pharmacists
年齡:參考人員主要集中于26~45歲。各年齡段參考人員《藥事法規》平均分數在50 歲之前及50 歲以后隨年齡增加分別呈升高和下降趨勢(見圖2)。

圖2 各年齡段參考人員《藥事法規》成績分布Fig.2 Distribution of scores in Pharmacy Administration and Regulations among examinees of different ages
工作年限:參考人員工作年限主要集中于0~15年。參考人員工作年限為0~5 年時《藥事法規》分數較高;6~10 年時出現較明顯下降,但此后直至40 年,總體呈升高趨勢;超過40年后該科分數出現顯著下降。當工作年限達到35 年之后,參考人員群體分數分布存在正偏斜,可能是由于較高分數的參考人員數量較多,較低分數的數量較少,且由于樣本量相對較少,數據代表性不夠,引入了一定的偏差,導致該年齡段的分數表現與其他年齡段不同(見圖3)。

圖3 各工作年限參考人員《藥事法規》成績分布Fig.3 Distribution of scores in Pharmacy Administration and Regulations among examinees of different years of work experience
專業:參考人員的專業主要為藥學類和醫學類。中藥學類、中醫學類專業的參考人員的平均分數均高于醫學類專業(見圖4)。

圖4 各專業參考人員《藥事法規》成績分布Fig.4 Distribution of scores in Pharmacy Administration and Regulations among examinees of different majors
單位類型:參考人員主要分布于藥品零售、批發企業。藥品零售、藥品批發、藥品生產企業參考人員的平均分數較醫療衛生機構和其他領域高(見圖5)。

圖5 各單位類型參考人員《藥事》平均成績Fig.5 Distribution of scores in Pharmacy Administration and Regulations among examinees from different types of employing organizations
職稱:多數參考人員無職稱,與多數執業藥師在藥品零售企業不參與職稱評定有關。《藥事法規》平均分數與參考人員的職稱呈正相關;但無職稱參考人員群體平均分數也較高(見圖6)。

圖6 不同職稱參考人員《藥事法規》成績分布Fig.6 Distribution of scores in Pharmacy Administration and Regulations among examinees with different professional titles
培養方式:全日制教育和非全日制教育參考人數相當,但接受全日制教育參考人員的平均分數相對更高(見圖7)。

圖7 不同培養方式參考人員《藥事法規》成績分布Fig.7 Distribution of scores in Pharmacy Administration and Regulations among examinees with different training methods
學歷:參考人員以大學專科為主,大學本科次之。該科目平均分數與學歷基本呈正相關(見圖8)。

圖8 不同學歷參考人員《藥事法規》成績分布Fig.8 Average scores of Pharmacy Administration and Regulations in examinees with different educational backgrounds
方差分析:《藥學(一)》《藥學(二)》《藥學綜合》考試成績、年齡、工作年限、專業、單位類型、職稱、培養方式、學歷分組數據F值較大且均大于F臨界,PR(>F)顯著性水平接近于0,提示上述因素對《藥事法規》分數的影響均顯著(見表3)。

表3 方差分析結果Tab.3 Results of analysis of variance
《藥事法規》是一門綜合性學科,覆蓋藥品的生產、經營和使用環節,內容包括藥品監管法律法規和政策,以考試作為導向,要求參考人員掌握必備的藥事管理與法規知識,合法合規地從事執業活動。學好藥事管理與法規不但需要參考人員認真學習和記憶,還需結合工作經驗加以理解掌握,故系統的藥學教育及適當的工作年限有助于對藥事管理與法規知識融會貫通。
4 個考試科目之間關系密切,各科目知識點相輔相成,故參考人員備考時需系統學習,尤其是未接受過藥學全日制教育的參考人員更需全面綜合掌握各科目的內容。《藥事法規》分數與參考人員的工作年限和年齡相關。隨著工作年限和年齡的增加,工作經驗不斷豐富,考生對藥事管理與法規的理解程度逐步加深。但工作年限滿40 年或年齡滿50 歲以后,多數人員即將退休,備考積極性減弱,其學習能力和記憶力均有所下降。此外,執業藥師職業資格藥學類考試對不同學歷、專業有0~5年的工作年限要求。由于藥學本專業、學歷較高的參考人員比相關專業、學歷較低的人員工作年限要求相應縮短,故工作年限為0~5 年的人員組成專業度和學歷均更高,其《藥事法規》分數也更高。
中藥學類和中醫學類參考人員通常是在已具備中藥類執業藥師資格的情況下報名藥學類執業藥師職業資格考試,其中取得藥學類執業藥師職業資格的,可注冊為藥學與中藥學類執業藥師(俗稱“雙證執業藥師”)[9]。這類參考人員已參加過《藥事法規》科目考試并取得合格成績,故平均分數比藥學類和醫學類參考人員略高。
醫學類和在醫療機構工作的參考人員由于工作崗位需求、專業背景等方面原因,對《藥事法規》知識的掌握較薄弱,故在備考執業藥師職業資格藥學類考試時要重點加強該方面學習,從而進一步提高藥學服務水平[10]。
根據本研究中所得結果,建議開展考前輔導的施教機構根據受教育人群的特點因材施教,針對不同學歷、專業等人員進行課程設計,從而提升課程的針對性和有效性[11]。各藥學專業高等學校要按實際需要開展藥學專業教育,進一步加強藥學服務型人才培養,從而為公眾提供更優質的藥學服務。