[Genabel-commits] r1089 - pkg/GenABEL/R

[noreply at r-forge.r-project.org]

Author: nd_0001
Date: 2013-01-31 10:48:48 +0100 (Thu, 31 Jan 2013)
New Revision: 1089

Modified:
   pkg/GenABEL/R/estlambda.R
Log:
minor improvement in estlambda

Modified: pkg/GenABEL/R/estlambda.R
===================================================================
--- pkg/GenABEL/R/estlambda.R	2013-01-27 20:32:44 UTC (rev 1088)
+++ pkg/GenABEL/R/estlambda.R	2013-01-31 09:48:48 UTC (rev 1089)
@@ -62,7 +62,7 @@
 	}
 	data <- sort(data)
 	ppoi <- ppoints(data)
-	ppoi <- sort(qchisq(1-ppoi,df))
+	ppoi <- sort(qchisq(ppoi,df=df,lower.tail=FALSE))
 	data <- data[1:ntp]
 	ppoi <- ppoi[1:ntp]
 #	s <- summary(lm(data~offset(ppoi)))$coeff