swh:1:snp:ffdd0a7d2c8ea15ad41d45b3b178f668bd942287
Raw File
Tip revision: 65dce028095b06ee5ed19673e3045548876c41d3 authored by Derek Young on 30 December 2007, 00:00:00 UTC
version 0.3.3
Tip revision: 65dce02
flaremixEM.R
flaremixEM<-function(y, x, lambda = NULL, beta = NULL, sigma = NULL, alpha = NULL, nu=NULL,
		epsilon = 1e-04, maxit = 10000, verb = FALSE, restart=50){
j=1
while(j<=length(nu)){
temp=try(try.flare(y=y, x=x, lambda = lambda, beta = beta, sigma = sigma, alpha = alpha, nu=nu[j],
		epsilon = epsilon, maxit = maxit, verb = verb, restart=restart),silent=TRUE)

if(class(temp)=="try-error") j=j+1 else j=2^100
}



if(j==(length(nu)+1)) stop(paste("Too many attempts. Select a different barrier constant.","\n"))


temp

}

back to top