# how do I convert a vcf file for a GWAS study (using R package vcfR)

I have a vcf file with individual level genetic data, that I read using the R package vcfR

gen <- read.vcfR('test.vcf')

Scanning file to determine attributes.
File attributes:
meta lines: 9
variant count: 152
column count: 3798
All meta lines processed.
gt matrix initialized.
Character matrix gt created.
Character matrix gt rows: 152
Character matrix gt cols: 3798
skip: 0
nrows: 153
row_num: 0
Processed variant: 153
All variants processed


How do I convert this file to a data frame that I can use for a GWAS or Mendelian Randomization study?

I was told that this should be a data frame where each column is a genetic variant and it take values 0, 1, or 2, depending on the number of "mutations".

I never worked with genetic data before and do not know much about genetics.

• Welcome to Bioinformatics.StackExchange! GWAS is a broad area, and there are numerous tools you could use. The format of your data will depend on which tool you are using, and hopefully this will be clearly documented in the tool's manual. If on the other hand you're trying to implement your own tool, you'll have to be much more specific about what you want (or are required) to accomplish. – Daniel Standage Mar 26 at 16:35
• Welcome to the site! Could you please explain what have you tried? For instance, what happens if you use as(gen, "data.frame")? – llrs Mar 26 at 16:57
• I found that I can use extract.gt, then transform the matrix. Entries in this matrix are "0/0", "0/1", "1/0", and "1/1". I now have to recode these to 0, 1, 1, 2. – spore234 Mar 27 at 7:59