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<-6 > 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="")) + } 51 1 124.976 24.56 67.119 48.824 67.196 48.87 67.019 49.066 1217.15 51 2 124.976 24.56 29.419 12.627 29.507 12.676 29.411 12.823 869.24 51 3 124.976 24.56 23.774 8.305 24.094 8.415 23.742 8.559 817.45 51 4 124.976 24.56 20.033 7.595 20.517 7.722 19.585 9.05 779.97 51 5 124.976 24.56 17.854 7.316 19.807 7.614 17.022 11.297 788.54 51 6 124.976 24.56 16.302 7.082 18.854 7.162 14.905 12.32 791.99 51 7 124.976 24.56 14.114 7.531 16.553 7.687 13.073 12.103 797.63 51 8 124.976 24.56 14.258 7.096 16.653 7.009 11.161 15.242 792.12 52 1 125.85 28.331 39.483 14.959 39.532 14.987 39.488 14.981 963.67 52 2 125.85 28.331 28.984 7.171 29.234 7.28 28.96 7.746 864.33 52 3 125.85 28.331 24.646 5.095 25.598 5.232 24.483 5.427 848.36 52 4 125.85 28.331 21.949 5.112 24.122 5.123 21.559 8.01 849.68 52 5 125.85 28.331 19.876 5.594 21.834 5.662 18.66 11.361 840.09 52 6 125.85 28.331 19.825 5.319 21.777 5.299 16.175 13.96 835.76 52 7 125.85 28.331 15.769 6.762 19.081 6.649 13.641 16.449 849.29 52 8 125.85 28.331 13.828 7.617 17.743 7.242 11.712 17.753 853.3 53 1 127.705 25.81 24.408 2.067 24.416 2.067 24.404 2.088 732.43 53 2 127.705 25.81 22.202 2.2 23.557 2.169 21.95 4.613 739.68 53 3 127.705 25.81 20.399 2.502 22.165 2.48 19.536 6.759 742.71 53 4 127.705 25.81 18.341 2.953 20.932 2.851 17.326 9.019 747.44 53 5 127.705 25.81 17.102 3.271 19.992 3.086 15.136 11.171 748.88 53 6 127.705 25.81 15.561 3.85 17.755 3.693 12.882 13.109 749.29 53 7 127.705 25.81 13.962 4.543 16.504 4.356 11.052 14.839 751.87 53 8 127.705 25.81 15.613 3.746 17.642 3.643 9.352 16.533 746.73 54 1 130.901 26.985 32.77 7.545 32.794 7.555 32.771 7.634 873.79 54 2 130.901 26.985 25.478 2.998 25.757 3.03 25.335 3.591 797.1 54 3 130.901 26.985 23.25 3.042 24.693 3.079 22.539 6.407 799.75 54 4 130.901 26.985 21.068 3.243 23.395 3.249 19.895 8.709 804.59 54 5 130.901 26.985 18.405 3.838 21.346 3.741 17.381 10.759 815.76 54 6 130.901 26.985 19.397 3.413 21.649 3.34 15.096 12.746 801.71 54 7 130.901 26.985 16.416 4.266 19.265 4.05 12.971 14.831 808.56 54 8 130.901 26.985 14.933 4.81 17.276 4.693 10.936 16.67 812.18 55 1 131.701 24.503 57.4 31.623 57.524 31.717 57.381 31.83 1143.11 55 2 131.701 24.503 42.533 21.157 43.266 21.699 42.629 21.88 1055.39 55 3 131.701 24.503 29.351 9.993 30.133 10.434 29.224 10.402 936.5 55 4 131.701 24.503 23.601 7.576 24.275 7.816 23.37 8.052 877.98 55 5 131.701 24.503 20.656 5.961 22.144 6.318 19.326 8.029 866.62 55 6 131.701 24.503 18.438 5.754 20.112 6.087 17.008 8.903 860.04 55 7 131.701 24.503 17.717 5.765 19.372 5.801 14.752 10.874 858.75 55 8 131.701 24.503 17.101 5.798 19.047 5.821 12.907 13.473 857.54 56 1 126.607 24.803 44.215 18.466 44.286 18.516 44.236 18.505 1018.33 56 2 126.607 24.803 30.422 8.336 30.897 8.598 30.558 8.928 891.8 56 3 126.607 24.803 24.633 5.032 25.253 5.238 24.619 5.922 838.79 56 4 126.607 24.803 21.135 4.578 22.712 5.033 21.059 6.435 829.66 56 5 126.607 24.803 19.909 4.532 21.715 4.635 17.984 9.576 820.49 56 6 126.607 24.803 16.57 5.47 19.24 5.552 15.662 11.067 829.78 56 7 126.607 24.803 16.376 5.167 18.736 5.165 13.577 12.838 823.89 56 8 126.607 24.803 15.189 5.388 17.56 5.323 11.781 15.585 827.59 57 1 123.018 23.755 24.173 2.858 24.182 2.86 24.167 2.85 727.52 57 2 123.018 23.755 21.118 2.13 21.441 2.136 20.949 3.283 708.51 57 3 123.018 23.755 20.111 2.087 21.214 2.07 18.898 5.257 711.88 57 4 123.018 23.755 18.057 2.375 20.274 2.317 16.867 7.382 723.25 57 5 123.018 23.755 15.937 2.83 19.042 2.605 14.809 9.383 735.31 57 6 123.018 23.755 15.954 2.748 18.633 2.593 12.791 11.286 729.27 57 7 123.018 23.755 15.661 2.791 17.93 2.68 10.96 13.196 729.36 57 8 123.018 23.755 15.751 2.757 17.908 2.636 9.281 14.744 728 58 1 126.352 24.042 24.692 3.073 24.704 3.075 24.687 3.121 737.94 58 2 126.352 24.042 22.105 2.5 23.501 2.502 21.895 5.468 739.28 58 3 126.352 24.042 19.79 2.705 21.265 2.678 19.464 5.695 735.61 58 4 126.352 24.042 18.204 3.072 19.857 3.053 17.254 8.097 731.12 58 5 126.352 24.042 16.056 3.738 18.332 3.654 14.931 10.422 740.71 58 6 126.352 24.042 15.073 3.967 17.66 3.875 12.995 11.674 740.23 58 7 126.352 24.042 13.323 4.663 16.322 4.449 11.247 13.504 754.11 58 8 126.352 24.042 12.163 5.15 14.848 4.941 9.785 14.977 756.99 59 1 123.066 26.144 40.206 16.764 40.26 16.802 40.18 16.68 971.55 59 2 123.066 26.144 24.144 3.21 24.277 3.237 24.114 3.338 774.37 59 3 123.066 26.144 22.508 3.188 23.892 3.201 21.407 6.113 778.83 59 4 123.066 26.144 21.082 3.245 23.366 3.192 19.024 8.588 785.12 59 5 123.066 26.144 18.867 3.653 21.551 3.616 16.587 10.984 790.71 59 6 123.066 26.144 18.336 3.651 20.837 3.588 14.507 13.091 791.57 59 7 123.066 26.144 18.958 3.451 20.951 3.442 12.602 14.462 790.05 59 8 123.066 26.144 14.982 4.716 17.623 4.567 10.682 16.38 798.93 60 1 131.74 25.156 62.935 34.038 63.087 34.149 62.89 34.177 1186.37 60 2 131.74 25.156 43.3 17.51 44.064 18.019 43.314 17.621 1062.73 60 3 131.74 25.156 29.867 9.47 30.414 9.772 29.331 10.126 936.89 60 4 131.74 25.156 23.376 5.835 24.155 6.118 22.872 7.1 862.75 60 5 131.74 25.156 20.636 4.909 22.556 5.185 19.854 7.613 864.87 60 6 131.74 25.156 18.968 4.775 21.29 4.766 16.915 9.607 854.53 60 7 131.74 25.156 17.579 4.977 19.612 4.879 14.549 12.073 853.4 60 8 131.74 25.156 16.24 5.236 18.826 5.125 12.07 14.776 857.63 > > > >