Well you can just create your own, fake ones. So, if your genotyped animals are in file marker.geno.clean with this aspect
64000670990546 12012020210211
45214790003 12111120211101
45214790004 22111020110102
you can simply create a random number for each animal, and a file with "animal" and overall mean, using awk:
awk '{print 1,$1,rand()}' marker.geno.clean > phenotypes.txt
and then use the basic renumf90 parameter file in https://artadia.blogspot.com/2019/03/basic-renumf90-parameter-file.html.
After that, you edit the renf90.par file to put the required OPTIONS. For instance, these are the options to get lots of information, check heritability of gene content
OPTION msg 100
OPTION h2_gene_content
No comments:
Post a Comment