# NOT RUN {
x<-matrix(c(12,9,14,11,21,16,15,24,20,17,14,11),nrow=4)
y<-matrix(c(11,21,6),nrow=1)
a<-mn(x)
b<-std(x)
a1<-mn(y)
b1<-std(y)
lam<-0.5
ma<-memG(a,b,x)
na<-nonmemS(ma,lam)
ia<-imemSFS(ma,na)
ra<-rmemSFS(ma,na,ia)
mb<-memG(a1,b1,y)
nb<-nonmemS(mb,lam)
ib<-imemSFS(mb,nb)
rb<-rmemSFS(mb,nb,ib)
k<-1
simWWLWW9(ma,na,mb,nb,ia,ib,ra,rb,k)
#[1] 0.04488958 0.04334510 0.08280306 0.08280306
# }
Run the code above in your browser using DataLab