Hierarchical linear model of MPA effects

Project information

Project members
brice.semmens

Sharing
Public


Description
These are a series of Matlab routines to carry out a hierarchical linear model of fish population trajectories that includes a linear adjustment for a "reserve effect" at sites designated as MPAs. %The model is set up as follows: %FIRST LEVEL------------------------------------------------------------------------------------------------------- %MODEL:***** y~ intercept + abundance_index * (S_g + R_t +R_nt) + error *** %where: %S_g is the slope for a given species in a given geographic region (drawn %from a normally distributed species specific hyperparameter) %R_t is a species specific additive adjustment to the slope that is drawn %from a normally distributed hyper parameter representing the reserve effect on harvested species. Activated by a dummy variable. %R_nt is the species specific additive adjustment to the slope that is drawn from a normally distributed hyper parameter %representing the reserve effect on non-harvested species. Activated by a dummy variable. %both intercept and error terms are nuisance parameters and do not need hyperparameters. I suppose you could put in an %error hyper parameter with an informative prior just to keep the model in check. %SECOND LEVEL (hyper parameters all normally distributed)----------------------------------------------------------- %S_g_hat -- These are slopes for a given species from which S_g's are drawn %R_t_hat -- This is the reserve effect for targeted species from which R_t's across targeted species are drawn %R_nt_hat -- This is the reserve effect for non-targed species from which %R_nt's across non-targeted species are drawn %--------------------------------------------------------------------------------------------------------------------------


Files

Hierarchical linear model of MPA effects /
  • Download Hierarchical linear model of MPA effects as zip archive Download
Name Status Rev Size Updated Downloads Brief Description Actions
Item is folder data_generation
Stable
- 0 items 01/03/2008 12:43 PST 6 6 d
(none)
Item is folder single hyper-reserve effect
Stable
- 2 items 02/05/2008 15:16 PST 89 13 d
(none)
hierarch_like3.m
Caution
1 2.8 KB 01/03/2008 12:41 PST 40 40 f
(none)
hierarch_linear2.m
Caution
1 9.6 KB 01/03/2008 12:40 PST 39 39 f
(none)
hierarch_linear_selected.m
Caution
1 9.6 KB 01/03/2008 12:42 PST 28 28 f
(none)
output_allspps.txt
Stable
1 32.2 KB 01/03/2008 12:43 PST 49 49 f
(none)
output_selected_species.txt
Stable
1 15.8 KB 01/03/2008 12:43 PST 75 75 f
(none)

Keywords



Downloads

43 project downloads
326 individual file downloads

Sculpin 0.2 | xhtml | problems or comments? | report bugs