[adegenet-forum] xval warning

LAURA NICOLE WOODINGS 17869067 at students.latrobe.edu.au
Sun Nov 1 03:22:08 CET 2015


Hi,

I am trying to use xval to cross-validate the number of PCs to retain for a DAPC using the below command and I am getting the following error msg:

xval_nopopbp75_cnt75_cl96_maf5to1 <- xvalDapc(nopop_bp75_cnt75_cl96_maf5to1NoNA at tab, nopop_clust$grp)

There were 50 or more warnings (use warnings() to see the first 50)
> warnings()
Warning messages:
1: In if (result == "overall") { ... :
  the condition has length > 1 and only the first element will be used
2: In if (result == "groupMean") { ... :
  the condition has length > 1 and only the first element will be used

Using find.clusters 2 groups are found but 1 group has 83 individuals and the other group has 4. Is the discrepancy in sample size causing the problem with xval? And if so how can I determine the optimal number of PC's to keep? I have attached my input file and the script that I use.

Also when conduction DAPC should you remove the NA's in the data? The DAPC seems to work when the Nas are not removed, but when I run it with them removed I get different results.

Thanks for your help!

Cheers,
Laura

PhD Candidate
Department of Ecology, Environment and Evolution| La Trobe University | Bundoora 3086 | Australia

m: +61 408 642 006 | e: 17869067 at students.latrobe.edu.au
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/adegenet-forum/attachments/20151101/1495d997/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: neutral_nopop_75bp_cnt75_cl96_2-16_ml80._maf0.05-0.01_missing0.8_hwe0.001.stru
Type: application/octet-stream
Size: 186036 bytes
Desc: neutral_nopop_75bp_cnt75_cl96_2-16_ml80._maf0.05-0.01_missing0.8_hwe0.001.stru
URL: <http://lists.r-forge.r-project.org/pipermail/adegenet-forum/attachments/20151101/1495d997/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Xval.R
Type: application/octet-stream
Size: 598 bytes
Desc: Xval.R
URL: <http://lists.r-forge.r-project.org/pipermail/adegenet-forum/attachments/20151101/1495d997/attachment-0003.obj>


More information about the adegenet-forum mailing list