similar to: Summary: Faster way for weighted matching

For each element in w I want to find a good match (subscript number) of an element in x. x and w can be long. Instead of just finding the closest match I want to use weighted multinomial sampling (which I've already figured out once I have the probabilities) where the probabilities come from the tricube function of absolute differences between donor and target values, but normalized to sum

Hi Johannes, I came up with the following (assuming that if step[i] is independent of step[i-1] then it is also independent of -step[i-1]): Cheers Joe > # your (unbounded) random walk > k <- cumsum(c(0,sample(c(-1,1), 1000, rep = T))) > > #shift it to positive, to allow following calculations k <- k - min(k) > # set your bound (remember min(k) is now zero) wid <- 10

Hi everyone, I would like to identify the case by groups that is just bigger that avg plus sd. For example, using species as group and petal.wid as my variable in the iris data. What's the better way to doit? creating a function? So,the question is to identify the single element of each species that is just larger than a cut-off point (i.e. larger than mean + sd) I made this, but I can not

Hi Everyone, I would like to add errorbars to a lattice plot. I already have a function that adds error bars to a plot (originally written by George Gilchrist) which takes as arguments the coordinates of the point and then the size of the error bar. I would like to integrate it in my lattice code but I seem to be stumped. I can put the error bar function inside my code as if it were a panel

Ubuntu 8.04 Wine 1.1.10 The sonic mega collection contains almost every genesis sonic game, along with some game gear sonics. The game installed flawlessly, and the no-cd patch worked also. However, upon loading, it shows a faint shadow of the famous SEGA startup screen and then locks up and must be force killed. Changing it to various compat modes did nothing. Normal Output: Code: name at

Hi, I have encountered this error when attempting a One-way Repeated-measure ANOVA with my data. I have read the "Anova in car: SSPE apparently deficient rank" thread by I'm not sure the within-subject interaction has more degrees of freedom than subjects in my case. I have prepared the following testing script: rm(list = ls())

Prasad wrote: > I wrote a function in R which uses tcltk package .... essentially I wanted > to give within that function, a widget with 2 radiobuttons to choose > between plotting Precip and Temperature plots. After the user has chosen > one of the radiobuttons there is another widget that asking him to identify > outliers. However, I am having a lot of problems...what R does is

Hello, I wrote a function in R which uses tcltk package .... essentially I wanted to give within that function, a widget with 2 radiobuttons to choose between plotting Precip and Temperature plots. After the user has chosen one of the radiobuttons there is another widget that asking him to identify outliers. However, I am having a lot of problems...what R does is evaluate the whole function

Peter Thank you for the information. I accidentally deleted Ken's post without having read it. Ken' s thought is great but as you said awful to implement I thought that capture.output would come in handy some time when I first saw it on an unrelated reply. Just thought :- the latex listings package may have alternatives If I remember correctly it has wrapping and other goodies but I

PR #7690 points out that if the confidence intervals (+/-1.58 IQR/sqrt(n)) in a boxplot with notch=TRUE are larger than the hinges -- which is most likely to happen for small n and asymmetric distributions -- the resulting plot is ugly, e.g.: set.seed(1001) npts <- 5 X <- rnorm(2*npts,rep(3:4,each=npts),sd=1) f <- factor(rep(1:2,each=npts)) boxplot(X~f) boxplot(X~f,notch=TRUE) I can

Hi all, When compiling some programs using the Mips backend i'm getting this assert message on lib/CodeGen/LiveInterval.cpp:227: "Range is not entirely in interval!" I don't know yet if it's something that is missing on the backend code or why the range to be removed it outside the interval, does anyone have any clue? A more detailed output is attached. The program i tried

Dear useRs, Is there a way to 'properly' format %d when plotting more than one page on png()? 'Properly' means to me with leading 0s, so that the PNGs become easy to navigate in a file/image browser. Lacking a better solution I ended up using the code below, but would much prefer something like png("test_%d.png",bg="white",width=1000,height=700) where %d

Here is a test case: extern double floor(double); extern double floor_(double); double x = 1.5; double y, y_; void foo() { double y = floor(x); double y_ = floor_(x); } If I compile this for Mips16, it calls the proper helper function for floor_ but not for floor, because the signature for floor in callee info is wrong. Args[0] = void RetTy = void

Both of these bugs are not a real harm, however, they have been annoying me for too long ... ;-) 1) print(2^30, digits = 12) #- exponential form; unnecessarily! formatC(2^30, digits = 12) #- shows you what you'd want above ## S-plus is okay here; note that the problem also affects ## paste(.) & format(.) : options(digits=10) paste(2^(4*1:8)) S-plus gives [1] "16"

My destination is :/Example(folder) / ITD(folder) / LW(folder) ,user(lw) of group(ITD) can enter and do everything(rwx) in LW(folder),the following is ACL settings On Freebsd: ----------------------- /Example(folder) : # owner: root # group: wheel user::rwx group::r-x other::r-x --------------------- /Example(folder) / ITD(folder): # owner: root # group: wheel user::rwx group::r-x group:ITD:r-x

Does it make a difference if you give the "-target" option to clang? $ clang -target mips-linux-gnu ch8_3.cpp -o ch8_3.bc -emit-llvm -c The .s file generated this way looks quite different from the one in your email. On Tue, Feb 19, 2013 at 5:06 PM, Jonathan <gamma_chen at yahoo.com.tw> wrote: > I didn't have Mips board. I compile as the commands and check the asm >

I'm just passing this patch on to the list as I know Brian is busy and I had proded him to file the bug report. http://llvm.org/bugs/show_bug.cgi?id=14721 "binutils on OpenBSD is quite old, which causes the build on mips64el to fail on MipsJITInfo.cpp - it complains of illegal operands. The attachment patches it so that it will build. All the patch does is replace the register aliases

I didn't have Mips board. I compile as the commands and check the asm output as below. 1. Question: The distance of caller arg[4] and arg[5] is 4 bytes. But the the callee get every arg[] by 8 bytes offset (arg_ptr1+8 or arg_ptr2+8). I assume the #BB#4 and #BB#5 are the arg_ptr which is the pointer to access the stack arguments. 2. Question: Stack memory 28($sp) has no initial value. If

i have this code typedef long long L; typedef integer I; void test2(auto L p0, auto L p1, auto L p2, auto L p4, ...) { va_list select; va_start (select, p4); report( va_arg(select,L) ); report( va_arg(select,I) ); report( va_arg(select,L) );

Could you file a bug? On Wed, Feb 27, 2013 at 2:47 AM, Stavropoulos Nikos < n.stavropoulos at think-silicon.com> wrote: > i have this code > > typedef long long L; > typedef integer I; > void test2(auto L p0, auto L p1, auto L p2, auto L p4, ...) > { > va_list select; > > va_start (select, p4); > > report( va_arg(select,L) ); > report(