diff command.R @ 0:913ea6991ee4 draft

Uploaded
author sebastian-boegel
date Thu, 20 Dec 2012 10:37:01 -0500
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/command.R	Thu Dec 20 10:37:01 2012 -0500
@@ -0,0 +1,24 @@
+args <- commandArgs(trailingOnly = TRUE)
+
+x<-unlist(strsplit(args[2],split=","))
+x<-as.numeric(x)
+args[3]
+paril<-1-pnorm(as.numeric(args[1]),mean(x),sd(x))
+poutlier<-pbinom(0,length(x),paril)
+1-poutlier
+
+x<-unlist(strsplit(args[5],split=","))
+x<-as.numeric(x)
+args[6]
+#1-exp(log(pnorm(as.numeric(args[4]),mean(x),sd(x)))*length(x))
+paril<-1-pnorm(as.numeric(args[4]),mean(x),sd(x))
+poutlier<-pbinom(0,length(x),paril)
+1-poutlier
+
+x<-unlist(strsplit(args[8],split=","))
+x<-as.numeric(x)
+args[9]
+#1-exp(log(pnorm(as.numeric(args[7]),mean(x),sd(x)))*length(x))
+paril<-1-pnorm(as.numeric(args[7]),mean(x),sd(x))
+poutlier<-pbinom(0,length(x),paril)
+1-poutlier