# 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)
ha<-hmemIFS(ma,na)
mb<-memG(a1,b1,y)
nb<-nonmemS(mb,lam)
hb<-hmemIFS(mb,nb)
k<-1
p<-2
simL(ma,na,mb,nb,ha,hb,p,k)
#[1] 0.7022635 0.6896045 0.8890488 0.8890488
# }
Run the code above in your browser using DataLab