> i=grep('f toN',rownames(x)) > x1=(unique(x[i,1,drop=F])) > rn=c(as.numeric(gsub('.+g1 t','',rownames(x1))),91:99) > for (j in 1:length(rn)) { + cat('\npsi(',rn[j],'interval=',intervals[rn[j]],'):\n'); + print(getpsi(rn[j])) + } psi( 1 interval= 1 ): f m N R M f 0.6621744 0.0000000 0.1689128 0.1689128 0.0000000 m 0.0000000 0.6665206 0.0000000 0.0000000 0.3334794 N 0.0000000 0.0000000 0.2922841 0.7077159 0.0000000 R 0.0000000 0.0000000 0.1662733 0.8337267 0.0000000 M 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000 psi( 3 interval= 1 ): f m N R M f 0.7113236 0.000000 0.1443382 0.1443382 0.000000 m 0.0000000 0.701836 0.0000000 0.0000000 0.298164 N 0.0000000 0.000000 0.4044058 0.5955942 0.000000 R 0.0000000 0.000000 0.1789342 0.8210658 0.000000 M 0.0000000 0.000000 0.0000000 0.0000000 1.000000 psi( 5 interval= 1 ): f m N R M f 0.8220748 0.000000 0.0889626 0.0889626 0.000000 m 0.0000000 0.874558 0.0000000 0.0000000 0.125442 N 0.0000000 0.000000 0.7032824 0.2967176 0.000000 R 0.0000000 0.000000 0.2837450 0.7162550 0.000000 M 0.0000000 0.000000 0.0000000 0.0000000 1.000000 psi( 7 interval= 3 ): f m N R M f 0.8785354 0.0000000 0.1214646 0.0000000 0.0000000 m 0.0000000 0.9188044 0.0000000 0.0000000 0.0811956 N 0.0000000 0.0000000 1.0000000 0.0000000 0.0000000 R 0.0000000 0.0000000 0.1382449 0.8617551 0.0000000 M 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000 psi( 91 interval= 5 ): f m N R M f 0.8785354 0.0000000 0.1214646 0.0000000 0.0000000 m 0.0000000 0.9188044 0.0000000 0.0000000 0.0811956 N 0.0000000 0.0000000 1.0000000 0.0000000 0.0000000 R 0.0000000 0.0000000 0.1382449 0.8617551 0.0000000 M 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000 psi( 92 interval= 5 ): f m N R M f 0.8785354 0.0000000 0.1214646 0.0000000 0.0000000 m 0.0000000 0.9188044 0.0000000 0.0000000 0.0811956 N 0.0000000 0.0000000 1.0000000 0.0000000 0.0000000 R 0.0000000 0.0000000 0.1382449 0.8617551 0.0000000 M 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000 psi( 93 interval= 5 ): f m N R M f 0.8785354 0.0000000 0.1214646 0.0000000 0.0000000 m 0.0000000 0.9188044 0.0000000 0.0000000 0.0811956 N 0.0000000 0.0000000 1.0000000 0.0000000 0.0000000 R 0.0000000 0.0000000 0.1382449 0.8617551 0.0000000 M 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000 psi( 94 interval= 1 ): f m N R M f 0.8785354 0.0000000 0.1214646 0.0000000 0.0000000 m 0.0000000 0.9188044 0.0000000 0.0000000 0.0811956 N 0.0000000 0.0000000 1.0000000 0.0000000 0.0000000 R 0.0000000 0.0000000 0.1382449 0.8617551 0.0000000 M 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000 psi( 95 interval= 2 ): f m N R M f 0.8785354 0.0000000 0.1214646 0.0000000 0.0000000 m 0.0000000 0.9188044 0.0000000 0.0000000 0.0811956 N 0.0000000 0.0000000 1.0000000 0.0000000 0.0000000 R 0.0000000 0.0000000 0.1382449 0.8617551 0.0000000 M 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000 psi( 96 interval= 1 ): f m N R M f 0.8785354 0.0000000 0.1214646 0.0000000 0.0000000 m 0.0000000 0.9188044 0.0000000 0.0000000 0.0811956 N 0.0000000 0.0000000 1.0000000 0.0000000 0.0000000 R 0.0000000 0.0000000 0.1382449 0.8617551 0.0000000 M 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000 psi( 97 interval= 2 ): f m N R M f 0.6621744 0.0000000 0.1689128 0.1689128 0.0000000 m 0.0000000 0.6665206 0.0000000 0.0000000 0.3334794 N 0.0000000 0.0000000 0.2922841 0.7077159 0.0000000 R 0.0000000 0.0000000 0.1662733 0.8337267 0.0000000 M 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000 psi( 98 interval= 2 ): f m N R M f 1 0 0 0 0 m 0 1 0 0 0 N 0 0 1 0 0 R 0 0 0 1 0 M 0 0 0 0 1 psi( 99 interval= 2 ): f m N R M f 1 0 0 0 0 m 0 1 0 0 0 N 0 0 1 0 0 R 0 0 0 1 0 M 0 0 0 0 1 > sink()