R version 2.10.1 (2009-12-14) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > library(coda) Loading required package: lattice > > S<-11000 ; burn<-1000 ; odens<-1 > > R0<-6 > KS<-7 > seeds<- 1:10 + (KS-1)*10 > MS<-NULL > > for(seed in seeds) + { + set.seed(seed) + source("datagen.r") # generate data + for(R in 1:8) # for each presumed rank + { + source("hbayes.r") + mse<-c(mean(Y^2), # size of matrix + mean((Y-M0)^2), # MSE of Y + mean((Y-M.ps/nss)^2) , # RSS of mean + mean((M0-M.ps/nss)^2), # MSE of of mean + mean((Y-M.mode)^2) , # RSS of mode + mean((M0-M.mode)^2), # MSE of mode + mean((Y-M.ols)^2), # RSS of ols + mean((M0-M.ols)^2)) # MSE of ols + cat(seed,R,round(100*mse,3)," ",round(DIC,2),"\n") + MS<-rbind(MS,c(seed,R,mse,DIC,eSS)) + } + dput(MS,paste("hbayes.SS",KS,".R0",R0,".results",sep="")) + } 61 1 131.584 24.767 65.916 40.837 66.154 41.028 65.898 41.121 1209.71 61 2 131.584 24.767 41.431 16.719 41.949 17.082 41.334 16.751 1037.67 61 3 131.584 24.767 25.344 7.798 25.649 7.847 25.327 8.453 853.63 61 4 131.584 24.767 20.397 5.12 21.569 5.351 20.293 6.119 814.8 61 5 131.584 24.767 17.792 5.291 20.152 5.368 17.257 8.881 810.84 61 6 131.584 24.767 15.585 6.009 17.778 5.746 14.918 10.96 798.7 61 7 131.584 24.767 15.051 5.968 17.059 5.817 13.166 12.505 795.68 61 8 131.584 24.767 12.966 6.786 15.766 6.515 10.915 14.806 807.32 62 1 128.67 25.583 36.813 10.215 36.846 10.231 36.814 10.243 929.87 62 2 128.67 25.583 28.346 6.481 28.722 6.651 28.45 7.501 857.17 62 3 128.67 25.583 23.657 4.554 24.45 4.798 23.604 6.231 819.94 62 4 128.67 25.583 20.192 4.912 21.791 5.19 20.33 8.359 820.06 62 5 128.67 25.583 18.389 4.957 20.896 4.99 17.433 10.53 811.02 62 6 128.67 25.583 17.106 5.102 19.238 5.165 15.105 12.132 808.98 62 7 128.67 25.583 14.962 5.85 17.493 5.863 13.077 13.861 813.99 62 8 128.67 25.583 15.839 5.147 18.105 5.029 11.119 15.829 808.11 63 1 124.478 27.722 76.292 58.388 77.465 59.552 77.052 60.013 1290.39 63 2 124.478 27.722 44.135 23.681 44.462 23.918 44.102 23.765 1065.9 63 3 124.478 27.722 32.518 13.335 34.085 14.251 33.019 14.278 997.33 63 4 124.478 27.722 24.112 9.7 24.709 9.928 23.658 10.354 882.79 63 5 124.478 27.722 20.018 9.428 20.683 9.568 19.622 11.323 846.43 63 6 124.478 27.722 17.764 9.168 19.733 9.513 17.065 13.122 849.74 63 7 124.478 27.722 15.454 9.532 17.434 9.621 14.565 13.939 851.91 63 8 124.478 27.722 14.918 9.364 17.19 9.41 11.735 16.79 847.39 64 1 130.81 24.18 56.453 34.195 56.605 34.321 56.456 34.412 1135.58 64 2 130.81 24.18 32.418 13.333 32.702 13.528 32.264 13.213 917.27 64 3 130.81 24.18 24.405 7.493 24.741 7.589 24.051 7.595 835.31 64 4 130.81 24.18 20.657 5.979 21.366 6.091 20.515 7.163 811.75 64 5 130.81 24.18 18.702 5.539 20.774 5.932 17.72 9.572 806.11 64 6 130.81 24.18 17.138 5.719 19.53 5.631 15.227 11.724 813.43 64 7 130.81 24.18 16.182 5.776 18.362 5.689 13.125 13.11 813.11 64 8 130.81 24.18 15.222 5.939 17.761 5.919 11.495 14.071 817.25 65 1 121.992 25.799 43.065 18.919 43.132 18.969 43.069 18.888 1005.21 65 2 121.992 25.799 27.682 6 28.016 6.017 27.759 6.394 841.15 65 3 121.992 25.799 22.899 3.092 23.613 2.987 22.619 4.438 791.81 65 4 121.992 25.799 21.717 3.113 23.137 3.139 20.053 7.169 797.23 65 5 121.992 25.799 20.168 3.306 22.213 3.293 17.773 9.262 803.84 65 6 121.992 25.799 19.554 3.299 21.573 3.234 15.669 11.298 807.05 65 7 121.992 25.799 15.642 4.302 18.851 4.055 13.503 12.949 830.3 65 8 121.992 25.799 17.705 3.601 19.951 3.501 11.453 15.111 813.88 66 1 136.673 28.596 34.75 7.079 34.772 7.087 34.754 7.174 902.06 66 2 136.673 28.596 27.87 3.819 28.446 4.025 27.822 4.288 853.1 66 3 136.673 28.596 24.465 3.317 26.047 3.234 24.032 6.111 841.23 66 4 136.673 28.596 22.297 3.418 24.996 3.343 21.493 8.759 850.91 66 5 136.673 28.596 20.067 3.801 23.496 3.532 19.056 10.728 862.08 66 6 136.673 28.596 21.238 3.421 23.321 3.303 16.728 13.14 840.31 66 7 136.673 28.596 16.178 4.949 20.375 4.511 14.263 15.313 878.75 66 8 136.673 28.596 18.384 4.009 21.118 3.831 12.093 16.811 853.66 67 1 122.212 24.266 40.658 19.627 40.724 19.676 40.68 19.669 977.94 67 2 122.212 24.266 23.745 5.168 23.961 5.167 23.78 5.543 767.37 67 3 122.212 24.266 20.884 3.647 21.904 3.591 20.844 4.777 764.85 67 4 122.212 24.266 19.004 3.67 21.045 3.625 18.353 7.38 768.22 67 5 122.212 24.266 17.316 3.945 19.514 3.873 16.066 8.894 770.2 67 6 122.212 24.266 15.91 4.235 18.36 4.167 13.797 11.808 773.11 67 7 122.212 24.266 13.697 5.039 16.671 4.769 12.093 13.197 783.22 67 8 122.212 24.266 13.303 5.071 15.948 4.905 10.02 14.916 780.54 68 1 131.374 26.064 32.061 6.547 32.081 6.555 32.064 6.604 863.5 68 2 131.374 26.064 25.473 3.539 25.764 3.605 25.313 4.281 800.38 68 3 131.374 26.064 22.002 3.451 23.061 3.482 22.047 5.777 788.56 68 4 131.374 26.064 19.69 3.778 21.611 3.782 19.153 7.842 792.79 68 5 131.374 26.064 18.752 3.715 20.786 3.633 16.545 10.405 791.44 68 6 131.374 26.064 15.849 4.76 18.788 4.659 14.525 12.227 809 68 7 131.374 26.064 17.09 4.038 19.323 3.974 12.712 13.859 793.48 68 8 131.374 26.064 15.422 4.53 18.069 4.372 10.947 15.595 798.89 69 1 121.432 21.069 41.798 24.951 41.866 25.009 41.811 25.071 991.14 69 2 121.432 21.069 27.753 11.647 27.949 11.773 27.778 11.806 844.12 69 3 121.432 21.069 21.361 6.226 21.624 6.351 21.237 6.401 766.98 69 4 121.432 21.069 19.11 5.401 20.298 5.706 18.537 6.653 764.57 69 5 121.432 21.069 17.508 5.13 19.235 5.361 16.344 8.18 769.67 69 6 121.432 21.069 15.648 5.099 18.157 5.282 14.007 9.316 775.49 69 7 121.432 21.069 14.937 5.054 17.393 5.166 12.179 10.684 776.39 69 8 121.432 21.069 13.373 5.189 15.978 5.108 10.365 12.869 784.8 70 1 130.293 24.009 38.3 18.282 38.341 18.31 38.295 18.487 948.65 70 2 130.293 24.009 26.403 7.14 26.611 7.21 26.415 7.47 818.57 70 3 130.293 24.009 21.499 5.262 22.224 5.484 21.394 6.223 782.75 70 4 130.293 24.009 19.138 5.214 21.438 5.001 18.447 8.493 778.79 70 5 130.293 24.009 16.254 6.181 18.857 6.179 15.849 11.116 776.91 70 6 130.293 24.009 16.169 5.636 17.989 5.724 13.599 12.393 769.03 70 7 130.293 24.009 14.404 6.17 16.512 6.304 11.749 14.153 768.41 70 8 130.293 24.009 14.174 6.113 16.224 6.278 9.942 15.469 764.42 > > > >