Displaying 2 results from an estimated 2 matches for "logposterior".
2007 Jun 01
0
Metropolis code help
...nsaid,DF$diuretic,DF$diuretic*DF$nsaid)
y<- DF$Var3
lF1=plogis(X%*%as.vector(beta),log.p=TRUE)
lF2=plogis(-X%*%as.vector(beta),log.p=TRUE)
sum(y*lF1+(1-y)*lF2)
}
# Use a uniform prior for p
logprior <- function(beta,y,X) 0
# The log posterior is the sum. It's the target of our MCMC run
logposterior <- function(beta,y,X) logprior(beta,y,X)+logitll(beta,y,X)
start <- c(0,0,0,0)
sims <- Metropolis(logposterior, start, 10000, sd=1)
se_sims <- apply(sims, 2, sd)
sims <- Metropolis(logposterior, start, 10000,sd=se_sims)
cbind(rbind(mean(sims[1001:10000,1]),mean(sims[1001:10000,2]...
2009 Jan 13
3
problem whit Geneland
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 <-