[adegenet-forum] extract.PLINKmap
Jombart, Thibaut
t.jombart at imperial.ac.uk
Mon Jan 28 16:13:51 CET 2013
Hello,
Can you please send a toy dataset and the code to reproduce the error?
Cheers
Thibaut
________________________________________
From: adegenet-forum-bounces at lists.r-forge.r-project.org [adegenet-forum-bounces at lists.r-forge.r-project.org] on behalf of H.Papoli Yazdi [h.papoli.yazdi at student.rug.nl]
Sent: 21 January 2013 15:27
To: adegenet-forum at r-forge.wu-wien.ac.at
Subject: [adegenet-forum] extract.PLINKmap
Hello,
I am encountering a problem in importing my SNP position data for a genlight object in adegenet.
I upload my genlight object using "read.PLINK" function. My .map file has the following format which is the standard PLINK format (chromosome, snp ID, genetic distance and physical position)
1 F0100010 0 7.29827
1 F0100020 0 127.006
1 F0100030 0 11.4989
I used "extract.PLINKmap" function to add the chromosome and position information to my genlight object, but it gives me the error message as follows:
"Warning message:
In matrix (txt, ncol=4, byrow=TRUE):
the length of the data [44706] is not a divisor or multiple of the number of lines [11177]"
Therefore, it does not load the position information. The problem is that it does not transform my .map file into a matrix of 4 columns and 44706 rows but it reduces the number of rows to 11177.
I would appreciate any suggestion very much.
Homa
More information about the adegenet-forum
mailing list