search for: axion

Hi all, I've a vector with entries, which are all of the same type, e.g. string: k <- c("bb", "bb", "bb", "aa", "cc", "cc") and want to create a second vector containing the number of each entry in k in the same order as in k, i.e. c(3, 1, 2) or: k <- c(5,5,5,5,2,2,4) => c(4,2,1) thanks

Hi all, I have a data frame with "clustered" rows as follows: Cu1 x1 y1 z1 ... Cu1 x2 y2 z2 ... Cu1 x3 y3 z3 ... # end of first cluster Cu1 Cu2 x4 y4 z4 ... Cu2 x5 y5 z5 Cu2 ... # end of second cluster Cu2 Cu3 ... ... "cluster"-size is 3 in the example above (rows making up a cluster are always consecutive). Is there any faster way to sample n clusters (with

...2. Re: Running rtest - how to/ help (indian scorpio) > 3. Re: how to read .sps (SPSS file extension)? (Eik Vettorazzi) > 4. Re: how to read .sps (SPSS file extension)? (Eik Vettorazzi) > 5. Re: plotting density for truncated distribution (Chris Andrews) > 6. construct a vector (axionator) > 7. Re: construct a vector (Marc Schwartz) > 8. Re: construct a vector (Richard.Cotton@hsl.gov.uk) > 9. Re: multiple imputation with fit.mult.impute in Hmisc - how > to replace NA with imputed value? (Frank E Harrell Jr) > 10. Re: construct a vector (Wacek Kusniercz...

Hi, I want to draw a contour plot of the following function: z = y*x + epsilon, where x ~ N(y, 1) and epsilon ~ N(0, sigma) with sigma fixed (e.g. 1) But didnt manage to feed "contour" with the right input. Thanks for your help. Armin

...t this is more of a client problem than a Samba/UNIX one but I can't prove it. BTW, the print job which causes the problem always comes out OK. Has anyone seen this phenomenon before and if so, has anyone got a fix/workaround for it? Any information would be appreciated. Ken Hamer khamer@axion.bt.co.uk

Ken Hamer (B83 131) <K.Hamer@axion.bt.co.uk> wrote > BUT...... After printing a document, Samba appears to hang up the connection, > although another seems to start pretty quickly. ... > server.netbios- client.1065 8557 0 8760 0 CLOSE_WAIT > > The connection will stay in this state until...

Hi, I would like to call R's lbfgsb function from my C/C++ code by including R_ext/Applic.h and linking against libR. Currently, I am allocating memory for x (and the other input arrays for lbfgsb) in my C/C++ code via malloc/new. However, this gives a segmentation fault when executing the program. I tried to allocate x via PROTECT(x = NEW_NUMERIC(n)); x_p = NUMERIC_POINTER(x);. This compiles

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

...stion is, is this a safe option to use? I get the impression that it is, but I'm not too familiar with NT's locking and need some reassurance before trusting our data to it. Thanks! Andrew Return-Path: <phil.hardern@bt.com> Delivered-To: samba@lists.samba.org Received: from gollum.axion.bt.co.uk (gollum.axion.bt.co.uk [132.146.17.41]) by lists.samba.org (Postfix) with ESMTP id E16B0463C for <samba@lists.samba.org>; Fri, 15 Jun 2001 05:53:30 -0700 (PDT) Received: from cclmsent02.lon.bt.com by gollum (local) with ESMTP; Fri, 15 Jun 2001 13:58:10 +0100 Received: by cclmse...

Hi all, I want to optimize a function fn using optim. This function fn calls the .C function in it with a function name and arguments given to fn, i.e. something like this pseudocode: fn <- function(par, "some params for fn") {... .C("Cfunction", ...) ...} optim(par, fn, "some params for fn") In my case this doesn't work, but is it possible in general and I

Hi all, does anybody know about R implementations for quantile regression for longitudinal data? I am just aware of a very basic version of R. Koenker's approach using fixed effects. Thanks in advance Armin

Hi, I wanted to use R's lbfgsb method for minimization from C. Unfortunately, my toy examples always crashes (segmentation fault). What's wrong with it? double eval(int n, double* par, void *ex) { double result = 0; for (int i=0; i<n; ++i) { result += par[i]*par[i]; } printf("result=%.2f\n", result); return result; } void grad(int n, double *par, double *gr,

Hi, I have to draw samples from an asymmetric-Laplace-Normal distribution: f(u|y, x, beta, phi, sigma, tau) \propto exp( - sum( ( abs(lo) + (2*tau-1)*lo )/(2*sigma) ) - 0.5/phi*u^2), where lo = (y - x*beta) and y=(y_1, ..., y_n), x=(x_1, ..., x_n) -- sorry for this huge formula -- A WinBUGS Gibbs sampler and the HI package arms sampler were used with the same initial data for all parameters. I