## Function translated using 'matlab.to.r()'
## Then manually adjusted to make work
## Author: Andrew Hooker
mf_all <- function(model_switch,xt_ind,x,a,bpop,d,sigma,docc,poped.db){
returnArgs <- switch(poped.db$settings$iFIMCalculationType+1,
mf3(model_switch,xt_ind,x,a,bpop,d,sigma,docc,poped.db), #Default (with no assumption that bpop and b are uncorrelated)
mf3(model_switch,xt_ind,x,a,bpop,d,sigma,docc,poped.db), #Reduced FIM
stop("Not yet implemented"), #Weighted models
stop("Not yet implemented"), #Loc models
mf3(model_switch,xt_ind,x,a,bpop,d,sigma,docc,poped.db), #Reduced FIM with derivative of SD sigma
mf3(model_switch,xt_ind,x,a,bpop,d,sigma,docc,poped.db), #FULL FIM parameterized with A,B,C matrices & derivative of variance
mf7(model_switch,xt_ind,x,a,bpop,d,sigma,docc,poped.db), #Calculate one model switch at a time, good for large matrices
mf3(model_switch,xt_ind,x,a,bpop,d,sigma,docc,poped.db)) #Reduced FIM parameterized with A,B,C matrices & derivative of variance
if(is.null(returnArgs)) stop(sprintf('Unknown FIM-calculation type'))
ret <- returnArgs[[1]]
poped.db <- returnArgs[[2]]
return(list( ret= ret,poped.db =poped.db ))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.