From guillot at agroparistech.fr  Wed Oct 10 16:35:41 2007
From: guillot at agroparistech.fr (guillot at agroparistech.fr)
Date: Wed, 10 Oct 2007 16:35:41 +0200 (CEST)
Subject: [Rd] package Geneland / Rgui under windows (PR#9964)
Message-ID: <20071010143541.CF35F28800A8@mail.pubhealth.ku.dk>

Full_Name: Gilles Guillot
Version: 2.6.0
OS: windows XP professional
Submission from: (NULL) (129.240.88.50)


This sequence of command does not work in the R gui 2.6.0:

library(Geneland)

set.seed(1)

data <- simdata(nindiv=200,
                coord.lim=c(0,1,0,1) ,
                number.nuclei=5 ,
                allele.numbers=rep(10,20),
                IBD=FALSE,
                npop=2,
                give.tess.grid=FALSE)

geno <- data$genotypes
coord <- t(data$coord.indiv)

path.mcmc <- paste(tempdir(),"/",sep="") 

set.seed(1)
mcmcFmodel(coordinates=coord,
           genotypes=geno,
           path.mcmc=path.mcmc,
           rate.max=10,
           delta.coord=0,
           npopmin=1,
           npopinit=5,
           npopmax=5,
           nb.nuclei.max=50,
           nit=500,
           thinning=1,
           freq.model="Dirichlet",
           varnpop=FALSE,
           spatial=TRUE)


the call to mcmcFmodel freezes R

The same sequence of command works in the R command line of R 2.6.0 
and also in the R GUI of 2.5.1



[Rd] package Geneland / Rgui under windows (PR#9964)
guillot at agroparistech.fr