search for: mcmcpack

I'm completely new to R and am trying to install an add-on package for one of our faculty members. I've had no problems with most of them, but am completely stuck trying to figure out why MCMCpack won't install.. any help is VERY much appreciated!!! here's what I get when I try to install the package: * Installing *source* package 'MCMCpack' ... checking for C++ compiler default output file name... a.out checking whether the C++ compiler works... yes checking whether we are...

I installed R-2.3.1 and ran make check without problem, but I am having trouble installing several packages using gcc (GCC) 3.2.3 20030502 (Red Hat Linux 3.2.3-53) on Red Hat Enterprise Linux AS release 3 (Taroon Update 7) Kernel 2.4.21-40.ELsmp on an x86_64 Below are the messages from MCMCpack. Perhaps my systems are messed up? Paul Gilbert > install.packages("MCMCpack") Warning in install.packages("MCMCpack") : argument 'lib' is missing: using /home/mfa/gilp/toolchain/R/asd05546/R-2.3.1/site-library --- Please select a CRAN mirror for use in this sessi...

Hi list, I tried to install MCMCpack to R-2.12.0, got the following error, CC -m64 -library=stlport4 -I/apps/sparcv9/R-2.12.0/lib/R/include -DSCYTHE_COMPILE_DIRECT -DSCYTHE_DEBUG=0 -DHAVE_TRUNC -DHAVE_IEEEFP_H -I/opt/csw/include -KPIC -g -c MCMCSVDreg.cc -o MCMCSVDreg.o "error.h", line 598: Error: The function "...

I'm having some problems in installing some source packages under linux. As an example, MCMCpack. An error is raised when linking: > install.packages("MCMCpack") [...] * Installing *source* package 'MCMCpack' ... checking for C++ compiler default output file name... a.out checking whether the C++ compiler works... yes checking whether we are cross compiling... no checking...

Hello Everyone Both the MCMCpack and the bayesm libraries allow us to make draws from the Inverse Wishart distribution. But I wanted to find out how exactly is the Inverse Wishart distribution parameterized in these libraries. The reason I ask is the following: Now its generally standard to express Inverse Wishart as IW(0.5 *...

Dear All I try to use your R package MCMCpack and I have encountered the following problem: The following code works fine: library(MCMCpack) line <- list(X = c(-2,-1,0,1,2), Y = c(1,3,3,3,5)) posterior1 <- MCMCregress(X~Y, data=line) summary(posterior1) But as long as I try the following lines library(MCMCpack) li...

Hi, I'm having issues with what I believe is a bug in the MCMCpack's MCMCmetrop1R function. I have code that basically looks like this: posterior.sampler <- function(data, prior.mu){ log.posterior <- function(theta) log.likelihood(data, theta) + log.prior(prior.mu, theta) post.samples <- MCMCmetrop1R(log.posterior, theta.init=prior.mu, burnin=1...

Hi there, I am running MCMCpack (MCMCirt1d model) on some files (26 items x about 800 – 1200 individuals). I have a problem. When I am working on “big” files, the R program crashes. More precisely I got the following Microsoft Warning: ******************************************* Microsoft Visual C++ Runtime Library...

Dear R users, Thanks in advance. I am Deb, Statistician at NSW Department of Commerce, Sydney. I am using R 2.9.1 on Windows XP. This has reference to the package “MCMCpack”. My objective is to select a better model using various alternatives. I have provided here an example code from MCMCpack.pdf. The matrix of Bayes Factors is: model1 model2 model3 model1 1.000 14.08 7.289 model2 0.071 1.00 0.518 model3 0.137 1.93 1.000 &gt...

Hi, I've been testing some models with the MCMCpack library. I can run the process and get a nice model "object". I can easily see the summary and even plot it. I can't seem to figure out how to: 1) Access the final coefficients in the model 2) Turn the coefficients into a model so I can then run predictions using them. A summary...

Hi, I've been having problems getting MCMCpack installed on my Solaris 10 x86 system. I've installed gcc from opencsw.org, and, for the most part, other R packages install nicely. >From what I can gather, the build fails due to a mismatch of include files. I *think* it wants to use the include files that are contained within the MCMCpa...

Hi R users, I am using the package MCMCpack and I would compute the Deviance Information Criterion for selecting models, instead of using bayes factors. I know that a similar question have been sent previously to the list by others, but I have found no response. Can somebody help me? best Giacomo

Hi everyone, I am running an ordinal probit using the Bayesian MCMCpack and I am getting an error saying "attempt for find suitable starting values failed" Here is my code: > posterior <- MCMCoprobit(y ~ x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8 +x9 + x10 + x11 + x12 +x13 , beta.start=c(-10, 0.05, 0.02, 0.04, 0.98, 0.61, -0.29, 0.91, -0.82, 1.34, 0.68, 0...

Compilation of MCMCpack under freebsd 6.2 i386 fails because of the following cryptic error: * Installing *source* package 'MCMCpack' ... checking for C++ compiler default output file name... a.out checking whether the C++ compiler works... yes checking whether we are cross compiling... no checking for suffix o...

[This email is either empty or too large to be displayed at this time]

R users- In moving from R 2.0.0 to R 2.1.0 in Windows, I have encountered a problem with the "riwish" command in the package "MCMCpack". I've searched the documentation and can't seem to figure it out. For example: Define a matrix: > lam <- matrix(c(.00233,-.00057,-.00057,.00190),2,2) and then use the riwish command to generate a random inverse-Wishart variate: > riwish(72,lam) In version 2.1.0, the r...

Hello, just to follow up a question from last week. Here what I've done so far (here an example): library(MCMCpack) Y=c(15,14,23,18,19,9,19,13) X1=c(0.2,0.6,0.45,0.27,0.6,0.14,0.1,0.52) X2a=c(17,22,21,18,19,25,8,19) X2b=c(22,22,29,34,19,26,17,22) X2 <- function()runif(length(X2a), X2a, X2b) model1 <- MCMCregress(Y~X1+X2()) summary(model1) but I am not sure if my X2-function is working in the MCMCpack...

...nd(1,x) > p <- ncol(c.df) > > # marginal log-prior of beta[] > logpriorfun <- function(beta, mu, gshape, grate) + { + logprior = -p*log(2) + log(gamma(p+gshape)) - log(gamma(gshape)) + + gshape*log(grate) - (p+gshape)* log(grate+sum(abs(beta))) + return(logprior) + } > require(MCMCpack) Loading required package: MCMCpack Loading required package: coda Loading required package: lattice Loading required package: MASS ## ## Markov Chain Monte Carlo Package (MCMCpack) ## Copyright (C) 2003-2007 Andrew D. Martin and Kevin M. Quinn ## ## Support provided by the U.S. National Science Fo...

Hello All, I'm trying to run a one-dimenional irt model using the packages MCMC and coda on a rather large set of roll-call voting data with many missing observations. Here's a sample of the code: Post10<- MCMCirt1d (Italy10, burnin = 1000, mcmc=50000, thin=100, verbose=TRUE, theta.constraints = list(V549=1, V443=-1)) The MCMCirt1d command seems to work fine, but when I try to

Hi, A question. When I run gelman.diag and gelman.plot with mcmc lists obtained from MCMCregress, the results are following. > post.R <- MCMCregress(Size~Age+Status, data = data, burnin = 5000, mcmc = 100000, + thin = 10, verbose = FALSE, beta.start = NA, sigma2.start = NA, + b0 = 0, B0 = 0, nu = 0.001, delta = 0.001) > post1.R <- MCMCregress(Size~Age+Status, data