Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*Unverified author*

R Software Module

rwasp_histogram.wasp

Title produced by software

Histogram

Date of computation

Tue, 14 Feb 2012 07:54:51 -0500

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/2012/Feb/14/t1329224139cggwyaoyt8ltazy.htm/, Retrieved Mon, 29 Apr 2024 14:42:40 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=162280, Retrieved Mon, 29 Apr 2024 14:42:40 +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

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

-     [Histogram] [] [2012-02-09 12:52:45] [af6b347fe17b8a57e076538e87e25d4e]
- R P     [Histogram] [Prijs zalm per kilo] [2012-02-14 12:54:51] [d41d8cd98f00b204e9800998ecf8427e] [Current]
-    D      [Histogram] [] [2012-02-14 14:30:29] [af6b347fe17b8a57e076538e87e25d4e] 
- RM D      [Kernel Density Estimation] [] [2012-02-14 14:37:45] [af6b347fe17b8a57e076538e87e25d4e] 

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	'Gertrude Mary Cox' @ cox.wessa.net

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 1 seconds \tabularnewline
R Server & 'Gertrude Mary Cox' @ cox.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=162280&T=0

[TABLE]
[ROW][C]Summary of computational 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]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gertrude Mary Cox' @ cox.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=162280&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=162280&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 computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	'Gertrude Mary Cox' @ cox.wessa.net

Frequency Table (Histogram)
Bins	Midpoint	Abs. Frequency	Rel. Frequency	Cumul. Rel. Freq.	Density
[17.4,17.5[	17.45	2	0.033333	0.033333	0.333333
[17.5,17.6[	17.55	1	0.016667	0.05	0.166667
[17.6,17.7[	17.65	4	0.066667	0.116667	0.666667
[17.7,17.8[	17.75	1	0.016667	0.133333	0.166667
[17.8,17.9[	17.85	2	0.033333	0.166667	0.333333
[17.9,18[	17.95	2	0.033333	0.2	0.333333
[18,18.1[	18.05	3	0.05	0.25	0.5
[18.1,18.2[	18.15	2	0.033333	0.283333	0.333333
[18.2,18.3[	18.25	2	0.033333	0.316667	0.333333
[18.3,18.4[	18.35	2	0.033333	0.35	0.333333
[18.4,18.5[	18.45	2	0.033333	0.383333	0.333333
[18.5,18.6[	18.55	0	0	0.383333	0
[18.6,18.7[	18.65	1	0.016667	0.4	0.166667
[18.7,18.8[	18.75	1	0.016667	0.416667	0.166667
[18.8,18.9[	18.85	3	0.05	0.466667	0.5
[18.9,19[	18.95	1	0.016667	0.483333	0.166667
[19,19.1[	19.05	1	0.016667	0.5	0.166667
[19.1,19.2[	19.15	1	0.016667	0.516667	0.166667
[19.2,19.3[	19.25	0	0	0.516667	0
[19.3,19.4[	19.35	3	0.05	0.566667	0.5
[19.4,19.5[	19.45	4	0.066667	0.633333	0.666667
[19.5,19.6[	19.55	4	0.066667	0.7	0.666667
[19.6,19.7[	19.65	2	0.033333	0.733333	0.333333
[19.7,19.8[	19.75	1	0.016667	0.75	0.166667
[19.8,19.9[	19.85	4	0.066667	0.816667	0.666667
[19.9,20[	19.95	1	0.016667	0.833333	0.166667
[20,20.1[	20.05	3	0.05	0.883333	0.5
[20.1,20.2[	20.15	0	0	0.883333	0
[20.2,20.3[	20.25	1	0.016667	0.9	0.166667
[20.3,20.4[	20.35	0	0	0.9	0
[20.4,20.5[	20.45	2	0.033333	0.933333	0.333333
[20.5,20.6[	20.55	2	0.033333	0.966667	0.333333
[20.6,20.7[	20.65	0	0	0.966667	0
[20.7,20.8[	20.75	1	0.016667	0.983333	0.166667
[20.8,20.9[	20.85	0	0	0.983333	0
[20.9,21[	20.95	0	0	0.983333	0
[21,21.1[	21.05	0	0	0.983333	0
[21.1,21.2]	21.15	1	0.016667	1	0.166667

\begin{tabular}{lllllllll}
\hline
Frequency Table (Histogram) \tabularnewline
Bins & Midpoint & Abs. Frequency & Rel. Frequency & Cumul. Rel. Freq. & Density \tabularnewline
[17.4,17.5[ & 17.45 & 2 & 0.033333 & 0.033333 & 0.333333 \tabularnewline
[17.5,17.6[ & 17.55 & 1 & 0.016667 & 0.05 & 0.166667 \tabularnewline
[17.6,17.7[ & 17.65 & 4 & 0.066667 & 0.116667 & 0.666667 \tabularnewline
[17.7,17.8[ & 17.75 & 1 & 0.016667 & 0.133333 & 0.166667 \tabularnewline
[17.8,17.9[ & 17.85 & 2 & 0.033333 & 0.166667 & 0.333333 \tabularnewline
[17.9,18[ & 17.95 & 2 & 0.033333 & 0.2 & 0.333333 \tabularnewline
[18,18.1[ & 18.05 & 3 & 0.05 & 0.25 & 0.5 \tabularnewline
[18.1,18.2[ & 18.15 & 2 & 0.033333 & 0.283333 & 0.333333 \tabularnewline
[18.2,18.3[ & 18.25 & 2 & 0.033333 & 0.316667 & 0.333333 \tabularnewline
[18.3,18.4[ & 18.35 & 2 & 0.033333 & 0.35 & 0.333333 \tabularnewline
[18.4,18.5[ & 18.45 & 2 & 0.033333 & 0.383333 & 0.333333 \tabularnewline
[18.5,18.6[ & 18.55 & 0 & 0 & 0.383333 & 0 \tabularnewline
[18.6,18.7[ & 18.65 & 1 & 0.016667 & 0.4 & 0.166667 \tabularnewline
[18.7,18.8[ & 18.75 & 1 & 0.016667 & 0.416667 & 0.166667 \tabularnewline
[18.8,18.9[ & 18.85 & 3 & 0.05 & 0.466667 & 0.5 \tabularnewline
[18.9,19[ & 18.95 & 1 & 0.016667 & 0.483333 & 0.166667 \tabularnewline
[19,19.1[ & 19.05 & 1 & 0.016667 & 0.5 & 0.166667 \tabularnewline
[19.1,19.2[ & 19.15 & 1 & 0.016667 & 0.516667 & 0.166667 \tabularnewline
[19.2,19.3[ & 19.25 & 0 & 0 & 0.516667 & 0 \tabularnewline
[19.3,19.4[ & 19.35 & 3 & 0.05 & 0.566667 & 0.5 \tabularnewline
[19.4,19.5[ & 19.45 & 4 & 0.066667 & 0.633333 & 0.666667 \tabularnewline
[19.5,19.6[ & 19.55 & 4 & 0.066667 & 0.7 & 0.666667 \tabularnewline
[19.6,19.7[ & 19.65 & 2 & 0.033333 & 0.733333 & 0.333333 \tabularnewline
[19.7,19.8[ & 19.75 & 1 & 0.016667 & 0.75 & 0.166667 \tabularnewline
[19.8,19.9[ & 19.85 & 4 & 0.066667 & 0.816667 & 0.666667 \tabularnewline
[19.9,20[ & 19.95 & 1 & 0.016667 & 0.833333 & 0.166667 \tabularnewline
[20,20.1[ & 20.05 & 3 & 0.05 & 0.883333 & 0.5 \tabularnewline
[20.1,20.2[ & 20.15 & 0 & 0 & 0.883333 & 0 \tabularnewline
[20.2,20.3[ & 20.25 & 1 & 0.016667 & 0.9 & 0.166667 \tabularnewline
[20.3,20.4[ & 20.35 & 0 & 0 & 0.9 & 0 \tabularnewline
[20.4,20.5[ & 20.45 & 2 & 0.033333 & 0.933333 & 0.333333 \tabularnewline
[20.5,20.6[ & 20.55 & 2 & 0.033333 & 0.966667 & 0.333333 \tabularnewline
[20.6,20.7[ & 20.65 & 0 & 0 & 0.966667 & 0 \tabularnewline
[20.7,20.8[ & 20.75 & 1 & 0.016667 & 0.983333 & 0.166667 \tabularnewline
[20.8,20.9[ & 20.85 & 0 & 0 & 0.983333 & 0 \tabularnewline
[20.9,21[ & 20.95 & 0 & 0 & 0.983333 & 0 \tabularnewline
[21,21.1[ & 21.05 & 0 & 0 & 0.983333 & 0 \tabularnewline
[21.1,21.2] & 21.15 & 1 & 0.016667 & 1 & 0.166667 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=162280&T=1

[TABLE]
[ROW][C]Frequency Table (Histogram)[/C][/ROW]
[ROW][C]Bins[/C][C]Midpoint[/C][C]Abs. Frequency[/C][C]Rel. Frequency[/C][C]Cumul. Rel. Freq.[/C][C]Density[/C][/ROW]
[ROW][C][17.4,17.5[[/C][C]17.45[/C][C]2[/C][C]0.033333[/C][C]0.033333[/C][C]0.333333[/C][/ROW]
[ROW][C][17.5,17.6[[/C][C]17.55[/C][C]1[/C][C]0.016667[/C][C]0.05[/C][C]0.166667[/C][/ROW]
[ROW][C][17.6,17.7[[/C][C]17.65[/C][C]4[/C][C]0.066667[/C][C]0.116667[/C][C]0.666667[/C][/ROW]
[ROW][C][17.7,17.8[[/C][C]17.75[/C][C]1[/C][C]0.016667[/C][C]0.133333[/C][C]0.166667[/C][/ROW]
[ROW][C][17.8,17.9[[/C][C]17.85[/C][C]2[/C][C]0.033333[/C][C]0.166667[/C][C]0.333333[/C][/ROW]
[ROW][C][17.9,18[[/C][C]17.95[/C][C]2[/C][C]0.033333[/C][C]0.2[/C][C]0.333333[/C][/ROW]
[ROW][C][18,18.1[[/C][C]18.05[/C][C]3[/C][C]0.05[/C][C]0.25[/C][C]0.5[/C][/ROW]
[ROW][C][18.1,18.2[[/C][C]18.15[/C][C]2[/C][C]0.033333[/C][C]0.283333[/C][C]0.333333[/C][/ROW]
[ROW][C][18.2,18.3[[/C][C]18.25[/C][C]2[/C][C]0.033333[/C][C]0.316667[/C][C]0.333333[/C][/ROW]
[ROW][C][18.3,18.4[[/C][C]18.35[/C][C]2[/C][C]0.033333[/C][C]0.35[/C][C]0.333333[/C][/ROW]
[ROW][C][18.4,18.5[[/C][C]18.45[/C][C]2[/C][C]0.033333[/C][C]0.383333[/C][C]0.333333[/C][/ROW]
[ROW][C][18.5,18.6[[/C][C]18.55[/C][C]0[/C][C]0[/C][C]0.383333[/C][C]0[/C][/ROW]
[ROW][C][18.6,18.7[[/C][C]18.65[/C][C]1[/C][C]0.016667[/C][C]0.4[/C][C]0.166667[/C][/ROW]
[ROW][C][18.7,18.8[[/C][C]18.75[/C][C]1[/C][C]0.016667[/C][C]0.416667[/C][C]0.166667[/C][/ROW]
[ROW][C][18.8,18.9[[/C][C]18.85[/C][C]3[/C][C]0.05[/C][C]0.466667[/C][C]0.5[/C][/ROW]
[ROW][C][18.9,19[[/C][C]18.95[/C][C]1[/C][C]0.016667[/C][C]0.483333[/C][C]0.166667[/C][/ROW]
[ROW][C][19,19.1[[/C][C]19.05[/C][C]1[/C][C]0.016667[/C][C]0.5[/C][C]0.166667[/C][/ROW]
[ROW][C][19.1,19.2[[/C][C]19.15[/C][C]1[/C][C]0.016667[/C][C]0.516667[/C][C]0.166667[/C][/ROW]
[ROW][C][19.2,19.3[[/C][C]19.25[/C][C]0[/C][C]0[/C][C]0.516667[/C][C]0[/C][/ROW]
[ROW][C][19.3,19.4[[/C][C]19.35[/C][C]3[/C][C]0.05[/C][C]0.566667[/C][C]0.5[/C][/ROW]
[ROW][C][19.4,19.5[[/C][C]19.45[/C][C]4[/C][C]0.066667[/C][C]0.633333[/C][C]0.666667[/C][/ROW]
[ROW][C][19.5,19.6[[/C][C]19.55[/C][C]4[/C][C]0.066667[/C][C]0.7[/C][C]0.666667[/C][/ROW]
[ROW][C][19.6,19.7[[/C][C]19.65[/C][C]2[/C][C]0.033333[/C][C]0.733333[/C][C]0.333333[/C][/ROW]
[ROW][C][19.7,19.8[[/C][C]19.75[/C][C]1[/C][C]0.016667[/C][C]0.75[/C][C]0.166667[/C][/ROW]
[ROW][C][19.8,19.9[[/C][C]19.85[/C][C]4[/C][C]0.066667[/C][C]0.816667[/C][C]0.666667[/C][/ROW]
[ROW][C][19.9,20[[/C][C]19.95[/C][C]1[/C][C]0.016667[/C][C]0.833333[/C][C]0.166667[/C][/ROW]
[ROW][C][20,20.1[[/C][C]20.05[/C][C]3[/C][C]0.05[/C][C]0.883333[/C][C]0.5[/C][/ROW]
[ROW][C][20.1,20.2[[/C][C]20.15[/C][C]0[/C][C]0[/C][C]0.883333[/C][C]0[/C][/ROW]
[ROW][C][20.2,20.3[[/C][C]20.25[/C][C]1[/C][C]0.016667[/C][C]0.9[/C][C]0.166667[/C][/ROW]
[ROW][C][20.3,20.4[[/C][C]20.35[/C][C]0[/C][C]0[/C][C]0.9[/C][C]0[/C][/ROW]
[ROW][C][20.4,20.5[[/C][C]20.45[/C][C]2[/C][C]0.033333[/C][C]0.933333[/C][C]0.333333[/C][/ROW]
[ROW][C][20.5,20.6[[/C][C]20.55[/C][C]2[/C][C]0.033333[/C][C]0.966667[/C][C]0.333333[/C][/ROW]
[ROW][C][20.6,20.7[[/C][C]20.65[/C][C]0[/C][C]0[/C][C]0.966667[/C][C]0[/C][/ROW]
[ROW][C][20.7,20.8[[/C][C]20.75[/C][C]1[/C][C]0.016667[/C][C]0.983333[/C][C]0.166667[/C][/ROW]
[ROW][C][20.8,20.9[[/C][C]20.85[/C][C]0[/C][C]0[/C][C]0.983333[/C][C]0[/C][/ROW]
[ROW][C][20.9,21[[/C][C]20.95[/C][C]0[/C][C]0[/C][C]0.983333[/C][C]0[/C][/ROW]
[ROW][C][21,21.1[[/C][C]21.05[/C][C]0[/C][C]0[/C][C]0.983333[/C][C]0[/C][/ROW]
[ROW][C][21.1,21.2][/C][C]21.15[/C][C]1[/C][C]0.016667[/C][C]1[/C][C]0.166667[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=162280&T=1

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

As an alternative you can also use a QR Code:

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

Frequency Table (Histogram)
Bins	Midpoint	Abs. Frequency	Rel. Frequency	Cumul. Rel. Freq.	Density
[17.4,17.5[	17.45	2	0.033333	0.033333	0.333333
[17.5,17.6[	17.55	1	0.016667	0.05	0.166667
[17.6,17.7[	17.65	4	0.066667	0.116667	0.666667
[17.7,17.8[	17.75	1	0.016667	0.133333	0.166667
[17.8,17.9[	17.85	2	0.033333	0.166667	0.333333
[17.9,18[	17.95	2	0.033333	0.2	0.333333
[18,18.1[	18.05	3	0.05	0.25	0.5
[18.1,18.2[	18.15	2	0.033333	0.283333	0.333333
[18.2,18.3[	18.25	2	0.033333	0.316667	0.333333
[18.3,18.4[	18.35	2	0.033333	0.35	0.333333
[18.4,18.5[	18.45	2	0.033333	0.383333	0.333333
[18.5,18.6[	18.55	0	0	0.383333	0
[18.6,18.7[	18.65	1	0.016667	0.4	0.166667
[18.7,18.8[	18.75	1	0.016667	0.416667	0.166667
[18.8,18.9[	18.85	3	0.05	0.466667	0.5
[18.9,19[	18.95	1	0.016667	0.483333	0.166667
[19,19.1[	19.05	1	0.016667	0.5	0.166667
[19.1,19.2[	19.15	1	0.016667	0.516667	0.166667
[19.2,19.3[	19.25	0	0	0.516667	0
[19.3,19.4[	19.35	3	0.05	0.566667	0.5
[19.4,19.5[	19.45	4	0.066667	0.633333	0.666667
[19.5,19.6[	19.55	4	0.066667	0.7	0.666667
[19.6,19.7[	19.65	2	0.033333	0.733333	0.333333
[19.7,19.8[	19.75	1	0.016667	0.75	0.166667
[19.8,19.9[	19.85	4	0.066667	0.816667	0.666667
[19.9,20[	19.95	1	0.016667	0.833333	0.166667
[20,20.1[	20.05	3	0.05	0.883333	0.5
[20.1,20.2[	20.15	0	0	0.883333	0
[20.2,20.3[	20.25	1	0.016667	0.9	0.166667
[20.3,20.4[	20.35	0	0	0.9	0
[20.4,20.5[	20.45	2	0.033333	0.933333	0.333333
[20.5,20.6[	20.55	2	0.033333	0.966667	0.333333
[20.6,20.7[	20.65	0	0	0.966667	0
[20.7,20.8[	20.75	1	0.016667	0.983333	0.166667
[20.8,20.9[	20.85	0	0	0.983333	0
[20.9,21[	20.95	0	0	0.983333	0
[21,21.1[	21.05	0	0	0.983333	0
[21.1,21.2]	21.15	1	0.016667	1	0.166667

Figure 1

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = 30 ; par2 = grey ; par3 = FALSE ; par4 = Unknown ;

Parameters (R input):

par1 = 30 ; par2 = grey ; par3 = FALSE ; par4 = Unknown ;

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

par1 <- as.numeric(par1)
if (par3 == 'TRUE') par3 <- TRUE
if (par3 == 'FALSE') par3 <- FALSE
if (par4 == 'Unknown') par1 <- as.numeric(par1)
if (par4 == 'Interval/Ratio') par1 <- as.numeric(par1)
if (par4 == '3-point Likert') par1 <- c(1:3 - 0.5, 3.5)
if (par4 == '4-point Likert') par1 <- c(1:4 - 0.5, 4.5)
if (par4 == '5-point Likert') par1 <- c(1:5 - 0.5, 5.5)
if (par4 == '6-point Likert') par1 <- c(1:6 - 0.5, 6.5)
if (par4 == '7-point Likert') par1 <- c(1:7 - 0.5, 7.5)
if (par4 == '8-point Likert') par1 <- c(1:8 - 0.5, 8.5)
if (par4 == '9-point Likert') par1 <- c(1:9 - 0.5, 9.5)
if (par4 == '10-point Likert') par1 <- c(1:10 - 0.5, 10.5)
bitmap(file='test1.png')
if(is.numeric(x[1])) {
if (is.na(par1)) {
myhist<-hist(x,col=par2,main=main,xlab=xlab,right=par3)
} else {
if (par1 < 0) par1 <- 3
if (par1 > 50) par1 <- 50
myhist<-hist(x,breaks=par1,col=par2,main=main,xlab=xlab,right=par3)
}
} else {
plot(mytab <- table(x),col=par2,main='Frequency Plot',xlab=xlab,ylab='Absolute Frequency')
}
dev.off()
if(is.numeric(x[1])) {
myhist
n <- length(x)
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,hyperlink('histogram.htm','Frequency Table (Histogram)',''),6,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Bins',header=TRUE)
a<-table.element(a,'Midpoint',header=TRUE)
a<-table.element(a,'Abs. Frequency',header=TRUE)
a<-table.element(a,'Rel. Frequency',header=TRUE)
a<-table.element(a,'Cumul. Rel. Freq.',header=TRUE)
a<-table.element(a,'Density',header=TRUE)
a<-table.row.end(a)
crf <- 0
if (par3 == FALSE) mybracket <- '[' else mybracket <- ']'
mynumrows <- (length(myhist$breaks)-1)
for (i in 1:mynumrows) {
a<-table.row.start(a)
if (i == 1)
dum <- paste('[',myhist$breaks[i],sep='')
else
dum <- paste(mybracket,myhist$breaks[i],sep='')
dum <- paste(dum,myhist$breaks[i+1],sep=',')
if (i==mynumrows)
dum <- paste(dum,']',sep='')
else
dum <- paste(dum,mybracket,sep='')
a<-table.element(a,dum,header=TRUE)
a<-table.element(a,myhist$mids[i])
a<-table.element(a,myhist$counts[i])
rf <- myhist$counts[i]/n
crf <- crf + rf
a<-table.element(a,round(rf,6))
a<-table.element(a,round(crf,6))
a<-table.element(a,round(myhist$density[i],6))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
} else {
mytab
reltab <- mytab / sum(mytab)
n <- length(mytab)
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Frequency Table (Categorical Data)',3,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Category',header=TRUE)
a<-table.element(a,'Abs. Frequency',header=TRUE)
a<-table.element(a,'Rel. Frequency',header=TRUE)
a<-table.row.end(a)
for (i in 1:n) {
a<-table.row.start(a)
a<-table.element(a,labels(mytab)$x[i],header=TRUE)
a<-table.element(a,mytab[i])
a<-table.element(a,round(reltab[i],4))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.tab')
}

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code