-
Notifications
You must be signed in to change notification settings - Fork 0
/
Data preparation.R
22 lines (19 loc) · 1 KB
/
Data preparation.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
#R script
conda activate r-env
Library(rMVP)
library(rMVP)
MVP.Data(fileNum="numeric.txt",
filePhe="trait.txt",
fileMap="map.txt",
#sep.num="\t",
#sep.map="\t",
#sep.phe="\t",
fileKin=FALSE,
filePC=FALSE,
#priority="memory",
#maxLine=10000,
out="mvp"
)
#check your geno.map file, to have the header in caps - SNP CHROM and POS
#else it reads it as another marker and throws the error markers in genotype and map don't match, which is extremely weird if you ask me
#DO NOT name your files with numbers eg mvp50. The MVP command cannot find these files and gives an error of