For genomic estimated breeding value(EBV) and Genome wide association studies(GWAS) using JWAS and XSim. The workflow includes
- Simulating SNP data using XSim. Files generated: pedigree.txt, phenotype.txt, genotype.txt, map.txt and ID.txt
- Statistical Analysis for EBV and GWAS using JWAS
- Includes XSimPreProcess package for automating XSim -> JWAS analysis.
- To use GWAS function for WPPA in JWAS package run
Pkg.add(PackageSpec(name=“JWAS”, rev=“master”))
and then import JWAS.misc byusing JWAS.misc
Test-Split.jl
splits your input files into training and testing sets for Single Step Analysis on JWAS
push!(LOAD_PATH, "/Path/To/My/Module/")
to add XSimPreProcess Module to julia. Requires this push for module use everytime julia is reopned- Copy/paste the command in the file
~/.julia/config/startup.jl
will extend LOAD_PATH on every Julia startup.
Usage of XSimPreProcess module for SSBR-analysis in this folder
ARGS = 5
ARG[1] = Input Genotype Path ARG[2] = Processed Data Path ARG[3] = Analysis Data Path ARG[4] = Number of QTL to set for analysis ARG[5] = Number of Markers to include for Phenotype Simulation
- Make 3 Directories - 1 for input data, 1 for processed data, 1 for analysis
- Files inside Input Data - Genotype file, XSim-Map-File
- Files inside Processed Data - Processed genotype, phenotype, pedigree, mapfile, rowIDs and other files as per pre-processing script
- Files inside Analysis Data - Output files of JWAS run
- Run command as follows
julia Phenotype-SSBR-Simulation.jl .\Directory1// .\Directory2// .\Directory3// 5 1