[adegenet-forum] Problem converting structure file to genind (length of 'dimnames' [2] not equal to array extent)
Jombart, Thibaut
t.jombart at imperial.ac.uk
Wed Jul 29 11:27:13 CEST 2015
Hi Charlotte,
the conversion looks OK to me. In genind, data are recoded as alleles counts for each [locus.allele]. So if tb109 is "-9" (=missing), then all the alleles for tb109 are 'NA'.
If you want to double check that conversion worked, it is probably best to compare your input to genind2df(x) where 'x' is your converted genind.
Best
Thibaut
________________________________
From: Charlotte Hurry [charlotte.hurry at griffithuni.edu.au]
Sent: 29 July 2015 08:51
To: Jombart, Thibaut
Cc: adegenet-forum at lists.r-forge.r-project.org
Subject: Re: [adegenet-forum] Problem converting structure file to genind (length of 'dimnames' [2] not equal to array extent)
Hello again
I hope you don't mind but I'd like to double check my file conversion to genind. I understand that missing values are transformed to NA but I don't understand why I have NA values at all loci for some individuals. I have read the information online regarding Genind objects but nothing that explains the conversion of allele data to the values in the Genind object.
For instance my genind object is as follows:
> tab(obj[loc=1])
tb109.118 tb109.120 tb109.114 tb109.116 tb109.112 tb109.122 tb109.108
GM4_283 1 1 0 0 0 0 0
GM336 0 0 1 1 0 0 0
GM337 NA NA NA NA NA NA NA
E276 2 0 0 0 0 0 0
E277 1 0 0 1 0 0 0
E278 0 0 2 0 0 0 0
Whereas my original data file is this:
tb109 tb117 tb129 tb130 tb136 tb187 tb188
GM4_283 1 118 108 132 125 122 168 159
GM4_283 1 120 112 138 129 122 176 168
GM336 1 114 112 128 123 124 168 -9
GM336 1 116 112 136 123 128 176 -9
GM337 1 -9 110 126 123 126 168 159
GM337 1 -9 112 126 123 126 176 174
E276 2 118 108 -9 -9 122 166 165
E276 2 118 112 -9 -9 122 176 190
E277 2 116 108 126 119 124 168 159
E277 2 118 112 138 129 128 176 174
What is the mechanism that forces all values of GM336 to be NA values in the Genind, and is this correct?
Thank you once again.
Charlotte
On 29 July 2015 at 01:57, Jombart, Thibaut <t.jombart at imperial.ac.uk<mailto:t.jombart at imperial.ac.uk>> wrote:
Hi Charlotte,
thanks for a detailed report.
STRUCTURE files are a bit of a pain because of all the possible options, so in general if you can get them in a different format (e.g. GENETIX, Fstat) it will make your life easier.
Here, your instruction seems okay given the sample of the file, but the file itself seems wrong. STRUCTURE's doc suggests you should remove 'Ind' and 'Pop' from the first line.
http://pritchardlab.stanford.edu/structure_software/release_versions/v2.3.4/structure_doc.pdf
It works on the sample provided.
Was this file generated by STRUCTURE directly?
If the STRUCTURE standards have changed, please fill in an issue and we'll adapt read.structure:
https://github.com/thibautjombart/adegenet/issues
Best
Thibaut
________________________________
From: adegenet-forum-bounces at lists.r-forge.r-project.org<mailto:adegenet-forum-bounces at lists.r-forge.r-project.org> [adegenet-forum-bounces at lists.r-forge.r-project.org<mailto:adegenet-forum-bounces at lists.r-forge.r-project.org>] on behalf of Charlotte Hurry [charlotte.hurry at griffithuni.edu.au<mailto:charlotte.hurry at griffithuni.edu.au>]
Sent: 24 July 2015 05:08
To: adegenet-forum at lists.r-forge.r-project.org<mailto:adegenet-forum at lists.r-forge.r-project.org>
Subject: [adegenet-forum] Problem converting structure file to genind (length of 'dimnames' [2] not equal to array extent)
Hello
I have been trying to convert my structure file to a genind object but I am having some issues. I am getting this error and I can't find an appropriate fix. I have looked through forum history and found a similar issue but the fix suggested didn't work for me.
my file looks like this:
Ind Pop tb109 tb117 tb129 tb130 tb136 tb187 tb188
GM4_283 1 118 108 132 125 122 168 159
GM4_283 1 120 112 138 129 122 176 168
GM336 1 114 112 128 123 124 168 0
GM336 1 116 112 136 123 128 176 0
GM337 1 0 110 126 123 126 168 159
GM337 1 0 112 126 123 126 176 174
E276 2 118 108 0 0 122 166 165
E276 2 118 112 0 0 122 176 190
E277 2 116 108 126 119 124 168 159
E277 2 118 112 138 129 128 176 174
E278 2 114 108 124 121 122 168 162
E278 2 114 112 126 125 128 176 180
I used the script
> tems<-read.structure("TBtest1.str",row.marknames=1,onerowperind=FALSE, n.ind=190, n.loc=7, col.lab=1, col.pop=2, ask=TRUE)
Which brings up the following message, and I press RETURN:
Which other optional columns should be read (press 'return' when done)? 1:
Converting data from a STRUCTURE .stru file to a genind object...
Then get this error:
Error in `colnames<-`(`*tmp*`, value = c("Ind", "Pop", "tb109", "tb117", :
length of 'dimnames' [2] not equal to array extent
Has anybody got a fix to get my original STRUCTURE file to work as a Genind object.?
As an aside, on the forum I found a work around under the title "Trouble converting to genid object", which converts a read.table document into genind which I followed and I got this error:
Warning message:
In .local(.Object, ...) : NAs introduced by coercion
and this file:
> head(as.matrix(obj1))
V1 V2 V3 V4 V5 V6 V7 V8 V9
001 NA NA NA NA NA NA NA NA NA
002 NA 1 118 108 132 125 122 168 159
003 NA 1 120 112 138 129 122 176 168
004 NA 1 114 112 128 123 124 168 0
005 NA 1 116 112 136 123 128 176 0
006 NA 1 0 110 126 123 126 168 159
So i would like a different solution to my problem if possible.
Many thanks
Charlotte
Charlotte
--
Charlotte Hurry
PhD Candidate
Australian Rivers Institute
Griffith University
07 37356655
http://www.griffith.edu.au/environment-planning-architecture/australian-rivers-institute/hdr-students/charlotte-hurry
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/adegenet-forum/attachments/20150729/133cb17d/attachment.html>
More information about the adegenet-forum
mailing list