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:
- have 5 random starts
- pick best runs them
- 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
Post a Comment