Entering edit mode
5.0 years ago
chiti.agarwal
▴
10
Hi
I am running GAPIT in R for GWAS. When I use external kinship and CV file the program shows error
Error in if ((n_gd + n_cv) > n_ran) { : argument is of length zero
In addition: There were 19 warnings (use warnings() to see them)
If i dont use these files it works. But I want to use these external files.
The code i used for GWAS
myY <- read.table("pythiumread_gwas.txt", head = TRUE)
myG <- read.table("Ca1-8_429_20042017_filtered_chiti_maf_0.05_missing_0.8_noCA_nocontig_noscaf_noLD.hmp.txt" , head = FALSE)
myCV <- read.table("mdp_population_structure.txt", head = TRUE)
myKI <- read.table("KI.txt", head = FALSE)
myGAPIT <- GAPIT(
Y=myY,
G=myG,
CV=myCV,
KI=myKI,
model=c("GLM","MLM","MLMM","FarmCPU")
)
I would be grateful if you can tell me why I am getting this error.
Hi,
You should have just edited your previous post. Please do not open multiple posts for the same question.