User Tools

Site Tools


readme.thrgibbs1

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.thrgibbs1 [2013/07/17 16:23] shogoreadme.thrgibbs1 [2016/10/17 19:52] – [Options] shogo
Line 2: Line 2:
  
 =====Summary===== =====Summary=====
-Gibbs sampler for threshold-linear mixed models. The original program (THRGIBBSF90) was written by DeukHwan Lee in 2001 based on GIBBS2F90. Rewritten by Shogo Tsuruta in 2004. See PREGSF90 with genotypes (SNP) for options.+Gibbs sampler for threshold-linear mixed models. The original program (THRGIBBSF90) was written by DeukHwan Lee in 2001 based on GIBBS2F90 and formulas by Rob Tempelman. Rewritten by Shogo Tsuruta in 2004. See PREGSF90 with genotypes (SNP) for options.
 \\ \\
 THRGIBBS1F90 implements Gibbs sampler for mixed threshold-linear models involving multiple categorical and linear variables. Thresholds and variances can be estimated or assumed. Another version of thrgibbs1f90b for binary responses is available. THRGIBBS1F90 implements Gibbs sampler for mixed threshold-linear models involving multiple categorical and linear variables. Thresholds and variances can be estimated or assumed. Another version of thrgibbs1f90b for binary responses is available.
 +\\
 +See PREGSF90 with genotypes (SNP) for options.
  
 =====Parameters===== =====Parameters=====
Line 22: Line 24:
 OPTION fixed_var all 1 2 3 OPTION fixed_var all 1 2 3
 </file> </file>
-Store all samples for solutions in "all_solutions" and posterior means and SD for 1, 2, and 3 effects in "final_solutions".+Store all samples for solutions in "all_solutions" and posterior means and SD for 1, 2, and 3 effects in "final_solutions", assuming that (co)variances in the parameter file are known.
 <file> <file>
 OPTION fixed_var mean OPTION fixed_var mean
 </file> </file>
-Only posterior means and SD for solutions are calculated for all effects in "final_solutions".+Only posterior means and SD for solutions are calculated for all effects in "final_solutions", assuming that (co)variances in the parameter file are known.
 <file> <file>
 OPTION fixed_var mean 1 2 3 OPTION fixed_var mean 1 2 3
-Only posterior means and SD for solutions are calculated for effects 1, 2, and 3 in "final_solutions".+Only posterior means and SD for solutions are calculated for effects 1, 2, and 3 in "final_solutions", assuming that (co)variances in the parameter file are known.
 </file> </file>
 <file> <file>
 OPTION solution all OPTION solution all
 </file> </file>
-Store all samples for solutions in "all_solutions" and posterior means and SD for all effects.+Store all samples for solutions in "all_solutions" and posterior means and SD for all effects. The file "all_solutions" could be very large.
 <file> <file>
 OPTION solution all 1 2 3 OPTION solution all 1 2 3
 </file> </file>
-Store all samples for solutions in "all_solutions" and posterior means and SD for 1, 2, and 3 effects.+Store all samples for solutions in "all_solutions" and posterior means and SD for 1, 2, and 3 effects. The file "all_solutions" could be very large.
 <file> <file>
 OPTION solution mean OPTION solution mean
 </file> </file>
-Only posterior means and SD for solutions are calculated for all effects in "final_solutions".+Only posterior means and SD for solutions are calculated for all effects in "final_solutions" while sampling (co)variances.
 <file> <file>
 OPTION solution mean 1 2 3 OPTION solution mean 1 2 3
 </file> </file>
-Only posterior means and SD for solutions are calculated for effects 1, 2, and 3 in "final_solutions".+Only posterior means and SD for solutions are calculated for effects 1, 2, and 3 in "final_solutions" while sampling (co)variances.
 <file> <file>
 OPTION cont 10000 OPTION cont 10000
 </file> </file>
-"10000" is the number of samples run previously. The user can restart the program from the last run.+"10000" is the number of samples run previously. The user can restart the program from the last run. This option requires "binary_final_solutions", "gibbs_samples", and "fort.99" files.
 <file> <file>
 OPTION prior 5 2 -1 5  OPTION prior 5 2 -1 5 
Line 68: Line 70:
 Set the fixed thresholds. No need to set 0 for binary traits. Set the fixed thresholds. No need to set 0 for binary traits.
 <file> <file>
-OPTION residual 1.+OPTION residual 1
 </file> </file>
 Set the residual variance = 1. Set the residual variance = 1.
 <file> <file>
-OPTION censored xx.+OPTION pos_def
 </file> </file>
-Negative values of the last category in the data set indicate censored records. "xx" determines the lower and upper limit of the category + xx when sampling from the distribution.+Specify checking pos-def for fixed effects. 
 +<file> 
 +OPTION censored 1 0 
 +</file> 
 +Negative values for the categorical trait in the data set indicate censored records. "1 0" determines that the first categorical trait is censored and the second uncensored. 
 +<file>
 OPTION SNP_file snp OPTION SNP_file snp
 </file> </file>
 Specify the SNP file name to use genotype data. Specify the SNP file name to use genotype data.
- 
readme.thrgibbs1.txt · Last modified: 2024/03/25 18:22 by 127.0.0.1

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki