From 16f5b1d3c0136413e92b2be0f20d365e92e9cd1c Mon Sep 17 00:00:00 2001
From: jranke <jranke@edb9625f-4e0d-4859-8d74-9fd3b1da38cb>
Date: Tue, 18 May 2010 12:58:38 +0000
Subject: Much more complete version that was just submitted to CRAN.

git-svn-id: svn+ssh://svn.r-forge.r-project.org/svnroot/kinfit/pkg/mkin@9 edb9625f-4e0d-4859-8d74-9fd3b1da38cb
---
 R/mkinerrmin.R | 14 ++++++++++++++
 1 file changed, 14 insertions(+)
 create mode 100644 R/mkinerrmin.R

(limited to 'R/mkinerrmin.R')

diff --git a/R/mkinerrmin.R b/R/mkinerrmin.R
new file mode 100644
index 00000000..7a922eb4
--- /dev/null
+++ b/R/mkinerrmin.R
@@ -0,0 +1,14 @@
+mkinerrmin <- function(errdata, n.parms, alpha = 0.05)
+{
+  means.mean <- mean(errdata$value_mean, na.rm=TRUE)
+
+	df = length(errdata$value_mean) - n.parms
+  
+	f <- function(err)
+	{
+		(sum((errdata$value_mean - errdata$value_pred)^2/((err * means.mean)^2)) - 
+		 qchisq(1 - alpha,df))^2
+	}
+	err.min <- optimize(f, c(0.01,0.9))$minimum
+	return(list(err.min = err.min, n.optim = n.parms, df = df))
+}
-- 
cgit v1.2.1