File: LokSkala.R

package info (click to toggle)
jsxgraph 1.11.1%2Bdfsg1-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 24,984 kB
  • sloc: javascript: 81,688; xml: 5,869; java: 1,072; php: 281; makefile: 189; python: 174; cpp: 76; sh: 12
file content (19 lines) | stat: -rw-r--r-- 407 bytes parent folder | download | duplicates (8)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
x <- commandArgs(TRUE)
x <- as.numeric(unlist(strsplit(x, ";")))

MW <- mean(x)
SD <- sd(x)
Med <- median(x)
Mad <- mad(x)

## Radius = 0.6
A1 <- 1.7756*Mad^2
A2 <- 1.4322*Mad^2
a2 <- -0.5034*Mad
b <- 1.8390*Mad
u <- A1*(x-Med)/Mad^2
v <- A2*(((x-Med)/Mad)^2-1)/Mad - a2
w <- pmin(1, b/sqrt(u^2 + v^2))
res <- c(Med, Mad) + c(mean(u*w), mean(v*w))

paste(round(c(MW, SD, Med, Mad, res), 3), collapse = ";")