How to solve the problem when I get error message "column 1 is not a genotype object" in LDheatmap in R?
3
1
Entering edit mode
7.1 years ago
vyoouh ▴ 20

I created the files similar to en example given in http://cran.r-project.org/web/packages/LDheatmap/vignettes/LDheatmap.pdf. Which are CEUDist (only one column of distance) and CEUSNP (genotypic data, with the markers name). However, I obtained the error saying  "column 1 is not a genotype object"

CEUSNP

snp1    snp2    snp3    snp4    snp5    snp6
G/G    A/A    C/C    G/G    G/G    G/G
A/A    G/G    C/C    G/G    A/A    G/G
G/G    G/G    C/C    G/G    G/G    G/G
G/G    G/G    C/C    G/G    A/A    G/G
G/G    G/G    C/C    T/T    A/A    G/G
G/G    G/G    A/A    T/T    G/G    G/G
G/G    G/G    A/A    G/G    G/G    G/G
A/A    G/G    C/C    T/T    A/A    G/G
A/A    G/G    C/C    G/G    G/G    G/G
A/A    G/G    <NA>    G/G    A/A    G/G
A/A    G/G    A/A    G/G    A/A    G/G
G/G    G/G    A/A    G/G    G/G    G/G
G/G    G/G    C/C    G/G    A/A    G/G

 

and 

CEUDist 

123456
123524
123457
123525
123458
123526

Could you tell me what I did wrong in this case?

Also, can I use CC instead of C/C and NA instead of <NA> in the genotypic file?

It will be great if any of you can share an example files with me.

 

Thank you in advance.

Jane

LDheatmap • 5.6k views
ADD COMMENT
0
Entering edit mode

Thank you for your answer! I have met the same question, but I dont know how to use the fuction genotype, could you give the examples? I appreciate. Deletion should be reserved to errors or inappropriate content.

ADD REPLY
1
Entering edit mode
7.1 years ago

The vignette for that package is terrible. You need to load the genetics package and then use genotype() when creating the data.frame. In short, you need two columns for each individual, one column per allele. That genotype function will then return a properly formatted genotype object. You'll do that for each SNP and save the results in a data.frame. They really should have given an example in the vignette.

ADD COMMENT
0
Entering edit mode

Devon Ryan, Hi Devon Ryan. I am trying to add genotype object to my SNP data set. I am unable to use the codes described in the subsequent comment as my map file is not generated from the vcf tools. I used tassel generated VCF files instead where the V2 column contains the alleles (either form separated by /). The error in R says "Error in expectedGenotypes(alleles = alleles, haplotype = TRUE) : Zero length alleles vector." when i use the following codes to construct the dataframe

class(SNPdata[,1])

for(i in 1:num){SNPdata[,i]<-as.genotype(SNPdata[,i])}

Error in expectedGenotypes(alleles = alleles, haplotype = TRUE) : Zero length alleles vector.

Can you kindly help me with the same. Thank you in advance.

ADD REPLY
0
Entering edit mode

Devon Ryan , I was wondering if you had the chance to look at my request pertaining to problem in plotting LD plot using LDheatmap package. Your kind help is greatly appreciated. If you need more clarity on the specific issue with my problem, please do let me know. Many thanks for your help.

Bornali

ADD REPLY
1
Entering edit mode
4.8 years ago
jiaobingke ▴ 20

I came across the same question as you do. I convert the data.frame of my own data into a genotype object using the Devon Ryan's suggestion and It works! There is a as.genotype() function in genetics package.I use this function to convert the normal data.frame to a genotype object. The example is

num<-ncol(SNPdata)        # SNPdata is a data.frame which looks like the CEUSNP above.
colnames(SNPdata)<-map$V2 # Name the SNPdata columns with the 2nd column of map file(produced by vcftools).

class(SNPdata[,1]) # You will get "factor".

for(i in 1:num){
 SNPdata[,i]<-as.genotype(SNPdata[,i]) # convert the columns of the original data frame into  genotype objects
}

class(SNPdata[,1])  #Here you can see the class becomes "genotype factor".Then you can put the SNPdata into LDheatmap().

Thank you and Devon Ryan!

ADD COMMENT
0
Entering edit mode

Glad it works :)

ADD REPLY
0
Entering edit mode
3.8 years ago
ATCG ▴ 320

how can one create a data.frame which looks like the CEUSNP from a vcf file?.

ADD COMMENT

Login before adding your answer.

Traffic: 2355 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6