yairf / SDS

R code to compute the Singleton Density Score (SDS)
28 stars 9 forks source link

issue about individual with no singleton #3

Open WeiCSong opened 3 years ago

WeiCSong commented 3 years ago

Hi Dr. Field, I'm running SDS on a simulation dataset where a small number of individuals did not have any singleton (the simulated chromosome was relatively short). Running compute_SDS returned an empty .txt, and I found that when reading in s-file,read.table collapsed all empty lines and broke the order of individuals. My question is, does singleton number=0 significantly violates the basic assumption of SDS? If not, can I modify the compute_SDS to restore empty lines when reading s_file? Thanks in advance for your help.