# HG changeset patch # User cropgeeks # Date 1524218318 14400 # Node ID 12db48e80a98894e4b109fe3cd88931c2b54d437 # Parent 5dab264e38c2e65f62bb0cca1bc50a844b193957 Uploaded diff -r 5dab264e38c2 -r 12db48e80a98 dartseq_seeduk_1c.R --- 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.