<paulandpen <at> optusnet.com.au> writes:
> 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????
>
str() is your friend...
-----
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)
}
> str(out)
List of 6
$ Sigmadraw : bayesm.var [1:5, 1:100] 7.37 7.04 5.56 5.00 4.50 ...
..- attr(*, "class")= chr [1:3] "bayesm.var" "bayesm.mat" "mcmc"
..- attr(*, "mcpar")= num [1:3] 1 5 1
$ mudraw : bayesm.mat [1:5, 1:10] 6.21 6.54 6.49 6.56 6.59 ...
..- attr(*, "class")= chr [1:2] "bayesm.mat" "mcmc"
..- attr(*, "mcpar")= num [1:3] 1 5 1
$ taudraw : bayesm.mat [1:5, 1:1811] 0.1862 -1.4988 -1.0348 0.0754 1.0258
...
..- attr(*, "class")= chr [1:2] "bayesm.mat" "mcmc"
..- attr(*, "mcpar")= num [1:3] 1 5 1
$ sigmadraw : bayesm.mat [1:5, 1:1811] 1.11 1.29 1.47 1.2 1.95 0.72 0.63 0.69
1.2 1.14 ...
..- attr(*, "class")= chr [1:2] "bayesm.mat" "mcmc"
..- attr(*, "mcpar")= num [1:3] 1 5 1
$ Lambdadraw: bayesm.mat [1:5, 1:4] 3.35 2.73 2.58 2.49 2.73 ...
..- attr(*, "class")= chr [1:2] "bayesm.mat" "mcmc"
..- attr(*, "mcpar")= num [1:3] 1 5 1
$ edraw :Classes 'bayesm.mat', 'mcmc' atomic [1:5] 0.000 0.001 0.001 0.001
0.001
.. ..- attr(*, "mcpar")= num [1:3] 1 5 1
>
Looks like you want sigmadraw (lowercase!)? So try
out$sigmadraw
Dieter
______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
