Download Statistical Analysis and Histogramming Reference

7.7. Fitting with finite Monte Carlo statistics
13
72.841
0.98350
1
2
3
4
5
6
7
8
9
10
11
12
155
40.320
-32.156
29.461
-20.018
7.7912
-18.482
-7.0469
-9.0691
4.4675
6.4382
-3.9874
3.1189
0
20
60
80
100
50
110
10
130
90
140
30
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
REGRESSOR STANDARD DEVIATION
CONFIDENCE INTERVAL
1
0.70650
[ 39.144
, 41.495
]
2
1.0509
[ -33.904
, -30.407
]
3
0.85147
[ 28.044
, 30.877
]
4
0.90853
[ -21.530
, -18.506
]
5
0.82336
[ 6.4213
, 9.1611
]
6
1.0259
[ -20.189
, -16.775
]
7
0.78401
[ -8.3514
, -5.7425
]
8
1.0731
[ -10.854
, -7.2836
]
9
0.68036
[ 3.3355
, 5.5995
]
10
0.83660
[ 5.0463
, 7.8301
]
11
0.85716
[ -5.4135
, -2.5612
]
12
1.1372
[ 1.2267
, 5.0110
]
DOUBLE PRECISION FUNCTION FPARAM (X)
DOUBLE PRECISION COEFF,P,P0,P1,P2,HELEFT,HBASFT
DIMENSION X(1),COEFF(12),IBASFT( 1,12)
DATA COEFF/ 0.40319615E+02,-0.32155589E+02, 0.29460772E+02,
+-0.20017895E+02, 0.77912196E+01,-0.18481896E+02,
+-0.70469122E+01,-0.90690550E+01, 0.44674803E+01,
+ 0.64381900E+01,-0.39873663E+01, 0.31188760E+01
+/
DATA IBASFT/ 0, 20, 60, 80,100, 50,110, 10,130, 90,140, 30
+/
FPARAM=0.
DO 25 K=1,12
P=1.
DO 15 I=1, 1
NUM=IBASFT(I,K)/10
ITYP=IBASFT(I,K)-NUM*10
IF (NUM.NE.0) THEN
IF (ITYP.EQ.0) THEN
P0=1.
P1=2*X (I)-1.
DO 10 J=2,NUM
P2=2*(2*X (I)-1.)*P1-P0
P0=P1
10 P1=P2
P=P*P1
END IF
IF (ITYP.EQ.1) P=P*HELEFT(NUM,X (I))
IF (ITYP.EQ.2) THEN
P=HBASFT(NUM,X )
GOTO 20
END IF
END IF
15 CONTINUE
20 FPARAM=FPARAM+COEFF(K)*P
25 CONTINUE
RETURN
END