### Using Model Averaging in R Mark

Posted:

**Tue Nov 22, 2022 5:33 pm**Hello all,

I am running CJS mark-resight models for an overwinter survival study I am doing on White-throated Sparrows. I am looking at any potential effects of Bird Age (juv vs adult), Sex, or Morph (tan-stripe vs white-stripe) on Phi/p within the species, and any interactive effects of the covariates on survival and resight probability. Total, I have three data sets from two different study sites and seasons. Because each data set has multiple equally competitive models (anywhere from 2 to 9), I was going to use model averaging to get my real estimates for Phi/p. This is where I run into two issues.

1) The model.averaging output from R Mark has 440 or 880 entries depending on whether I run Phi/p averaging separately or together.

2) The output table does not seem to give a clear indication of what real estimate is associated with which covariate or group. The most it gives me is something like "Phi g000 c1 a3 t4" if I run Phi and p separate. If I run them together, it doesn't give me any sort of label for the estimates.

So my first question, is there any way to cut down on the number of entries it gives me to make it more readable? I've tried cutting down my model.collect table to just include the competitive models, but that didn't make a difference.

Second question, is there a way to get the output labeled in a way that's easier to interpret what the effects of the covariates are? I am not having any issue with that if I look at the models individually. For example, one of my top models for one data set has an interactive effect of Bird Age and Morph on Phi. Looking at the outputs from the model, I can see that juvenile tan-stripes have the highest survival estimate, followed by adult white-stripes, adult tan-stripes, and juvenile white-stripes. But other competitive models include effects from Morph alone, Sex alone, and interactive effects of Sex and Morph. Once I do model averaging, it becomes unclear what the effects are and what the differences are between the covariates/groups, as well as what the real estimates are.

Am I just misunderstanding how model averaging is used or works within R Mark? Please let me know if you have any clarifying questions, and thank you for your help and time.

-Nate

I am running CJS mark-resight models for an overwinter survival study I am doing on White-throated Sparrows. I am looking at any potential effects of Bird Age (juv vs adult), Sex, or Morph (tan-stripe vs white-stripe) on Phi/p within the species, and any interactive effects of the covariates on survival and resight probability. Total, I have three data sets from two different study sites and seasons. Because each data set has multiple equally competitive models (anywhere from 2 to 9), I was going to use model averaging to get my real estimates for Phi/p. This is where I run into two issues.

1) The model.averaging output from R Mark has 440 or 880 entries depending on whether I run Phi/p averaging separately or together.

2) The output table does not seem to give a clear indication of what real estimate is associated with which covariate or group. The most it gives me is something like "Phi g000 c1 a3 t4" if I run Phi and p separate. If I run them together, it doesn't give me any sort of label for the estimates.

So my first question, is there any way to cut down on the number of entries it gives me to make it more readable? I've tried cutting down my model.collect table to just include the competitive models, but that didn't make a difference.

Second question, is there a way to get the output labeled in a way that's easier to interpret what the effects of the covariates are? I am not having any issue with that if I look at the models individually. For example, one of my top models for one data set has an interactive effect of Bird Age and Morph on Phi. Looking at the outputs from the model, I can see that juvenile tan-stripes have the highest survival estimate, followed by adult white-stripes, adult tan-stripes, and juvenile white-stripes. But other competitive models include effects from Morph alone, Sex alone, and interactive effects of Sex and Morph. Once I do model averaging, it becomes unclear what the effects are and what the differences are between the covariates/groups, as well as what the real estimates are.

Am I just misunderstanding how model averaging is used or works within R Mark? Please let me know if you have any clarifying questions, and thank you for your help and time.

-Nate