Re: [R] Procrustes

2015-02-19 Thread Pierrick Bruneau
Hi Tara, Providing a simple example script that reproduces your case and using it to support your question would increase your chances to obtain an answer. Best, Pierrick On Thu, Feb 19, 2015 at 5:43 PM, Tara Dirilgen tara.diril...@ucdconnect.ie wrote: I have been using R to calculate the

Re: [R] A good way to debug a c++ library embedded to an R code

2015-02-02 Thread Pierrick Bruneau
There's also a section about this in Writing R extensions: http://cran.r-project.org/doc/manuals/r-release/R-exts.html#Debugging-compiled-code Pierrick On Mon, Feb 2, 2015 at 10:42 PM, Qiang Kou q...@umail.iu.edu wrote: If you are familiar with GDB, you can just start R by R -d gdb. Dirk gave

[R] Difference between eigs() and eigen()

2014-12-29 Thread Pierrick Bruneau
Dear R users and contributors, I recently observed a difference between the outputs of the classic eigen() function, and the Arnoldi variant eigs() that extracts only the few first eigenpairs. Here is some sample code illustrating the problem: library(rARPACK) library(speccalt) set.seed(1) #

[R] Sorting method used by order()

2014-12-02 Thread Pierrick Bruneau
Hi all, In ?order, the sorting method used appears as the method argument to sort.list(), but I cannot make out which is used by default when calling order(), and how to parametrize it. Does someone have a clue there? Thanks by advance, Pierrick [[alternative HTML version deleted]]

Re: [R] Sorting method used by order()

2014-12-02 Thread Pierrick Bruneau
method is radix, AFAIK. Cheers Petr -Original Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Pierrick Bruneau Sent: Tuesday, December 02, 2014 12:33 PM To: r-h...@stat.math.ethz.ch Subject: [R] Sorting method used by order() Hi all, In ?order

Re: [R] Why does sin(pi) not return 0?

2013-09-26 Thread Pierrick Bruneau
Just to add a small note, sin(pi) is below machine precision : .Machine$double.eps [1] 2.220446e-16 (see ?.Machine for exact definition) - if myval is bounded to be positive, one safe way of testing equality to 0 would then be myval .Machine$double.eps On Thu, Sep 26, 2013 at 2:48 PM,

Re: [R] problem subsetting data.frame in R version 2.15.2 for Windows

2013-03-21 Thread Pierrick Bruneau
Hi Borja, You may issue: attach(data) which results in adding your column names to the search path of R for name resolving. Pierrick Bruneau CRP Gabriel Lippmann On Wed, Mar 20, 2013 at 11:17 PM, Borja . borjalato...@outlook.com wrote: Good day. I create a data frame like

Re: [R] problem subsetting data.frame in R version 2.15.2 for Windows

2013-03-21 Thread Pierrick Bruneau
...@gmail.com wrote: On Mar 21, 2013, at 7:39, Pierrick Bruneau pbrun...@gmail.com wrote: Hi Borja, You may issue: attach(data) No -- bad idea -- dangerous -- confusing statefulness, etc. (See explanations in the archives as to why) which results in adding your column names

Re: [R] help on extracting values from a matrix

2013-03-20 Thread Pierrick Bruneau
Hi, Well, simply : d[a24,2] should do the trick. Pierrick Bruneau CRP Gabriel Lippmann On Wed, Mar 20, 2013 at 12:53 PM, Andras Farkas motyoc...@yahoo.com wrote: 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

Re: [R] highlight overlapping region of two densities

2013-03-20 Thread Pierrick Bruneau
inelegant distribution tails... Maybe somebody will have a more elegant suggestion ? Pierrick Bruneau CRP Gabriel Lippmann On Wed, Mar 20, 2013 at 1:32 PM, andrija djurovic djandr...@gmail.comwrote: Hi all. I would like to highlight overlapping regions of two densities and I could not find a way

[R] Unexpected behaviour of apply()

2013-03-08 Thread Pierrick Bruneau
, Pierrick Bruneau http://www.bruneau44.com [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

Re: [R] Unexpected behaviour of apply()

2013-03-08 Thread Pierrick Bruneau
) ? Thanks by advance for your help, Pierrick Bruneau http://www.bruneau44.com [[alternative HTML version deleted]] __** R-help@r-project.org mailing list https://stat.ethz.ch/mailman/**listinfo/r-helphttps://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Using eigen() for extracting only few major eigenpairs

2013-02-05 Thread Pierrick Bruneau
, but with accounting for my specifics, it is not worth the effort. Thanks again, Pierrick Bruneau CRP Gabriel Lippmann On Fri, Feb 1, 2013 at 10:50 AM, Pierrick Bruneau pbrun...@gmail.comwrote: Thanks a lot, I'll test this package very soon. As it seems general purpose (ie not specifically fitted to the square

Re: [R] Using eigen() for extracting only few major eigenpairs

2013-02-01 Thread Pierrick Bruneau
Thanks a lot, I'll test this package very soon. As it seems general purpose (ie not specifically fitted to the square symmetric context), I hope the advantage over the standard routine for symmetric matrices remains significant. Pierrick Bruneau CRP Gabriel Lippmann On Thursday, January 31, 2013

[R] Using eigen() for extracting only few major eigenpairs

2013-01-31 Thread Pierrick Bruneau
library (and contribute the package myself afterwards) ? Thanks by advance for your help, Pierrick Bruneau CRP Gabriel Lippmann [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Using eigen() for extracting only few major eigenpairs

2013-01-31 Thread Pierrick Bruneau
Thanks a lot ! I'm going to try this very soon. As the package seems general purpose (ie not dedicated to real symmetric matrices), I hope it remains more efficient than the eigen function parametrized for symmetric matrices. Pierrick Bruneau CRP Gabriel Lippmann On Thu, Jan 31, 2013 at 2:51 PM

[R] Rd2pdf freeze

2012-11-04 Thread Pierrick Bruneau
help. Pierrick Bruneau Research Associate CRP Gabriel Lippmann (Luxemburg) [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] Rd2pdf freeze

2012-11-04 Thread Pierrick Bruneau
I naively thought it would be able to examine the compressed file, Thanks for your answer :) On Sun, Nov 4, 2012 at 2:55 PM, Duncan Murdoch murdoch.dun...@gmail.comwrote: On 12-11-04 8:50 AM, Pierrick Bruneau wrote: Hi everyone, From the currently available version of the package VBmix, I

[R] Filling points in a trellis object

2012-10-10 Thread Pierrick Bruneau
(theplot.xyplot.points.panel.1.1, 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

[R] impossible to fill point glyphs in a lattice plot

2012-10-10 Thread Pierrick Bruneau
missing ? Thanks by advance for your help, Pierrick Bruneau Research Fellow CRP Gabriel Lippmann [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] Problem with grid.rect

2012-09-27 Thread Pierrick Bruneau
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