library(mvmeta)
y = matrix(0,8,2)
S = inputcov(matrix(runif(8*2,0.1,2),8,2,dimnames=list(NULL,c("V1","V2"))),cor=c(0.5))
Psi = inputcov(1:2,cor=0.3)
mydat=mvmetaSim(y,S,Psi)
y1=mydat[,1]
y2=mydat[,2]
s1=sqrt(S[,1])
s2=sqrt(S[,3])
data=data.frame(y1, s1, y2, s2)
galaxy(data, type = "continuous", method="galaxy.cl", k=2, L=1, estimator="R0", maxiter=150)Run the code above in your browser using DataLab