search

jgx65 / quantinemo

a swiss knife to simulate complex demographic and genetic scenarios
https://www2.unil.ch/popgen/softwares/quantinemo/
8 stars 0 forks source link

quanti_save_phenotype does not report the correct patch #26

Closed jgx65 closed 4 years ago

jgx65 commented 4 years ago

the patch reported with quanti_save_phenotype set to 1 is always patch 1.

frederic-michaud commented 4 years ago

Hi Jérôme,

Sorry. I don't get what the problem is exactly. Only individual from patch 1 are saved? Can you provide an input file as an example? Thanks!

Frédéric

jgx65 commented 4 years ago

Hi Frederic, sorry, I wasn't clear. All individuals phenotypes are reported, but they are all reported to belong to patch 1.

The first column of output files pheno_g200_r*.phe is made of 1s only instead of the corresponding patch number. see the example ini file below

# Simulations
generations 200
replicates 10
folder "Island"
random_per_replicate 1

# Metapopulation
patch_number 12
patch_capacity 50
dispersal_rate 0.02

# Quantitative trait
quanti_nb_trait 1
quanti_loci 10
quanti_all  2
quanti_mutation_model  2
quanti_mutation_rate 1e-5
quanti_allelic_var 0.005
quanti_environmental_model 0
quanti_heritability 0
quanti_environmental_proportion 1
quanti_va_model 1

quanti_save_phenotype 1
quanti_phenot_logtime 200
quanti_phenot_filename "pheno"

Cheers

frederic-michaud commented 4 years ago

Hello Jérôme,

Do you have a way to compile the code to see if the last commit 790ef30 solved your problem?

Cheers, Frédéric

jgx65 commented 4 years ago

Hello Frederic, this one is solved, many thanks!

Jérôme Goudet Dept. Ecology & Evolution Biophore, UNIL-Sorge CH-1015 Lausanne http://www.unil.ch/dee/goudet-group mail:jerome.goudet@unil.ch Tel:+41 (0)21 692 4242

From: frederic-michaud notifications@github.com Sent: Wednesday, 13 May 2020 21:05 To: jgx65/quantinemo quantinemo@noreply.github.com Cc: Jérôme Goudet jerome.goudet@unil.ch; Author author@noreply.github.com Subject: Re: [jgx65/quantinemo] quanti_save_phenotype does not report the correct patch (#26)

Hello Jérôme,

Do you have a way to compile the code to see if the last commit 790ef30https://github.com/jgx65/quantinemo/commit/790ef3028f47bfdd1d55cedab8869ca95cd77d65 solved your problem?

Cheers, Frédéric

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/jgx65/quantinemo/issues/26#issuecomment-628186147, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ABWOPWWDIAQDFFESOUYDJGLRRLVMXANCNFSM4M7U7GTA.