diff dartseq_seeduk_1c.R @ 22:12db48e80a98 draft

Uploaded
author cropgeeks
date Fri, 20 Apr 2018 05:58:38 -0400
parents ef0611a290ea
children
line wrap: on
line diff
--- a/dartseq_seeduk_1c.R	Fri Apr 20 05:35:52 2018 -0400
+++ b/dartseq_seeduk_1c.R	Fri Apr 20 05:58:38 2018 -0400
@@ -65,7 +65,7 @@
 #[1] "Completed"
 #
 #
-gl_call_rate <-  gl.filter.callrate(gl,method = 'loc', t=args[2])
+gl_call_rate <-  gl.filter.callrate(gl,method = 'loc', t=as.double(args[2]))
 #Reporting for a genlight object
 #Note: Missing values most commonly arise from restriction site mutation.
 #
@@ -77,7 +77,7 @@
 #  No. of individuals: 376 
 #  No. of populations:  0 
 #
-gl_rep <- gl.filter.repavg(gl_call_rate,t=args[3])
+gl_rep <- gl.filter.repavg(gl_call_rate,t=as.double(args[3]))
 #Reporting for a genlight object
 #Note: RepAvg is a DArT statistic reporting reproducibility averaged across alleles for each locus. 
 
@@ -108,7 +108,7 @@
 #Individuals with less than 60% missing values = 375 [99.7%]
 #[1] "Completed"
 
-gl_final <-  gl.filter.callrate(gl_rep,method = 'ind', t=args[4])
+gl_final <-  gl.filter.callrate(gl_rep,method = 'ind', t=as.double(args[4]))
 #Reporting for a genlight object
 #Note: Missing values most commonly arise from restriction site mutation.