similar to: Re: Bayesian stats

Hi, I'm using R 1001 for Windows NT and the sm library. I'm trying to create plots for my data set like Bowman and Azzalini have in Figure 1.8 (p. 9) of their book for my data (i.e. a contour plot for each group in my data set and its all plotted on 1 plot). The problem I'm having is that R is not drawing closed contour lines for each group. Sometimes it does; other times it

Is there an easy way, for MCLUST in particular, to force its graphical output to all black and not use colors to show the clusters/groups? It shows the groups quite nicely by using triangles, crosses, etc. Unfortunately, my printer doesn't convert reds and yellows too nicely and the result is a indistinct figure on the page. (Yes, I have access to a color printer but my intended publication

Hi all. Are there any R functions around that do quick logistic regression with a Gaussian prior distribution on the coefficients? I just want posterior mode, not MCMC. (I'm using it as a step within an iterative imputation algorithm.) This isn't hard to do: each step of a glm iteration simply linearizes the derivative of the log-likelihood, and, at this point, essentially no

Dear Friends, According to Gelman et al (2003), "...Bayesian P-values are defined as the probability that the replicated data could be more extreme than the observed data, as measured by the test quantity p=pr[T(y_rep,tetha) >= T(y,tetha)|y]..." where p=Bayesian P-value, T=test statistics, y_rep=data from replicated experiment, y=data from original experiment, tetha=the function

Dear friends, Over the past weeks, I have been asking a lot of questions about how to use R in Bayesian analysis. I am brand new to R, but I am very pleased with it. I started with winbugs but I found winbugs to be a limited software, not bad but has several limitations. By contrast, R allows the analyst to tackle any problem with a huge set of tools for any kind of analysis. I love R. In

Hi, I am trying to create what are called "spider diagrams" in the geochemical literature using R. A spider diagram is basically a plot of the atomic number versus the concentration on a log scale. Lines are drawn from each atomic number for each sample. Right now, my data frame looks like: SAMPLE SITE V3 LA.NASC LAATNUM CE.NASC CEATMNUM 1 1A:001

On Fri, May 31, 2019 at 01:28:55AM -0400, Scott Kostyshak wrote: > The attached patch adds some sanity checks to the "type" argument of > quantile(). Output from the following commands show the change of > behavior with the current patch: > > vec <- 1:10 > quantile(vec, type = c(1, 2)) > quantile(vec, type = 10) > quantile(vec, type = "aaa")

I have downloaded the Umacs (Universal Markov chain sampler) and submitted the following sample code from Kerman and Gelman.   s <-Sampler( J=8, sigma.y  =c(15,10,16,11,9,11,10,18),           y  =c(28, 8,-3,7,-1,1,18,12),      theta =Gibbs(theta.update,theta.init),           V =Gibbs(V.update,mu.init),         mu =Gibbs(mu.update,mu.init),         tau =Gibbs(tau.update,tau.init),       

I am new to R (ex SAS user) , and I cannot merge two data frames without getting duplicated rows in the results. How to avoid this happening without using the unique() function? 1. First data frame is called "tmv" with 6 variables and 239 rows: > tmv[1:10,] temps nom prenom sexe dist style 1 01:59:36 Cyr Steve H 45 free 2 02:09:55 Gosselin

Is there a way to mark a ticket as a potential regression in the bug tracker? I think the following issue is a regression: https://bugs.r-project.org/bugzilla/show_bug.cgi?id=17684 I've just tested (2020-06-12 r78687) and what I believe to be a regression is still there. I don't think the bug has bitten many people, so I don't think it is critical, but often it is helpful to mark

Hello Aki, Timo, according to git this fix should be in 2.2.27, which I'm running, so I guess this isn't it or something else is missing. See: https://bugs.debian.org/cgi-bin/bugreport.cgi?bug=859700 Regards, Christian On Thu, 6 Apr 2017 15:37:33 +0900 Christian Balzer wrote: > Hello, > > On Thu, 6 Apr 2017 09:24:23 +0300 Aki Tuomi wrote: > > > On 06.04.2017 07:02,

Hello, this is on Debian Jessie box, dovecot 2.2.27 from backports, imap-hibernate obviously enabled. I've been seeing a few of these since starting this cluster (see previous mail), they all follow the same pattern, a user who accesses their mailbox with both POP3 and IMAP deletes mails with POP3 and the IMAP (imap-hibernate really) is getting confused and upset about this: --- Apr 6

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

Dear Friends, I am trying to use the gelman-rubin convergence test. I generated a matrix samp[10,000x86] with the gibbs sampler. the test requires the creation of "mcmc" objects. Since I don't know how to define samp as a "mcmc" object, I tried to create one mcmc object by means of the mcmc() function. With this function I tried to create a mcmc object dul from samp but I

Hello, On Thu, 6 Apr 2017 20:04:38 +0900 Christian Balzer wrote: > Hello Aki, Timo, > > according to git this fix should be in 2.2.27, which I'm running, so I > guess this isn't it or something else is missing. > See: > https://bugs.debian.org/cgi-bin/bugreport.cgi?bug=859700 > Any update on this? Christian > Regards, > > Christian > > On Thu, 6

Thanks for the info. I do have one further question for you. On your servers that are currently handling 50k IMAP sessions, how many users does that correspond to? Since many users will have multiple IMAP sessions on multiple devices, I?d like to hear about some real-world numbers that could be used for budgeting a new project like mine. Also, do you use Dovecot IMAP proxies in front of your

Hello, On Thu, 6 Apr 2017 09:24:23 +0300 Aki Tuomi wrote: > On 06.04.2017 07:02, Christian Balzer wrote: > > > > Hello, > > > > this is on Debian Jessie box, dovecot 2.2.27 from backports, > > imap-hibernate obviously enabled. > > > > I've been seeing a few of these since starting this cluster (see previous > > mail), they all follow the

On 06.04.2017 07:02, Christian Balzer wrote: > > Hello, > > this is on Debian Jessie box, dovecot 2.2.27 from backports, > imap-hibernate obviously enabled. > > I've been seeing a few of these since starting this cluster (see previous > mail), they all follow the same pattern, a user who accesses their mailbox > with both POP3 and IMAP deletes mails with POP3 and

thanks for your help happy to say that the performance dramatically improved after i use the high performance settings from here http://wiki.dovecot.org/LoginProcess grep Login: /var/log/mail.log.1 |wc -l with the mail.log being of a typical, busy day. 412992 i also picked up the imap and pop3 connections during peak hours [root at ns1 domains]# doveadm who | awk

Hi Scott, This question is appropriate for the r-help mailing list, but probably off-topic here on r-devel. Best, Ista On Tue, May 1, 2018 at 2:57 PM, Scott Kostyshak <skostyshak at ufl.edu> wrote: > I have very little knowledge about file encodings and would like to > learn more. > > I've read the following pages to learn more: > >