Mercurial > repos > sebastian-boegel > seq2hla
comparison command.R @ 0:913ea6991ee4 draft
Uploaded
author | sebastian-boegel |
---|---|
date | Thu, 20 Dec 2012 10:37:01 -0500 |
parents | |
children |
comparison
equal
deleted
inserted
replaced
-1:000000000000 | 0:913ea6991ee4 |
---|---|
1 args <- commandArgs(trailingOnly = TRUE) | |
2 | |
3 x<-unlist(strsplit(args[2],split=",")) | |
4 x<-as.numeric(x) | |
5 args[3] | |
6 paril<-1-pnorm(as.numeric(args[1]),mean(x),sd(x)) | |
7 poutlier<-pbinom(0,length(x),paril) | |
8 1-poutlier | |
9 | |
10 x<-unlist(strsplit(args[5],split=",")) | |
11 x<-as.numeric(x) | |
12 args[6] | |
13 #1-exp(log(pnorm(as.numeric(args[4]),mean(x),sd(x)))*length(x)) | |
14 paril<-1-pnorm(as.numeric(args[4]),mean(x),sd(x)) | |
15 poutlier<-pbinom(0,length(x),paril) | |
16 1-poutlier | |
17 | |
18 x<-unlist(strsplit(args[8],split=",")) | |
19 x<-as.numeric(x) | |
20 args[9] | |
21 #1-exp(log(pnorm(as.numeric(args[7]),mean(x),sd(x)))*length(x)) | |
22 paril<-1-pnorm(as.numeric(args[7]),mean(x),sd(x)) | |
23 poutlier<-pbinom(0,length(x),paril) | |
24 1-poutlier |