<div><br></div>Hi Yurii and Lennart,<div>Thanks for sending these thoughtful replies to my inquiry.</div><div><br></div><div>In reply to Yurii's questions:</div><div>- I would like changes that I make to be public and incorporated into one of the branches of the source code.</div>
<div>- I would also like to follow the "normal" path of submitting a patch for evaluation. I don't necessarily need commit-access to the repository.</div><div>- Regarding the version that I should start with:</div>
<div>  -- my natural inclination would be to start with the version in which the survival analysis is working (<span style="background-color:rgb(255,255,255);color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px">1.3.0-fixed or coxph-enabled 0.2.0) and then turn my attention to what isn't working in 0.2.0.  I guess what I do depends on the project's desires and the ease to which I can actually get the statistical tests working (with whatever file input) in 0.2.0. I haven't seen the 0.2.0 code (locally, I have </span><font color="#222222" face="arial, sans-serif">v.0.1-9e (2011.05.15) and I haven't even attempted to run survival analysis).</font></div>
<div><span style="background-color:rgb(255,255,255);color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="background-color:rgb(255,255,255);color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px">In reply to </span><span style="background-color:rgb(255,255,255);color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px">Lennart's queries:</span></div>
<div><span style="background-color:rgb(255,255,255);color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px">- I do want to focus on the implementation of the weighted Cox regression.</span></div><div><font color="#222222" face="arial, sans-serif">- I am familiar (in principle) with subversion and have scanned the first few chapters of the "Version Control with Subversion" book.  We have svn version 1.4.2 installed on the Broad servers; I'll probably install the latest version locally.</font></div>
<div><span style="background-color:rgb(255,255,255);color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px"><br></span></div><div><font color="#222222" face="arial, sans-serif">I don't want to promise that I can fix the filevector problem-- but perhaps the problem will become clear while working with the Cox regression module.</font></div>
<div><font color="#222222" face="arial, sans-serif"><br></font></div><div><font color="#222222" face="arial, sans-serif">~Alisa</font></div><div><span style="background-color:rgb(255,255,255);color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px"><br>
</span></div><div><br><div class="gmail_quote">On Wed, Jul 25, 2012 at 9:52 AM, Yury Aulchenko <span dir="ltr"><<a href="mailto:yurii.aulchenko@gmail.com" target="_blank">yurii.aulchenko@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Dear Lennart, Alisa, all,<div><div><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<div><div><div><br></div></div></div></div></span></div></span></div></span></div></span></div></span></span></div><div><div class="im"><div>On Jul 25, 2012, at 10:44 AM, L.C. Karssen wrote:</div><br><blockquote type="cite">
<div>In short, the status of Cox regression in ProbABEL is the following:<br>- In PA 0.2.0 Cox regression is disabled by default because it doesn't<br>work correctly (but it can be enabled easily). Looking back it seems<br>
that it broke with the integration of the filevector/DatABEL code,<br>somewhere between the 0.1-3 and 0.1-9 releases.<br>- For those who want functional Cox regression (although not with<br>filevector data) I suggest to use PA 0.1-3. Recently we've implemented a<br>
bugfix in that code base and we are almost ready to make that version<br>available to the public.<br><br>So, at the moment we have two code bases. Ideally we would like to get<br>the filevector problem solved in the PA 0.2.0 branch (the bugfix from<br>
0.1-3 has also been implemented there). If you think you are up to it<br>you are more than welcome to see if you can fix the filevector problem<br>in 0.2.0!<br><br>If you want to focus more on the implementation of the weighted Cox<br>
regression I guess it's easier to base your work on the 0.1-3-fixes code<br>where you can easily make test runs. Once it works you can forward-port<br>your improvements to the 0.2.0 branch.<br><br>Yurii, do you agree?<br>
<br></div></blockquote><div><br></div></div>This is very good summary of where we stand now.</div><div><br></div><div>The big decision is whether to </div><div><br></div><div>a) start with 0.2.0 and figure out why Cox does not work there (and then proceed to inclusion of robust and weighted). This may take time - we do not quite have a clue why it does not work. </div>
<div><br></div><div>or </div><div><br></div><div>b) start directly with 1.3.0-fixed (which is due to be released soon) and incorporate robust and weighted,  I think, must be relatively straightforward because actually the code is there, it is only that we do not pass right arguments and data for the code to do the job</div>
<div><br></div><div>A bit of problem with "b" is of cause that we will be getting more and more into that branch (which, because it can not use DatABEL/filevector format, is very RAM-demanding, while 0.2.0 can operate with <1Gb RAM). On the other side, if the thing which really matters is to get this functionality ASAP, and RAM is no concern, (b) is clearly the way to go. As Lennart mentioned, other possibility is to do (b) then go for (a) - though merging branches could be a nightmare.</div>
<div><br></div><div>Under any scenario you decide upon, we will be happy to help in planning and implementation (but see below).</div><div><br></div><div>Alisa, to make things clear from the start - do you plan to release you modifications of the code into "open", e.g. incorporate these modifications into one of the mainstream branches available for other people as part of ProbABEL releases? The reason I ask is that there are examples that people make "private" versions which they never release into open and keep for their private use. Of cause under latter scenario we have less motivation to actively help. </div>
<div><br></div><div><div class="im"><blockquote type="cite"><div><br>On the technical side of things: do you know how to work with<br>Subversion? That would make life easier, because then I could set up a<br>separate branch for you to work with. So that, once your feature is<br>
tested and works as expected, we can easily integrate the changes in the<br>main branch(es). (If all this sounds like gibberish submitting a patch<br>also works :-).)<br><br></div></blockquote><div><br></div></div>Just to comment: the "normal" procedure would of cause be that people first submit a patch and then can request for commit-access to repository. On the other side, Alisa already participated in development of ProbABEL. And on the final note, the decision of granting commit-access is the decision of current active developers (meaning - for ProbABEL - Lennart, and, in lesser degree, Maarten and me). I will be happy to go along with any decision worked out by Lennart and Alisa :)</div>
<div><br></div><div>best wishes,</div><div>Yurii<br><div><br></div><div><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<div><div><div>-------------------------------------------------------</div><div><div><div style="word-wrap:break-word"><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<div><div>Yurii S. Aulchenko, PhD, Dr. Habil</div><div>Director,YuriiA consulting</div><div>Rotterdam, The Netherlands</div><div>yurii [dot] aulchenko [monkey] gmail [dot] com</div></div><div><br></div></div></span></div>
</span></div></div></div></div></div></div></span></div></span></div></span></div></span></div></span></span></div><div><div class="h5"><blockquote type="cite"><div><br>Let us know what you think and if we can be of any help!<br>
<br>Looking forward to your contributions!<br><br>Lennart Karssen.<br><br><br><br>On 24/07/12 17:26, Alisa Manning wrote:<br><blockquote type="cite"><br></blockquote><blockquote type="cite">Hello--<br></blockquote><blockquote type="cite">
In the past, I've been involved in code contributions to ProbABEL<br></blockquote><blockquote type="cite">(extracting covariance between regression estimates for snp and snp*E<br></blockquote><blockquote type="cite">for interaction tests).<br>
</blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I've set aside some time this summer to work on implementing a weighted<br></blockquote><blockquote type="cite">Cox regression model in ProbABEL -- the only change that's required in<br>
</blockquote><blockquote type="cite">the code is to implement robust standard error estimation in the cox<br></blockquote><blockquote type="cite">proportional hazard functionality.<br></blockquote><blockquote type="cite">
<br></blockquote><blockquote type="cite">(Aside: I see that this is a feature request)<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I noticed that the Cox regression model is under active development in<br>
</blockquote><blockquote type="cite">ProbABEL version 0.2.0.  The purpose of this email is to check in with<br></blockquote><blockquote type="cite">the group regarding my plan:<br></blockquote><blockquote type="cite"><br>
</blockquote><blockquote type="cite">1) check out some latest developers version of the code<br></blockquote><blockquote type="cite">2) play around with the coxph code (if it isn't fixed)<br></blockquote><blockquote type="cite">
3) evaluate and test in a statistical fashion (I am a biostatistician<br></blockquote><blockquote type="cite">after all!)<br></blockquote><blockquote type="cite">4) implement robust standard error estimation<br></blockquote>
<blockquote type="cite">5) submit patch for review <br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Comments and ideas are welcome.<br></blockquote><blockquote type="cite"><br></blockquote>
<blockquote type="cite">Thanks!<br></blockquote><blockquote type="cite">Alisa<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">-- <br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">
Alisa Knodle Manning, Ph.D.<br></blockquote><blockquote type="cite">Postdoctoral Research Fellow<br></blockquote><blockquote type="cite">Broad Institute<br></blockquote><blockquote type="cite">617-714-7662<br></blockquote>
<blockquote type="cite"><a href="mailto:amanning@broadinstitute.org" target="_blank">amanning@broadinstitute.org</a> <<a href="mailto:amanning@broadinstitute.org" target="_blank">mailto:amanning@broadinstitute.org</a>><br>
</blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">_______________________________________________<br></blockquote>
<blockquote type="cite">genabel-devel mailing list<br></blockquote><blockquote type="cite"><a href="mailto:genabel-devel@lists.r-forge.r-project.org" target="_blank">genabel-devel@lists.r-forge.r-project.org</a><br></blockquote>
<blockquote type="cite"><a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/genabel-devel" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/genabel-devel</a><br></blockquote><blockquote type="cite">
<br></blockquote><br>-- <br>-----------------------------------------------<br>dr. L.C. Karssen<br>Erasmus MC<br>Department of Epidemiology<br>Room Ee-2224<br><br>Postbus 2040<br>3000 CA Rotterdam<br>The Netherlands<br><br>
phone: <a href="tel:%2B31-10-7044217" value="+31107044217" target="_blank">+31-10-7044217</a><br>fax: <a href="tel:%2B31-10-7044657" value="+31107044657" target="_blank">+31-10-7044657</a><br>e-mail: <a href="mailto:l.karssen@erasmusmc.nl" target="_blank">l.karssen@erasmusmc.nl</a><br>
GPG key ID: 0E1D39E3<br>-----------------------------------------------<br><br><br><br>_______________________________________________<br>genabel-devel mailing list<br><a href="mailto:genabel-devel@lists.r-forge.r-project.org" target="_blank">genabel-devel@lists.r-forge.r-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/genabel-devel" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/genabel-devel</a></div></blockquote></div></div></div><br></div>
</div></blockquote></div><br><br clear="all"><div><br></div>-- <br><br>Alisa Knodle Manning, Ph.D.<br>Postdoctoral Research Fellow<br>Broad Institute<br>617-714-7662<br><a href="mailto:amanning@broadinstitute.org" target="_blank">amanning@broadinstitute.org</a><br>
<br>
</div>