R/EMVSbest.R

Defines functions EMVSbest

Documented in EMVSbest

EMVSbest<-function(result){

  if(result$independent == F){
        posts=result$log_g_function
        posts[is.finite(posts)==F]=NaN
	which<-which.max(posts)
        logpost<-posts[which]
        print("Best Model Found")
        list(log_g_function=logpost,indices=(1:dim(result$betas)[2])[result$prob_inclusion[which,]>0.5])
  } else {
  print("Log posterior not available for independent prior case")
}
}

Try the EMVS package in your browser

Any scripts or data that you put into this service are public.

EMVS documentation built on Oct. 13, 2021, 5:09 p.m.