[R] problem whit Geneland

[Christoph Heibl]

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:    heibl at 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 at 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.