User Tools

Site Tools


readme.gibbs2

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
Next revisionBoth sides next revision
readme.gibbs2 [2012/08/17 15:06] – [Running] ignacioreadme.gibbs2 [2016/09/30 19:29] – [Running] andres
Line 1: Line 1:
 ====== GIBBS2F90 ====== ====== GIBBS2F90 ======
-\\ 
-Ignacy Misztal, University of Georgia, 3/7/2001 
  
 ===== Summary ===== ===== Summary =====
 +Ignacy Misztal, University of Georgia, 3/7/2001
 +\\
 Gibbs2f90 is a modification of gibbs1f90 to sample correlated effects jointly. This modification results in better mixing for models with correlated effects such as maternal or random-regression models. Memory requirements and CPU time per round are somewhat higher than in gibbs1f90. Gibbs2f90 is a modification of gibbs1f90 to sample correlated effects jointly. This modification results in better mixing for models with correlated effects such as maternal or random-regression models. Memory requirements and CPU time per round are somewhat higher than in gibbs1f90.
 +\\
 +See PREGSF90 with genotypes (SNP) for options.
  
 ===== Running ===== ===== Running =====
 +You can run as 
 +
 +<file>gibbs2f90 parameter_file --rounds 1000000 --burnin 10000 - -thin 10 --thinprint 10 </file>
 +
 +in which case you give the number of rounds, the burn-in, the storage of samples in the file ''gibbs_samples'', and the interval between printouts in the screen. Otherwise you can simply type ''gibbs1f90'' and you will be asked the parameter file and the questions below:
 +
 +
 number of samples and length of burn-in? number of samples and length of burn-in?
  
Line 40: Line 49:
 OPTION prior 5 2 -1 5  OPTION prior 5 2 -1 5 
 </file> </file>
-2 is the degree of belief for the 5th effect, and 5 is the degree of belief for the residual.\\ 
 The (co)variance  priors are specified in the parameter file.\\ The (co)variance  priors are specified in the parameter file.\\
 Degree of belief for all random effects should be specified using the following structure:\\ Degree of belief for all random effects should be specified using the following structure:\\
-OPTION prior eff1 db1    eff2 db2 ...    effn dbn    -1 dbres'\\ +OPTION prior eff1 db1    eff2 db2 ...    effn dbn    -1 dbres\\ 
-effx correspond to the effect number and dbx to the degree of belief for this random effect.\\ +effx correspond to the effect number and dbx to the degree of belief for this random effect-1 corresponds to the degree of belief of the residual variance.\\ 
--1 correspond to the degree of belief of the residual variance.+In this example 2 is the degree of belief for the 5th effect, and 5 is the degree of belief for the residual.\\
 <file> <file>
 OPTION seed 123 321 OPTION seed 123 321
readme.gibbs2.txt · Last modified: 2024/03/25 18:22 by 127.0.0.1

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki