Get 10% of triobra1 and calculate REAP /projects/relatedness/plink-1.07-x86_64/plink –file trio1_bra.phasing –thin 0.2 –noweb –geno 0 –recode –out trio1_bra.phasing.20 /projects/relatedness/plink-1.07-x86_64/plink –file trio1_bra.phasing –thin 0.1 –noweb –geno 0 –recode –out trio1_bra.phasing.10 Of these, 956761 snps are new, 5637 already exist Cálculo só com os brasileiros merge3.ped
Month: February 2013
Get 10% of triobra1 and calculate REAP /projects/relatedness/plink-1.07-x86_64/plink –file trio1_bra.phasing –thin 0.2 –noweb –geno 0 –recode –out trio1_bra.phasing.20 /projects/relatedness/plink-1.07-x86_64/plink –file trio1_bra.phasing –thin 0.1 –noweb –geno 0 –recode –out trio1_bra.phasing.10 Of these, 956761 snps are new, 5637 already exist Cálculo só com os brasileiros merge3.ped
Report 210213
Vcf::throw(‘Vcf4_1=HASH(0x1ce37c8)’, ‘The command “tabix -l /projects/1000genomes/integrated_call_s…’) called at /lgc/programs/vcftools_0.1.10/lib/perl5/site_perl/Vcf.pm line 2630
VcfReader::get_chromosomes(‘Vcf4_1=HASH(0x1ce37c8)’) called at /lgc/programs/vcftools_0.1.10/bin/vcf-merge line 195
main::init_cols(‘HASH(0x18d5c18)’, ‘Vcf4_1=HASH(0x18dcd30)’) called at /lgc/programs/vcftools_0.1.10/bin/vcf-merge line 277
main::merge_vcf_files(‘HASH(0x18d5c18)’) called at /lgc/programs/vcftools_0.1.10/bin/vcf-merge line 12
Traceback (most recent call last):
File “/lgc/scripts/reap.py”, line 51, in