Differences

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

--- readme.thrgibbs1 [2012/05/29 13:19] – shogo
+++ readme.thrgibbs1 [2012/05/29 17:07] – shogo
@@ Line 1: / Line 1: @@
 ====== THRGIBBS1F90 ======
 \\
-Gibbs sampler for threshold-linear mixed models
+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.
 =====Summary=====
-Based on THRGIBBSF90 written by DeukHwan Lee at the University of Georgia, 2001 based on GIBBS2F90. Rewritten by Shogo Tsuruta, 2004. Last modified 09/05/2007
-=====Functionality=====
 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.
@@ Line 55: / Line 50: @@
 OPTION cont 10000
 </file>
-"10000" is the number of samples run previously.
+"10000" is the number of samples run previously. The user can restart the program from the last run.
-The user can restart the program from the last run.
 <file>
 OPTION prior 5
@@ Line 68: / Line 62: @@
 OPTION thresholds 0.0 1.0 2.0
 </file>
-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>
-OPTION residual 1
+OPTION residual 1.
 </file>
-Set the residual variance = 1
+Set the residual variance = 1.
 <file>
-OPTION censored xx
+OPTION censored xx.
 </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.