search

ness-lab / mrt

0 stars 2 forks source link

Parse LDhat output to obtain rho #8

Closed James-S-Santangelo closed 4 years ago

James-S-Santangelo commented 4 years ago

Depends on successful completion of #7

The LDhat output for all simulations will have to be parse to get rho, which will likely be easiest using Python.

Reminder: The rho value needs to be divided by the sequence length (currently fixed at 1e8), which can then be divided by r to get Ne

It would be best if this script iterates over LDhat output files, calculates rho and Ne, and writes these values to a CSV with separate rows for each parameter combination (which should be maintained in filenames). We can then join this CSV with the theta CSV generated by this script for analysis in R.

James-S-Santangelo commented 4 years ago

Closing since no longer using LDhat