Re: blm: Extracting simulated values of parameters (blm.fit, betaMatrix, scaleMatrix)

["Southworth, Harry" <Harry.Southworth@astrazeneca.com>]

The object returned by blm is an S4 class object.

 

mod = blm( y ~ x, data )

getSlots( mod )

 

The slot called "model" is a list, the first element of which also seems to have an S4 class.

 

getSlots( mod@model[[ 1 ]] )

 

The slot called "betaMatix" appears to contain the simulated values for the intercepts and slopes, and the slot called "scaleMatrix" contains the simulated values of the scale parameter.

 

m = mod@model[[ 1 ]]@betaMatrix

 

hist( m[ , 1 ] ) # intercept

 

s = mod@model[[ 1 ]]@scaleMatrix

 

hist( s )