Re: [R] Maximum likelihood with analytical Hessian and

2014-12-18 Thread Xavier Robin
carefully. > > Generally the work in computing the Hessian ( # obs * (# parameters)^2 > in size) is not worth the effort, but there are problems for which it > does make a lot of sense. > > JN > > On 14-12-18 06:00 AM, r-help-requ...@r-project.org wrote: >> Message: 1

[R] Maximum likelihood with analytical Hessian and

2014-12-17 Thread Xavier Robin
Dear list, I have an optimization problem that I would like to solve by Maximum Likelihood. I have analytical functions for the first and second derivatives of my parameters. In addition, some parameters are constrained between 0 and 1, while some others can vary freely between -Inf and +Inf. I a

[R] attach 'name' argument ignored with a file?

2011-11-22 Thread Xavier Robin
Dear useRs & experRts, I have the feeling that the 'name' argument to the attach function is ignored when 'what' is a file name. Here is an example: > save(letters, file="letters.RData") > letters.env <- attach("letters.RData", name="letters") > search() > letters.env The name on the search path

Re: [R] Equivalent of win.print for Linux

2011-10-28 Thread Xavier Robin
Thank you, it was exactly what I was looking for! Regards, Xavier Le 28. 10. 11 15:45, Prof Brian Ripley a écrit : > See the help for postscript ... especially the 'Printing' section. > > > On Fri, 28 Oct 2011, Xavier Robin wrote: > >> Hi! >> >>

Re: [R] Equivalent of win.print for Linux

2011-10-28 Thread Xavier Robin
inter=blah) Now the question is, is there a printer device equivalent to win.print for Linux? Thanks, Xavier On Fri, Oct 28, 2011 at 15:23, Raphael Saldanha wrote: > I have never used, but take a look on ?dev.print > > On Fri, Oct 28, 2011 at 11:16 AM, Xavier Robin wrote: > >> H

[R] Equivalent of win.print for Linux

2011-10-28 Thread Xavier Robin
Hi! In Windows the "win.print" function allows plotting directly to a printer (or copying an open device to the printer). This is very convenient to quickly print a plot once it looks good. Is there an equivalent function under Linux? For example through CUPS, IPP, LPD or other ? Obviously with

[R] [R-pkgs] pROC 1.4.3: compare two ROC curves in R

2011-04-01 Thread Xavier Robin
find more information in our paper [4] and on pROC website: http://www.expasy.org/tools/pROC/ Hope you'll find it useful! Xavier Robin -- References: [1] DeLong ER, DeLong DM, Clarke-Pearson DL (1988) Comparing the areas under two or more correlated receiver operating characteristic curves:

Re: [R] Integration of two lines

2011-01-26 Thread Xavier Robin
Hans W Borchers wrote : > First define a function from those points: > > fx <- approxfun(x, f_x) > fy <- approxfun(y, f_y) > f <- function(x) abs(fx(x)-fy(x)) > > and now you can apply integrate() or trapz(): > > xx <- sort(c(x, y)) > yy <- f(xx) > trapz(xx, yy) > > tra

Re: [R] Integration of two lines

2011-01-25 Thread Xavier Robin
Le 25.01.2011 15:23, Rmh a écrit : > g <- function(x) abs(f1(x)-f2(x)) > > now you have one function and you can integrate it. Thank you Rich. Unfortunately I have no f1 and f2 functions, only a set of observed points on two lines - and no idea about the underlying distribution to create a funct

[R] Integration of two lines

2011-01-25 Thread Xavier Robin
Hello, I need to integrate the absolute difference between two lines measured on different points. # For example : x <- seq(0, 1, 1/100) f_x <- runif(101) + x y <- seq(0, 1, 1/23) f_y <- runif(24) + (1 - y) plot(x, f_x, type="l") lines(y, f_y) Then I would like to compute Integral( | f_x - f_y

[R] Propagate vector attributes to data frame

2008-10-29 Thread Xavier Robin
o extract them. attributes(myDataFrame$myColum) is NULL and there is nothing in attributes(myDataFrame) that looks like a possible myAttr. Is there a way so that attributes of numeric vectors could be kept in the data.frame? Thanks in advance, Xavier Robin -- Xavier Robin Biomedical Proteomics

Re: [R] xy.coords in text

2008-10-24 Thread Xavier Robin
options include using the grconvertX and grconvertY functions. > > Hope this helps, > > -- > Gregory (Greg) L. Snow Ph.D. > Statistical Data Center > Intermountain Healthcare > [EMAIL PROTECTED] > 801.408.8111 > > >> -Original Message- >> Fr

[R] xy.coords in text

2008-10-23 Thread Xavier Robin
es that > 'y' may be missing since 'xy.coords(x,y)' is used for construction of the > coordinates. Can someone explain me this difference? And optionally how to plot text in the "topleft" part of the plot without using legend? Thanks, Xavier -- Xavier Robin Bi

Re: [R] Rpart and case weights: working with functions

2008-06-03 Thread Xavier Robin
S Ellison a écrit : As a work-round, try something lik test.function <- function (formula.str, data) { weights <- rep(.1, 8) rpart(as.formula(formula.str), data, weights) } which works if formula.str is a string instead of a formula object. Ok, I think I got the point. The only problem

[R] Rpart and case weights: working with functions

2008-06-03 Thread Xavier Robin
do to get it work (i.e. using the copy of weights defined inside the function)? Thanks in advance for your help. Xavier Robin -- Xavier Robin Biomedical Proteomics Research Group (BPRG) Department of Structural Biology and Bioinformatics (DBSB) Geneva University Medical Center (CMU) 1, rue Mich