Mar 20, 2019, 12:32:55 PM3/20/19
to R/qtl discussion
I'm having trouble with read.cross() for type "4way", reading from a 'csvr' file. I get the error
1: In summary.cross(cross) :
The genetic maps should all be matrices with two rows.
Currently, I have my data formatted so that;
Col1 = markername
Col2 = chromosome id
Col3 = male map position
Col4 = female map position
But clearly this is not what read.cross() expects. I am unable to find any examples of a 'correctly formatted' 4way cross anywhere, and write.cross() will not allow the export of 4way crosses to csv format, so I can't see the correct way to format this so that read.cross() can read it as a matrix with two rows.
Any help would be appreciated.
test <- read.cross('csvr','.', 'test.txt', genotypes=NULL, estimate.map=FALSE, crosstype="4way")