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
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월달... 이후의 실제의 국내 추세를 계속 관찰 주시를 해봐야 할 것으로 사료됩니다.
또한 많은 돌발변수들이 많은 관계로 어느 누구도 정확한 예측을 하기는 힘들 것으로 생각됩니다.
한편 바이러스가 퍼지면 감염경로를 추적하여 방어체계를 구축하더라도 일반적으로 바이러스는 무생물 처럼 비리온 상태로 존재하다가
그 지역에 전파된 이후 지속적으로 병을 재발 시키게 될 가능성도 있으며...
특히 변종 바이러스가 나타나는 경우에 기존 약의 효과는 장담 할 수 없기 때문에...
더더욱 감염확산방지체계를 초기에 좀더 철저하게 관리 해야만 한다고 보입니다.
((( 치료약 )))
현재 완벽한 약은 아니나...
미국의 Z맵, 일본의 아비간, 캐나다 TKM-에볼라 3종 및 라미부딘(제픽스: 에이즈 치료제, B형 간염치료제) 1종 등이 거론 되고 있는 것으로 보입니다.
' 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
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 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
암튼....
현재 저 또한 제일 걱정이 되는 것은 맨처음 계산된 이론적인 전파속도 보다 훨씬 빠르게 만일 우리나라 도심에 급속도로 감염확산시...
현재 어느정도 효과가 있는 약을 외국 선진국으로 부터 그리 쉽게 물량을 충분하게 구하기가 쉽지 않을 것으로 보아...
또한 감염증상 발현 후 적절한 치료약으로 처방하지 못하면 통상 7일~9일 안에 사망한다는 사실과
현재 인터넷 정보에 의하면 불안정한 약으로 보이나...
미국의 Z맵, 일본의 아비간, 캐나다 TKM-에볼라 3종 및
라이베리아 현지의 의사 고비 로건이 에볼라 환자에게 투여하여 효과를 크게 보았다는 라미부딘 1종이 더 있는 것으로 확인중이니
그나마 당행이기는 하지만 얼마나 국내에서 에볼라환자 발생시 적절한 약의 보급이 이루어져 에볼라 감염확산을 막을 수 있을지는
그 지역의 감염확산 정도에 따라서 많은 차이가 있을 것으로 생각이 됩니다.
암튼...암튼... 저의 단순한 계산 결과와 같이 불행한 결과가 되지 않기를 간절히 바래볼 뿐 입니다.
라미부딘(제픽스: 에이즈 치료제, B형 간염치료제)로써 내성이 크다지만 에볼라 치료에 효과가 좋았다고 한다면...
혹시 라미부딘 유사약인 내성이 낮은 바라크루드와 세비보...등도 에볼라 치료에 효과가 있을지 매우 궁금함...
왜냐하면 국제적으로 인정된 약의 보급이 늦어져서 바로 죽어야할 판이라면 라미부딘 에이즈 치료제나, B형 간염약이라도 먹어야 할지 모르잖아요...
암튼...암튼...
한편... 국제적으로 실제 감염수치 결과값을 단순 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
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
2차함수로 계산 결과는...
Y = A + B*X + C*X^2
A = 733.0668
B = -1162.761
C = 198.8368
STD.DEV. = 987.4241
Prod.Corr.Coef. = .781975
2차 방정식에 따른 결과값을 얻기 위해서 다음과 같이 프로그램을 돌려보면...
10 A = 733.0668: B = -1162.761: C = 198.8368
11 OPEN "infect99.txt" FOR OUTPUT AS #1
20 FOR I = 1 TO 33
30 E = A + B * I + C * I * I
40 PRINT I; " Day", "Ebolavirus Infection Number = "; USING "##############"; E
41 PRINT #1, I, USING "##############"; E
50 NEXT I
51 CLOSE #1
11 OPEN "infect99.txt" FOR OUTPUT AS #1
20 FOR I = 1 TO 33
30 E = A + B * I + C * I * I
40 PRINT I; " Day", "Ebolavirus Infection Number = "; USING "##############"; E
41 PRINT #1, I, USING "##############"; E
50 NEXT I
51 CLOSE #1
다음과 같은 2차 방정식 추정 결과를 얻습니다.
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
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
3차함수로 계산 결과는...
Y = A + B*X + C*X^2 + D*X^3
A = -41.41331
B = 861.1438
C = -354.3622
D = 37.11634
STD.DEV. = 319.3323
Prod.Corr.Coef. = .781975
3차 방정식에 따른 결과값을 얻기 위해서 다음과 같이 프로그램을 돌려보면...
10 A = -41.41331: B = 861.1438: C = -354.3622: D = 37.11634
11 OPEN "infect3d.txt" FOR OUTPUT AS #1
20 FOR I = 1 TO 33
30 E = A + B * I + C * I * I + D * I * I * I
40 PRINT I; " Day", "Ebolavirus Infection Number = "; USING "##############"; E
41 PRINT #1, I, USING "##############"; E
50 NEXT I
51 CLOSE #1
11 OPEN "infect3d.txt" FOR OUTPUT AS #1
20 FOR I = 1 TO 33
30 E = A + B * I + C * I * I + D * I * I * I
40 PRINT I; " Day", "Ebolavirus Infection Number = "; USING "##############"; E
41 PRINT #1, I, USING "##############"; E
50 NEXT I
51 CLOSE #1
다음과 같은 3차 방정식 추정 결과를 얻습니다.
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 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차 방정식 예상 추정값이 좀더 정확할 수 있을 것으로 보입니다.
이론값과 실제값의 편차가 상당히 크게 나타나고 있으니 나름 참고해서 보시면 좋겠습니다.
글쓰기