"cmlmomco" <- function(f,para) {
if(! are.par.valid(para)) return()
Qlo <- par2qua(0,para,paracheck=FALSE)
Mu <- vector(mode="numeric", length=length(f))
for(i in 1:length(f)) {
if(f[i] == 1) { Mu[i] <- par2qua(1,para,paracheck=FALSE); next }
tmp <- NULL
"afunc" <- function(p) {
return(par2qua(p,para,paracheck=FALSE))
}
try(tmp <- integrate(afunc, f[i], 1))
Mu[i] <- ifelse(is.null(tmp), NA, tmp$value/(1-f[i]))
}
return(Mu)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.