search for: axionator

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 Kusnierczyk)...

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

Dear Folks, I've got a generally well-behaved Solaris 2.5.1 machine running a generally well-behaved version of 1.9.17p4 running as a daemon. Everything works beautifully including browsing with the help of a central WINS database on an NT box and printing using a group of HP printers with JetDirect cards. BUT...... After printing a document, Samba appears to hang up the connection,

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 the user tries to access > something on the network drive

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]

\\Pettit Connect As: Password: On the server side, I didn't get far because the domain name come up NULL so it failed on that test. Below are two logs log98.txt and log2000.txt that show full debug info. I might just have a stupid configuration or installation mistakes. I am sure that my basic configuration must work since it is a very typical setup. Any help will be highly appreciated

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