Mercurial > repos > elixir-it > vinyl_optimizer

diff GENEO_VINYL.R @ 3:221db2eb3c8e draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author elixir-it
date Wed, 22 Jul 2020 19:23:14 +0000
parents 6e4eb4856874
children
line wrap: on
line diff
--- a/GENEO_VINYL.R	Wed Jul 22 19:20:30 2020 +0000
+++ b/GENEO_VINYL.R	Wed Jul 22 19:23:14 2020 +0000
@@ -52,7 +52,7 @@
 
         	F=fisher.test(m,alternative="greater")
         	Fpv=F$p.value
-        	Fodds=F$estimate
+		Fodds=F$estimate
 		localScore=0.5*-log10(Fpv)+0.3*Fodds-0.2*posT #+0.175*posR-0.125*posT
 		if (localScore>score)
 		{
@@ -68,11 +68,14 @@
 	}
 	Command=system("rm *.ofile *.ovcfile *.osummary",intern=FALSE)
 	outV=paste(round(x,digits=2),collapse="\t");
+	if (surv1>60)
+	{
+		surv1=62;
+	}
 	cat(paste(outV,surv1,totR,surv2,totT,pval,rat,score,"\n",sep="\t"),file=ofile,append=T);
 	return(score*-1);
 }
 library(genalg)
 
-G=rbga(stringMin=minsV,stringMax=maxV,popSize=20,iters=25,evalFunc=evalVINYL)
+G=rbga(stringMin=minsV,stringMax=maxV,popSize=50,iters=50,evalFunc=evalVINYL)
 #cat(summary(G),file=ofile,append=T)	
-