c1=data.frame(v1=rnorm(30,0),v2=rnorm(30,0),v3=factor(round(runif(30))+1))
c2=data.frame(v1=rnorm(30,2),v2=rnorm(30,4),v3=factor(round(runif(30))+4))
data=rbind(c1,c2)
res=mspec(data, k = 2)
summary(res)
plot(res)
Run the code above in your browser using DataLab