To implement, first copy all the files and save them under the same directory. The functions needed for implementing the variance component (VC) score test are stored in "hap.vc.fun.r". The main function is "vc.score.fun", and whose arguments are explained inside the function. Please read and run the commands in "example.r" line by line, and check if the same results (or similar result for resampling method) with what listed in the file could be obtained. =========================================================================== (1) example.r : the example file to show how to use the code posted here (2) hap.vc.fun.r: main functions for implementing VC-score test (3) mygeno: an example data matrix (genotypic data in Schaid et al. 2002 format; use "help(haplo.score)" to see the specific format). (4) myy.qt: an example vector of quantitative traits (5) myy.bt: an example vector of binary traits (6) myx.adj: an example matrix of non-genetic covariances =========================================================================== Please email me (jytzeng@stat.ncsu.edu) when running into any questions. Thanks for being interested in this method! Jung-Ying July 5, 2007