Displaying 3 results from an estimated 3 matches for "mcmc1".
Did you mean:
mcmc
2007 Aug 07
1
bayesm - question about 'rscaleUsage function'
...ode 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........
2011 Aug 04
0
Problems with Z in rhierMnlRwMixture using bayesm
...],nd=NULL,Xd=NULL,INT=TRUE,base=1)
lgtdata[[ind]]=list(y=y,X=X,id=id[i]); ind=ind+1
}
nlgt=length(lgtdata)
Z=NULL
nlgt=length(lgtdata)
for(i in 1:nlgt) {
Z=rbind(Z,Demog[Demog[,1]==lgtdata[[i]]$id,3:9])
}
Z=log(Z)
Z[,1]=Z[,1]-mean(Z[,1])
Z[,2]=Z[,2]-mean(Z[,2])
keep=5
R=5000
mcmc1=list(keep=keep,R=R)
out=rhierMnlRwMixture(Data=list(p=p,lgtdata=lgtdata,Z=Z),Prior=list(ncomp=1),Mcmc=mcmc1)
Many thanks for all your help and time.
Regards,
Elcin
--
View this message in context: http://r.789695.n4.nabble.com/Problems-with-Z-in-rhierMnlRwMixture-using-bayesm-tp3717887p37178...
2007 Mar 30
0
problem using mcmcsamp() with glmer models containing interaction terms in fixed effects
...I(sqrt(counts))~as.factor(visit)|region, data=data1)
# two lmer models
lmer1 <- lmer (counts~visit+visit:zone + (1|region), data=data1,
family=poisson(link="log"))
lmer2 <- lmer (counts~visit+visit:zone + (visit+visit:zone|region),
data=data1, family=poisson(link="log"))
mcmc1 <- mcmcsamp(lmer1, n=10000) # works
mcmc2 <- mcmcsamp(lmer2, n=10000) # doesn't work
[[alternative HTML version deleted]]