# Need help in coding- Monte Carlo simulation for assessing cure against virus

11 views (last 30 days)
SethC on 30 Mar 2020
Answered: Hiro Yoshino on 30 Mar 2020
Hi everyone,
Am working on a thesis and need to seek help in coding (for Monte Carlo). Basically, have had some Matlab experience before but very basic...
Problem:
1. Each viral strain contains 7 - 11 different Proteins (A, B, C... K)
2. In a controlled space, the space &time dimensions are X and T
3. Proteins occur in variable amounts of quantities q (qA, qB... qK)
4. There are 3 types of viral strains:
• Type A = {A B C D E F G J} (contains J) - probability of occuring = 10%
• Type B = {A B C D E F G K} (contains K) - probability of occuring = 5%
• Type C = {A B C D E F G H} (common strain) - probability of occuring = 85%
5. Three chains of vaccines were developed to test within the same controlled space
• Type B1 = {A B C D E F G J} (can handle common strain and J)
• Type B2 = {A B C D E F G K} (Can handle common strain and K)
• Type B3 = {A B C D E F G H} (to handle common strain)
• Each Type's occurence is a function of space and time, X and T
• Each type starts off with quantity Q (QA, QB... QK) and depletes as it attends to a viral strain:
• After 1 week (or 168 hours), the quantity of vaccine will be boosted back to original levels
Example:
Type B1 = f(X ,T)
QB1 (X,T0) = (QA QB ..... QK)
QB1 (X,T1) = (QA-qA,1 QB-qB,1 ....)
QB1 (X, 169) = (QA QB .... QK) [Replenished to full]
6. Need to find out how many of these viral strains in the controlled sample can be destroyed by just these Types of vaccines (B1 / B2 / B3).
• Success case happens when quantity of vaccine > quantity of virus (Example: AND(Qa > qa, Qb > qb, Qc>qc...)
Would like to seek help to code/ simulate the above such that it is possible to know how many of the viral strains can be successfully destroyed in a one year period (365 X 24 hours)
thanks