File tree 2 files changed +6
-6
lines changed
2 files changed +6
-6
lines changed Original file line number Diff line number Diff line change 1
1
# !/usr/bin/perl -w
2
- # @(#)pique-input.pl 2021-05-26 A.J.Travis and A.Douglas
2
+ # @(#)pique-input.pl 2022-01-27 A.J.Travis and A.Douglas
3
3
4
4
#
5
5
# PIQUE - Parallel Identification Of QTL's using EMMAX
@@ -386,7 +386,7 @@ sub recode {
386
386
# convert .vcf to .ped + .map
387
387
if ( $type eq " vcf" ) {
388
388
must_exist(" $input .vcf" );
389
- plink(" --recode -- vcf $input .vcf --out $input " );
389
+ plink(" --vcf $input .vcf --recode ped --out $input " );
390
390
}
391
391
392
392
# convert .tped + .tfam to .ped + .map
@@ -399,7 +399,7 @@ sub recode {
399
399
# recode .ped + .map
400
400
must_exist(" $input .ped" );
401
401
must_exist(" $input .map" );
402
- plink(" --recode12 --file $opt_i --out $output " );
402
+ plink(" --recode 12 --file $opt_i --out $output " );
403
403
404
404
# check output files
405
405
must_exist(" $output .ped" );
@@ -919,7 +919,7 @@ sub transpose {
919
919
if ($opt_v ) {
920
920
print " Transpose $prefix ...\n " ;
921
921
}
922
- plink(" --recode -- transpose --file $prefix --out $prefix " );
922
+ plink(" --recode transpose --file $prefix --out $prefix " );
923
923
return ;
924
924
}
925
925
Original file line number Diff line number Diff line change 1
- # @(#)Makefile 2021-05-26 A.J.Travis and A.Douglas
1
+ # @(#)Makefile 2022-01-27 A.J.Travis and A.Douglas
2
2
3
3
#
4
4
# PIQUE - Parallel Identification Of QTL's using EMMAX
@@ -70,7 +70,7 @@ test_vcf: test_vcf.vcf
70
70
$(PIQUE-RUN ) -i $@ _GWAS -k IBS -c $@ _GWAS.covar -n $(THREADS )
71
71
72
72
test_vcf.vcf : $(IPREFIX ) .ped $(IPREFIX ) .map
73
- $(PLINK ) --file $(IPREFIX ) --recode- vcf --out $(basename $@ )
73
+ $(PLINK ) --file $(IPREFIX ) --recode vcf --out $(basename $@ )
74
74
75
75
ld : input
76
76
$(PLINK ) --file $(IPREFIX ) _recode12 --r --out $(OPREFIX ) _r
You can’t perform that action at this time.
0 commit comments