[R-sig-genetics] Query regarding DAPC using ADEGENT in R

[Bhuller, Ravneet]

Hello Thibaut,

The genind object is of size 8.9 Mb.

Yesterday, when I tried find.clusters in the Rmarkdown document, it did not ask for any input. I wanted to know how many components to retain to run DAPC.

However, today, when I tried find.clusters in R (without Rmarkdown), it worked absolutely fine.

But another problem is that pop function is giving an error. Attached is the screenshot. My data is DNA multiple-sequence-alignment FASTA files of 357 genes, which are concatenated. Each file has same 220 bacterial strains. And my meta data includes 220 strains distributed into 19 serovars (integers). Is my meta data in the correct required form?

[cid:943D7881-6E60-444E-8038-5754883AF5B2 at dlink.com]

Kind regards,

Rav




On 1 Sep 2016, at 16:17, Thibaut Jombart <thibautjombart at gmail.com<mailto:thibautjombart at gmail.com>> wrote:

Hi,

what is the size of the dataset? Are you sure it did not display a graph of
cumulative variance and waits for your input? As you specified it, the
function is interactive and will wait on your input before giving results.

Best
Thibaut


--
Dr Thibaut Jombart
Lecturer, Department of Infectious Disease Epidemiology
Imperial College London
https://sites.google.com/site/thibautjombart/
https://github.com/thibautjombart
Twitter: @TeebzR <https://twitter.com/TeebzR>

On 1 September 2016 at 00:33, Graham Reynolds <greynold at unca.edu> wrote:

Hi Rav, if you have multiple cores on your CPU you can distribute the task
among them, which really speeds it up. That might at least let you see
whether it is a problem with the dataset.

grp <-find.clusters(app, max.n.clust=n, parallel = "multicore", ncpus = 4)

You might need to install the Parallel package

Cheers


R. Graham Reynolds, Ph.D.
Assistant Professor of Biology
University of North Carolina Asheville
www.caribbeanboas.org



*From:* R-sig-genetics <r-sig-genetics-bounces at r-project.org> on behalf of
Bhuller, Ravneet <ravneet.bhuller13 at imperial.ac.uk>

*Sent:* Wednesday, August 31, 2016 5:39 AM
*To:* r-sig-genetics at r-project.org
*Subject:* [R-sig-genetics] Query regarding DAPC using ADEGENT in R

Dear All,

I am trying to do DAPC on my genind object (app) using ADEGENET in R. But
when I run the following command:

grp <- find.clusters(app)

it is processing continuously for more than 2 hours but not producing any
output. Any guidance will be very much appreciated.

Regards,

Rav

       [[alternative HTML version deleted]]

_______________________________________________
R-sig-genetics mailing list
R-sig-genetics at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-genetics
R-sig-genetics Info Page - stat.ethz.ch
<https://stat.ethz.ch/mailman/listinfo/r-sig-genetics>
stat.ethz.ch
The R-sig-genetics mailing list is devoted to friendly, insightful,
exciting discussions on genetic data analysis using the R software.
Covered
topics include:



       [[alternative HTML version deleted]]

_______________________________________________
R-sig-genetics mailing list
R-sig-genetics at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-genetics


[[alternative HTML version deleted]]

_______________________________________________
R-sig-genetics mailing list
R-sig-genetics at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-genetics

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-genetics/attachments/20160902/3f957198/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: DAPC ERROR.png
Type: image/png
Size: 84722 bytes
Desc: DAPC ERROR.png
URL: <https://stat.ethz.ch/pipermail/r-sig-genetics/attachments/20160902/3f957198/attachment-0001.png>