From fff1fc581da5b4ff935ebd4d7ded02f750472fdc Mon Sep 17 00:00:00 2001 From: jranke Date: Tue, 27 Mar 2012 01:03:18 +0000 Subject: Start of the transition to fitting transformed parameters. Many things are broken now (see TODO list) git-svn-id: svn+ssh://svn.r-forge.r-project.org/svnroot/kinfit/pkg/mkin@20 edb9625f-4e0d-4859-8d74-9fd3b1da38cb --- R/mkinstart.R | 23 ----------------------- 1 file changed, 23 deletions(-) delete mode 100644 R/mkinstart.R (limited to 'R/mkinstart.R') diff --git a/R/mkinstart.R b/R/mkinstart.R deleted file mode 100644 index 46e5bed..0000000 --- a/R/mkinstart.R +++ /dev/null @@ -1,23 +0,0 @@ -mkinstart <- function(model, data, mode = "auto") -{ - if (class(model) != "mkinmod") stop("The first argument must be a model of class mkinmod") - names <- model$parms - observed <- names(model$map) - if(!all(observed %in% levels(data$name))) stop("The data must contain the observed variables used in the model") - for (obs in observed) - { - tmp <- subset(data, name == obs) - max <- tmp[which.max(tmp$value), ] - type = names(model$map[[obs]])[[1]] - kinmodel <- ifelse(type == "SFORB", "DFOP", type) - tmp.longdata <- subset(data, name == obs & time >= max$time) - tmp.widedata <- mkin_long_to_wide(tmp.longdata, outtime = "t") - names(tmp.widedata) <- c("t", "parent") - tmp.fit <- kinfit( - kindata = tmp.widedata, - kinmodels = kinmodel, - parent.0.user = max$value) - if(class(tmp.fit[[kinmodel]]) == "try-error") stop(paste("Automatic generation of starting parameters failed\nkinfit failed to find a", kinmodel, "fit for", obs)) - tmp.results <- kinresults(tmp.fit) - } -} -- cgit v1.2.1