[Rd] extending package with function calling an Objective Caml program

2007-06-21 Thread Tamara Steijger
Hallo, we are trying to extend the R package multcompView in agreement with the author Hans-Peter Piepho. The function multcompLetters implements so far a heuristic. We would like to add a function that implements an exact algorithm and returns a provable optimum result. This algorithm

Re: [Rd] extending package with function calling an Objective Caml program

2007-06-21 Thread Prof Brian Ripley
On Wed, 20 Jun 2007, Tamara Steijger wrote: Hallo, we are trying to extend the R package multcompView in agreement with the author Hans-Peter Piepho. The function multcompLetters implements so far a heuristic. We would like to add a function that implements an exact algorithm and returns a

[Rd] Basic Question on error(_(...

2007-06-21 Thread statmobile
I'm sorry to bother this list with such trivial questions, but I'm trying to take Prof. Ripley's advice in porting some Lapack wrappers into my own code, because they are not public. I'm specifically choosing programs from R-2.5.0/src/modules/lapack/Lapack.c. I can't seem to understand the

Re: [Rd] Basic Question on error(_(...

2007-06-21 Thread Peter Dalgaard
[EMAIL PROTECTED] wrote: I'm sorry to bother this list with such trivial questions, but I'm trying to take Prof. Ripley's advice in porting some Lapack wrappers into my own code, because they are not public. I'm specifically choosing programs from R-2.5.0/src/modules/lapack/Lapack.c. I

Re: [Rd] Basic Question on error(_(...

2007-06-21 Thread Prof Brian Ripley
On Thu, 21 Jun 2007, Peter Dalgaard wrote: [EMAIL PROTECTED] wrote: I'm sorry to bother this list with such trivial questions, but I'm trying to take Prof. Ripley's advice in porting some Lapack wrappers into my own code, because they are not public. I'm specifically choosing programs from

Re: [Rd] segfault during cbind

2007-06-21 Thread Martin Morgan
To amend, the traceback is only Traceback: 1: cbind(probe, gb) aborting ... Oops. Martin Martin Morgan [EMAIL PROTECTED] writes: The following code results in a seg fault. sessionInfo() R version 2.6.0 Under development (unstable) (2007-06-21 r42013) x86_64-unknown-linux-gnu

Re: [Rd] segfault during cbind

2007-06-21 Thread Prof Brian Ripley
I think it needs to be case LISTSXP: PROTECT(u = coerceVector(u, mode)); k = LENGTH(u); idx = (!isMatrix(u)) ? rows : k; for (i = 0; i idx; i++) SET_VECTOR_ELT(result, n++,

Re: [Rd] segfault during cbind

2007-06-21 Thread Martin Morgan
Yes, that seems to do the trick, for both seg fault and valgrind! Thanks very much. Martin Prof Brian Ripley [EMAIL PROTECTED] writes: I think it needs to be case LISTSXP: PROTECT(u = coerceVector(u, mode)); k = LENGTH(u);

Re: [Rd] segfault during cbind

2007-06-21 Thread Prof Brian Ripley
On Thu, 21 Jun 2007, Martin Morgan wrote: Yes, that seems to do the trick, for both seg fault and valgrind! Thanks very much. Thanks: this and a similar one elsewhere seem long-standing bugs in bind.c Will commit shortly. Martin Prof Brian Ripley [EMAIL PROTECTED] writes: I think it