# NOT RUN {
#When data set y consist of only one row use k=1
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)
mb<-memG(a1,b1,y)
nb<-nonmemS(mb,lam)
p<-2
t<-2
k<-1
simBA(ma,na,mb,nb,p,t,k)
#0.7072291 0.6947466 0.8919850 0.8919850
#When data set y having more than one rows
#use k = the number of rows of data set y
x<-matrix(c(12,9,14,11,21,16,15,24,20,17,14,11),nrow=4)
y<-matrix(c(11,24,21,12,6,11),nrow=2)
a<-mn(x)
b<-std(x)
a1<-mn(y)
b1<-std(y)
lam<-0.5
ma<-memG(a,b,x)
na<-nonmemS(ma,lam)
mb<-memG(a1,b1,y)
nb<-nonmemS(mb,lam)
p<-2
t<-2
sim<-c()
for(k in 1:nrow(y)){sim<-rbind(sim,simBA(ma,na,mb,nb,p,t,k))}
sim
# [,1] [,2] [,3] [,4]
#[1,] 0.7072291 0.6947466 0.8919850 0.8919850
#[2,] 0.9410582 0.9843247 0.7380007 0.7380007
# }
Run the code above in your browser using DataLab