Skip to content

Model notation in methods draft #39

@christophe-rd

Description

@christophe-rd

@lizzieinvancouver as you requested, starting an issue with the equation text:

$$\ln(\text{ringwidth})_i \sim \text{Normal}(\mu, \sigma_y)$$

$$ \mu = \alpha + \alpha_{\text{species}} + \alpha_{\text{site}} + \alpha_{\text{treeid}} + \alpha_{\text{year}} + \beta_{\text{species}} \times \text{season predictor} $$ $$ \alpha \sim \text{Normal}(1, 3), \quad $$ $$\alpha_{\text{species}} \sim \text{Normal}(0, 6), \quad $$ $$\alpha_{\text{site}} \sim \text{Normal}(0, \sigma_{\alpha,\text{site}})$$ $$\alpha_{\text{treeid}} \sim \text{Normal}(0, \sigma_{\alpha,\text{treeid}}), \quad $$ $$\beta_{\text{species}} \sim \text{Normal}(0, 0.5) $$ $$ \sigma_{\alpha,\text{site}} \sim \text{Normal}(0, 1), \quad$$ $$\sigma_{\alpha,\text{treeid}} \sim \text{Normal}(0, 1), \quad $$ $$\sigma_y \sim \text{Normal}+(0, 3)$$

I wasn't sure if you wanted me to also add the text that describes the equation. In case you did, here it is (I didn't include the model diagnostics and package details):

where $i$ is each observation of the dataset. Intercepts on species ($\alpha_{species}$) and year ($\alpha_{year}$) and slope on species ($\beta_{species}$) were not pooled (fixed effects), and we partially pooled (random effects) site ($\alpha_{site}$) and treeid ($\alpha_{treeid}$) to account for the similarity of the measurements within these groups.

Also, you pointed to $\sigma_{y}$ with "missing its normal". I'm not sure what you mean by that...

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions