diff options
| author | jranke <jranke@edb9625f-4e0d-4859-8d74-9fd3b1da38cb> | 2010-05-18 12:58:38 +0000 |
|---|---|---|
| committer | jranke <jranke@edb9625f-4e0d-4859-8d74-9fd3b1da38cb> | 2010-05-18 12:58:38 +0000 |
| commit | 16f5b1d3c0136413e92b2be0f20d365e92e9cd1c (patch) | |
| tree | edab7b27ba4253b9dd45520e504c54851f65f26f /R/mkinerrmin.R | |
| parent | 30cbb4092f6d2d3beff5800603374a0d009ad770 (diff) | |
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
Diffstat (limited to 'R/mkinerrmin.R')
| -rw-r--r-- | R/mkinerrmin.R | 14 |
1 files changed, 14 insertions, 0 deletions
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))
+}
|