Mercurial > repos > sebastian-boegel > seq2hla
view command.R @ 1:155b796033b6 draft default tip
Uploaded
author | sebastian-boegel |
---|---|
date | Fri, 21 Dec 2012 03:46:15 -0500 |
parents | 913ea6991ee4 |
children |
line wrap: on
line source
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