Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*Unverified author*

R Software Module

rwasp_hierarchicalclustering.wasp

Title produced by software

Hierarchical Clustering

Date of computation

Sun, 04 Nov 2007 10:02:58 -0700

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2007/Nov/04/tx0mjk7v5g436zd1194195693.htm/, Retrieved Sun, 05 May 2024 17:11:19 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=364, Retrieved Sun, 05 May 2024 17:11:19 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

Estimated Impact

222

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)

-       [Hierarchical Clustering] [WS 5 Various EDA ...] [2007-11-04 17:02:58] [4bd8a0043457404de73994ae0e323922] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

7,4	6,4	8,7
7,2	6,3	8,5
7	6	8,2
6,6	5,3	8,3
6,4	5,2	8
6,4	5,1	8,1
6,8	5,4	8,7
7,3	5,8	9,3
7	5,6	8,9
7	5,6	8,8
6,7	5,4	8,4
6,7	5,4	8,4
6,3	5,5	7,3
6,2	5,5	7,2
6	5,3	7
6,3	5,7	7
6,2	5,6	6,9
6,1	5,5	6,9
6,2	5,6	7,1
6,6	5,9	7,5
6,6	6	7,4
7,8	7	8,9
7,4	6,6	8,3
7,4	6,6	8,3
7,5	6,3	9
7,4	6,3	8,9
7,4	6,3	8,8
7	6,3	7,8
6,9	6,2	7,8
6,9	6,2	7,8
7,6	6,3	9,2
7,7	6,4	9,3
7,6	6,4	9,2
8,2	7,8	8,6
8	7,7	8,5
8,1	7,7	8,5
8,3	7,7	9
8,2	7,7	9
8,1	7,6	8,8
7,7	7,5	8
7,6	7,4	7,9
7,7	7,4	8,1
8,2	7,5	9,3
8,4	7,6	9,4
8,4	7,6	9,4
8,6	8,1	9,3
8,4	7,8	9
8,5	8	9,1
8,7	7,9	9,7
8,7	7,9	9,7
8,6	7,8	9,6
7,4	6,7	8,3
7,3	6,6	8,2
7,4	6,6	8,4
9	7,7	10,6
9,2	7,9	10,9
9,2	8	10,9
8,5	7,7	9,6
8,3	7,5	9,3
8,3	7,6	9,3
8,6	7,8	9,6
8,6	7,8	9,5
8,5	7,7	9,5
8,1	7,4	9
8,1	7,5	8,9
8	7,2	9
8,6	7,5	10,1
8,7	7,6	10,2
8,7	7,6	10,2
8,6	7,8	9,5
8,4	7,7	9,3
8,4	7,7	9,3
8,7	8,2	9,4
8,7	8,2	9,3
8,5	8,1	9,1
8,3	7,8	9
8,3	7,8	8,9
8,3	7,7	9
8,1	6,7	9,8
8,2	6,7	10
8,1	6,7	9,8
8,1	7,2	9,4
7,9	6,9	9
7,7	6,8	8,9
8,1	7,2	9,3
8	7,1	9,1
7,7	6,9	8,8
7,8	6,9	8,9
7,6	6,7	8,7
7,4	6,5	8,6
7,7	6,6	9,1
7,9	6,6	9,3
7,6	6,5	8,9

Summary of compuational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

\begin{tabular}{lllllllll}
\hline
Summary of compuational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 2 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=364&T=0

[TABLE]
[ROW][C]Summary of compuational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]2 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=364&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=364&T=0

As an alternative you can also use a QR Code:

The GUIDs for individual cells are displayed in the table below:

Summary of compuational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

Summary of Dendrogram
Label	Height
1	0
2	0
3	0
4	0
5	0
6	0
7	0
8	0
9	0
10	0
11	0
12	0.0999999999999996
13	0.0999999999999996
14	0.0999999999999996
15	0.0999999999999996
16	0.0999999999999996
17	0.0999999999999996
18	0.0999999999999996
19	0.0999999999999996
20	0.0999999999999996
21	0.0999999999999996
22	0.0999999999999996
23	0.100000000000001
24	0.127614237491539
25	0.127614237491541
26	0.133333333333333
27	0.133333333333335
28	0.137377344785322
29	0.141421356237309
30	0.141421356237309
31	0.141421356237309
32	0.141421356237309
33	0.141421356237309
34	0.141421356237310
35	0.141421356237310
36	0.141421356237310
37	0.141421356237310
38	0.173205080756888
39	0.176417624662798
40	0.176417624662798
41	0.1764176246628
42	0.188561808316412
43	0.199999999999999
44	0.210018769324858
45	0.217400800258808
46	0.217719700768865
47	0.230940107675849
48	0.230940107675851
49	0.235702260395515
50	0.244948974278317
51	0.245436036690136
52	0.254165767818701
53	0.279458180291988
54	0.282842712474619
55	0.286168582658929
56	0.298142396999971
57	0.301763809020503
58	0.352859547920897
59	0.397606774342517
60	0.438685717649667
61	0.44684678446316
62	0.463087612001388
63	0.474034131450302
64	0.531864212827577
65	0.533693132821101
66	0.549237892498694
67	0.554234759029406
68	0.590009844047375
69	0.59481699012736
70	0.660920148281909
71	0.66111711741648
72	0.891825466289151
73	0.957332835878085
74	0.957506881651033
75	1.32373256223693
76	1.43916639772776
77	1.7635439904068
78	1.88870440693586
79	2.09246075267243
80	2.09391601590158
81	2.41833042900632
82	2.60806974096245
83	3.04547285938668
84	3.47791412803025
85	4.54834978713410
86	5.02997035987643
87	5.59184354295347
88	6.44758256224646
89	8.39491989814402
90	9.18402332816775
91	27.8558742683003
92	53.6444792355086

\begin{tabular}{lllllllll}
\hline
Summary of Dendrogram \tabularnewline
Label & Height \tabularnewline
1 & 0 \tabularnewline
2 & 0 \tabularnewline
3 & 0 \tabularnewline
4 & 0 \tabularnewline
5 & 0 \tabularnewline
6 & 0 \tabularnewline
7 & 0 \tabularnewline
8 & 0 \tabularnewline
9 & 0 \tabularnewline
10 & 0 \tabularnewline
11 & 0 \tabularnewline
12 & 0.0999999999999996 \tabularnewline
13 & 0.0999999999999996 \tabularnewline
14 & 0.0999999999999996 \tabularnewline
15 & 0.0999999999999996 \tabularnewline
16 & 0.0999999999999996 \tabularnewline
17 & 0.0999999999999996 \tabularnewline
18 & 0.0999999999999996 \tabularnewline
19 & 0.0999999999999996 \tabularnewline
20 & 0.0999999999999996 \tabularnewline
21 & 0.0999999999999996 \tabularnewline
22 & 0.0999999999999996 \tabularnewline
23 & 0.100000000000001 \tabularnewline
24 & 0.127614237491539 \tabularnewline
25 & 0.127614237491541 \tabularnewline
26 & 0.133333333333333 \tabularnewline
27 & 0.133333333333335 \tabularnewline
28 & 0.137377344785322 \tabularnewline
29 & 0.141421356237309 \tabularnewline
30 & 0.141421356237309 \tabularnewline
31 & 0.141421356237309 \tabularnewline
32 & 0.141421356237309 \tabularnewline
33 & 0.141421356237309 \tabularnewline
34 & 0.141421356237310 \tabularnewline
35 & 0.141421356237310 \tabularnewline
36 & 0.141421356237310 \tabularnewline
37 & 0.141421356237310 \tabularnewline
38 & 0.173205080756888 \tabularnewline
39 & 0.176417624662798 \tabularnewline
40 & 0.176417624662798 \tabularnewline
41 & 0.1764176246628 \tabularnewline
42 & 0.188561808316412 \tabularnewline
43 & 0.199999999999999 \tabularnewline
44 & 0.210018769324858 \tabularnewline
45 & 0.217400800258808 \tabularnewline
46 & 0.217719700768865 \tabularnewline
47 & 0.230940107675849 \tabularnewline
48 & 0.230940107675851 \tabularnewline
49 & 0.235702260395515 \tabularnewline
50 & 0.244948974278317 \tabularnewline
51 & 0.245436036690136 \tabularnewline
52 & 0.254165767818701 \tabularnewline
53 & 0.279458180291988 \tabularnewline
54 & 0.282842712474619 \tabularnewline
55 & 0.286168582658929 \tabularnewline
56 & 0.298142396999971 \tabularnewline
57 & 0.301763809020503 \tabularnewline
58 & 0.352859547920897 \tabularnewline
59 & 0.397606774342517 \tabularnewline
60 & 0.438685717649667 \tabularnewline
61 & 0.44684678446316 \tabularnewline
62 & 0.463087612001388 \tabularnewline
63 & 0.474034131450302 \tabularnewline
64 & 0.531864212827577 \tabularnewline
65 & 0.533693132821101 \tabularnewline
66 & 0.549237892498694 \tabularnewline
67 & 0.554234759029406 \tabularnewline
68 & 0.590009844047375 \tabularnewline
69 & 0.59481699012736 \tabularnewline
70 & 0.660920148281909 \tabularnewline
71 & 0.66111711741648 \tabularnewline
72 & 0.891825466289151 \tabularnewline
73 & 0.957332835878085 \tabularnewline
74 & 0.957506881651033 \tabularnewline
75 & 1.32373256223693 \tabularnewline
76 & 1.43916639772776 \tabularnewline
77 & 1.7635439904068 \tabularnewline
78 & 1.88870440693586 \tabularnewline
79 & 2.09246075267243 \tabularnewline
80 & 2.09391601590158 \tabularnewline
81 & 2.41833042900632 \tabularnewline
82 & 2.60806974096245 \tabularnewline
83 & 3.04547285938668 \tabularnewline
84 & 3.47791412803025 \tabularnewline
85 & 4.54834978713410 \tabularnewline
86 & 5.02997035987643 \tabularnewline
87 & 5.59184354295347 \tabularnewline
88 & 6.44758256224646 \tabularnewline
89 & 8.39491989814402 \tabularnewline
90 & 9.18402332816775 \tabularnewline
91 & 27.8558742683003 \tabularnewline
92 & 53.6444792355086 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=364&T=1

[TABLE]
[ROW][C]Summary of Dendrogram[/C][/ROW]
[ROW][C]Label[/C][C]Height[/C][/ROW]
[ROW][C]1[/C][C]0[/C][/ROW]
[ROW][C]2[/C][C]0[/C][/ROW]
[ROW][C]3[/C][C]0[/C][/ROW]
[ROW][C]4[/C][C]0[/C][/ROW]
[ROW][C]5[/C][C]0[/C][/ROW]
[ROW][C]6[/C][C]0[/C][/ROW]
[ROW][C]7[/C][C]0[/C][/ROW]
[ROW][C]8[/C][C]0[/C][/ROW]
[ROW][C]9[/C][C]0[/C][/ROW]
[ROW][C]10[/C][C]0[/C][/ROW]
[ROW][C]11[/C][C]0[/C][/ROW]
[ROW][C]12[/C][C]0.0999999999999996[/C][/ROW]
[ROW][C]13[/C][C]0.0999999999999996[/C][/ROW]
[ROW][C]14[/C][C]0.0999999999999996[/C][/ROW]
[ROW][C]15[/C][C]0.0999999999999996[/C][/ROW]
[ROW][C]16[/C][C]0.0999999999999996[/C][/ROW]
[ROW][C]17[/C][C]0.0999999999999996[/C][/ROW]
[ROW][C]18[/C][C]0.0999999999999996[/C][/ROW]
[ROW][C]19[/C][C]0.0999999999999996[/C][/ROW]
[ROW][C]20[/C][C]0.0999999999999996[/C][/ROW]
[ROW][C]21[/C][C]0.0999999999999996[/C][/ROW]
[ROW][C]22[/C][C]0.0999999999999996[/C][/ROW]
[ROW][C]23[/C][C]0.100000000000001[/C][/ROW]
[ROW][C]24[/C][C]0.127614237491539[/C][/ROW]
[ROW][C]25[/C][C]0.127614237491541[/C][/ROW]
[ROW][C]26[/C][C]0.133333333333333[/C][/ROW]
[ROW][C]27[/C][C]0.133333333333335[/C][/ROW]
[ROW][C]28[/C][C]0.137377344785322[/C][/ROW]
[ROW][C]29[/C][C]0.141421356237309[/C][/ROW]
[ROW][C]30[/C][C]0.141421356237309[/C][/ROW]
[ROW][C]31[/C][C]0.141421356237309[/C][/ROW]
[ROW][C]32[/C][C]0.141421356237309[/C][/ROW]
[ROW][C]33[/C][C]0.141421356237309[/C][/ROW]
[ROW][C]34[/C][C]0.141421356237310[/C][/ROW]
[ROW][C]35[/C][C]0.141421356237310[/C][/ROW]
[ROW][C]36[/C][C]0.141421356237310[/C][/ROW]
[ROW][C]37[/C][C]0.141421356237310[/C][/ROW]
[ROW][C]38[/C][C]0.173205080756888[/C][/ROW]
[ROW][C]39[/C][C]0.176417624662798[/C][/ROW]
[ROW][C]40[/C][C]0.176417624662798[/C][/ROW]
[ROW][C]41[/C][C]0.1764176246628[/C][/ROW]
[ROW][C]42[/C][C]0.188561808316412[/C][/ROW]
[ROW][C]43[/C][C]0.199999999999999[/C][/ROW]
[ROW][C]44[/C][C]0.210018769324858[/C][/ROW]
[ROW][C]45[/C][C]0.217400800258808[/C][/ROW]
[ROW][C]46[/C][C]0.217719700768865[/C][/ROW]
[ROW][C]47[/C][C]0.230940107675849[/C][/ROW]
[ROW][C]48[/C][C]0.230940107675851[/C][/ROW]
[ROW][C]49[/C][C]0.235702260395515[/C][/ROW]
[ROW][C]50[/C][C]0.244948974278317[/C][/ROW]
[ROW][C]51[/C][C]0.245436036690136[/C][/ROW]
[ROW][C]52[/C][C]0.254165767818701[/C][/ROW]
[ROW][C]53[/C][C]0.279458180291988[/C][/ROW]
[ROW][C]54[/C][C]0.282842712474619[/C][/ROW]
[ROW][C]55[/C][C]0.286168582658929[/C][/ROW]
[ROW][C]56[/C][C]0.298142396999971[/C][/ROW]
[ROW][C]57[/C][C]0.301763809020503[/C][/ROW]
[ROW][C]58[/C][C]0.352859547920897[/C][/ROW]
[ROW][C]59[/C][C]0.397606774342517[/C][/ROW]
[ROW][C]60[/C][C]0.438685717649667[/C][/ROW]
[ROW][C]61[/C][C]0.44684678446316[/C][/ROW]
[ROW][C]62[/C][C]0.463087612001388[/C][/ROW]
[ROW][C]63[/C][C]0.474034131450302[/C][/ROW]
[ROW][C]64[/C][C]0.531864212827577[/C][/ROW]
[ROW][C]65[/C][C]0.533693132821101[/C][/ROW]
[ROW][C]66[/C][C]0.549237892498694[/C][/ROW]
[ROW][C]67[/C][C]0.554234759029406[/C][/ROW]
[ROW][C]68[/C][C]0.590009844047375[/C][/ROW]
[ROW][C]69[/C][C]0.59481699012736[/C][/ROW]
[ROW][C]70[/C][C]0.660920148281909[/C][/ROW]
[ROW][C]71[/C][C]0.66111711741648[/C][/ROW]
[ROW][C]72[/C][C]0.891825466289151[/C][/ROW]
[ROW][C]73[/C][C]0.957332835878085[/C][/ROW]
[ROW][C]74[/C][C]0.957506881651033[/C][/ROW]
[ROW][C]75[/C][C]1.32373256223693[/C][/ROW]
[ROW][C]76[/C][C]1.43916639772776[/C][/ROW]
[ROW][C]77[/C][C]1.7635439904068[/C][/ROW]
[ROW][C]78[/C][C]1.88870440693586[/C][/ROW]
[ROW][C]79[/C][C]2.09246075267243[/C][/ROW]
[ROW][C]80[/C][C]2.09391601590158[/C][/ROW]
[ROW][C]81[/C][C]2.41833042900632[/C][/ROW]
[ROW][C]82[/C][C]2.60806974096245[/C][/ROW]
[ROW][C]83[/C][C]3.04547285938668[/C][/ROW]
[ROW][C]84[/C][C]3.47791412803025[/C][/ROW]
[ROW][C]85[/C][C]4.54834978713410[/C][/ROW]
[ROW][C]86[/C][C]5.02997035987643[/C][/ROW]
[ROW][C]87[/C][C]5.59184354295347[/C][/ROW]
[ROW][C]88[/C][C]6.44758256224646[/C][/ROW]
[ROW][C]89[/C][C]8.39491989814402[/C][/ROW]
[ROW][C]90[/C][C]9.18402332816775[/C][/ROW]
[ROW][C]91[/C][C]27.8558742683003[/C][/ROW]
[ROW][C]92[/C][C]53.6444792355086[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=364&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=364&T=1

As an alternative you can also use a QR Code:

The GUIDs for individual cells are displayed in the table below:

Summary of Dendrogram
Label	Height
1	0
2	0
3	0
4	0
5	0
6	0
7	0
8	0
9	0
10	0
11	0
12	0.0999999999999996
13	0.0999999999999996
14	0.0999999999999996
15	0.0999999999999996
16	0.0999999999999996
17	0.0999999999999996
18	0.0999999999999996
19	0.0999999999999996
20	0.0999999999999996
21	0.0999999999999996
22	0.0999999999999996
23	0.100000000000001
24	0.127614237491539
25	0.127614237491541
26	0.133333333333333
27	0.133333333333335
28	0.137377344785322
29	0.141421356237309
30	0.141421356237309
31	0.141421356237309
32	0.141421356237309
33	0.141421356237309
34	0.141421356237310
35	0.141421356237310
36	0.141421356237310
37	0.141421356237310
38	0.173205080756888
39	0.176417624662798
40	0.176417624662798
41	0.1764176246628
42	0.188561808316412
43	0.199999999999999
44	0.210018769324858
45	0.217400800258808
46	0.217719700768865
47	0.230940107675849
48	0.230940107675851
49	0.235702260395515
50	0.244948974278317
51	0.245436036690136
52	0.254165767818701
53	0.279458180291988
54	0.282842712474619
55	0.286168582658929
56	0.298142396999971
57	0.301763809020503
58	0.352859547920897
59	0.397606774342517
60	0.438685717649667
61	0.44684678446316
62	0.463087612001388
63	0.474034131450302
64	0.531864212827577
65	0.533693132821101
66	0.549237892498694
67	0.554234759029406
68	0.590009844047375
69	0.59481699012736
70	0.660920148281909
71	0.66111711741648
72	0.891825466289151
73	0.957332835878085
74	0.957506881651033
75	1.32373256223693
76	1.43916639772776
77	1.7635439904068
78	1.88870440693586
79	2.09246075267243
80	2.09391601590158
81	2.41833042900632
82	2.60806974096245
83	3.04547285938668
84	3.47791412803025
85	4.54834978713410
86	5.02997035987643
87	5.59184354295347
88	6.44758256224646
89	8.39491989814402
90	9.18402332816775
91	27.8558742683003
92	53.6444792355086

Figure 1

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = ward ; par2 = ALL ; par3 = FALSE ; par4 = FALSE ;

Parameters (R input):

par1 = ward ; par2 = ALL ; par3 = FALSE ; par4 = FALSE ;

R code (references can be found in the software module):

par3 <- as.logical(par3)
par4 <- as.logical(par4)
if (par3 == 'TRUE'){
dum = xlab
xlab = ylab
ylab = dum
}
x <- t(y)
hc <- hclust(dist(x),method=par1)
d <- as.dendrogram(hc)
str(d)
mysub <- paste('Method: ',par1)
bitmap(file='test1.png')
if (par4 == 'TRUE'){
plot(d,main=main,ylab=ylab,xlab=xlab,horiz=par3, nodePar=list(pch = c(1,NA), cex=0.8, lab.cex = 0.8),type='t',center=T, sub=mysub)
} else {
plot(d,main=main,ylab=ylab,xlab=xlab,horiz=par3, nodePar=list(pch = c(1,NA), cex=0.8, lab.cex = 0.8), sub=mysub)
}
dev.off()
if (par2 != 'ALL'){
if (par3 == 'TRUE'){
ylab = 'cluster'
} else {
xlab = 'cluster'
}
par2 <- as.numeric(par2)
memb <- cutree(hc, k = par2)
cent <- NULL
for(k in 1:par2){
cent <- rbind(cent, colMeans(x[memb == k, , drop = FALSE]))
}
hc1 <- hclust(dist(cent),method=par1, members = table(memb))
de <- as.dendrogram(hc1)
bitmap(file='test2.png')
if (par4 == 'TRUE'){
plot(de,main=main,ylab=ylab,xlab=xlab,horiz=par3, nodePar=list(pch = c(1,NA), cex=0.8, lab.cex = 0.8),type='t',center=T, sub=mysub)
} else {
plot(de,main=main,ylab=ylab,xlab=xlab,horiz=par3, nodePar=list(pch = c(1,NA), cex=0.8, lab.cex = 0.8), sub=mysub)
}
dev.off()
str(de)
}
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Summary of Dendrogram',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Label',header=TRUE)
a<-table.element(a,'Height',header=TRUE)
a<-table.row.end(a)
num <- length(x[,1])-1
for (i in 1:num)
{
a<-table.row.start(a)
a<-table.element(a,hc$labels[i])
a<-table.element(a,hc$height[i])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.tab')
if (par2 != 'ALL'){
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Summary of Cut Dendrogram',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Label',header=TRUE)
a<-table.element(a,'Height',header=TRUE)
a<-table.row.end(a)
num <- par2-1
for (i in 1:num)
{
a<-table.row.start(a)
a<-table.element(a,i)
a<-table.element(a,hc1$height[i])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')
}

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code