N <- 100 x <- seq(from=0,to=1,length=N) x <- x[c(-1,-N)] y <- x X <- c(40,40) Y <- c(40,5000) beta.x <- dbeta(x,X[1]+1,X[2]+1,log=TRUE) beta.y <- dbeta(y,Y[1]+1,Y[2]+1,log=TRUE) #prev.data <- c(560,1000) prev.data <- X+Y p <- x beta.p <- d…