This is a test version of Biostars. For the public version, visit https://www.biostars.org.
explanation of code
observed_statistics_np = np.array(observed_statistics) 
allDifs = []
goodIndexes = []
for i in range(num_replicates):
    d = np.linalg.norm(observed_statistics_np - np.array(sim_stats[i]))
    allDifs.append(d)
    if( d < 50):
        goodIndexes.append(i)

allDifs
len(goodIndexes)

Can anyone please help me understand what this code does ?

np.linalg.norm statistics
• 618 views • link
updated by rpolicastro rpolicastro • written by desptsp • desptsp

0 answers

No answers yet.

Log in to answer this question.

More posts like this