Mercurial > repos > sebastian-boegel > seq2hla
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