3
function preprocess () {
4
echo ('require(moments)\n');
7
function calculate () {
8
vars = "substitute (" + str_replace ("\n", "), substitute (", trim (getValue ("x"))) + ")";
9
if (getValue ("narm")) narm = ", na.rm=TRUE";
10
else narm = ", na.rm=FALSE"
13
echo ('objects <- list (' + vars + ')\n');
14
echo ('results <- data.frame (\'Variable Name\'=rep (NA, length (objects)), check.names=FALSE)\n');
15
echo ('for (i in 1:length (objects)) {\n');
16
echo (' var <- eval (objects[[i]], envir=globalenv ())\n');
17
echo (' results[i, \'Variable Name\'] <- rk.get.description (objects[[i]], is.substitute=TRUE)\n');
20
if (getValue ("skewness")) {
21
echo (' results[i, \'Skewness\'] <- skewness (var' + narm + ')\n');
23
if (getValue ("kurtosis")) {
24
echo (' results[i, \'Kurtosis\'] <- kurtosis (var' + narm + ')\n');
25
echo (' results[i, \'Excess Kurtosis\'] <- results[i, \'Kurtosis\'] - 3\n');
27
if (getValue ("geary")) {
28
echo (' results[i, \'Geary Kurtosis\'] <- geary (var' + narm + ')\n');
31
if (getValue ("length")) {
33
echo (' results[i, \'Length\'] <- length (var)\n');
34
echo (' results[i, \'NAs\'] <- sum (is.na(var))\n');
38
function printout () {
39
echo ('rk.header ("Skewness and Kurtosis")\n');
40
echo ('rk.results (results)\n');