similar to: No subject

Is there any way of getting just the LaTex output which R CMD Rd2dvi generates on the way to produces a .dvi file, which could then be incoorporated into another document? Thanks, Jonathan -- Dr. Jonathan Myles e-mail:jonathan.myles at mrc-bsu.cam.ac.uk MRC Biostatistics Unit Tel. 01223 330372 Institute of Public Health FAX 01223 330388 University Forvie Site

Hi Does anyone know of any R coding/functions for MCMC approaches? I am currently using BUGS but I wonder if the bazaar has produced anything? I think I am pushing BUGS to it's limit and possibly past it at the moment. John -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Thanks to everyone who helped with my query---Professors Riply, Bates, Lumley, and Tim Keitt. I said <snip> > My worries are > > a) That I've used calloc rather than S_alloc throughout. Professors Bates and Ripley both pointed out that calloc should be replaced by Calloc and free by Free, which, I'll do. ---- I said: > b) That I've used the following

Hi all, Thanks to Prof. Ripley, Prof. Gentleman, and Simon Wood (did I miss anyone?). The problem seemed to have gone away. Everyone suggested using some malloc debugger (such as Electric Fence). All I did was following half of what BDR suggested below, i.e., changing all the S_alloc() calls to Calloc() and Free(). I didn't get to try efence, and the problem seems to have disappeared! As

Hello, thank you very much for your previous answers about the C++ code. I am interested in the application of the Gibbs Sampler in the IRT models, so in the function MCMCirt1d and MCMCirtkd. I've found the C++ source codes, as you suggested, but I cannot find anything about the Gibbs Sampler. All the files are for the Metropolis algorithm. Maybe I am not able to read them very well, by the

Hello, I am a chemical student and I make use of 'lme/lmer function' to handle experiments in split-plot structures. I know about the mcmcsamp and I think that it's very promissory. I would like knowing "the concept behind" of the mcmcsamp function. I do not want the C code of the MCMCSAMP function. I would like to get the "pseudo-algorithm" to understanding that

Dear R users, I'm trying to carry out monte carlo integration of a posterior density function which is the product of a normal and a gamma distribution. The problem I have is that the density function always returns 0. How can I solve this problem? Here is my code #generate data x1 <- runif(100, min = -10, max = 10) y <- 2 * x1^2 + rnorm(100) # # # # # # # # Model 0 # # # # # # #

I am trying to get our database interface (PADI) working with R. The code does some memory allocation and for Splus there is an ifdef which controls whether malloc or S_alloc is called. I did nm R.binary | grep alloc to see if this was supported and I find there are some choices: [2490] | 446036| 272|FUNC |GLOB |0 |7 |R_alloc [806] | 693956| 312|FUNC |LOCL |0 |7

Tried to look for package GLLMGibbs-(package for fitting mixed models by Gibbs samppling contributed by Myles' and Clayton) in CRAN but not found it, Does anybody know if it is still available? thanks Oarabile

An embedded and charset-unspecified text was scrubbed... Name: n?o dispon?vel URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110405/06486a8a/attachment.pl>

I am trying to write a C function to create a vector of integers that can be used by the R calling function. I do not know the size of the vector in the R calling function. (Well, actually, I have an upper limit on the size, but that is so large that R cannot allocate it. What I'm doing in the function is to do a sieving procedure, and the result will be small enough to fit into my

Where are the archives of r-help kept, please? Jonathan -- Dr. Jonathan Myles e-mail:jonathan.myles at mrc-bsu.cam.ac.uk MRC Biostatistics Unit Tel. 01223 330371 Institute of Public Health FAX 01223 330388 University Forvie Site Robinson Way CAMBRIDGE CB2 2SR -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list --

Is there any way of running R in ESS using -n and -v arguments other than the default? Thanks Jonathan -- Dr. Jonathan Myles e-mail:jonathan.myles at mrc-bsu.cam.ac.uk MRC Biostatistics Unit Tel. 01223 330371 Institute of Public Health FAX 01223 330388 University Forvie Site Robinson Way CAMBRIDGE CB2 2SR

Dear R-users (and developers), I am looking for an efficient framework to carry out parameter estimations based on MCMC (optionally with specified priors). My goal is as follow: * take ANY R-function returning a likelihood-value (this function may itself call external programmes or other code!) * run a sampler that covers the multidimensional parameter space (thus creating a posterior

Dear R-devel, Last week I got several responses to my question about debugging segfaults in my code (original post below). After I changed the S_alloc() calls to Calloc()/Free(), the symptom was gone, but I was told to keep looking. So I did: o Switched to Calloc/Free. Electric Fence did not find any problem. o Put assert(index < bound); assert(index >=0); everywhere in the C routine

I'm using R for calculations and plots in a statistical paper. Is there a canonical reference to use in the bibliography? (Or is it not common to acknowledge computer programs in journal papers?) -- Bj?rn-Helge Mevik <bhm at math.uio.no> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Does anyone know what exactly S_alloc() does? (Yes, it's porting time again ...) Thanks, -k =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To:

Dear all, I'm currently experiencing big troubles with my C code called by .Call in a R function. I know this may not be the right place for such things but these errors are driving me crazy. And I hope other people may have already experienced these problems so that they could give me good suggestions. Sorry if I'm completely out of topic. Well let's go... I'm writing a

I would like to estimate a 95% highest density area for a multivariate parameter space (In the context of anova). Unfortunately I have only experience with univariate kernel density estimation, which is remarkebly easier :) Using Gibbs, i have sampled from a posterior distirbution of an Anova model with k means (mu) and 1 common residual variance (s2). The means are independent of eachother, but

Hi, I am trying to figure out the observed acceptance rate and M, using generalised rejection sampling to generate a sample from the posterior distribution for p. I have been told my code doesn't work because I need to "take the log of the expression for M, evaluate it and then exponentiate the result." This is because R is unable to calculate high powers such as 545.501. As