aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorPjotr Prins2015-03-17 10:39:58 +0300
committerPjotr Prins2015-03-17 10:39:58 +0300
commit7c958796142e9fe0cc3155d10581d68ae2bbc5df (patch)
tree4949fc749c45f509e83b2f30a325516f85b1823a
parent7268ed46c9bce6fea58b55723c325426bdeeb617 (diff)
downloadgenenetwork2-7c958796142e9fe0cc3155d10581d68ae2bbc5df.tar.gz
lmm.py: bring to original state
-rw-r--r--wqflask/wqflask/my_pylmm/pyLMM/lmm.py4
1 files changed, 0 insertions, 4 deletions
diff --git a/wqflask/wqflask/my_pylmm/pyLMM/lmm.py b/wqflask/wqflask/my_pylmm/pyLMM/lmm.py
index c42f9fb7..994b6be7 100644
--- a/wqflask/wqflask/my_pylmm/pyLMM/lmm.py
+++ b/wqflask/wqflask/my_pylmm/pyLMM/lmm.py
@@ -688,15 +688,11 @@ class LMM:
optimum.
"""
- print("H=",H)
- print("X=",X)
- print("REML=",REML)
n = len(self.LLs)
HOpt = []
for i in range(1,n-2):
if self.LLs[i-1] < self.LLs[i] and self.LLs[i] > self.LLs[i+1]:
HOpt.append(optimize.brent(self.LL_brent,args=(X,REML),brack=(H[i-1],H[i+1])))
- print("HOpt=",HOpt)
if np.isnan(HOpt[-1][0]):
HOpt[-1][0] = [self.LLs[i-1]]