Variable selection for multilevel phylogenetic model
Open
Nobody has claimed this yet.
enhancement
- Dominant language
- R
- Stars
- 114
- Forks
- 31
- PR merge metrics
- No merged PRs in 30d
Description
Hi,
I am trying to use projpred for variable selection on a phylogenetic multilevel model fit with brms.
variable_response ~ x1 + x2 + (1|gr(species, cov = A))
However, I get the following error:
n model.frame.default(data = list(. variable_response.1 = c(0.580296874164707, :
invalid type (list) for variable 'gr(species, cov = A)
In addition: Warning message:
Some Pareto k diagnostic values are slightly high. See help('pareto-k-diagnostic') for details.
Thanks in advance,
Contributor guide
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start by reproducing the reported projpred workflow with the brms formula containing gr(species, cov = A), then inspect how that term is handled before the model.frame.default failure. Done means variable selection can run for this multilevel phylogenetic model without the reported invalid-type error; the Pareto k warning is also part of the reproduction context.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- r
- Domain
- data, machine-learning
- Issue type
- Bug
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 35/100