#matrix from Waller example for fsIndeterminacy from fungible
Lambda <- matrix(c(.8, 0,
.7, 0,
.6, 0,
0, .5,
0, .4,
0, .3), 6, 2, byrow=TRUE)
fsi(Lambda) #should be .7675 and .3837
A <- matrix( c(.8, 0,
.8, 0,
0, .6,
0, .6), 4,2, byrow=TRUE)
fsi(A) #should be .781 .529 (from Nicewander Appendix 1)
Run the code above in your browser using DataLab