similar to: Efficiency: whole vector and for loop methods

Hello, I am doing mcmc=10000 simulations from a posterior distribution of the parameters of a mixture of K=6 normal densities. I have mcmc by K matrices simMeans, simVars and simWeights containing the simulation output: one row for each simulation, one column for each normal component of the mixture. One thing I would like to do is a plot of the posterior predictive density. In order to do that

mtex[3]<-"I need\'s a \"double quote\" with no backslash" mtex[3] [1] "I need's a \"double quote\" with no backslash" so how is it done? Thanks in advance, graham lawrence _________________________________________________________________ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list --

Hi Graham, Your guess is right. The event loop stuff is handled differently on Windows than on Unix-like OSes. There is a hook, called tcl_do that is hard-coded into the R Windows GUI. It is called whenever the console is idle. RGtk2 is able to hook into that to process GTK+ events. The problem is that the two packages replace each other's handler. So in this case tcltk is blocking GTK+ from

Dear Team, Don't understand output below. Please advise. Windows98 and rw1041(patched) from 1/25/02. Graham Lawrence C:\>PATH=C:\bin;C:\mingw32;C:\Perl;C:\HtmlWork;C:\R\rw1041\bin;C:\Tcl;C:\Tcl\bin C:\>Rcmd INSTALL C:\R\lattice_0.4-0.tar.gz Bad command or file name C:\>Rcmd C:\>Rcmd Install Bad command or file name C:\>dir C:\R\rw1041\bin Volume in drive C has no

Hi, Is there any documentation to explain R's error diagnostics, e.g. messages like >source("../src/sources/routine.R") Read 8 records Read 2 records Error in "[<-"(*tmp*, !test, value = rep(no, length = length(ans))[!test]) : incompatible types The nature of the problem incompatible types, is clear; but I'm more interested in its location, and feel

I've been using R for several months now, in fact it was my intro to "real" computer programming. It used to frustrate me no end, but lately I've noticed my reaction has changed. Nowadays its much more like "wow, is that cool, or what?". So I guess I've gotten over the hump of the learning curve. Comparing what I now do with R to my previous hodge-podge of

Hi all, I'm curious if anyone has coded an Expectation-Maximization algorithm that could help me model some circular data I have. I'd like to model it as a mixture of uniform and Von Mises centered on 0, so the only free parameters is the mixing proportion and the kappa of the Von Mises. I couldn't find anything in the contributed packages that seemed to suit this purpose. Any

Dear R-help, I have a series of data.frames (i1,i2,...,in) all containing the same number of items, but dissimilar content. I would like to instal them in another data.frame, say index, so that I can access their items with index[i,j]. No matter how I try to set index up, its subframes cannot be indexed, because they all have the row number of 1. What am I doing wrong? (I get the same

> ssh has to set NONBLOCK otherwise it can, well, block - there's > no way for ssh to know a priori how much data it can write to a fd. I don't know anything about how ssh is structured, but I think it must be a bit more complicated than that. Ssh only sets O_NONBLOCK on an fd if isatty(fd) returns false, so it's able to function with blocking input and output if the relevant

Hello, I would like to please ask if someone would explain how r reads characters and numbers differently. Using read.csv, I had a matrix that resembled the following, only with many more ids and data: ID Visit variable 2 1 5 2 1 3 2 3 4 2 41 1 2 42 34 2 5 54 2 9 1 2 10 3 2 12 5 5 1 54 5 2 9 5 3 3 5 41 54 5 41 2 5 5 235 5 9 4 5 10 2 5 12

Dear R-help, Why does the second grep reject when the first grep accepts? Why does regexpr accept what grep rejects? >eqtn [1] "(1.2*A%*%B)/2" >if(grep("A\\%\\*\\%B",eqtn)==1)erind<-1 >if(grep("A\\%\\*\\%A",eqtn)==1)erind<-1 Error in if (grep("A\\%\\*\\%A", eqtn) == 1) erind <- 1 : missing value where logical needed

Dear list, A while ago, Vadim asked opinions on improving efficiency of sample() with prob, e.g. sample with replacement with weight. ( https://stat.ethz.ch/pipermail/r-devel/2004-September/030844.html ) He did not post what he ended up with this problem though. I am having exactly the same problem. I need to sample with replacement from a population of size 10 million with fitness values for

Chris, nice segue to Swift ! :-), but... The question is what should LLVM do with UB in general, saying that we are going to change one specific idiom from undefined to defined glosses over the real question: why should we ever optimize / delete any UB at all ? This “depressing and faintly terrifying thing” as you call it, should be viewed not as an opportunity for optimization, but

Dear R-help, I would like to be able to run the following code sequence as a source routine. If I paste it into R via the clipboard it works as expected, but if I source the code instead then the last 3 statements fail. I've also tried writing to the file in place of the sink sequence, but that also hits a snag. R 1.6, w98e2, dfr is a data frame containing the content of Julian

Dear R-Help, Searching for a function that will copy the value of an R object to the windows clipboard led me to "tkclipboard.append(...)". Will this function do what I seek? Or if not, does anyone know of such an R function? If tkclipboard.append(...) is what I need, please explain how I can set TCL_LIBRARY. With R 1.6 on W98 I got the following R dialog. (If it is relevant,

Hi all, The survivals functions can be tested by the Log-rank test and others, for example the Gehan-Breslow. The graham breslow work with the alpha values. But I don't know how is the Gehan-Breslow test with R. Somebody know a type function?.. or other suggestions? Any help will be really appreciated Jos? Bustos Marine Biologist Master Apllied Stat Program University of Concepci?n

Hi rsync team, I thought that rsync would try to overlap computing and IO on both machines. I'm rsyncing a large tree (18G) and am keeping an eye on that. Suddenly my side completely stopped. No IO visible, no CPU time spent. The otherside was doing 100% CPU. Then the other side started to do disk IO. Then suddenly the activities moved over to my side, and I saw things moving again in the

Hi, Here's the updated RFC for representing scalable vector types and associated constants in IR. I added a section to address questions that came up on the recent patch review. -Graham =================================================== Supporting Scalable Vector Architectures in LLVM IR =================================================== ========== Background ========== *ARMv8-A

Hi, Are there any objections to going ahead with this? If not, we'll try to get the patches reviewed and committed after the 7.0 branch occurs. -Graham > On 2 Jul 2018, at 10:53, Graham Hunter <Graham.Hunter at arm.com> wrote: > > Hi, > > I've updated the RFC slightly based on the discussion within the thread, reposted below. Let me know if I've missed

Hi Graham, Thanks for your kind explanation. There was internal discussion about it. If possible, can you let me know the Clang/LLVM CodeGen patches for the vector type on phabricator please? I would like to check what kinds of the restrictions the type causes on Clang/LLVM. Thanks, JinGu Kang ________________________________ From: Graham Hunter <Graham.Hunter at arm.com> Sent: 28 May