search for: bruneau

...at[3,1] <- FALSE apply(mat, 1, indexes) ---- I would expect a 10-cell list with integer(0) in each cell - instead I get integer(0), which wrecks my further process. Is there a simple way to get the result I expect (and the only consistent one, IMHO) ? Thanks by advance for your help, Pierrick Bruneau http://www.bruneau44.com [[alternative HTML version deleted]]

...omposition when only few eigenpairs are needed : did somebody already care about this problem (through a contributed package for example) ? Or do I have to directly try to mess around with the LAPACK library (and contribute the package myself afterwards) ? Thanks by advance for your help, Pierrick Bruneau CRP Gabriel Lippmann [[alternative HTML version deleted]]

...PS [bup@localhost]: connect failed: No such host do_notify: ntype 0x0005 (COMMBAD) Communications with UPS bup@localhost lost I found that if I put the -p it work but i'm not sure of what that's implying. Any clues of how to fix that problem ? Where to look? Thanks in advance -- Mathieu Bruneau aka ROunofF === GPG keys available @ http://rounoff.darktech.org -- No virus found in this outgoing message. Checked by AVG Free Edition. Version: 7.1.371 / Virus Database: 267.14.1/206 - Release Date: 2005-12-16

...", "center")), vp=vp1) (I know I could do this in a simpler way, but let's suppose I want to go through the pushviewport / viewport def / grid.rect with vp argument sequence) I must be missing something obvious ! What am I doing wrong ? Thanks by advance for your help. Pierrick Bruneau CRP Gabriel Lippmann [[alternative HTML version deleted]]

...LSE, : VBmix_0.2.9.tar.gz:30932: unexpected '}' Warning in parse_Rd("VBmix_0.2.9.tar.gz", encoding = "unknown", fragment = FALSE, : Rd files are encoded as regular ASCII text, I cannot get what is going wrong... Any idea ? Thanks by advance for your help. Pierrick Bruneau Research Associate CRP Gabriel Lippmann (Luxemburg) [[alternative HTML version deleted]]

...ter vector. > x=c(2.00,1.20,5.00,6.56) > y= as.character(x) > y [1] "2" "1.2" "5" "6.56" What I want is : [1] "2.00" "1.20" "5.00" "6.56" Does someone know how to do this please ? Benoit Bruneau

With the following code : I would like to plot 4 points, and have the circle and diamond shapes filled with grey. What am I missing ? Thanks by advance for your help, Pierrick Bruneau Research Fellow CRP Gabriel Lippmann -- View this message in context: http://r.789695.n4.nabble.com/Filling-points-in-a-trellis-object-tp4645679.html Sent from the R help mailing list archive at Nabble.com.

Dear R contributors, Say I want to debug some C code invoked through .Call() - say "varbayes" in the VBmix package. following the instructions in "Writing R Extensions", I perform the following actions : R -d gdb run library(VBmix) CTRL+C break varbayes signal 0 mod <- varbayes(as.matrix(iris)[,1:4], 2) The breakpoint is indeed activated, seemingly at the correct position

Dear R contributors, Considering the following sample C code, that illustrates two possible uses of a Cholesky decomp for inverting a matrix, equally valid at least in theory: SEXP test() { int d = 2; int info = 0; double mat[4] = {2.5, 0.4, 0.4, 1.7}; double id[4] = {1.0, 0.0, 0.0, 1.0}; double lmat[3]; F77_CALL(dpotrf)("L", &d, mat, &d, &info); lmat[0] = mat[0]; lmat[1]

Dear R contributors, I'm running into trouble when trying to allocate some large (but in theory viable) vector in the context of C code bound to R through .Call(). Here is some sample code summarizing the problem: SEXP test() { int size = 10000000; double largevec[size]; memset(largevec, 0, size*sizeof(double)); return(R_NilValue); } If size if small enough (up to 10^6), everything is

On Thu, 2014-12-11 at 14:00 +0100, Pierrick Bruneau wrote: > Dear R contributors, > > Say I want to debug some C code invoked through .Call() - say > "varbayes" in the VBmix package. following the instructions in > "Writing R Extensions", I perform the following actions : > > R -d gdb > run > library(...

Hi Pierrick, You're storing largevec on the stack, which is probably causing a stack overflow. Allocate largvec on the heap with malloc or one of the R memory allocation routines instead and it should work fine. Karl On Thu, Dec 18, 2014 at 12:00 AM, Pierrick Bruneau <pbruneau at gmail.com> wrote: > > Dear R contributors, > > I'm running into trouble when trying to allocate some large (but in > theory viable) vector in the context of C code bound to R through > .Call(). Here is some sample code summarizing the problem: > > SEXP...

This isn't the help list for LAPACK, but as far as I can tell, dppsv expects a symmetric matrix input compacted as triangular, not a Choleski decomposed one. So try assigning lmat before the call to dpotrf. -pd > On 20 Dec 2014, at 22:06 , Pierrick Bruneau <pbruneau at gmail.com> wrote: > > Dear R contributors, > > Considering the following sample C code, that illustrates two possible > uses of a Cholesky decomp for inverting a matrix, equally valid at > least in theory: > > SEXP test() { > > int d = 2; > in...

Dear R-Users, I'm trying to setup a personal repository for a few packages I'm working on. I am on R-Forge but I still need to have various versions of my package that R-Forge does not build (for R 2.8.1 for example). So I followed the instructions in this document: Hhttp:// cran.r-project.org/doc/manuals/R-admin.html#Setting-up-a-package-repository and used this function as

...t;, gp=gpar(lwd=c(2,2,2,2), cex=c(2,3,3,2), col=1, fill=c("transparent", "grey", "grey", "transparent"))) I would like to plot 4 points, and have the circle and diamond shapes filledwith grey. What am I missing ? Thanks by advance for your help, Pierrick Bruneau Research Fellow CRP Gabriel Lippmann [[alternative HTML version deleted]]

Dear All,   any thoughts on how I can do the following:   let us say we have:   a <-c(2,4,16,28,48) b <-c(10,4,2,0.4,0.03) d <-cbind(a,b)   what I would like to do is to extract values of column b in the matrix based on the corresponding values of column a. For example: I would like to extract all b values that have a corresponding a value that is less than 24 into a numeric vector, so

Good day. I create a data frame like this: > data <- data.frame(a=1:10,b=11:20,c=21:30) I can subset this data.frame by saying: > data[data$a>7,] and I get this result a b c 8 8 18 28 9 9 19 29 10 10 20 30 I understand I should get the same result by saying

Hi all. I would like to highlight overlapping regions of two densities and I could not find a way to do it. Here is the sample code: myd <- c(2,4,5, 4,3,2,2,3,3,3,2,3,3,4,2,4,3,3,3,2,2.5, 2, 3,3, 2.3, 3, 3, 2, 3) myd1 <- myd-2 plot(range(density(myd)$x, density(myd1)$x), range(density(myd)$y, density(myd1)$y), type = "n") lines(density(myd), col=1, lwd=4)