Unknown Parent Groups and PyPedal

I had a visiting scientist in the lab the last couple of weeks and had to get a manuscript finished. But here we are again. Here’s the deal — no unknown parent groups in version 2.0. It turns out to be harder to do well that I first thought, so rather than do a poor job I decided to punt it. However, I’m not going to leave you hanging, so the 2.0 release will come with only two three small fixes:

  • Change the processing of ASD pedigrees to use hashlib instead of the md5 module. This may result in user-selectable hash functions in the event that you have trouble with collisions.
  • For some reason, the logfile is not being created properly. That needs fixed.
  • The NewPedigree::save() method is going to get a new argument so that you can specify the format of the output file.

So, okay, less talking and more coding.

Posted in Articles | Tagged | Leave a comment

So, what’s the deal with PyPedal?

I’m glad you asked that. If you’ve seen the website recently then you know that there’s been no development activity for quite a long time. What can I say? I’ve been pretty busy on other things. The important point is that PyPedal is not abandonware.

On September 1, 2010, I uploaded Release Candidate 9 and let it go into the wild. In addition to adding GENES 1.20 input and output, the changes include some minor documentation updates and a few small code fixes.

In Release Candidate 10, which will be released this week, all I’ve done is add an option to pyp_nrm.fast_a_matrix() to add coefficients of inbreeding read from the pedigree file to the diagonal elements for founders when forming the relationship matrix. I need to update the documentation and check the files once more because I think that there are still some stray hard-coded debugging messages which need to be removed.

The only other thing I’m thinking of adding is support for unknown parent groups, and possibly adding methods to the NewPedigree class that support adding (merging) and subtracting (splitting) pedigrees, but I think that’s going to get punted to the 2.1 development branch because it’s a much harder problem than it may sound on the surface.

I still need bug reports and feature requests. I’ll do my best to fix bugs promptly. I can’t  make any promises about feature requests. My thinking right now is that I need to focus on polishing things. The demographics (pyp_demog), graphics (pyp_graphics) and reporting (pyp_reports) modules all could use a lot of fine-tuning to make them more useful. The documentation always needs work. I’m still concerned about the implementation of the algorithms in pyp_metrics for calculating effective ancestor numbers. They’re much less fun that implementing new features, but probably necessary. And then there’s the whole issue of dense SNP data…

Posted in Articles | Tagged | Leave a comment

ElleryPedia

Can’t get enough of Thing One’s wit and wisdom? Visit ElleryPedia for more!

Posted in Articles | Tagged | Leave a comment

About to board MSY to BWI.

Posted on by john.b.cole | Leave a comment

At BWI to start the odyssey.

Posted on by john.b.cole | Leave a comment

Tweet for Reward Points Day 2! I want to win 500 RP from @mafiawars #mafiawars! http://zyn.ga/AF

Posted on by john.b.cole | Leave a comment

I’m #reading Firehouse by David Halberstam http://bit.ly/dcKF3L

Posted on by john.b.cole | Leave a comment

I’m #reading The Family by Jeff Sharlet http://bit.ly/aRrKn6

Posted on by john.b.cole | Leave a comment

I’m #reading The Family by Jeff Sharlet http://bit.ly/aRrKn6

Posted on by john.b.cole | Leave a comment

I’m #reading The Catcher in the Rye by J.D. Salinger http://bit.ly/cDWPhe

Posted on by john.b.cole | Leave a comment