x값 : 확산 개월수    ( 0 값 기준은 2013년 12월 시점 )

                                                                                                                                y값 : 확산 감염자수 (명)

2013/12  ... 1
2014/01  ... ?
2014/02  ... ?
2014/03/22 ... 49
2014/04/21 ... 242
2014/05/23 ... 270
2014/06/22 ... 599
2014/07/23 ... 1201
2014/08/20 ... 2615
2014/09/23 ... 6574
2014/10/23 ... 10141

... ?

... ?

... ?

 

중에서 9개의 데이터만 사용하였으며 10개월째 실제 감염수치와 이론상의 감염수치를 같게 고려한 감염전파속도로 볼 때에

지난번 월 단위로 2배 증가하는 결과 보다 추세 분석 결과 좀더 가파른 2.5154배 정도가 되고 있음을 확인하였음.

본 결과는 아프리카에서 집중적으로 전파되는 정도와 앞으로의 아프리카에서 세계로 전파되는 정도를 모두 함께 나타내고 있으므로...

실제 우리나라 서울 도심에 에볼라바이러스가 확산된다면 이보다 월씬 빠른 결과를 나타낼 것 으로 추정을 하고 있습니다.

왜냐하면 감염확산이 집중된 아프리카보다도 우리나라 서울과 수원과 같은 도시의 경우 주거밀집도가 높으며 교통수단과 화폐수단이

아프리카보다 발달이 되어있는 관계로 훨씬 감염확산이 빠르리라 추정되고 있습니다.

저의 짧은 소견인지는 모르나...

혹시 우리나라에 에볼라 바이러스가 들어온다면 아래의 월별 감염전파 속도 보다 훨씬 커질 것으로 예상을 하고 있습니다만...

좀더 적확한 결과 및 11월달 12월달... 이후의 실제의 국내 추세를 계속 관찰 주시를 해봐야 할 것으로 사료됩니다.

또한 많은 돌발변수들이 많은 관계로 어느 누구도 정확한 예측을 하기는 힘들 것으로 생각됩니다.

한편 바이러스가 퍼지면 감염경로를 추적하여 방어체계를 구축하더라도 일반적으로 바이러스는 무생물 처럼 비리온 상태로 존재하다가

그 지역에 전파된 이후 지속적으로 병을 재발 시키게 될 가능성도 있으며...

특히 변종 바이러스가 나타나는 경우에 기존 약의 효과는 장담 할 수 없기 때문에...

더더욱 감염확산방지체계를 초기에 좀더 철저하게 관리 해야만 한다고 보입니다.

 

' QBASIC SOURCE FILE

10 E = 1
11 OPEN "Infect.Txt" FOR OUTPUT AS #1
20      FOR i = 1 TO 33
30          E = E * 2.5154
40          PRINT i; " Day", "Ebolavirus Infection Number = "; USING "##############"; E
41          PRINT #1, i, USING "##############"; E
50      NEXT i
51 CLOSE #1

 

RUN 결과

월                        예상 감염자수

 1                         3
 2                         6
 3                        16
 4                        40
 5                       101
 6                       253
 7                       637
 8                      1603
 9                      4031
 10                    10141
 11                    25508
 12                    64163
 13                   161396
 14                   405974
 15                  1021188
 16                  2568696                                                서울인구수  =  약 1000만 ?
 17                  6461298
 18                 16252747
 19                 40882160
 20                102834984
 21                258671104
 22                650661248
 23               1636673280
 24               4116887808                                              세계인구수 ? 육박
 25              10355618816                                             
 26              26048522240
 27              65522450432
 28             164815159296
 29             414576050176
 30            1042824560640
 31            2623120867328
 32            6598198034432
 33           16597107343360

 

 

암튼...암튼...

한편... 국제적으로 실제 감염수치 결과값을 단순 2차 함수와 3차 함수에 근접시켜 본 결과는 모두 표준편차가 심하지만 각각 다음과 같은 결과를 나타내고 있습니다.

' 사용된 프로그램 소스파일

10 CLS : KEY OFF
15 LOCATE 5, 17: PRINT "*********************************************"
20 LOCATE 6, 17: PRINT "*                                           *"
25 LOCATE 7, 17: PRINT "*        Least Squares fitting Program      *"
30 LOCATE 8, 17: PRINT "*                                           *"
35 LOCATE 9, 17: PRINT "*              ( Version 1.00 )             *"
40 LOCATE 10, 17: PRINT "*                                           *"
45 LOCATE 11, 17: PRINT "*    (C) Copyright 1991  by  $$$$ $$$ $$    *"
50 LOCATE 12, 17: PRINT "*                                           *"
55 LOCATE 13, 17: PRINT "*********************************************"
60 LOCATE 22, 17: PRINT "Press any Key ....."
65 A$ = INKEY$: IF A$ = "" THEN 65
100 'SELECT MENU
110 CLS
120 COLOR 0, 7: LOCATE 3, 26: PRINT " SELECTION ": COLOR 7, 0
130 LOCATE 6, 26: PRINT "1. FIRST DEGREE EQUATION "
140 LOCATE 8, 26: PRINT "2. SECOND DEGREE EQUATION "
150 LOCATE 10, 26: PRINT "3. THIRD DEGREE EQUATION "
160 LOCATE 12, 26: PRINT "4. QUIT "
170 LOCATE 22, 26: INPUT "SELECT NUMBER ? ", NO$
180 ON VAL(NO$) GOTO 210, 810, 1480, 2260
190 GOTO 100
200 END
210 ' first degree
220 CLS : PRINT
230 PRINT "<<<  least squares fitting for first degree equation  >>>"
240 PRINT
250 PRINT "Y  =  A  +  B * X "
260 PRINT STRING$(57, "_")
270 PRINT
280 INPUT "data number  "; N
290 PRINT
300 DIM X(N), Y(N)
310 FOR K = 1 TO N
320     PRINT "point ("; K; ","; K; ") ";
330     INPUT X(K), Y(K)
340 NEXT K
350 FOR J = 1 TO N
360     U1 = U1 + X(J)
370     U2 = U2 + X(J) * X(J)
380     V1 = V1 + Y(J)
390     U1V1 = U1V1 + X(J) * Y(J)
400 NEXT J
410 FOR Z = 1 TO N
420     DX = DX + (X(Z) - (U1 / N)) ^ 2
430     DY = DY + (Y(Z) - (V1 / N)) ^ 2
440     DXY = DXY + (X(Z) - (U1 / N)) * (Y(Z) - (V1 / N))
450 NEXT Z
460 R = DXY / SQR(DX * DY)
470 B = (U1 * V1 - N * U1V1) / (U1 * U1 - N * U2)
480 A = (V1 - B * U1) / N
490 CLS : PRINT
500 COLOR 0, 7: PRINT " <<<  result  >>> ": COLOR 7, 0
510 PRINT
520 PRINT "  Y  =  A  +  B * X "
530 PRINT
540 PRINT "  A  = "; A
550 PRINT "  B  = "; B
560 PRINT STRING$(70, "-")
570 PRINT "  X "; TAB(15); " Y ";
580 PRINT "measured"; TAB(30);
590 PRINT " Y calculated"
600 PRINT
610 ICHK = 88
620 FOR I = 1 TO N
630     Y = A + B * X(I)
640     DEV = DEV + (Y(I) - Y) * (Y(I) - Y)
650     PRINT " "; X(I);
660     PRINT TAB(15); Y(I);
670     PRINT TAB(30); Y
680     ICHK = I MOD 5
690     IF ICHK = 0 AND I < N THEN PRINT STRING$(70, "-"): GOSUB 2280
700 NEXT I
710 PRINT STRING$(70, "-")
720 DEV = DEV / (N - 1)
730 SD = SQR(DEV)
740 PRINT "  STD.DEV.  = "; SD
750 PRINT "  Prod.Corr.Coef.  = "; R ' (?)
760 PRINT : PRINT : PRINT "< END > : Press any Key .....";
770 A$ = INKEY$: IF A$ = "" THEN 770
780 U1 = 0: U2 = 0: V1 = 0: U1V1 = 0: DX = 0: DY = 0: DXY = 0: DEV = 0
790 ERASE X, Y
800 GOTO 100
810 ' second degree
820 CLS : PRINT
830 PRINT "<<<  least squares fitting for second degree equation  >>>"
840 PRINT
850 PRINT "Y  =  A  +  B * X  +  C * X ^ 2 "
860 PRINT STRING$(58, "_")
870 PRINT
880 INPUT "data number  "; N
890 PRINT
900 DIM X(N), Y(N)
910 FOR K = 1 TO N
920     PRINT "point ("; K; ","; K; ") ";
930     INPUT X(K), Y(K)
940 NEXT K
950 FOR J = 1 TO N
960     U1 = U1 + X(J)
970     U2 = U2 + X(J) * X(J)
980     U3 = U3 + X(J) * X(J) * X(J)
990    U4 = U4 + X(J) * X(J) * X(J) * X(J)
1000    V1 = V1 + Y(J)
1010    U1V1 = U1V1 + X(J) * Y(J)
1020    U2V1 = U2V1 + X(J) * X(J) * Y(J)
1030 NEXT J
1040 FOR Z = 1 TO N
1050     DX = DX + (X(Z) - (U1 / N)) ^ 2
1060     DY = DY + (Y(Z) - (V1 / N)) ^ 2
1070     DXY = DXY + (X(Z) - (U1 / N)) * (Y(Z) - (V1 / N))
1080 NEXT Z
1090 R = DXY / SQR(DX * DY)
1100 L1 = ((U2 * V1 - N * U2V1) * (U1 * U1 - N * U2) - (U1 * V1 - N * U1V1) * (U1 * U2 - N * U3))
1110 M1 = ((U2 * U2 - N * U4) * (U1 * U1 - N * U2) - (U2 * U1 - N * U3) * (U1 * U2 - N * U3))
1120 C = L1 / M1
1130 B = ((U1 * V1 - N * U1V1) - C * (U1 * U2 - N * U3)) / (U1 * U1 - N * U2)
1140 A = (V1 - (B * U1 + C * U2)) / N
1150 CLS : PRINT
1160 COLOR 0, 7: PRINT " <<<  result  >>> ": COLOR 7, 0
1170 PRINT
1180 PRINT "  Y  =  A  +  B * X  +  C * X ^ 2 "
1190 PRINT
1200 PRINT "  A  = "; A
1210 PRINT "  B  = "; B
1220 PRINT "  C  = "; C
1230 PRINT STRING$(70, "-")
1240 PRINT "  X "; TAB(15); " Y ";
1250 PRINT "measured"; TAB(30);
1260 PRINT " Y calculated"
1270 PRINT
1280 ICHK = 88
1290 FOR I = 1 TO N
1300     Y = A + B * X(I) + C * X(I) ^ 2
1310     DEV = DEV + (Y(I) - Y) * (Y(I) - Y)
1320     PRINT " "; X(I);
1330     PRINT TAB(15); Y(I);
1340     PRINT TAB(30); Y
1350     ICHK = I MOD 5
1360     IF ICHK = 0 AND I < N THEN PRINT STRING$(70, "-"): GOSUB 2370
1370 NEXT I
1380 PRINT STRING$(70, "-")
1390 DEV = DEV / (N - 1)
1400 SD = SQR(DEV)
1410 PRINT "  STD.DEV.  = "; SD
1420 PRINT "  Prod.Corr.Coef.  = "; R ' (?)
1430 PRINT : PRINT : PRINT "< END > : Press any Key ..... ";
1440 A$ = INKEY$: IF A$ = "" THEN 1440
1450 U1 = 0: U2 = 0: U3 = 0: U4 = 0: V1 = 0: U1V1 = 0: U2V1 = 0: DX = 0: DY = 0: DXY = 0: DEV = 0
1460 ERASE X, Y
1470 GOTO 100
1480 ' third degree
1490 CLS : PRINT
1500 PRINT "<<<  least squares fitting for third degree equation  >>>"
1510 PRINT
1520 PRINT "Y  =  A  +  B * X  +  C * X ^ 2  +  D * X ^ 3 "
1530 PRINT STRING$(57, "_")
1540 PRINT
1550 INPUT "data number  "; N
1560 PRINT
1570 DIM X(N), Y(N)
1580 FOR K = 1 TO N
1590     PRINT "point ("; K; ","; K; ") ";
1600     INPUT X(K), Y(K)
1610 NEXT K
1620 FOR J = 1 TO N
1630     U1 = U1 + X(J)
1640     U2 = U2 + X(J) * X(J)
1650     U3 = U3 + X(J) * X(J) * X(J)
1660     U4 = U4 + X(J) * X(J) * X(J) * X(J)
1670     U5 = U5 + X(J) * X(J) * X(J) * X(J) * X(J)
1680     U6 = U6 + X(J) * X(J) * X(J) * X(J) * X(J) * X(J)
1690     V1 = V1 + Y(J)
1700     U1V1 = U1V1 + X(J) * Y(J)
1710     U2V1 = U2V1 + X(J) * X(J) * Y(J)
1720     U3V1 = U3V1 + X(J) * X(J) * X(J) * Y(J)
1730 NEXT J
1740 L1 = ((U1 * V1 - N * U1V1) * (U1 * U2 - N * U3) - (U2 * V1 - N * U2V1) * (U1 * U1 - N * U2))
1750 M1 = ((U1 * U2 - N * U3) * (U1 * U2 - N * U3) - (U2 * U2 - N * U4) * (U1 * U1 - N * U2))
1760 N1 = ((U3 * U1 - N * U4) * (U1 * U2 - N * U3) - (U2 * U3 - N * U5) * (U1 * U1 - N * U2))
1770 L2 = ((U1 * V1 - N * U1V1) * (U3 * U1 - N * U4) - (U3 * V1 - N * U3V1) * (U1 * U1 - N * U2))
1780 M2 = ((U1 * U2 - N * U3) * (U3 * U1 - N * U4) - (U2 * U3 - N * U5) * (U1 * U1 - N * U2))
1790 N2 = ((U3 * U1 - N * U4) * (U3 * U1 - N * U4) - (U3 * U3 - N * U6) * (U1 * U1 - N * U2))
1800 FOR Z = 1 TO N
1810     DX = DX + (X(Z) - (U1 / N)) ^ 2
1820     DY = DY + (Y(Z) - (V1 / N)) ^ 2
1830     DXY = DXY + (X(Z) - (U1 / N)) * (Y(Z) - (V1 / N))
1840 NEXT Z
1850 R = DXY / SQR(DX * DY)
1860 D = (L2 * M1 - L1 * M2) / (M1 * N2 - M2 * N1)
1870 C = (L1 * N2 - L2 * N1) / (M1 * N2 - M2 * N1)
1880 B = ((U1 * V1 - N * U1V1) - C * (U1 * U2 - N * U3) - D * (U3 * U1 - N * U4)) / (U1 * U1 - N * U2)
1890 A = (V1 - B * U1 - C * U2 - D * U3) / N
1900 CLS : PRINT
1910 COLOR 0, 7: PRINT " <<<  result  >>> ": COLOR 7, 0
1920 PRINT
1930 PRINT "  Y  =  A  +  B * X  +  C * X ^ 2  + D * X ^ 3 "
1940 PRINT
1950 PRINT "  A  = "; A
1960 PRINT "  B  = "; B
1970 PRINT "  C  = "; C
1980 PRINT "  D  = "; D
1990 PRINT STRING$(70, "-")
2000 PRINT "  X "; TAB(15); " Y ";
2010 PRINT "measured"; TAB(30);
2020 PRINT " Y calculated"
2030 PRINT
2040 ICHK = 88
2050 FOR I = 1 TO N
2060     Y = A + B * X(I) + C * X(I) * X(I) + D * X(I) * X(I) * X(I)
2070     DEV = DEV + (Y(I) - Y) * (Y(I) - Y)
2080     PRINT " "; X(I);
2090     PRINT TAB(15); Y(I);
2100     PRINT TAB(30); Y
2110     ICHK = I MOD 5
2120     IF ICHK = 0 AND I < N THEN PRINT STRING$(70, "-"): GOSUB 2460
2130 NEXT I
2140 PRINT STRING$(70, "-")
2150 DEV = DEV / (N - 1)
2160 SD = SQR(DEV)
2170 PRINT "  STD.DEV.  = "; SD
2180 PRINT "  Prod.Corr.Coef.  = "; R ' (?)
2190 PRINT : PRINT : PRINT "< END > : Press any Key ..... ";
2200 A$ = INKEY$: IF A$ = "" THEN 2200
2210 U1 = 0: U2 = 0: U3 = 0: U4 = 0: U5 = 0: U6 = 0
2220 V1 = 0: U1V1 = 0: U2V1 = 0: U3V1 = 0: DX = 0: DY = 0: DXY = 0: DEV = 0
2230 ERASE X, Y
2240 GOTO 100
2250 BEEP: BEEP: GOTO 100
2260 ' quit
2270 CLS : KEY ON: SYSTEM
2280 ' first - space
2290 PRINT : PRINT "< NEXT > : Press any Key ..... ";
2300 A$ = INKEY$: IF A$ = "" THEN 2300
2310 IF ICHK = N THEN RETURN
2320 FOR CL = 11 TO 19
2330     LOCATE CL, 1: PRINT STRING$(70, " ")
2340 NEXT CL
2350 LOCATE 11, 1
2360 RETURN
2370 ' second - space
2380 PRINT : PRINT "< NEXT > : Press any Key ..... ";
2390 A$ = INKEY$: IF A$ = "" THEN 2390
2400 IF ICHK = N THEN RETURN
2410 FOR CL = 12 TO 20
2420     LOCATE CL, 1: PRINT STRING$(70, " ")
2430 NEXT CL
2440 LOCATE 12, 1
2450 RETURN
2460 ' third - space
2470 PRINT : PRINT "< NEXT > : Press any Key ..... ";
2480 A$ = INKEY$: IF A$ = "" THEN 2480
2490 IF ICHK = N THEN RETURN
2500 FOR CL = 13 TO 21
2510     LOCATE CL, 1: PRINT STRING$(70, " ")
2520 NEXT CL
2530 LOCATE 13, 1
2540 RETURN

 

 1                      -231
 2                      -797
 3                      -966
 4                      -737
 5                      -110
 6                       915
 7                      2337
 8                      4157
 9                      6374
 10                     8989
 11                    12002
 12                    15412
 13                    19221
 14                    23426
 15                    28030
 16                    33031
 17                    38430
 18                    44226
 19                    50421
 20                    57013
 21                    64002
 22                    71389
 23                    79174
 24                    87357
 25                    95937
 26                   104915
 27                   114291
 28                   124064
 29                   134235
 30                   144803
 31                   155770
 32                   167134
 33                   178895

 

 

 1                       502
 2                       560
 3                       355
 4                       109
 5                        45
 6                       386
 7                      1354
 8                      3172
 9                      6063
 10                    10250
 11                    15955
 12                    23401
 13                    32811
 14                    44407
 15                    58412
 16                    75049
 17                    94540
 18                   117108
 19                   142977
 20                   172367
 21                   205503
 22                   242607
 23                   283902
 24                   329610
 25                   379954
 26                   435156
 27                   495440
 28                   561029
 29                   632144
 30                   709008
 31                   791845
 32                   880877
 33                   976326

 

이 결과는 아프리카 현지에서 감염되는 수준으로만 국한지어 생각을 하기 보다는

              앞으로 있을 아프리카 현지에서 세계로 퍼져나가는 정도를 나타낸다는 의미로 보시는 것이 좋을 듯 사료 됩니다.

              2차 방정식 예상 추정값 보다는 3차 방정식 예상 추정값이 좀더 정확할 수 있을 것으로 보입니다.