This program is by Jean-Marc Baillard and is used here by permission.
This program is supplied without representation or warranty of any kind. Jean-Marc Baillard and The Museum of HP Calculators therefore assume no responsibility and shall have no liability, consequential or otherwise, of any kind arising from the use of this program material or any part thereof.
-This program calculates the median, quartiles, or any other
quantile, by linear interpolation.
-The data are supposed to be grouped like this:
xi | frequencies |
[ x1 ; x2 [ | n1 |
[ x2 ; x3 [ | n2 |
.......... | ....... |
[ xp; xp+1] | np |
Data Registers: All the following registers are to be initialized before executing "QTLE"
• R00: p • R01:
x1 • R02:
x2 .......... • Rpp: xp
• Rpp+1: xp+1
• Rpp+2: n1 • Rpp+3:
n2 ........ • R2p+2:
np
Flags: none
Subroutine: none
01 LBL "QTLE"
02 STO M
03 RCL 00
04 ENTER^
05 ST+ X
06 1
07 +
08 0
09 LBL 00
10 RCL IND Y
11 +
12 DSE Y
13 DSE Z
14 GTO 00
15 RCL M
16 *
17 0
18 LBL 01
19 ISG Z
20 CL X
21 RCL IND Z
22 +
23 X<Y?
24 GTO 01
25 ST- Y
26 X<> L
27 /
28 X<>Y
29 RCL 00
30 -
31 ENTER^
32 DSE X
33 X<> Z
34 RCL IND Y
35 RCL IND T
36 -
37 *
38 RCL IND Y
39 +
40 0
41 X<> M
42 SIGN
43 CLA
44 END
( 71 bytes / SIZE 2p+2 )
STACK | INPUTS | OUTPUTS |
Y | / | 1 |
X | r ( 0 < r < 1 ) | quantile |
L | / | r |
Example: Find the median Me , the first and last
quartile Q1 and Q3 , the first decile d1
and the third centile c3 of the following data:
x | n |
[ 10 ; 12 [ | 1 |
[ 12 ; 15 [ | 3 |
[ 15 ; 20 [ | 7 |
[ 20 ; 24 [ | 4 |
[ 24 ; 30 [ | 2 |
-There are 5 classes:
5 STO 00
then we store the 11 numbers 10 ; 12 ; 15 ; 20 ; 24 ;
30 ; 1 ; 3 ; 7 ; 4 ; 2 into registers R01 thru R11 respectively and
.5 XEQ "QTLE"
yields Me = Q2 = 18.2143
.25 R/S
-----
Q1 = 15.1786
.75 R/S
-----
Q3 = 21.75
.1
R/S
-----
d1 = 12.7
.03 R/S
-----
c3 = 11.02
Go back to the HP-41 software library
Go back to the general software library
Go
back to the main exhibit hall