Vincenzo,

there is no function 'mcmcFmodel' in Geneland. Perhaps you intend to use 'MCMC'?

See ?MCMC for a proper use.

Additionally you might profit from reading the Geneland Manual on
http://folk.uio.no/gillesg/Geneland/Geneland.html


HTH
Christoph


________________________________________________________

Christoph Heibl

Systematic Botany
Ludwig-Maximilians-Universität München
Menzinger Str. 67
D-80638 München
GERMANY

phone:     +49-(0)89-17861-251
e-mail:    he...@lmu.de

http://www.christophheibl.de/ch-home.html

SAVE PAPER - THINK BEFORE YOU PRINT



On Jan 13, 2009, at 11:20 AM, Vincenzo Landi wrote:

I do the these passages:
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)
But at the end R say to me that do had found the function "mcmcFmodel".
Any sugestion?
thanks

**********************************************************
**********************************************************
VincenzoLandi
AnimalScience Phd student

Faculty of Agricolture
Universityof Perugia (Italy).

Departamentof Genetica
Campus deRabanales
Universidadde Cordoba
Cordoba(Spain)

Contacts:

cell phone: +39 3395388713
Skype: vincenzozoot
Italy: +39 0755857109
Spain:+34 957218706







        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting- guide.html
and provide commented, minimal, self-contained, reproducible code.

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to