[R] optim evils

lso that the routine makes fine progress toward a minimum (as desired) but in the last iteration reverses course and returns a function value much larger than the starting value. What is going on here? All help is much appreciated. Michael Meyer [[alternative HTML version deleted]] __

[R] optim evils

It would take some effort to extract selfcontained code from the mass of code wherein this optimization is embedded. Moreover I would have to obtain permission from my employer to do so. This is not efficient. However some things are evident from the trace log which I have submitted: (a) L-BFGS-

[R] optim evils

ry to destill the code so that reproducible failure of L-BFGS-B occurs and post it here. Michael Meyer [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] error handling

Greetings,   Error handling does not seem to work as expected. When a handler for general conditions is provided in first position in the handler list it will be called regardless of what condition has been signalled.   See sample code, move  "condition = function(c){...}" to first position in h

[R] How to catch errors regarding the hessian in 'optim'

een to emerge from optim # do what you have to do if such is detected } Sorry if this answer seems trivial. Michael Meyer __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.

[R] optim evils

? Thanks, Michael Meyer __ 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 commented, minimal, self-contained, reproducible code.

[R] S4 base class

not, which would be a serious drawback to the R class mechanism. Thanks, Michael Meyer __ 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 a

[R] S4 base class

w do I refer to the base class object from the derived class? Many thanks  Michael Meyer __ 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 pr

[R] S4 base class

Quote By the way, your use of the syntax D::f and B::f suggests that you're thinking from a C++ point of view.  That's very likely to lead to frustration:  the S4 object system is very different from C++.  Methods don't belong to classes, they belong to generics. There is no such thing as D::f

[R] S4 base class

alling f(baseClassObject(this),...) The question is the following: How do I refer to the base class object from the derived class? Many thanks    Michael Meyer __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] S4 base class

@Martin Morgan, Duncan Murdoch: OK Thanks. I did not understand the callNextMethod. I will investigate this in detail. This is great! Thanks again,   Michael Meyer __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] triangular matrix inverse

Greetings,   is there an algorithm which computes the inverse of a triangular matrix while being aware of its triangular form?   I know about solve() but this is probably not efficient on a triangular matrix.   Thanks, Michael Meyer [[alternative HTML version deleted

[R] C interface

Greetings, I am trying to call simple C-code from R. I am on Windows XP with RTools installed. The C-function is #include #include #include #include // prevent name mangling extern "C" { SEXP __cdecl test(SEXP s){ SEXP result; PROTECT(result = NEW_NUMERIC(1)); double* ptr=NUMERIC_PO

[R] C Interface

Thanks for all replies. I will post the question on R-develop also since eventually I would like to compile more substantial C and C++ code into shared libraries. Michael [[alternative HTML version deleted]] __ R-help@r-project.org mailing lis

[R] OOP and passing by value

Greetings, I love the R system and am sincerely grateful for the great effort the product and contributors are delivering. My question is as follows: I am trying to use S4 style classes but cannot write functions that modify an object because paramter passing is by value. For example I want to d

Re: [R] OOP and passing by value

Thanks to all replies. I was able to get it running with the R.oo package. I hope this reply makes it to the proper thread. Michael [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo

[R] OO code organization

Greetings,   I am using S4 classes and the code is organized by putting each class into a separate source file in a separate folder:   In Folder base/base.R:defines class "base" and does   setGeneric("showSelf", function(this) standardGeneric("showSelf") ) setMethod("showSelf", signature(t

[R] OO code organization

Greetings,   I found the mistake. As this may be useful to others the details are as follows:   In the source file derived_j/derived_j.R of every derived class the source of the base source("..base/base.R") was included (as you would in C or C++). This file contains the statement setGener

[R] Get the filename from the given path

Do   splits <- strsplit(filePath,"file://%22,fixed=true)[[1/]] fileName <- splits[length(splits)]) [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting g

[R] plot curve on surface

: persp(x,y,z)   Now I would like to plot the curve z=z(x,x) in red on this surface. How do I do that?   Many thanks,   Michael Meyer __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting gui

[R] numerical integration

Greetings, Sorry, the last message was sent by mistake! Here it is again: I encounter a strange problem computing some numerical integrals on [0,oo). Define $$ M_j(x)=exp(-jax) $$ where $a=0.08$. We want to compute the $L^2([0,\infty))$-inner products $$ A_{ij}:=(M_i,M_j)=\int_0^\infty M_i(x)M_j(x

[R] numerical integrals

Greetings,   I encounter a strange problem computing some numerical integrals on [0,oo). Define $$ M_j(x)=exp(-jax) $$ where $a=0.08$. We want to compute the $L^2([0,\infty))$-inner products $$ A_{ij}:=(M_i,M_j)=\int_0^\infty M_i(x)M_j(x)dx $$ Analytically we have $$ A_{ij}=1/(a(i+j)). $$ In the co

[R] Rcpp, function signature

Greetings, I am an (very) grateful user of Rcpp. As such I defined a function // [[Rcpp::export]] NumericVector leftShift(NumericVector x){ for(int i=0;ihttps://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide comme

[R] Quantiles of sums of independent discrete random variables

ideas how I could approach this problem with the tools of R? Thanks in advance for all replies, Michael Meyer [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/ma

[R] Use of geometric mean .. in good data analysis

ic mean all asymptotic results no longer apply. Michael Meyer [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide h

[R] Quantiles of sums of independent discrete random variables

/10, FUN.VALUE=0 ))) g <- function(s) dnorm(s,mean=0,sd=0.01) cat("\n\n") for(i in 1:5){       cat("Doing convolution ",i,"\n")   g <- convolve(g,f)$density } cat("\nConvolutions finished, plotting density.") s <- seq(0,100,length.out=1024) matplot(s,g(s)

[R] DescTools::Quantile

quantiles are equal to 100! Is this function working or am I not using it correctly? Michael Meyer __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide