if (FALSE) {
##load car data
data(car)
## compute 5 runs of disjunctive model with 2 features
carem2<-plfm(maprule="disj",freq1=car$freq1,freqtot=car$freqtot,F=2,M=5)
## Compute a sample of the posterior distribution
## for the disjunctive model with two features
## use the posterior mode obtained with the previous plfm analysis
carbayes2<-bayesplfm(maprule="disj",freq1=car$freq1,freqtot=car$freqtot,F=2,
maxNiter=500,Nburnin=0,Nstep=100,Nchains=2,
start.bayes="fitted.plfm",fitted.plfm=carem2)
## compute a summary of the object generated by bayesplfm
summarycarbayes2<-summary(carbayes2)
}
Run the code above in your browser using DataLab