r - how to start glmulti with consensus -


given:

  • use of r/rstudio , library "glmulti"
  • we use "set seed" of 123465
  • we use noam ross data , tutorial glmulti (hyperlink)

except make 5 "sub-runs" use code:

fish.model_n <- glmulti(global.model,         level = 2,  #  @ main effects         method = "g",         crit="aicc",         report = true,         plotty = false,         popsize = 10,         mutrate = 0.025,         sexrate = 0.01,         imm = 0.33,         confsetsize = 3,         deltam = 10,         deltab = 0) 

and instead of "n" have values 1...5 indicating set was.

how combine them , use them more glmulti?

i want this:

  1. have 5 random starts
  2. pick best runs them
  3. use starting position larger run

thoughts:

  • it looks "consensus" finds better models list.
  • if put output of consensus input, 'glmulti' errors out saying "improper call of glmulti".

clarification (as requested):

  • consensus can take several previous runs , find "best" those.
  • a genetic algorithm builds on itself. "best" of previous iterations "mutated" or "cross-mated" create candidates evaluation on next iteration.
  • genetic models can sensitive initial conditions, , having multiple start points searches can speed search.
  • there doesn't seem clean way "inject" population, consensus or way. means can't take 10 small runs different random ic , mutation paths , find best, , use "good seed" start larger run.

it might add quite bit of extensibility 'glmulti' if "new seeds" fed "live gene pool" evaluation runs , if "current best" popped out without stopping run (much).


Comments

Popular posts from this blog

How to understand 2 main() functions after using uftrace to profile the C++ program? -

c# - Update a combobox from a presenter (MVP) -

How to put a lock and transaction on table using spring 4 or above using jdbcTemplate and annotations like @Transactional? -