search for: optimout

...ont)>0){ if(length(possAlleles)!=K+1) stop("invalid length for possible Allele Frequency values") if(is.null(names(possAlleles))) names(possAlleles)<-paste("AlleleFreq",1:length(possAlleles),sep="") return(list(pi=rep(NA,length=K+1),alleleSet=possAlleles,optimOut=NA,nBelowCutoff=nBelowCutoff,nZero=nZero,nMut=length(m))) } else return(list(pi=rep(NA,length=K+1),alleleSet=NA,optimOut=NA,nBelowCutoff=nBelowCutoff,nZero=nZero,nMut=length(m))) } ###Make 'weight matrix' of P(X=X[i]|X>0, P=p(k,s)) of size nStages x nCopies for each x -- ###an...

Dear Friends, I'm using the optim command to maximize a likelihood function. My optim command is as follows estim.out <- optim(beta, loglike, X=Xmain, Y=Y, hessian=T, method="BFGS", control=c(fnscale=-1, trace=1, REPORT=1)) Setting the report=1, gives me the likelihood function value (if i'm correct) at each step. The output from running this is as follows initial value