File: LokSkala.R

package info (click to toggle)
jsxgraph 0.99.5%2Bdfsg1-1
  • links: PTS
  • area: main
  • in suites: stretch
  • size: 25,808 kB
  • sloc: xml: 5,869; java: 1,072; python: 710; php: 363; makefile: 159; sh: 47
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 = ";")