Question: From Illumina Final report to PLINK input
4
gravatar for lisesandenbergh
5.0 years ago by
South Africa
lisesandenbergh40 wrote:
Received back the Final report (and accompanying files) for Illumina (ovine) SNP chip results and would now like to analyse the data in PLINK. In the past I used the PLINK plug-in for Genomestudio to create a .ped and .map files for PLINK. Unfortunately I no longer have access to Genomestudio, so I have been looking into using the .lgen format. I have been struggling to get the results from the Final report into the .lgen format using R. Is there a simpler way?
snp genome • 4.1k views
ADD COMMENTlink modified 5.0 years ago by njbernstein40 • written 5.0 years ago by lisesandenbergh40
2
gravatar for njbernstein
5.0 years ago by
njbernstein40
United States
njbernstein40 wrote:

I had to do this recently (C: Recoding lgen into Ped using plink)

​awk 'NR > 11 { next }  !seen[$2]++ {print $2,$2,0,0}' final_report.txt > Plate3_final_report.fam
awk 'NR > 11 { next }  {print $2,$2,$1,$3,$4}' Plate3_final_report.txt > Plate3_final_report.lgen
awk 'NR < 2 { next } {print $3,$2,0,$4}' SNP_Map.txt > Plate3_final_report.map
ADD COMMENTlink modified 13 months ago by _r_am32k • written 5.0 years ago by njbernstein40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1418 users visited in the last hour
_