11 months ago by
Seeing the error you post in the second image, the problem seems to be that you are not telling
read.table to treat the first line as the column names. To do that, you should use
TEST <- read.table(file = "C:/Users/lenovo/Desktop/TEST.qassoc", quote = "\"", comment.char = "", header=TRUE)
and the manhattan function should the work.
Another option would be to use karyoploteR to create the manhattan plots (it will NOT create the QQ plots, though). For that, you'll need to have a file with the lengths of the sheep chromosomes (a text file with a chromosome per line, and two columns: the name and the number of bases). If you have any problem getting it, ask me and I can help you get the data.
You can load the data and plot the manhattan plot with this code. The plot will be a bit different since karyoploteR plots the dot at the exact position of each SNP instead of plotting all of them with the same separation.
snps.data <- toGRanges("C:/Users/lenovo/Desktop/TEST.qassoc")
kp <- plotKaryotype(genome="sheep.chromosome.lengths.txt")
You can get more information on how to create manhattan plots even in custom genomes in the karyoploteR tutorial