<div dir="ltr">🤓👍<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, 3 Sep 2019 at 09:51, Anja Thormann <<a href="mailto:anja@ebi.ac.uk">anja@ebi.ac.uk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;">Hi Ramiro,<div><br></div><div>thank you for reporting this. There is indeed a bug in the code which passes the window size parameter to the LD calculation script. The window size used at the moment is set to 200kb instead of the requested 500kb. This will be corrected for our next release which is planned for the second half of September. </div><div><br></div><div>Best regards,</div><div>Anja<br><div><br><blockquote type="cite"><div>On 2 Sep 2019, at 16:39, Ramiro Magno <<a href="mailto:ramiro.magno@gmail.com" target="_blank">ramiro.magno@gmail.com</a>> wrote:</div><br class="gmail-m_6989626417606793240Apple-interchange-newline"><div><div dir="ltr">Hi devs,<div><br></div><div>When looking for proxy snps of rs7707921, the snp rs17243791 does not show up:</div><div><br></div><div><a href="https://rest.ensembl.org/ld/homo_sapiens/rs7707921/1000GENOMES:phase_3:CEU?window_size=500;d_prime=0;r2=0.8" target="_blank">https://rest.ensembl.org/ld/homo_sapiens/rs7707921/1000GENOMES:phase_3:CEU?window_size=500;d_prime=0;r2=0.8</a><br></div><div><br></div><div>But if I ask for their linkage disequilibrium explicitly:</div><div><br></div><div><a href="https://rest.ensembl.org/ld/homo_sapiens/pairwise/rs7707921/rs17243791?population_name=1000GENOMES:phase_3:CEU;d_prime=0;r2=0" target="_blank">https://rest.ensembl.org/ld/homo_sapiens/pairwise/rs7707921/rs17243791?population_name=1000GENOMES:phase_3:CEU;d_prime=0;r2=0</a><br></div><div><br></div><div>I can see that they are indeed in LD with an r^2 of 0.920 (above the 0.8 threshold used above).</div><div><br></div><div>The two snps are 215,400 bp away from each other, which is below the 250k up/downstream threshold, i.e. window_size=500.</div><div><br></div><div>Any clue?</div><div><br></div><div>Many thanks in advance, Ramiro Magno</div></div>
_______________________________________________<br>Dev mailing list    <a href="mailto:Dev@ensembl.org" target="_blank">Dev@ensembl.org</a><br>Posting guidelines and subscribe/unsubscribe info: <a href="https://lists.ensembl.org/mailman/listinfo/dev_ensembl.org" target="_blank">https://lists.ensembl.org/mailman/listinfo/dev_ensembl.org</a><br>Ensembl Blog: <a href="http://www.ensembl.info/" target="_blank">http://www.ensembl.info/</a><br></div></blockquote></div><br></div></div>_______________________________________________<br>
Dev mailing list    <a href="mailto:Dev@ensembl.org" target="_blank">Dev@ensembl.org</a><br>
Posting guidelines and subscribe/unsubscribe info: <a href="https://lists.ensembl.org/mailman/listinfo/dev_ensembl.org" rel="noreferrer" target="_blank">https://lists.ensembl.org/mailman/listinfo/dev_ensembl.org</a><br>
Ensembl Blog: <a href="http://www.ensembl.info/" rel="noreferrer" target="_blank">http://www.ensembl.info/</a><br>
</blockquote></div>