Displaying 2 results from an estimated 2 matches for "rmnlindepmetrop".
2006 Feb 08
0
bayesm, rmnlIndepMetrop
Hi,
I tried to use rmnlIndepMetrop (bayesm package) for my MNL model with 4
choice alternatives, 5 independent variables, 69 observations,
dim(X) [1] 276 5, nu=6. So I run such code:
if(nchar(Sys.getenv("LONG_TEST")) != 0) {R=2000} else {R=10}
set.seed(66)
df=read.table("X_metrop.dat",header=TRUE)
inp=a...
2006 Aug 02
1
RE
...s, but my question is when
I get the XD matrix as explained below the covariate matrix (X) which
initially had 3 columns now changed to a 9 rows and one column and two
additional variables X4 and X5 as explained belwo which I guess for the beta,
what is the role of these variables in the 'rmnlIndepMetrop'analysis.
example
x1<-runif(3,-1,1)
x2<-runif(3,0,1)
x3<-runif(3,10,50)
X<-cbind(x1,x2,x3)
X
x1 x2 x3
[1,] -0.9701396 0.4084203 41.31097
[2,] 0.3844539 0.4791997 36.85861
[3,] 0.2732056 0.5433642 13.14610
Xa<-cbind(X,-X)
XD<-createX(p=3,na=2,nd=NUL...