File: LokSkala.R

package info (click to toggle)
jsxgraph 1.5.0%2Bdfsg1-3
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 24,020 kB
  • sloc: javascript: 76,112; xml: 5,869; java: 1,072; python: 767; php: 181; makefile: 154; cpp: 76; sh: 9
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 = ";")