Data Cloning II

Project information

Project members
eric.ward e2holmes kwills

Sharing
Public


Description
Before you begin with this code, make sure you download the latest release of OpenBUGS(now 3.0.3) and it is installed to the default directory, "C:\Program Files\OpenBUGS". We've also found that it helps to have a cleaned up registry + defregmented hard drive before using OpenBUGS. The version of R that you're using should also be noted - R 2.5.1 and 2.7.0 work great with OpenBUGS, but there was a glitch in R 2.6.0 that prevented OpenBUGS from ever running through R. When using BUGS, I think that you can obviously use WinBUGS or OpenBUGS. While both programs should yield the same result, I had much more success with OpenBUGS, particularly in terms of the program not crashing during the burn-in (anyone familiar with 'trap' error messages knows what I'm talking about). There are 2 project files. The first (writeModel.r) is sourced by the second file (runDataCloning2.r). If you want to modify the bounds on parameters, you'll need to edit writeModel.r - otherwise, everything will be done behind the scenes. I've tried to predict many of the errors that might occur - for example, it doesn't make sense to have all sites belonging to the same group (m = 1) and an unconstrained covariance matrix. These files won't automate every possible MSSM that you want to fit, but they will do the majority of them. The R script files actually write the BUGS code for you, so you don't have to know any BUGS coding. You should however, be aware of the priors. Uniform priors are used on SDs and growth rates - to change these, look anywhere in the file for (~ dunif(a,b), where a and b are the limits). The code also includes the option of including an interaction matrix (multiple species) or density dependence matrix (multiple populations). The priors for all of these terms are bounded (0,1). Included are several examples of summarizing parameters, including the median, mean, and using density to estimate the mode. The actual MLEs are going to be the points associated with the best (lowest) density.


Files

Data Cloning II /
  • Download Data Cloning II as zip archive Download
Name Status Rev Size Updated Downloads Brief Description Actions
Item is folder old code fragments
Stable
- 3 items 10/28/2008 11:37 PDT 511 58 d
(none)
dat.csv
Stable
1 369 B 10/28/2008 11:37 PDT 142 142 f
Dummy data file
Data_Cloning_I_univariate.zip
Stable
1 1.2 KB 02/22/2009 17:30 PST 156 156 f
Zip file containing univariate code
runDataCloning2.r
Stable
1 5 KB 10/28/2008 11:37 PDT 160 160 f
This is the master file
writeModel.r
Stable
1 10.3 KB 10/28/2008 11:37 PDT 170 170 f
This is sourced; writes the BUGS file

Keywords



Downloads

100 project downloads
1139 individual file downloads

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