re: [BAD]<br><br>This is really strange -- the links work well for me! definitely we need to figure out what is the problem and how to solve that -- we need to share links to repository!<br><br>To start with -- Gulya, Lennart, et al., can you follow the links I have send? (see original mail below)<br>
<br>re: [GOOD]<br><br>These are interesting news! Can someone (Gulya?) run quick checks comparing the speed and accuracy of what Xia suggested versus standard polygenic?<br><br>Gulya: if you did not do this already, can you please subscribe to the genabel-devel list (<a href="http://lists.r-forge.r-project.org/mailman/listinfo/genabel-devel">http://lists.r-forge.r-project.org/mailman/listinfo/genabel-devel</a>) -- we want to keep the discussion open and documented, unless things discussed are too sensitive.<br>
<br>Yurii<br><br><div class="gmail_quote">On Wed, Feb 16, 2011 at 2:11 PM, Xia Shen <span dir="ltr"><<a href="mailto:xia.shen@lcb.uu.se">xia.shen@lcb.uu.se</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div style="word-wrap: break-word;">Hi Yurii,<div><br></div><div>Bad news and good news:</div><div><br></div><div>[BAD]</div><div>I cannot open the two links that you sent...</div><div><br></div><div>[GOOD]</div><div>One of my colleague happens to be playing with your 'polygenic' function, and has got stuck by some binary trait because the numerical algorithm didn't converge. However, treating the trait as normal works. To fit the same mixed model, I tried my 'hglm' package (published on the recent issue of The R Journal), it worked pretty well, and fitting binary, poisson, gamma, ... traits will all work.</div>
<div><br></div><div>The way to fit hglm using GenABEL results is straightforward:</div><div><br></div><div>1. create the IBS matrix G, has to be symmetric, just need to replicate the lower triangle of the IBS from GenABEL and overwrite the upper triangle</div>
<div>2. create L that satisfies LL' = G by</div><div> s <- svd(G)</div><div> L <- s$u%*%diag(sqrt(s$d))</div><div>3. call hglm like</div><div> hglm(y = phdata(...)$..., X = design matrix of fixed effects, Z = L, family = gaussian())</div>
<div> where NA has to be removed before calling, gaussian() can be changed to binomial(link = 'logit'), etc.</div><div><br>
The hglm function will return all the parameter estimates (fixed effects, variance components) and random effects predicts and residuals, etc. The algorithm is based on inter-connected GLMs (h-likelihood theory, the work by Lee & Nelder in the last decade), efficient and stable, without establishing the likelihood and calling numeric algorithm like FGLS.</div>
<div><br></div><div>I think the residuals from the mixed model is the most interesting stuff that will be passed into the next GWAS in GenABEL.</div><div><br></div><div>Hopefully this can serve some help :)</div><div><br>
</div><div>Cheers</div><div><br></div><div>Xia</div><div><br></div><div><br><div>
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: STHeiti; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre>
<font face="'Courier New'"><span style="font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre><font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre>
<font face="'Courier New'"><span style="white-space: normal;"><b><font size="3"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre>
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre>
<font face="'Courier New'"><b><font size="3"><span style="font-size: 11px;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre>
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre>
<font face="'Courier New'"><span style="white-space: normal;"><b><font size="3"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre>
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre>
<font face="'Courier New'"><b><font size="3"><span style="font-size: 11px;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre>
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre>
<font face="'Courier New'"><span style="white-space: normal;"><b><font size="3"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre>
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre>
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'" size="3"><span style="font-size: 11px; white-space: pre-wrap;"><b><span style="font-family: Monaco; font-weight: normal; white-space: normal; font-size: medium;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: STHeiti; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><div>
<span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre><font face="'Courier New'"><span style="font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre>
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre>
<font face="'Courier New'"><span style="white-space: normal;"><b><font size="3"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre>
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre>
<span style="font-size: 11px;"><span style="font-weight: normal; font-size: small;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre><span style="white-space: normal;"><span style="white-space: pre-wrap; font-size: small;"><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre>
<span style="white-space: normal;"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre>
<span style="white-space: normal;"><span style="white-space: pre-wrap; font-size: small;"><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre>
<span style="font-size: 11px;"><span style="font-weight: normal; font-size: small;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre><span style="white-space: normal;"><span style="white-space: pre-wrap; font-size: small;"><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre>
<span style="white-space: normal;"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre>
<span style="white-space: normal;"><span style="white-space: pre-wrap; font-size: small;"><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre>
<span style="font-size: 11px;"><span style="font-family: Monaco; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"></span><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap;"></span><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"></span><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"></span><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"></span><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"></span><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"></span><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"></span><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: STHeiti; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre>
</pre><span style="font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre></pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre>
</pre><span style="white-space: normal;"></span><span style="font-size: 11px; white-space: pre-wrap;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre>
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre>
</pre><span style="font-size: 11px;"></span><span style="font-family: Monaco; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"></span><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"></span><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"></span><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"></span><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"></span><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"></span><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: STHeiti; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"></span><div>
</div><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre></pre><span style="font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre>
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre>
</pre><span style="white-space: normal;"></span><span style="font-size: 11px; white-space: pre-wrap;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre>
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre>
</pre><span style="font-size: 11px;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre>
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre>
</pre><span style="white-space: normal;"></span><span style="font-size: 11px; white-space: pre-wrap;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre>
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre>
</pre><span style="font-size: 11px;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre>
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre>
</pre><span style="white-space: normal;"></span><span style="font-size: 11px; white-space: pre-wrap;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre>
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre>
</pre><span style="font-size: 11px;"> o \o/ _ o _| \ / |_ o_ \o/ o
/|\ | /\ _\o \o | o/ o/_ /\ | /|\
/ \ / \ | \ /) | ( \ /o\ / ) | (\ / | / \ / \
---------------------------------------------------------------</span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="font-size: 11px; white-space: pre-wrap;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre style="display: inline ! important;">
</pre><span style="font-size: 11px;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="font-size: 11px; white-space: pre-wrap;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre style="display: inline ! important;">
</pre><span style="font-size: 11px;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="font-size: 11px; white-space: pre-wrap;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre style="display: inline ! important;">
</pre><span style="font-size: 11px;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="font-size: 11px; white-space: pre-wrap;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre style="display: inline ! important;">
</pre><span style="font-size: 11px;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="font-size: 11px; white-space: pre-wrap;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre style="display: inline ! important;">
</pre><span style="font-size: 11px;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="font-size: 11px; white-space: pre-wrap;"></span><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre style="display: inline ! important;">
</pre><span style="font-size: 11px;">-------<br></span></span></pre></span></span></span></span></pre></span></span></span></span></pre></span></span></span></span></pre></span></span></span></pre></span></span></span></span></pre>
</span></span></span></span></pre></span></span></span></span></pre></span></span></span></pre></span></span></b></span></font></span></font></pre></span></b></span></span></font></b></span></font></pre></span></span></b></span></font></span></font></pre>
</span></b></span></font></pre></span></div></span></span></span></span></span></span></span></span></b></span></font></span></font></pre></span><font face="Monaco"><font face="'Courier New'" size="3"><b></b></font></font></span></b></span><b><font face="Monaco"><font face="'Courier New'" size="3"></font></font></b></font></span><font face="'Courier New'"><font face="Monaco"></font></font></font></pre>
</span></b></span><b><font face="Monaco"><font face="'Courier New'"></font></font></b></span><font face="Monaco"><font face="'Courier New'"></font></font></font></b></span><b><font size="3"><font face="Monaco"></font></font></b></font></pre>
</span><font face="'Courier New'"><b><font size="3"></font></b></font></span></b></span><b><font face="'Courier New'"><font size="3"></font></font></b><font face="'Courier New'"></font></font></span></font></pre>
</span></b></span><b><font face="Monaco"><font face="'Courier New'"></font></font></b></span><font face="Monaco"><font face="'Courier New'"></font></font></font></b></font></pre></span><font face="'Courier New'"><b><font size="3"></font></b></font></span></b></span><b><font face="'Courier New'"><font size="3"></font></font></b></font></span><font face="'Courier New'"><font face="'Courier New'"></font></font></font></pre>
</span></b></span><b><font face="Monaco"><font face="'Courier New'"></font></font></b></span><font face="Monaco"><font face="'Courier New'"></font></font></font></b></span><b><font size="3"><font face="Monaco"></font></font></b></font></pre>
</span><font face="'Courier New'"><b><font size="3"></font></b></font></span></b></span><b><font face="'Courier New'"><font size="3"></font></font></b><font face="'Courier New'"></font></font></span></font></pre>
</span></b></span><b><font face="Monaco"><font face="'Courier New'"></font></font></b></span><font face="Monaco"><font face="'Courier New'"></font></font></font></b></font></pre></span><font face="'Courier New'"><b><font size="3"></font></b></font></span></b></span><b><font face="'Courier New'"><font size="3"></font></font></b></font></span><font face="'Courier New'"><font face="'Courier New'"></font></font></font></pre>
</span></b></span><b><font face="Monaco"><font face="'Courier New'"></font></font></b></span><font face="Monaco"><font face="'Courier New'"></font></font></font></b></span><b><font size="3"><font face="Monaco"></font></font></b></font></pre>
</span><font face="'Courier New'"><b><font size="3"></font></b></font></span><font face="'Courier New'"><font size="3"> Xia Shen - PhD student in Bioinformatics, Uppsala University.</font></font></b></span><b><font face="'Courier New'"><font size="3"></font></font></b></font></span></font></pre>
</span><span><pre><font face="Monaco"><font face="'Courier New'"><b><font face="'Courier New'"><b><font face="'Courier New'" size="3"><span style="font-size: 11px;"> - Computational Genetics Group</span></font></b><span style="font-family: 'Courier New'; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<font face="'Courier New'"><span style="white-space: normal;"><b><font size="3"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span><pre style="display: inline ! important;">
<b><font face="'Courier New'" size="3"><span style="font-size: 11px;">: <a href="http://www.computationalgenetics.se" target="_blank">www.computationalgenetics.se</a></span></font></b></pre></span></span></b></span></font></span></font></pre>
</span></b></span></span></font></b></span></font></pre></span></span></b></span></font></span></font></pre></span></b></span></font></b></font></font></pre></span><span><pre><font face="Monaco"><font face="'Courier New'"><b><font face="'Courier New'"><b><font face="'Courier New'" size="3"><span style="font-size: 11px;"></span></font></b><span style="font-family: 'Courier New'; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<font face="'Courier New'"><span style="white-space: normal;"><b><font size="3"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span><pre style="display: inline ! important;">
<b><font face="'Courier New'" size="3"><span style="font-size: 11px;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<font face="'Courier New'"><span style="white-space: normal;"><b><font size="3"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span><pre style="display: inline ! important;">
<b><font face="'Courier New'" size="3"><span style="font-size: 11px;"> - Personal URL: </span></font></b></pre></span></span></b></span></font></span></font></pre></span></b></span></span></font></b></span></font></pre>
</span></span></b></span></font></span></font></pre></span></b></span></span></font></b></pre></span></span></b></span></font></span></font></pre></span></b></span></span></font></b></span></font></pre></span></span></b></span></font></span></font></pre>
</span></b></span><span style="font-family: 'Courier New'; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<font face="'Courier New'"><span style="white-space: normal;"><b><font size="3"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span><pre style="display: inline ! important;">
<b><font face="'Courier New'" size="3"><span style="font-size: 11px;"><a href="http://www.19850911.com" target="_blank">www.19850911.com</a></span></font></b></pre></span></span></b></span></font></span></font></pre>
</span></b></span></span></font></b></span></font></pre></span></span></b></span></font></span></font></pre></span></b></span></font></b></font></font></pre><pre><font face="Monaco"><font face="'Courier New'"><b><font face="'Courier New'"><b><font face="'Courier New'" size="3"><span style="font-size: 11px;"><br>
</span></font></b></font></b></font></font></pre></span></b></span><b><font face="Monaco"><font face="'Courier New'"><b><font face="'Courier New'"></font></b><font face="'Courier New'"><b><font size="3"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre>
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: STHeiti; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span><pre>
<b><span style="white-space: normal;"><pre><b><font face="'Courier New'" size="3"><span style="font-size: 11px;"> All judgements are, in their rationale, STATISTICS. - C.R.Rao</span></font></b></pre></span></b></pre>
</span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre><span style="font-family: 'Courier New',Courier,mono; font-size: small;"><b><font face="'Courier New'"><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre>
<b><font face="'Courier New'" size="3"><span style="font-size: 11px;">---------------------------------------------------------------</span></font></b><span style="font-family: 'Courier New'; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<font face="'Courier New'"><span style="white-space: normal;"><b><font size="3"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<font face="'Courier New'"><b><font size="3"><span style="font-size: 11px;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<font face="'Courier New'"><span style="white-space: normal;"><b><font size="3"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<font face="'Courier New'"><b><font size="3"><span style="font-size: 11px;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<font face="'Courier New'"><span style="white-space: normal;"><b><font size="3"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><b><font size="3"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Monaco; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="font-family: 'Courier New'; font-weight: bold; white-space: pre-wrap; font-size: small;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: STHeiti; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<span style="font-family: 'Courier New'; font-size: small;"><b><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<font face="Monaco"><span style="white-space: normal;"><font face="'Courier New'"><span style="white-space: pre-wrap; font-size: small;"><b><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<font face="'Courier New'"><span style="white-space: normal;"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<span style="white-space: normal;"><span style="white-space: pre-wrap; font-size: small;"><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<span style="font-size: 11px;"><span style="font-weight: normal; font-size: small;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<span style="white-space: normal;"><span style="white-space: pre-wrap; font-size: small;"><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<span style="white-space: normal;"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<span style="white-space: normal;"><span style="white-space: pre-wrap; font-size: small;"><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<span style="font-size: 11px;"><span style="font-weight: normal; font-size: small;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"><pre style="display: inline ! important;">
<span style="white-space: normal;"><span style="white-space: pre-wrap; font-size: small;"><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"><pre style="display: inline ! important;">
<span style="white-space: normal;"><span style="font-size: 11px; white-space: pre-wrap;"><span style="font-weight: normal; font-size: small;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif; font-weight: normal; white-space: normal; font-size: medium;"></span><pre style="display: inline ! important;">
</pre><span style="white-space: normal;"></span><span style="white-space: pre-wrap; font-size: small;"></span><span style="font-family: STHeiti; font-weight: normal; white-space: normal; font-size: medium;"></span><span style="font-family: Verdana,Arial,Helvetica,sans-serif;"></span><pre style="display: inline ! important;">
</pre><span style="font-size: 11px;">-------</span></span></span></pre></span></span></span></span></pre></span></span></span></pre></span></span></span></span></pre></span></span></span></span></pre></span></span></span></span></pre>
</span></span></span></pre></span></span></span></span></pre></span></span></span></span></font></pre></span></span></b></span></font></span></font></pre></span></b></span></pre></span></span></b></span></font></span></font></pre>
</span></span></span></font></b></font></span></font></pre></span></b></span><b><font face="Monaco"><font face="'Courier New'"><font size="3"></font></font></font></b></span><font face="Monaco"><font face="'Courier New'"><font size="3"></font></font></font></font></b></span><b><font size="3"><font face="Monaco"><font face="'Courier New'"></font></font></font></b></font></pre>
</span></span></b></span><b><font face="'Courier New'"><b><font size="3"><font face="Monaco"></font></font></b><font size="3"></font></font></b></font></span><font face="'Courier New'"><b><font face="'Courier New'"></font></b></font></font></pre>
</span><font face="Monaco"><font face="'Courier New'"><b></b></font></font></b></span></span><b><font face="Monaco"><font face="'Courier New'"></font></font></b><font face="Monaco"></font></font></b></font></pre>
</span></span><font face="'Courier New'"><b><font size="3"></font></b></font></b></span><b><font face="'Courier New'"><font size="3"></font></font></b></font></span><font face="'Courier New'"><b><font face="'Courier New'"></font></b></font></font></pre>
</span><font face="Monaco"><font face="'Courier New'"><b></b></font></font></b></span></span><b><font face="Monaco"><font face="'Courier New'"></font></font></b><font face="Monaco"></font></font></b></span></font></pre>
</span></span><font face="'Courier New'"><b><font size="3"></font></b></font></b></span><b><font face="'Courier New'"><font size="3"></font></font></b></font></span><font face="'Courier New'"><b><font face="'Courier New'"></font></b></font></font></pre>
</span><font face="Monaco"><font face="'Courier New'"><b></b></font></font></b></span></span><b><font face="Monaco"><font face="'Courier New'"></font></font></b><font face="Monaco"></font></font></b></font></pre>
</span></span><font face="'Courier New'"><b><font size="3"></font></b></font></b></span><b><font face="'Courier New'"><font size="3"></font></font></b></font></span><font face="'Courier New'"><b><font face="'Courier New'"></font></b></font></font></pre>
</span><font face="Monaco"><font face="'Courier New'"><b></b></font></font></b></span></span><b><font face="Monaco"><font face="'Courier New'"></font></font></b><font face="Monaco"></font></font></b></span></font></pre>
</span></span><font face="'Courier New'"><b><font size="3"></font></b></font></b></span><b><font face="'Courier New'"><font size="3"></font></font></b></font></span><font face="'Courier New'"><b><font face="'Courier New'"></font></b></font></font></pre>
</span><font face="Monaco"><font face="'Courier New'"><b></b></font></font></b></span><b><font face="Monaco"><font face="'Courier New'"></font></font></b><font face="Monaco"></font></pre></span></font></b></span><b><font face="'Courier New'"></font></b></pre>
</span><b></b></span></b></span></font></span></font></pre></span><font face="Monaco"><font face="'Courier New'"><b></b></font><b></b></font></span></span></font></b></font></font></font></b></font></pre></span></span><font face="'Courier New'"><b><font face="Monaco"><font face="'Courier New'"><font face="'Courier New'"><b><font size="3"><font face="Monaco"></font></font></b></font></font></font></b></font></span><font face="'Courier New'"><b><font face="Monaco"><font face="'Courier New'"><font face="'Courier New'"><font size="3"><font face="Monaco"></font></font></font></font></font></b></font></span><font face="'Courier New'"><b><font face="Monaco"><font face="'Courier New'"><font face="'Courier New'"><font size="3"></font></font></font></font></b></font></span></span></span><font face="'Courier New'"><b><font face="Monaco"><font face="'Courier New'"><font face="'Courier New'"></font></font></font></b><font face="Monaco"><font face="'Courier New'"></font></font></font></span><font face="'Courier New'"><font face="Monaco"></font></font></span><font face="'Courier New'"></font></span></span><font face="'Courier New'"></font></div>
<br></div><span><img src="cid:D01E5A5D-7F6F-465D-8A79-DB500ABD42CD@hgen.slu.se" width="120" height="120"></span>
</div><div><div></div><div class="h5">
<br><div><div>On Feb 10, 2011, at 3:35 PM, Yurii Aulchenko wrote:</div><br><blockquote type="cite"><div>Hi Gulya,<br><br>As I mentioned, I am also very bothered by the fact that 'polygenic'<br>works so slowly. I cc to Xia who mentioned he had some idea on how to<br>
make it working faster -- Xia, would be great if you could contribute<br>to the discussion! I also cc to GA devel list, so people know what is<br>going on.<br><br>Here I will describe how polygenic works currently (very far from optimal!)<br>
<br>'polygenic'<br><br><a href="https://r-forge.r-project.org/scm/viewvc.php/pkg/GenABEL/R/polygenic.R?view=markup&revision=623&root=genabel" target="_blank">https://r-forge.r-project.org/scm/viewvc.php/pkg/GenABEL/R/polygenic.R?view=markup&revision=623&root=genabel</a><br>
<br>is numerically maximizing likelihood function computed by procedure 'polylik'<br><br><a href="https://r-forge.r-project.org/scm/viewvc.php/pkg/GenABEL/R/polylik.R?view=markup&revision=545&root=genabel" target="_blank">https://r-forge.r-project.org/scm/viewvc.php/pkg/GenABEL/R/polylik.R?view=markup&revision=545&root=genabel</a><br>
<br>For maximization, standard R procedures 'nlm' or 'optim' can be used<br>-- see what happens in 'polygenic' code between lines 246-300.<br><br>Any ideas on how we can speed it up are very welcome!<br>
<br>Yurii<br></div></blockquote></div><br></div></div></blockquote></div><br>