Mathematical Study of SARS

... maybe solving equations is experience you need to gain understanding - but until you do understand, you're just solving equations. Richard P. Feynman

[HOME]

[Probability Distributions]

[Quantum Computation]


[Cryptography]

Severe Acute Respiratory Syndrome (SARS)

Typical pneumonia is usually caused by bacteria such as Streptococcus pneumoniae. Atypical pneumonia is usually caused by influenza virus, mycoplasma, chlamydia, adenovirus and other unknown micro-organisms. Common symptoms include fever, chills, cough, headache and general malaise.

Atypical pneumonia can be transmitted by respiratory droplets over a short distance of one metre, or by contact with a patient's respiratory secretions.

For more information, Hong Kong Government's Atypical pneumonia Web site

The outbreak of SARS in Block E, Amoy Gardens, Ngau Tau Kok was first reported on 26 March. I have done a simple mathematical study of this incident based on the following assumptions.

Cumulative Number of SARS Admission

Hong Kong 18 Districts Intra-spreading Dynamics

Calculated Probability

Symbol Assumption Value
Number of flat, N Number of floor = 26, Number of flat per floor = 10 260
Incubation period, d days 2 - 10 2
Number of lift, L 1 - 4 1
Fraction of coincident time, T Average time taken to and from the lift and flat = 5 minutes.
Consider one takes maximum 90 minutes from home to the workplace, the overlapping time at peak = 90 mminutes.
Each day the morning and evening is the peak period of lift usage.
10 / 180
= 1/18
Probability of infection by droplets, P As calculated by the formula 0.094
Probability of inflection by touching droplets in lift, t It is (1/No of floor) of touching the virus pad. After touching it, there are four cases that may cause infection: touching nose, eye, mouth, and his family member. 1 /26 * 1/4
= 1/104
Number of flat infected, I According to the newspaper, there were four infected families. Not counting the one who had member visited Wales Hospital, there were three infected. 3

Formula

For incubation 2, 3, 4, 5 days and no. of flat per floor 8-12, below is the calculated probability of infection by droplets.

Discussion - 30 March, 2003

Cumulative Number of SARS Admission

Logarithmic plot of the cumulative number of SARS admission from March 14 to May 22

May 23, 2003

On May 23, WHO removed Hong Kong and Guangdong from SARS travel advisory list.

Discussion

May 1, 2003

Ministry of Health of mainland China announces the number of confirmed and number of suspected of SARS daily. Department of Health of Hong Kong has followed this way of doing on April 25.

Discussion

April 24, 2003

Discussion

April 18, 2003

According to the government released data the total number of deaths up to April 18 is 69 and the total number of discharged patients is 322. The cumulative number of SARS patients is 1327. Department of Health calculates the death rate of the disease by 69 over 1327 equal to 0.052. In terms of percentage it is about 5%.

Discussion

April 11, 2003

On April 9, the Hong Kong Hospital Authority (HA) estimated that there would be about 1,800 to 3,000 SARS patients in Hong Kong at the end of April.

Discussion

April 3, 2003

Today local newspapers report government saying that in Wales Hospital during the SARS outbreak statistically every 25 persons contacted the patient had one infected atypical pneumonia, the probability is 4%.

For P = 0.04, (number of lift L, incubation period d)
      (L,d) = (1,4), (2,8), (3,12), (4,16)

The corresponding incubation periods matches with the Health Department estimation, i.e. 2-12 days.

Discussion

Epidemic Model

The linearity of the logarithmic plot of above suggests that

where k, C are constants, t is time, and I is the number of infected.

Consider, the rate of spreading is proportional to the number of infected and also the number of uninfected. Each infected person will contact fraction of people c in the incubation period. The differential equation, therefore, can be written as

where a is a constant, N is the number of people who are susceptible to infection, c is the fraction of the people that the patient contact in the incubation period. In the Amoy Gardens case, N is equal to the total number of residents of Block E.

The solution of this autonomous first order differential equation is

therefore, k = acN, is the slope of the logarithmic plot and C is the intercept of the t-axis.

Discussion

Hong Kong 18 Districts Intra-spreading Dynamics

Department of Health produces a list of the names of the buildings with SARS patients who have been hospitalized for treatment everyday. The list has been compiled from information provided by patients with confirmed diagnosis of SARS. If there are no further cases in 10 days, the names of the buildings will be taken off the list.

I have done a brief study of the spread of SARS in HK by using this list. The study is based on the idea of data clustering and dispersion of incidents. At the beginning, the data set was arranged in a normal distribution form, with the center at Kwun Tong and New Territories on the right and Kowloon and HK on the left. Date-sequential plotting of these data sets reveals both geographic and eventual patterns. The results may be used in the tracking and strategic planning of the control of spreading of the disease.

Data Preparation

Initial sorting and arrangement of the 18 Districts data (April 15)

Southern Central & Western Wan Chai Yau Tsim Mong Wong Tai Sin Sai Kung Kowloon City Sham Shui Po Eastern Kwun Tong Tai Po Shatin Kwai Tsing Tuen Mun Tsuen Wan North Yuen Long Islands
0 0 1 9 9 10 11 12 13 29 23 22 16 9 5 4 2 1

The value was peak at Kwun Tong. On the left are districts from Hong Kong & Kowloon except Sai Kung. On the right are districts from New Territories. Sai Kung is put on the HK & Kowloon side because there is a main road from Kowloon City to Wong Tai Sin to Sai Kung. Later we shall see this arrangement gives a meaningful result.

Note: I did not expect to see the distribution of Hong Kong, Kowloon, NT districts at the beginning of the data preparation. I just thinking of Amoy Gardens in Kwun Tong be the main cause of spreading should be placed at the central.

Results & Plotting (April 15 - May 22)

Discussion (May 2 - May 22) Discussion (April 25 - May 1) Discussion (May 2 - May 22) Discussion (April 15 - 24)

HongKong - China Inter-spreading Dynamics

Consider SARS carriers or patients arrived from China is a Poisson process. ...


Census-Taker Problem

A census-taker knocks on a door, and asks the woman inside how many children she has and how old they are.
Mother said, "I have three daughers, their ages are whole numbers, and the product of the ages is 36."
"That's not enough information," responded the census-taker.
"I'd tell you the sum of their ages, but you's still be stumped."
"I wish you'd tell me something more."
"Okay, my oldest daughter Annie likes dogs."

What are the ages of the three daughters?

Managing Risk