Hi all,

I have managed to get the r-scale usage algorithm to work, but I need to obtain the final results from this.  As I understand it, this code is designed to generate a matrix after processing and store it somewhere?????  

Here is the code.....

I get this part of the code, it all makes sense.

##
if(nchar(Sys.getenv("LONG_TEST")) != 0) {R=1000} else {R=5} 
{
data(customerSat)
surveydat = list(k=10,x=as.matrix(customerSat))
Mcmc1 = list(R=R)
set.seed(66)
out=rscaleUsage(Data=surveydat,Mcmc=Mcmc1)
summary(out$mudraw)
}

My question is how do I retrieve the results from this in a matrix format????

I want to extract the matrix in its complete form and save this as a file for further processing.

What I want to get is.....

            V1     V2      V3      V4
V1   1.47    0.99    1.14    1.56 
V2   1.67    0.93    1.35    1.21 

If I type 'out' see below in the GUI interface I get the following:

>out
       [,1288] [,1289] [,1290] [,1291] .......
[3,]    1.47    0.99    1.14    1.56    ........
[4,]    1.17    0.96    0.63    1.32    ........ 
[5,]    1.65    0.51    0.39    1.62    ....... 

>summary(out)
           Length Class      Mode   
Sigmadraw   500   bayesm.var numeric
mudraw       50   bayesm.mat numeric
taudraw    9055   bayesm.mat numeric
sigmadraw  9055   bayesm.mat numeric
Lambdadraw   20   bayesm.mat numeric
edraw         5   bayesm.mat numeric

Thanks Paul

