laninsky / genetic_diversity_diffs Goto Github PK

Dear Alana,
I am trying to run your script on my sequence data, I have short fragment (397 bp) of haploid sequence samples. It runs well at the beginning until I get an error message, please see a copy of my screen below:

[1] STEP ONE: Loading in all the variables
[1]
[1] An error message after this indicates your working directory is not valid
[1] Not to worry, your working directory IS valid! I've successfully set the working directory
[1]
[1] An error message after this indicates your file is not located in the directory you listed
[1] Not to worry, your file IS located in the directory! I'm pulling it into my memory to extract the parts we are interested in
[1]
[1] I have found the haplotype definition block in your arlequin file
Error in if (minhaplength == maxhaplength) { :
missing value where TRUE/FALSE needed

The script I have used is :
genetic_diversity_diffs("C:/Users/mbx/Desktop/pvalue","samples.arp",1000,"Y","Y")

My data file has been modified to meet the same criteria for the example file that found with the script, I have uploaded a text data file similar to the .arp one I have used.

kindest regards,

Ahmed

samples.txt