|
DataMuseum.dkPresents historical artifacts from the history of: MIKADOS |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about MIKADOS Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - metrics - download
Length: 5056 (0x13c0)
Notes: Mikados TextFile, Mikados_K
Names: »KORRELAT«
└─⟦4cd80c18f⟧ Bits:30003895/4xCOMAL3a.imd 4 * COMAL - SPC/1 - Disketter til bogen
└─⟦this⟧ »KORRELAT«
└─⟦bcf7b34c7⟧ Bits:30003895/4xCOMAL2.imd 4 * COMAL - SPC/1 - Disketter til bogen
└─⟦this⟧ »KORRELAT«
0010 // -------------------------- korrelation ---------------------- 0020 // Programmet finder en række deskriptorer for n talpar: 0030 // (x1,y1),(x2,y2), ... ,(xn,yn). 0040 // n angives i den første datasætning. 0050 // x1,x2, ... ,xn angives i de næste datasætninger 0060 // y1,y2, ... ,yn angives i de næste datasætninger 0070 // 0080 CLEAR 0090 READ N 0100 DIM X(N),Y(N),S$ OF 1 0110 SUM1:=0;SUM2:=0;SUM3:=0;SUM4:=0;SUM5:=0 0120 FOR I:=1 TO N DO 0130 READ X(I) 0140 NEXT I 0150 FOR I:=1 TO N DO 0160 READ Y(I) 0170 NEXT I 0180 FOR I:=1 TO N DO 0190 SUM1:=SUM1+X(I) 0200 SUM2:=SUM2+Y(I) 0210 SUM3:=SUM3+X(I)**2 0220 SUM4:=SUM4+Y(I)**2 0230 SUM5:=SUM5+X(I)*Y(I) 0240 NEXT I 0250 MX:=SUM1/N 0260 MY:=SUM2/N 0270 VX:=SUM3/N-MX*MX 0280 VY:=SUM4/N-MY*MY 0290 COVXY:=SUM5/N-MX*MY 0300 SX:=SQR(VX) 0310 SY:=SQR(VY) 0320 R:=COVXY/SX/SY 0330 A:=COVXY/VX 0340 B:=MY-A*MX 0350 PRINT USING "Observationssætet består af### talpar.":N 0360 PRINT 0370 PRINT "Observationssættet X:" 0380 FOR I:=1 TO N DO 0390 PRINT USING "+####, ":X(I), 0400 NEXT I 0410 PRINT 0420 PRINT 0430 PRINT "Observationssættet Y:" 0440 FOR I:=1 TO N DO 0450 PRINT USING "+####, ":Y(I), 0460 NEXT I 0470 PRINT 0480 PRINT 0490 PRINT "Observationssættene X og Y har følgende deskriptorer: " 0500 PRINT 0510 PRINT USING "M(X) =+####.###":MX 0520 PRINT USING "M(Y) =+####.###":MY 0530 PRINT USING "V(X) =+####.###":VX 0540 PRINT USING "V(Y) =+####.###":VY 0550 PRINT USING "S(X) =+####.###":SX 0560 PRINT USING "S(Y) =+####.###":SY 0570 PRINT USING "COV(X,Y) =+####.###":COVXY 0580 PRINT USING "R(X,Y) =+####.###":R 0590 PRINT 0600 PRINT "Y-regressionslinien y:=ax + b er fastlagt ved:" 0610 PRINT USING "a=+###.### b=+###.### ":A,B 0630 PRINT 0640 DATA 8 0650 DATA 5,7,8,10,8,11,9,7 0660 DATA 12,8,7,6,8,4,7,9 0670 PRINT USING "Liniens ligning: y = +###.###x +###.###":A,B 0690 CURSOR 70,24 0700 INPUT "Return!":S$