TWI812056B - 檢查藥物相互作用的方法和電子裝置 - Google Patents
檢查藥物相互作用的方法和電子裝置 Download PDFInfo
- Publication number
- TWI812056B TWI812056B TW111108870A TW111108870A TWI812056B TW I812056 B TWI812056 B TW I812056B TW 111108870 A TW111108870 A TW 111108870A TW 111108870 A TW111108870 A TW 111108870A TW I812056 B TWI812056 B TW I812056B
- Authority
- TW
- Taiwan
- Prior art keywords
- medication
- electronic device
- medication combination
- drug
- score
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 206010013710 Drug interaction Diseases 0.000 title claims abstract description 12
- 239000003814 drug Substances 0.000 claims abstract description 457
- 229940079593 drug Drugs 0.000 claims abstract description 457
- 239000000890 drug combination Substances 0.000 claims abstract description 37
- 230000004044 response Effects 0.000 claims abstract description 18
- 230000000875 corresponding effect Effects 0.000 claims description 129
- 239000013598 vector Substances 0.000 claims description 87
- 238000009826 distribution Methods 0.000 claims description 52
- 238000012216 screening Methods 0.000 claims description 14
- 230000002596 correlated effect Effects 0.000 claims description 8
- 238000010586 diagram Methods 0.000 description 11
- 230000003993 interaction Effects 0.000 description 6
- 238000002483 medication Methods 0.000 description 5
- 230000003321 amplification Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 206010067484 Adverse reaction Diseases 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000006838 adverse reaction Effects 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H70/00—ICT specially adapted for the handling or processing of medical references
- G16H70/40—ICT specially adapted for the handling or processing of medical references relating to drugs, e.g. their side effects or intended usage
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
- G16H10/60—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/10—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/70—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients
Abstract
提出一種檢查藥物相互作用的方法和電子裝置。方法包含:根據多個病歷產生第一用藥組合以及住院事件之間的第一勝算比、第二用藥組合以及住院事件之間的第二勝算比以及第三用藥組合以及住院事件之間的第三勝算比;根據第二勝算比產生對應於第一藥物的第一分數;根據第三勝算比產生對應於第二藥物的第二分數;以及響應於第一勝算比大於第一閾值,第一分數與第二分數的總和大於第二閾值,並且第一分數與第二分數的商數小於第三閾值,輸出第一用藥組合。
Description
本發明是有關於一種檢查藥物相互作用的方法和電子裝置。
病患時常需要在同一段時間內服用多種藥物。這些藥物之間可能會產生交互作用導致嚴重的不良反應,致使病患非預期住院。為了避免上述的情況發生,須對各種用藥組合進行相互作用的檢查。然而,用藥組合的數量龐大,逐一檢查各種用藥組合是非常沒有效率的。因此,如何提出一種能快速地檢查出具有高風險之用藥組合的方法,是本領域人員致力的目標之一。
本發明提供一種檢查藥物相互作用的方法和電子裝置,可輸出具有高風險的用藥組合以供使用者參考。
本發明的一種檢查藥物相互作用的方法,包含:取得多個病歷,其中多個病歷的至少其中之一指示服用第一用藥組合的病患是否發生住院事件;根據多個病歷產生用藥組合集合,其中用藥組合集合包含第一用藥組合、第二用藥組合以及第三用藥組合,其中第一用藥組合與第二用藥組合均包含第一藥物,並且第一用藥組合與第三用藥組合均包含第二藥物;根據多個病歷產生第一用藥組合以及住院事件之間的第一勝算比、第二用藥組合以及住院事件之間的第二勝算比以及第三用藥組合以及住院事件之間的第三勝算比;根據第二勝算比產生對應於第一藥物的第一分數,其中第一分數與第二勝算比負相關;根據第三勝算比產生對應於第二藥物的第二分數,其中第二分數與第三勝算比負相關,其中第一分數大於或等於第二分數;以及響應於第一勝算比大於第一閾值,第一分數與第二分數的總和大於第二閾值,並且第一分數與第二分數的商數小於第三閾值,輸出第一用藥組合。
在本發明的一實施例中,上述的根據第二勝算比產生對應於第一藥物的第一分數的步驟包含:響應於第二勝算比大於風險閾值而標記第二用藥組合;根據標記的第二用藥組合產生第三分數,其中第三分數等於用藥組合集合中包含第一藥物但不包含第二藥物且被標記的用藥組合的數量除以用藥組合集合中包含第一藥物但不包含第二藥物的用藥組合的數量;以及根據第三分數計算第一分數,其中第一分數與第三分數的總和等於一。
在本發明的一實施例中,上述的根據多個病歷產生用藥組合集合的步驟包含:執行篩選流程以產生第一獨有用藥組合集合,包含:根據多個病歷以及隱含狄利克雷分布模型產生包含第一主題向量的K個主題向量,其中K為第一主題數量,其中K個主題向量分別對應於K個主題,其中K個主題包含對應於第一主題向量的第一主題,其中第一主題向量包含所有用藥組合的機率分布;自具有最大機率的用藥組合開始,從第一主題向量中選出多個重要用藥組合以產生第一重要用藥組合集合;以及根據第一重要用藥組合集合決定第一獨有用藥組合集合;以及根據第一獨有用藥組合集合產生用藥組合集合。
在本發明的一實施例中,上述的K個主題包含第二主題,其中根據第一重要用藥組合集合決定第一獨有用藥組合集合的步驟包含:響應於第一重要用藥組合包含於與第一主題相對應的第一重要用藥組合集合以及與第二主題相對應的第二重要用藥組合集合,將第一重要用藥組合自第一重要用藥組合集合刪除以產生第一獨有用藥組合集合。
在本發明的一實施例中,上述的根據第一獨有用藥組合集合產生用藥組合集合的步驟包含:重覆執行多次篩選流程以產生包含第一獨有用藥組合集合的多個獨有用藥組合集合;響應於多個獨有用藥組合集合中的第一用藥組合的數量大於數量閾值,根據第一用藥組合以產生對應於第一主題的第一穩定用藥組合集合;以及根據第一穩定用藥組合集合產生用藥組合集合。
在本發明的一實施例中,上述的根據第一穩定用藥組合集合產生用藥組合集合的步驟包含:根據多個病歷以及隱含狄利克雷分布模型產生分別對應於多個病歷的多個病歷向量,其中多個病歷向量的每一者包含K個主題的機率分布;根據K個主題的機率分布判斷多個病歷中與第一主題相對應的病歷集合;計算病歷集合中的至少一病歷佔病歷集合的比率,其中至少一病歷指示第一穩定用藥組合集合中的至少一用藥組合;以及響應於比率大於比率閾值,根據第一穩定用藥組合集產生用藥組合集合,其中用藥組合集合包含第一穩定用藥組合集中的多個用藥組合。
在本發明的一實施例中,上述的病歷集合中的第一病歷對應於K個主題的第一機率分布,其中根據K個主題的機率分布判斷多個病歷中與第一主題相對應的病歷集合的步驟包含:響應於第一機率分布中的最大機率對應於第一主題,判斷第一病歷對應於第一主題。
在本發明的一實施例中,上述的方法更包含:根據多個病歷以及隱含狄利克雷分布模型產生對應於第一主題數量的第一指標以及對應於第二主題數量的第二指標;以及比較第一指標與第二指標以從第一主題數量與第二主題數量中選出第一主題數量以作為K。
在本發明的一實施例中,上述的產生對應於第一主題數量的第一指標的步驟包含:根據多個病歷、隱含狄利克雷分布模型以及第一主題數量產生K個主題向量;以及計算K個主題向量的所有2-組合的平均相似度以作為第一指標。
在本發明的一實施例中,上述的產生對應於第一主題數量的第一指標的步驟包含:根據多個病歷、隱含狄利克雷分布模型以及第一主題數量產生分別對應於多個病歷的多個病歷向量,其中多個病歷向量的每一者包含K個主題的機率分布;根據K個主題的機率分布判斷多個病歷中與第一主題相對應的至少一病歷;以及根據至少一病歷的個數與多個病歷的總數計算比率以作為第一指標。
在本發明的一實施例中,上述的根據K個主題的機率分布判斷多個病歷中與第一主題相對應的至少一病歷的步驟包含:自多個病歷向量取得對應於至少一病歷的K個主題的第一機率分布;以及響應於第一機率分布中的最大機率對應於第一主題且大於機率閾值,判斷至少一病歷對應於第一主題。
在本發明的一實施例中,上述的產生對應於第一主題數量的第一指標的步驟包含:根據多個病歷、隱含狄利克雷分布模型以及第一主題數量產生分別對應於多個病歷的多個病歷向量,其中多個病歷向量的每一者包含K個主題的機率分布;根據K個主題的機率分布將多個病歷分為K個群組,其中K個群組分別對應於K個主題;根據K個群組計算群組間距離的第一統計值;根據K個群組計算群組內距離的第二統計值;以及計算第一統計值與第二統計值的比率以作為第一指標。
在本發明的一實施例中,上述的根據K個群組計算群組間距離的第一統計值的步驟包含:計算K個主題向量之間的多個距離;以及將多個距離相加以求得第一統計值。
在本發明的一實施例中,上述的K個群組包含第一群組以及第二群組,其中根據K個群組計算群組內距離的第二統計值的步驟包含:計算第一群組中的多個元素之間的多個距離以產生對應於第一群組的第一群組內距離總和;以及將對應於第一群組的第一群組內距離總和與對應於第二群組的第二群組內距離總和相加以求得第二統計值。
本發明的一種檢查藥物相互作用的電子裝置,包含處理器以及收發器。處理器耦接收發器並且經配置以執行:通過收發器取得多個病歷,其中多個病歷的至少其中之一指示服用第一用藥組合的病患是否發生住院事件;根據多個病歷產生用藥組合集合,其中用藥組合集合包含第一用藥組合、第二用藥組合以及第三用藥組合,其中第一用藥組合與第二用藥組合均包含第一藥物,並且第一用藥組合與第三用藥組合均包含第二藥物;根據多個病歷產生第一用藥組合以及住院事件之間的第一勝算比、第二用藥組合以及住院事件之間的第二勝算比以及第三用藥組合以及住院事件之間的第三勝算比;根據第二勝算比產生對應於第一藥物的第一分數,其中第一分數與第二勝算比負相關;根據第三勝算比產生對應於第二藥物的第二分數,其中第二分數與第三勝算比負相關,其中第一分數大於或等於第二分數;以及響應於第一勝算比大於第一閾值,第一分數與第二分數的總和大於第二閾值,並且第一分數與第二分數的商數小於第三閾值,通過收發器輸出第一用藥組合。
基於上述,本發明可從眾多的用藥組合中篩選出具有高風險的用藥組合,並可確認所述用藥組合成為高風險的原因並非來自於用藥組合中的藥物本身,而是來自於藥物的交互作用。
為了使本發明之內容可以被更容易明瞭,以下特舉實施例作為本發明確實能夠據以實施的範例。另外,凡可能之處,在圖式及實施方式中使用相同標號的元件/構件/步驟,係代表相同或類似部件。
圖1根據本發明的一實施例繪示一種檢查藥物交互作用的電子裝置100的示意圖。電子裝置100可包含處理器110、儲存媒體120以及收發器130。
處理器110例如是中央處理單元(central processing unit,CPU),或是其他可程式化之一般用途或特殊用途的微控制單元(micro control unit,MCU)、微處理器(microprocessor)、數位信號處理器(digital signal processor,DSP)、可程式化控制器、特殊應用積體電路(application specific integrated circuit,ASIC)、圖形處理器(graphics processing unit,GPU)、影像訊號處理器(image signal processor,ISP)、影像處理單元(image processing unit,IPU)、算數邏輯單元(arithmetic logic unit,ALU)、複雜可程式邏輯裝置(complex programmable logic device,CPLD)、現場可程式化邏輯閘陣列(field programmable gate array,FPGA)或其他類似元件或上述元件的組合。處理器110可耦接至儲存媒體120以及收發器130,並且存取和執行儲存於儲存媒體120中的多個模組和各種應用程式。
儲存媒體120例如是任何型態的固定式或可移動式的隨機存取記憶體(random access memory,RAM)、唯讀記憶體(read-only memory,ROM)、快閃記憶體(flash memory)、硬碟(hard disk drive,HDD)、固態硬碟(solid state drive,SSD)或類似元件或上述元件的組合,而用於儲存可由處理器110執行的多個模組或各種應用程式。
收發器130以無線或有線的方式傳送及接收訊號。收發器130還可以執行例如低噪聲放大、阻抗匹配、混頻、向上或向下頻率轉換、濾波、放大以及類似的操作。
處理器110可通過收發器130取得分別對應於N位病患的N個病歷,其中N為正整數。病歷可指示病患曾經服用過的用藥組合,其中一用藥組合可包含兩種藥物。舉例來說,用藥組合(α,β)可包含藥物α以及藥物β。病歷還可指示病患是否發生非預期的住院事件。表1為N個病歷的示意圖。以病歷#1為例,病歷#1指示對應於病歷#1的病患曾經服用過藥物A和藥物B的組合(A,B)、藥物X和藥物Y的組合(X,Y)以及藥物U和藥物W的組合(U,W)。病歷#1還指示對應於病歷#1的病患曾經發生非預期的住院事件。也就是說,病患曾經服用過的用藥組合可能導致非預期的住院事件發生。以病歷#N為例,病歷#N指示對應於病歷#N的病患曾經服用過藥物C和藥物D的組合(C,D)以及藥物O和藥物P的組合(O,P)。病歷#N還指示對應於病歷#N的病患未曾發生非預期的住院事件。也就是說,病患曾經服用過的用藥組合並未導致非預期的住院事件發生。病歷上的用藥組合例如是以解剖學治療學及化學分類(anatomical therapeutic chemical,ATC)代碼的形式記載的。
表1
病歷 | 用藥組合 | 住院事件 |
#1 | (A,B),(X,Y),(U,W) | 是 |
#2 | (G,H) | 是 |
… | … | … |
#N | (C,D),(O,P) | 否 |
在本實施例中,假設N個病歷共記載了M個用藥組合,其中M為正整數。以表1為例,M個用藥組合至少包含用藥組合(A,B)、(X,Y)、(U,W)、(G,H)、(C,D)以及(O,P)。M個用藥組合即為N個病例中的所有用藥組合的聯集。
處理器110可根據N個病歷產生包含多個用藥組合的用藥組合集合,進而從用藥組合集合挑選出具有高風險交互作用的用藥組合以供使用者參考。首先,處理器110可使用隱含狄利克雷分布(Latent Dirichlet allocation,LDA)模型分析N個病歷。LDA模型的參數可包含主題數量K,其中K為正整數。K可決定LDA模型所產生的輸出與K個主題有關。處理器110可先決定最佳主題數量Kopt的值。
具體來說,LDA模型的輸出可關聯於主題(topic)和詞語(word)。處理器110可將N個病歷輸入至LDA模型以產生分別對應於K個主題(或稱為用藥模式)的K個主題向量。每一個主題向量可包含所有用藥組合(即:M個用藥組合)的機率分布。用藥組合即為LDA模型的詞語。主題向量包含所有詞語的機率分布(word distribution)。換句話說,主題向量可以是包含M個機率的向量,其中M個機率分別對應於M個用藥組合。
表2為K個主題向量的範例。以對應於主題#1的主題向量為例,主題向量可至少包含對應於用藥組合(A,B)的機率值「0.20」、對應於用藥組合(C,D)的機率值「0.05」以及對應於用藥組合(X,Y)的機率值「0.20」。主題向量中的所有元素(即:M個機率)的總和等於「1」。
表2
編號 | 用藥組合 | 主題#1 | 主題#2 | … | 主題#K |
#1 | (A,B) | 0.20 | 0.00 | … | 0.10 |
#2 | (C,D) | 0.05 | 0.10 | … | 0.20 |
… | … | … | … | … | … |
#M | (X,Y) | 0.20 | 0.15 | … | 0.05 |
另一方面,LDA模型還可根據N個病歷產生分別對應於N個病歷的N個病歷向量。每一個病歷向量可包含K個主題的機率分布(topic distribution)。換句話說,病歷向量可以是包含K個機率的向量,其中K個機率分別對應於K個主題。表3為N個病歷向量的範例。以對應於病歷#1的病歷向量為例,病歷向量可至少包含對應於主題#1的機率值「0.20」、對應於主題#2的機率值「0.00」以及對應於主題#K的機率值「0.10」。病歷向量中的所有元素(即:K個機率)的總和等於「1」。
表3
病歷 | 主題#1 | 主題#2 | … | 主題#K |
#1 | 0.20 | 0.00 | … | 0.10 |
#2 | 0.05 | 0.10 | … | 0.20 |
… | … | … | … | … |
#N | 0.00 | 0.05 | … | 0.15 |
處理器110可根據主題之間的相似度、具有偏重主題(即:具有偏重用藥模式)之病歷的比率以及分群效能指標等因素來決定最佳主題數量Kopt的值。
為了找出各種不同的用藥模式,主題與主題之間的差異越大越好。也就是說,主題與主題之間的相似度越低越好。在一實施例中,處理器110可計算K個主題向量的所有2-組合(2-combination,總共有
個)的平均相似度以作為決定最佳主題數量Kopt的指標。相似度例如是餘弦(cosine)相似度或雅卡爾(Jaccard)相似度,但本揭露不限於此。以表2為例,假設K等於「3」。處理器110可計算對應於主題#1的主題向量[0.20 0.05 … 0.20]與對應於主題#2的主題向量[0.00 0.10 … 0.15]之間的相似度、對應於主題#1的主題向量[0.20 0.05 … 0.20]與對應於主題#3的主題向量[0.10 0.20 … 0.05]之間的相似度以及對應於主題#2的主題向量[0.00 0.10 … 0.15]與對應於主題#3的主題向量[0.10 0.20 … 0.05]之間的相似度等三個相似度,並且計算三個相似度的平均以取得平均相似度,如表4所示。
表4
主題 | 相似度 |
(#1,#2) | 0.82 |
(#1,#3) | 0.58 |
(#2,#3) | 0.68 |
平均相似度: 0.69 |
圖2根據本發明的一實施例繪示主題之間的平均相似度與主題數量K的關係的示意圖。比較各個主題數量K所對應的平均相似度後可發現,當K值越大時,主題之間的平均相似度越小。因此,若以平均相似度作為決定最佳主題數量Kopt的指標,則處理器110可選擇較大的值以作為最佳主題數量Kopt。
為了使每一個不同的病歷(或病人)被歸類在具有代表性的用藥模式,具有偏重主題之病歷的比率越大越好。在一實施例中,處理器110可根據與特定主題相對應的病歷的個數與所有病歷的總數(即:N)計算比率以作為決定最佳主題數量Kopt的指標。具體來說,處理器110可根據病歷向量中的K個主題的機率分布判斷所述病歷向量是否偏重於特定的主題。若機率分布中的最大機率對應於一特定主題且所述最大機率大於機率閾值,則處理器110可判斷所述病歷向量(或病歷)偏重於所述特定主題。
表5
主題 | 病歷#1 | 病歷#2 | 病歷#3 | 病歷#4 | 病歷#5 |
#1 | 0.60 | 0.05 | 0.40 | 0.10 | 0.35 |
#2 | 0.30 | 0.70 | 0.25 | 0.20 | 0.30 |
#3 | 0.10 | 0.25 | 0.35 | 0.70 | 0.35 |
表5為多個病歷向量的範例,其中假設N等於「5」,K等於「3」,並且機率閾值等於「0.50」。以病歷#1為例,處理器110可響應於病歷#1中的最大機率「0.60」對應於主題#1且大於機率閾值「0.50」而判斷病歷#1偏重於主題#1。以病歷#3為例,處理器110可響應於病歷#3中的最大機率「0.40」小於或等於「0.50」而判斷病歷#3並未偏重於任何主題。依此類推,處理器110可根據表5的資料取得所有病歷的每一者所偏重的主題,如表6所示。
表6
病歷#1 | 病歷#2 | 病歷#3 | 病歷#4 | 病歷#5 | |
偏重主題 | 主題#1 | 主題#2 | 無 | 主題#3 | 無 |
在取得每一個病歷所偏重的主題後,處理器110可計算與特定主題相對應的病歷的個數與所有病歷的總數的比率以作為指標。以表6為例,與特定主題相對應的病歷(即:病歷#1、病歷#2和病歷#4)的個數等於「3」並且所有病歷的總數N等於「5」。處理器110可計算出比率「3/5」以作為用以決定最佳主題數量Kopt的指標。
圖3根據本發明的一實施例繪示具有偏重主題的病歷之比率與主題數量K的關係的示意圖。比較各個主題數量K所對應的比率後可發現,當K值越小時,具有偏重主題的病歷之比率越大。因此,若以具有偏重主題的病歷之比率作為決定最佳主題數量Kopt的指標,則處理器110可選擇較小的值以作為最佳主題數量Kopt。
在一實施例中,處理器110可以分群效能指標作為決定最佳主題數量Kopt的指標。首先,處理器110可將病歷分配至特定主題的群組。具體來說,處理器110可根據對應於病歷的病例向量中的K個主題的機率分布將所述病歷分配至K個群組的其中之一,其中K個群組分別對應於K個主題。處理器110例如可將病歷分配至對應於病歷例向量中的最大機率的主題之群組。以表5為例,處理器110可將病歷#1分配至對應於主題#1的群組,將病歷#2分配至對應於主題#2的群組,將病歷#3分配至對應於主題#1的群組,並且將病歷#4分配至對應於主題#3的群組。若病歷向量中存在複數個最大機率,則處理器110可根據預設規則或隨機地將病歷分配至對應於複數個最大機率的其中之一的主題。以表5為例,處理器110可根據預設規則或隨機地將病歷#5分配至對應於主題#1和主題#3的其中之一的群組。
在為N個病例分配好群組後,處理器110可根據K個群組計算對應於群組間距離的第一統計值以及對應於群組內距離的第二統計值。分群效能指標可等於第一統計值與第二統計值的比率。
第一統計值例如是K個主題向量之間的多個距離的總和。具體來說,第一統計值可為K個群組的所有2-組合的距離之總和。舉例來說,假設K等於「3」且K個群組包含群組#1、群組#2和群組#3。處理器110可計算群組#1與群組#2之間的距離、群組#1與群組#3之間的距離以及群組#2與群組#3之間的距離等三個距離,並且計算三個距離的總和以取得第一統計值。距離可根據主題向量之間的距離計算。舉例來說,對應於主題#1的群組#1和對應於主題#2的群組#2之間的距離可等於對應於主題#1的主題向量(例如:表2的[0.20 0.05 … 0.20])與對應於主題#2的主題向量(例如:表2的[0.00 0.10 … 0.15])等兩個主題向量之間的距離。群組之間的距離越大,代表分群的效能越好。因此,第一統計值可與分群效能指標呈正比。
第二統計值例如是分別對應於K個群組的K個群組內距離的總和。每一個群組的群組內距離可等於群組中的多個元素之間的多個距離的總和。更詳細來說,對應於群組的群組內距離可為所述群組中的所有元素的所有2-組合的距離的總和。以對應於主題#1的群組#1為例,假設群組#1包含病歷#1、病歷#2和病例#3等三個元素(即:N個病例中有三個病例與主題#1相對應),處理器110可計算病歷#1和病歷#2之間的距離、病歷#1和病歷#3之間的距離以及病歷#2和病歷#3之間的距離等三個距離,並且將三個距離相加以取得群組#1的群組內距離。
在一實施例中,處理器110可將病歷中記載的用藥組合向量化以計算病歷之間的距離。以表1的病歷#1和病歷#2為例,若處理器110欲計算病歷#1和病歷#2之間的距離,處理器110可將病歷#1的用藥組合「(A,B),(X,Y),(U,W)」轉換為一向量並可將病歷#2的用藥組合「(G,H)」轉換為另一向量。處理器110可計算兩個向量之間的距離以作為病歷#1和病歷#2之間的距離。群組內的元素之間的距離越小,代表分群的效能越好。因此,第二統計值可與分群效能指標呈反比。
圖4根據本發明的一實施例繪示分群效能指標與主題數量K的關係的示意圖。比較各個主題數量K所對應的比率後可發現,當K值越大時,分群效能指標越大。因此,若以分群效能指標作為決定最佳主題數量Kopt的指標,則處理器110可選擇較大的值以作為最佳主題數量Kopt。
處理器110可根據圖2、圖3和圖4來決定最佳主題數量Kopt。在決定好最佳主題數量Kopt後,處理器110可將最佳主題數量Kopt設為LDA模型參數中的主題數量K。接著,處理器110可根據主題數量K、LDA模型和N個病例產生分別對應於K個主題的K個主題向量(如表2的範例)和分別對應於N個病例的N個病例向量(如表3的範例)。
處理器110可執行篩選流程以為K個主題的每一者產生對應的獨有用藥組合集合。具體來說,處理器110可根據主題向量中所有用藥組合(即:M個用藥組合)的機率分布選出一或多個重要用藥組合以產生對應於所述主題向量的重要用藥組合集合。在一實施例中,處理器110可使用手肘法(elbow method)以自具有最大機率的用藥組合開始,從主題向量中選出多個重要用藥組合以產生重要用藥組合集合。
表7為重要用藥組合集合的範例。以主題#1為例,處理器110可根據對應於主題#1的主題向量將用藥組合依據機率的大小排列。具有越大機率的用藥組合排在越前面。接著,處理器110可利用手肘法找出經過排列之用藥組合的拐點(inflection point),並且選出排列在拐點之前的用藥組合以作為重要用藥組合。圖5根據本發明的一實施例繪示利用手肘法找出重要用藥組合的示意圖。處理器110可將特定主題(例如:主題#1)中的用藥組合依據機率的大小排列以畫出曲線50。若拐點為第四個用藥組合(即:具有第四大機率的用藥組合),處理器110可從M個用藥組合中選出排列在前四個的用藥組合(即:具有前四大機率的用藥組合)以產生重要用藥組合集合。主題#1的主題向量中,用藥組合(A,B)具有最大機率,用藥組合(C,D)具有次大機率,用藥組合(E,F)具有第三大機率,並且用藥組合(G,H)具有第四大機率。若用藥組合(G,H)對應於拐點,則處理器110可挑選排列在用藥組合(G,H)之前的用藥組合以產生重要用藥組合集合。
表7
主題 | 重要用藥組合集合 |
#1 | (A,B),(C,D),(E,F),(G,H) |
#2 | (G,H),(I,J),(X,Y),(S;W) |
… | … |
#K | (S;W),(O,P),(Q,R) |
在取得分別對應於K個主題的K個重要用藥組合集合後,處理器110可根據K個重要用藥組合集合產生分別對應於K個主題的K個獨有用藥組合集合。具體來說,處理器110可將重複出現於不同重要用藥組合集合的用藥組合刪除以產生獨有用藥組合集合。以表7為例,處理器110可響應於用藥組合(G,H)包含於與主題#1相對應的重要用藥組合集合且包含於與主題#2相對應的重要用藥組合而將用藥組合(G,H)自主題#1的重要用藥組合集合中刪除,並且將用藥組合(G,H)自主題#2的重要用藥組合集合中刪除,藉以產生獨有用藥組合集合,如表8所示。表8可為處理器110執行第一次篩選流程所產生的結果。
表8
主題 | 獨有用藥組合集合 |
#1 | (A,B),(C,D),(E,F), (G,H) |
#2 | (G,H),(I,J),(X,Y), (S;W) |
… | … |
#K | (S;W),(O,P),(Q,R) |
由於LDA演算法具有機率性,故前述的篩選流程中選出的重要且獨有的用藥組合可能僅是偶然出現在這次的篩選流程中。為了確保篩選出的用藥組合是穩定的,處理器110可重複執行多次篩選流程。具體來說,處理器110可執行多次篩選流程以產生多個獨有用藥組合集合。處理器110可響應於多個獨有用藥組合集合中的特定用藥組合的數量大於數量閾值而根據所述特定用藥組合產生穩定用藥組合集合,其中所述穩定用藥組合集合與所述特定用藥組合對應於相同的主題。
以對應於主題#1的獨有用藥組合集合中的用藥組合(A,B)為例,假設處理器110執行10次篩選流程且數量閾值為「6」。若用藥組合(A,B)出現的篩選流程如表9所示,則處理器110可響應於用藥組合(A,B)在10次篩選流程中出現的次數(即:7次)大於數量閾值而判斷用藥組合(A,B)是穩定的。據此,處理器110可根據用藥組合(A,B)產生對應於主題#1的穩定用藥組合集合,其中對應於主題#1的穩定用藥組合集合例如是表8中的用藥組合(A,B)、(C,D)和(E,F)。
表9
篩選流程 | #1 | #2 | #3 | #4 | #5 | #6 | #7 | #8 | #9 | #10 |
(A,B)存在? | 是 | 是 | 是 | 否 | 是 | 否 | 是 | 是 | 否 | 是 |
表10為各個主題的穩定用藥組合的範例。在產生分別對應於K個主題的K個穩定用藥組合集合後,處理器110可驗證各個穩定用藥組合集合是否符合足夠多人的用藥模式。具體來說,處理器110可根據病歷向量中K個主題的機率分布中的最大機率判斷病例所對應的主題。若病歷的病歷向量中的最大機率對應於特定主題,則處理器110可判斷所述病例對應於所述特定主題。若病歷向量中包含複數個最大機率,則處理器110可根據預設規則或隨機地判斷對應於複數個最大機率的其中之一的主題與病歷相對應。在完成判斷後,每個主題可對應於包含至少一病例的病歷集合。舉例來說,若N個病例中共有10個病例對應於主題#1,則代表主題#1所對應的病歷集合包含10個病例。
表10
主題 | 穩定用藥組合集合 |
#1 | (A,B),(C,D),(E,F) |
#2 | (I,J),(X,Y) |
… | … |
#K | (O,P),(Q,R) |
處理器110可計算病歷集合中的至少一病例佔病歷集合的比率,其中所述至少一病例指示穩定用藥組合集合中的至少一用藥組合。若比率大於比率閾值,代表符合穩定用藥組合集合所對應之主題(即:用藥模式)的病歷的樣本數足夠。據此,處理器110可根據穩定用藥組合集合產生最終的用藥組合集合。舉例來說,假設比率閾值為50%。若在N個病例中,有60%的病歷包含主題#1所對應的穩定用藥組合集合中的至少一用藥組合,代表符合主題#1之用藥模式的樣本(即:病歷)足夠。因此,處理器110可根據主題#1所對應的穩定用藥組合集合產生最終的用藥組合集合。表11為主題#1所對應的病歷集合的範例。相對來說,若在N個病例中,僅有40%的病歷包含主題#1所對應的穩定用藥組合集合中的至少一用藥組合,代表符合主題#1之用藥模式的樣本不足夠。因此,處理器110可不根據主題#1所對應的穩定用藥組合集合產生最終的用藥組合集合。
假設主題#1對應的病歷集合至少包含病歷#10、#11和#12(即:病歷#10、#11和#12的偏重主題為主題#1)。參考表10和表11,由於病歷#10記載的用藥組合(C,D)出現在主題#1的穩定用藥組合集合中,故處理器110可判斷病歷#10為上述至少一病例的其中之一。由於病歷#11記載的用藥組合(E,F)出現在主題#1的穩定用藥組合集合中,故處理器110可判斷病歷#11為上述至少一病例的其中之一。由於病歷#12記載的用藥組合(X,Y)並未出現在主題#1的穩定用藥組合集合中,故處理器110可判斷病歷#12非為上述至少一病例的其中之一。
表11
主題#1 | |||||
病歷 | … | #10 | #11 | #12 | ... |
用藥組合 | … | (C,D) | (E,F) | (X,Y) | … |
處理器110可根據上述的步驟檢查K個穩定用藥組合集合中的每一者是否符合足夠多人的用藥模式。若穩定用藥組合集合符合足夠多人的用藥模式,則處理器110可保留所述穩定用藥組合集合。若穩定用藥組合集合不符合足夠多人的用藥模式,則處理器110可刪除所述穩定用藥組合集合。據此,處理器110可從分別對應於K個主題的K個穩定用藥組合集合選出k個穩定用藥組合集合,其中k為小於或等於K的正整數。處理器110可取k個穩定用藥組合集合的聯集以取得最終的用藥組合集合,其中所述用藥組合集合可包含多個用藥組合。最終的用藥組合集合中的每一個用藥組合具有高重要性、高獨特性和高穩定性等特性,且符合大量病人的用藥模式。
在取得最終的用藥組合集合後,處理器110可為用藥組合集合中的用藥組合標記風險等級。具體來說,處理器110可根據N個病例為特定的用藥組合計算勝算比(odds ratio,OR),如方程式(1)和表12的混淆矩陣所示,其中
代表N個病例中曾經服用用藥組合且曾發生住院事件的病歷的數量,
代表N個病例中曾經服用用藥組合但未曾發生住院事件的病歷的數量,
代表N個病例中未曾服用用藥組合但曾發生住院事件的病歷的數量,並且
代表N個病例中未曾服用用藥組合且未曾發生住院事件的病歷的數量。
…(1)
表12
N個病例 | 曾發生住院事件 | 未曾發生住院事件 |
曾經服用該用藥組合 | ||
未曾服用該用藥組合 |
在計算完用藥組合集合中的每一個用藥組合的勝算比後,處理器110可根據用藥組合的勝算比為用藥組合標記風險等級。若用藥組合的勝算比大於風險閾值,代表所述用藥組合非常可能是導致住院事件發生的原因。據此,處理器110可將所述用藥組合標記為高風險。相對來說,若用藥組合的勝算比小於或等於風險閾值,代表所述用藥組合與住院事件發生較不相關。據此,處理器110可將所述用藥組合標記為低風險。表13為用藥組合的風險等級標記之範例。假設風險閾值為「1.3」,則處理器110可將勝算比大於「1.3」的用藥組合標記為高風險,並且將勝算比小於或等於「1.3」的用藥組合標記為低風險。
表13
用藥組合 | 勝算比 | 標記 |
(A,C) | 1.5 | 高風險 |
(A,D) | 0.9 | 低風險 |
(A,E) | 1.1 | 低風險 |
(A,F) | 2.7 | 高風險 |
(B,G) | 2.7 | 高風險 |
(B,H) | 0.8 | 低風險 |
(B,I) | 6.1 | 高風險 |
(B,J) | 3.0 | 高風險 |
在取得用藥組合集合中的每一個用藥組合的標記後,處理器110可根據標記的用藥組合的產生對應於所述用藥組合中的一藥物的風險組合分數(risk combination fraction,RCF)。以用藥組合(α,β)為例,為了確認用藥組合(α,β)中的藥物α與用藥組合集合中的其他藥物(即:除了藥物β以外的藥物)的組合是安全的,處理器110可根據方程式(2)計算對應於藥物α的風險組合分數RCF(或稱為「第三分數」),其中
為用藥組合集合中包含藥物α但不包含藥物β的用藥組合的數量,並且
為用藥組合集合中包含藥物α但不包含藥物β且被標記為高風險的用藥組合的數量。
…(2)
在取得藥物α的RCF後,處理器110可根據方程式(3)計算藥物α的正常組合分數(normal combination fraction,NCF)(或稱為「第一分數」或「第二分數」)。藥物α的正常組合分數越高,代表藥物α與除了藥物β的其他藥物之用藥組合的風險越低。藥物α的NCF可與用藥組合中包含藥物α但不包含藥物β的用藥組合的勝算比負相關。以表13為例,藥物A的NCF可與用藥組合(A,C)、(A,D)、(A,E)或(A,F)的勝算比負相關。
…(3)
以表13的藥物A為例,假設表13包含了用藥組合集合中除了用藥組合(A,B)以外的所有用藥組合,其中包含藥物A且被標記為高風險的用藥組合共有兩個,分別為用藥組合(A,C)和(A,F)。據此,處理器110可根據方程式(2)和(3)計算出藥物A的RCF等於「0.5」且藥物A的NCF等於「0.5」。以表13的藥物B為例,表13中包含藥物B且被標記為高風險的用藥組合共有三個,分別為用藥組合(B,G)、(B,I)和(B,J)。據此,處理器110可根據方程式(2)和(3)計算出藥物B的RCF等於「0.75」且藥物B的NCF等於「0.25」。
假設藥物α的NCF大於或等於藥物β的NCF,在取得藥物α的NCF和藥物β的NCF後,處理器110可計算藥物α的NCF和藥物β的NCF的商數(或比值)Q(α,β),如方程式(4)所示,其中NCF(α)為藥物α的NCF且NCF(β)為藥物β的NCF。Q(α,β)的數值大於或等於1,且該數值越低代表藥物α與除了藥物β的其他藥物之用藥組合的風險以及藥物β與除了藥物α的其他藥物之用藥組合的風險越相近。
…(4)
處理器110可根據以下三個條件判斷特定用藥組合是否具有高風險的相互作用。以用藥組合(α,β)為例,若用藥組合(α,β)的勝算比大於第一閾值,藥物α的NCF(α)與藥物β的NCF(β)的總和(或稱為NCF總和)大於第二閾值,並且藥物α的NCF(α)與藥物β的NCF(β)的商數(或稱為NCF商數)小於第三閾值,則處理器110可判斷用藥組合(α,β)具有高風險的交互作用。處理器110可通過收發器130輸出用藥組合(α,β)以供使用者參考。
表14為多個用藥組合之勝算比和NCF的相關參數。假設第一閾值為「2.0」,第二閾值為「1.2」,且第三閾值為「1.8」。由於用藥組合(E,F)的勝算比大於第一閾值,NCF總和大於第二閾值,且NCF商數小於第三閾值,故處理器110可判斷用藥組合(E,F)完全符合三個條件。據此,處理器110可輸出用藥組合(E,F)。由於用藥組合(G,H)的勝算比小於第一閾值,故處理器110可判斷用藥組合(G,H)並未完全符合三個條件。據此,處理器110可不輸出用藥組合(G,H)。由於用藥組合(A,B)的NCF總和小於第二閾值或NCF商數大於第三閾值,故處理器110可判斷用藥組合(A,B)並未完全符合三個條件。據此,處理器110可不輸出用藥組合(A,B)。由於用藥組合(C,D)的NCF商數大於第三閾值,故處理器110可判斷用藥組合(C,D)並未完全符合三個條件。據此,處理器110可不輸出用藥組合(C,D)。
表14
用藥組合 (α,β) | 勝算比 | NCF (α,其他) | NCF (β,其他) | NCF 總和 | NCF 商數 |
(E,F) | 2.38 | 0.695 | 0.652 | 1.347 | 1.065 |
(G,H) | 1.02 | 0.714 | 0.640 | 1.354 | 1.116 |
(A,B) | 8.53 | 0.500 | 0.250 | 0.750 | 2.000 |
(C,D) | 2.18 | 1.001 | 0.500 | 1.501 | 2.002 |
圖6根據本發明的一實施例繪示一種檢查藥物交互作用的方法的流程圖,其中所述方法可由如圖1所示的電子裝置100實施。在步驟S601中,取得多個病歷,其中多個病歷的至少其中之一指示服用第一用藥組合的病患是否發生住院事件。在步驟S602中,根據多個病歷產生用藥組合集合,其中用藥組合集合包含第一用藥組合、第二用藥組合以及第三用藥組合,其中第一用藥組合與第二用藥組合均包含第一藥物,並且第一用藥組合與第三用藥組合均包含第二藥物。在步驟S603中,根據多個病歷產生第一用藥組合以及住院事件之間的第一勝算比、第二用藥組合以及住院事件之間的第二勝算比以及第三用藥組合以及住院事件之間的第三勝算比。在步驟S604中,根據第二勝算比產生對應於第一藥物的第一分數,其中第一分數與第二勝算比負相關。在步驟S605中,根據第三勝算比產生對應於第二藥物的第二分數,其中第二分數與第三勝算比負相關,其中第一分數大於或等於第二分數。在步驟S606中,響應於第一勝算比大於第一閾值,第一分數與第二分數的總和大於第二閾值,並且第一分數與第二分數的商數小於第三閾值,輸出第一用藥組合。
綜上所述,本發明的電子裝置可通過隱含狄利克雷分布模型對多個病歷進行分析,基於用藥模式的相似度、病人偏重用藥模式或分群效能指標的因素選出用以對用藥組合群組的最佳數量。在根據最佳數量將各個用藥組合群組至特定的用藥模式後,電子裝置可根據用藥組合的重要性、獨特性、穩定性和樣本數量等因素篩選出最具有代表性的多個用藥組合。若特定用藥組合中的兩種安全藥物容易產生不良的交互作用,則輸出特定用藥組合的資訊供使用者參考。
100:電子裝置
110:處理器
120:儲存媒體
130:收發器
50:曲線
S601、S602、S603、S604、S605、S606:步驟
圖1根據本發明的一實施例繪示一種檢查藥物交互作用的電子裝置的示意圖。
圖2根據本發明的一實施例繪示主題之間的平均相似度與主題數量K的關係的示意圖。
圖3根據本發明的一實施例繪示具有偏重主題的病歷之比率與主題數量K的關係的示意圖。
圖4根據本發明的一實施例繪示分群效能指標與主題數量K的關係的示意圖。
圖5根據本發明的一實施例繪示利用手肘法找出重要用藥組合的示意圖。
圖6根據本發明的一實施例繪示一種檢查藥物交互作用的方法的流程圖。
S601、S602、S603、S604、S605、S606:步驟
Claims (14)
- 一種檢查藥物相互作用的方法,適用於電子裝置,包括:由所述電子裝置取得多個病歷,其中所述多個病歷的至少其中之一指示服用第一用藥組合的病患是否發生住院事件;由所述電子裝置根據所述多個病歷產生用藥組合集合,其中所述用藥組合集合包括所述第一用藥組合、第二用藥組合以及第三用藥組合,其中所述第一用藥組合與所述第二用藥組合均包括第一藥物,並且所述第一用藥組合與所述第三用藥組合均包括第二藥物;由所述電子裝置根據所述多個病歷產生所述第一用藥組合以及所述住院事件之間的第一勝算比、所述第二用藥組合以及所述住院事件之間的第二勝算比以及所述第三用藥組合以及所述住院事件之間的第三勝算比;由所述電子裝置根據所述第二勝算比產生對應於所述第一藥物的第一分數,其中所述第一分數與所述第二勝算比負相關;由所述電子裝置根據所述第三勝算比產生對應於所述第二藥物的第二分數,其中所述第二分數與所述第三勝算比負相關,其中所述第一分數大於或等於所述第二分數;以及由所述電子裝置響應於所述第一勝算比大於第一閾值,所述第一分數與所述第二分數的總和大於第二閾值,並且所述第一分數與所述第二分數的商數小於第三閾值,輸出所述第一用藥組合, 其中由所述電子裝置根據所述多個病歷產生所述用藥組合集合的步驟包括:由所述電子裝置執行篩選流程以產生第一獨有用藥組合集合,包括:由所述電子裝置根據所述多個病歷以及隱含狄利克雷分布模型產生包括第一主題向量的K個主題向量,其中K為第一主題數量,其中所述K個主題向量分別對應於K個主題,其中所述K個主題包括對應於所述第一主題向量的第一主題,其中所述第一主題向量包括所有用藥組合的機率分布;由所述電子裝置自具有最大機率的用藥組合開始,從所述第一主題向量中選出多個重要用藥組合以產生第一重要用藥組合集合;以及由所述電子裝置根據所述第一重要用藥組合集合決定所述第一獨有用藥組合集合;以及由所述電子裝置根據所述第一獨有用藥組合集合產生所述用藥組合集合。
- 如請求項1所述的方法,其中由所述電子裝置根據所述第二勝算比產生對應於所述第一藥物的所述第一分數的步驟包括:由所述電子裝置響應於所述第二勝算比大於風險閾值而標記所述第二用藥組合;由所述電子裝置根據標記的所述第二用藥組合產生第三分數, 其中所述第三分數等於所述用藥組合集合中包括所述第一藥物但不包含所述第二藥物且被標記的用藥組合的數量除以所述用藥組合集合中包括所述第一藥物但不包括所述第二藥物的用藥組合的數量;以及由所述電子裝置根據所述第三分數計算所述第一分數,其中所述第一分數與所述第三分數的總和等於一。
- 如請求項1所述的方法,其中所述K個主題包括第二主題,其中由所述電子裝置根據所述第一重要用藥組合集合決定所述第一獨有用藥組合集合的步驟包括:由所述電子裝置響應於第一重要用藥組合包含於與所述第一主題相對應的所述第一重要用藥組合集合以及與所述第二主題相對應的第二重要用藥組合集合,將所述第一重要用藥組合自所述第一重要用藥組合集合刪除以產生所述第一獨有用藥組合集合。
- 如請求項1所述的方法,其中由所述電子裝置根據所述第一獨有用藥組合集合產生所述用藥組合集合的步驟包括:由所述電子裝置重覆執行多次所述篩選流程以產生包括所述第一獨有用藥組合集合的多個獨有用藥組合集合;由所述電子裝置響應於所述多個獨有用藥組合集合中的所述第一用藥組合的數量大於數量閾值,根據所述第一用藥組合以產生對應於所述第一主題的第一穩定用藥組合集合;以及由所述電子裝置根據所述第一穩定用藥組合集合產生所述用藥組合集合。
- 如請求項4所述的方法,其中由所述電子裝置根據所述第一穩定用藥組合集合產生所述用藥組合集合的步驟包括:由所述電子裝置根據所述多個病歷以及所述隱含狄利克雷分布模型產生分別對應於所述多個病歷的多個病歷向量,其中所述多個病歷向量的每一者包括所述K個主題的機率分布;由所述電子裝置根據所述K個主題的所述機率分布判斷所述多個病歷中與所述第一主題相對應的病歷集合;由所述電子裝置計算所述病歷集合中的至少一病歷佔所述病歷集合的比率,其中所述至少一病歷指示所述第一穩定用藥組合集合中的至少一用藥組合;以及由所述電子裝置響應於所述比率大於比率閾值,根據所述第一穩定用藥組合集產生所述用藥組合集合,其中所述用藥組合集合包括所述第一穩定用藥組合集中的多個用藥組合。
- 如請求項5所述的方法,其中所述病歷集合中的第一病歷對應於所述K個主題的第一機率分布,其中由所述電子裝置根據所述K個主題的所述機率分布判斷所述多個病歷中與所述第一主題相對應的所述病歷集合的步驟包括:由所述電子裝置響應於所述第一機率分布中的最大機率對應於所述第一主題,判斷所述第一病歷對應於所述第一主題。
- 如請求項1所述的方法,更包括:由所述電子裝置根據所述多個病歷以及所述隱含狄利克雷分布模型產生對應於所述第一主題數量的第一指標以及對應於第二 主題數量的第二指標;以及由所述電子裝置比較所述第一指標與所述第二指標以從所述第一主題數量與所述第二主題數量中選出所述第一主題數量以作為K。
- 如請求項7所述的方法,其中由所述電子裝置產生對應於所述第一主題數量的所述第一指標的步驟包括:由所述電子裝置根據所述多個病歷、所述隱含狄利克雷分布模型以及所述第一主題數量產生所述K個主題向量;以及由所述電子裝置計算所述K個主題向量的所有2-組合的平均相似度以作為所述第一指標。
- 如請求項7所述的方法,其中由所述電子裝置產生對應於所述第一主題數量的所述第一指標的步驟包括:由所述電子裝置根據所述多個病歷、所述隱含狄利克雷分布模型以及所述第一主題數量產生分別對應於所述多個病歷的多個病歷向量,其中所述多個病歷向量的每一者包括所述K個主題的機率分布;由所述電子裝置根據所述K個主題的所述機率分布判斷所述多個病歷中與所述第一主題相對應的至少一病歷;以及由所述電子裝置根據所述至少一病歷的個數與所述多個病歷的總數計算比率以作為所述第一指標。
- 如請求項9所述的方法,其中由所述電子裝置根據所述K個主題的所述機率分布判斷所述多個病歷中與所述第一主題相對應的所述至少一病歷的步驟包括:由所述電子裝置自所述多個病歷向量取得對應於所述至少一病歷的所述K個主題的第一機率分布;以及由所述電子裝置響應於所述第一機率分布中的最大機率對應於所述第一主題且大於機率閾值,判斷所述至少一病歷對應於所述第一主題。
- 如請求項7所述的方法,其中由所述電子裝置產生對應於所述第一主題數量的所述第一指標的步驟包括:由所述電子裝置根據所述多個病歷、所述隱含狄利克雷分布模型以及所述第一主題數量產生分別對應於所述多個病歷的多個病歷向量,其中所述多個病歷向量的每一者包括所述K個主題的機率分布;由所述電子裝置根據所述K個主題的所述機率分布將所述多個病歷分為K個群組,其中所述K個群組分別對應於所述K個主題;由所述電子裝置根據所述K個群組計算群組間距離的第一統計值;由所述電子裝置根據所述K個群組計算群組內距離的第二統計值;以及由所述電子裝置計算所述第一統計值與所述第二統計值的比 率以作為所述第一指標。
- 如請求項11所述的方法,其中由所述電子裝置根據所述K個群組計算所述群組間距離的所述第一統計值的步驟包括:由所述電子裝置計算所述K個主題向量之間的多個距離;以及由所述電子裝置將所述多個距離相加以求得所述第一統計值。
- 如請求項11所述的方法,其中所述K個群組包括第一群組以及第二群組,其中由所述電子裝置根據所述K個群組計算所述群組內距離的所述第二統計值的步驟包括:由所述電子裝置計算所述第一群組中的多個元素之間的多個距離以產生對應於所述第一群組的第一群組內距離總和;以及由所述電子裝置將對應於所述第一群組的所述第一群組內距離總和與對應於所述第二群組的第二群組內距離總和相加以求得所述第二統計值。
- 一種檢查藥物相互作用的電子裝置,包括:收發器;以及處理器,耦接所述收發器並且經配置以執行:通過所述收發器取得多個病歷,其中所述多個病歷的至少其中之一指示服用第一用藥組合的病患是否發生住院事件;根據所述多個病歷產生用藥組合集合,其中所述用藥組合集 合包括所述第一用藥組合、第二用藥組合以及第三用藥組合,其中所述第一用藥組合與所述第二用藥組合均包括第一藥物,並且所述第一用藥組合與所述第三用藥組合均包括第二藥物;根據所述多個病歷產生所述第一用藥組合以及所述住院事件之間的第一勝算比、所述第二用藥組合以及所述住院事件之間的第二勝算比以及所述第三用藥組合以及所述住院事件之間的第三勝算比;根據所述第二勝算比產生對應於所述第一藥物的第一分數,其中所述第一分數與所述第二勝算比負相關;根據所述第三勝算比產生對應於所述第二藥物的第二分數,其中所述第二分數與所述第三勝算比負相關,其中所述第一分數大於或等於所述第二分數;以及響應於所述第一勝算比大於第一閾值,所述第一分數與所述第二分數的總和大於第二閾值,並且所述第一分數與所述第二分數的商數小於第三閾值,通過所述收發器輸出所述第一用藥組合,其中根據所述多個病歷產生所述用藥組合集合的步驟包括:執行篩選流程以產生第一獨有用藥組合集合,包括:根據所述多個病歷以及隱含狄利克雷分布模型產生包括第一主題向量的K個主題向量,其中K為第一主題數量,其中所述K個主題向量分別對應於K個主題,其中所述K個主題包括對應於所述第一主題向量的第一主題,其中所述 第一主題向量包括所有用藥組合的機率分布;自具有最大機率的用藥組合開始,從所述第一主題向量中選出多個重要用藥組合以產生第一重要用藥組合集合;以及根據所述第一重要用藥組合集合決定所述第一獨有用藥組合集合;以及根據所述第一獨有用藥組合集合產生所述用藥組合集合。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW111108870A TWI812056B (zh) | 2022-03-10 | 2022-03-10 | 檢查藥物相互作用的方法和電子裝置 |
CN202210665946.4A CN116779094A (zh) | 2022-03-10 | 2022-06-14 | 检查药物相互作用的方法和电子装置 |
US17/842,809 US20230290520A1 (en) | 2022-03-10 | 2022-06-17 | Method and electronic device of checking drug interaction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW111108870A TWI812056B (zh) | 2022-03-10 | 2022-03-10 | 檢查藥物相互作用的方法和電子裝置 |
Publications (2)
Publication Number | Publication Date |
---|---|
TWI812056B true TWI812056B (zh) | 2023-08-11 |
TW202336779A TW202336779A (zh) | 2023-09-16 |
Family
ID=87932276
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW111108870A TWI812056B (zh) | 2022-03-10 | 2022-03-10 | 檢查藥物相互作用的方法和電子裝置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20230290520A1 (zh) |
CN (1) | CN116779094A (zh) |
TW (1) | TWI812056B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104765947A (zh) * | 2015-03-02 | 2015-07-08 | 大连理工大学 | 一种面向大数据的潜在药物不良反应数据挖掘方法 |
TW201640389A (zh) * | 2015-02-17 | 2016-11-16 | 陳沛隆 | 藥物不良反應風險評估方法及其裝置 |
US20180166175A1 (en) * | 2015-05-22 | 2018-06-14 | Georgetown University | Discovery and analysis of drug-related side effects |
CN108831563A (zh) * | 2018-03-29 | 2018-11-16 | 南京邮电大学 | 一种用于判别药品不良反应信号分类检测的决策方法 |
CN110391028A (zh) * | 2010-08-13 | 2019-10-29 | 智能医学公司 | 用于生产个性化定制药物产品的系统和方法 |
CN105940114B (zh) * | 2013-08-19 | 2020-08-28 | 塞弗欧米公司 | 药物选择的计算机可读介质及系统 |
CN112309585A (zh) * | 2020-08-26 | 2021-02-02 | 国家药品监督管理局药品评价中心(国家药品不良反应监测中心) | 一种不良反应信号检测方法和装置 |
-
2022
- 2022-03-10 TW TW111108870A patent/TWI812056B/zh active
- 2022-06-14 CN CN202210665946.4A patent/CN116779094A/zh active Pending
- 2022-06-17 US US17/842,809 patent/US20230290520A1/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110391028A (zh) * | 2010-08-13 | 2019-10-29 | 智能医学公司 | 用于生产个性化定制药物产品的系统和方法 |
CN105940114B (zh) * | 2013-08-19 | 2020-08-28 | 塞弗欧米公司 | 药物选择的计算机可读介质及系统 |
TW201640389A (zh) * | 2015-02-17 | 2016-11-16 | 陳沛隆 | 藥物不良反應風險評估方法及其裝置 |
CN104765947A (zh) * | 2015-03-02 | 2015-07-08 | 大连理工大学 | 一种面向大数据的潜在药物不良反应数据挖掘方法 |
US20180166175A1 (en) * | 2015-05-22 | 2018-06-14 | Georgetown University | Discovery and analysis of drug-related side effects |
CN108831563A (zh) * | 2018-03-29 | 2018-11-16 | 南京邮电大学 | 一种用于判别药品不良反应信号分类检测的决策方法 |
CN112309585A (zh) * | 2020-08-26 | 2021-02-02 | 国家药品监督管理局药品评价中心(国家药品不良反应监测中心) | 一种不良反应信号检测方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN116779094A (zh) | 2023-09-19 |
TW202336779A (zh) | 2023-09-16 |
US20230290520A1 (en) | 2023-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9147167B2 (en) | Similarity analysis with tri-point data arbitration | |
JP6272555B2 (ja) | 計算機システム、及び、情報処理方法 | |
TW202029079A (zh) | 異常群體識別方法及裝置 | |
US11157697B2 (en) | Method and apparatus for extracting diagnosis object from medical document | |
JP2020537787A5 (zh) | ||
CN108346474A (zh) | 基于单词的类内分布与类间分布的电子病历特征选择方法 | |
Cheng et al. | Identification of homogeneous and heterogeneous variables in pooled cohort studies | |
TWI812056B (zh) | 檢查藥物相互作用的方法和電子裝置 | |
WO2022227171A1 (zh) | 关键信息提取方法、装置、电子设备及介质 | |
US20160283862A1 (en) | Multi-distance similarity analysis with tri-point arbitration | |
Wang et al. | Hybrid density-and partition-based clustering algorithm for data with mixed-type variables | |
CN111177356B (zh) | 一种酸碱指标医疗大数据分析方法及系统 | |
CN109767819B (zh) | 病历的分组方法、装置、存储介质和电子设备 | |
TWI790479B (zh) | 生理狀態評估方法與生理狀態評估裝置 | |
JP6705763B2 (ja) | 生成装置、生成方法および生成プログラム | |
CN107273430B (zh) | 一种数据存储方法和装置 | |
CN111667023B (zh) | 获取目标类别的文章的方法和装置 | |
Jung et al. | A machine learning method for selection of genetic variants to increase prediction accuracy of type 2 diabetes mellitus using sequencing data | |
CN110033862B (zh) | 一种基于加权有向图的中医量化诊断系统及存储介质 | |
Gada | Disease Prediction System using Machine Learning | |
CN110265151A (zh) | 一种基于ehr中异构时态数据的学习方法 | |
Rayan et al. | An Ensemble Model for Early Sepsis Prediction using Clinical Records from the Intensive Care Unit (ICU). | |
TWI776020B (zh) | 鎖定用戶興趣之裝置、方法及其電腦程式產品 | |
CN115312182B (zh) | 一种预测疫苗接种后不良反应转为严重不良反应风险的模型 | |
CN112906386B (zh) | 一种确定文本特征的方法及装置 |