annotate command.R @ 1:155b796033b6 draft default tip

Uploaded
author sebastian-boegel
date Fri, 21 Dec 2012 03:46:15 -0500
parents 913ea6991ee4
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
1 args <- commandArgs(trailingOnly = TRUE)
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
2
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
3 x<-unlist(strsplit(args[2],split=","))
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
4 x<-as.numeric(x)
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
5 args[3]
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
6 paril<-1-pnorm(as.numeric(args[1]),mean(x),sd(x))
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
7 poutlier<-pbinom(0,length(x),paril)
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
8 1-poutlier
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
9
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
10 x<-unlist(strsplit(args[5],split=","))
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
11 x<-as.numeric(x)
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
12 args[6]
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
13 #1-exp(log(pnorm(as.numeric(args[4]),mean(x),sd(x)))*length(x))
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
14 paril<-1-pnorm(as.numeric(args[4]),mean(x),sd(x))
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
15 poutlier<-pbinom(0,length(x),paril)
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
16 1-poutlier
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
17
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
18 x<-unlist(strsplit(args[8],split=","))
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
19 x<-as.numeric(x)
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
20 args[9]
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
21 #1-exp(log(pnorm(as.numeric(args[7]),mean(x),sd(x)))*length(x))
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
22 paril<-1-pnorm(as.numeric(args[7]),mean(x),sd(x))
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
23 poutlier<-pbinom(0,length(x),paril)
913ea6991ee4 Uploaded
sebastian-boegel
parents:
diff changeset
24 1-poutlier