diff bin/calc_pval.R @ 1:adc0f7765d85 draft

planemo upload

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/bin/calc_pval.R	Thu Sep 07 15:06:58 2017 -0400
@@ -0,0 +1,29 @@
+args = commandArgs();
+input_file=args[4];
+av=as.numeric(args[5]);
+sd=as.numeric(args[6]);
+nu=as.numeric(args[7]);
+
+
+data=read.table(input_file);
+
+col=4;
+x=data[,col];
+pup=pnorm(x, mean=av, sd=sd, lower.tail = FALSE);
+adpup=p.adjust(pup,method="fdr");
+pdn=pnorm(x, mean=av, sd=sd, lower.tail = TRUE);
+adpdn=p.adjust(pdn,method="fdr");
+
+p=pup;
+
+size = length(x);
+
+for (i in 1:size)
+{
+	if(x[i]<0) 
+	{
+		p[i]=pdn[i];
+	}
+}
+write(p, file=paste(input_file,"p",sep="."),ncolumns=1);
+