Re: [R] non-linear integer optimization?

2010-09-23 Thread Hans W Borchers
darckeen darckeen at hotmail.com writes: This is an example of the type of problem, and how i'm currently using optim() to solve it. I would classify this as an integer programming (IP) problem, it is not even non-linear, as there are no continuous variables. Your approach with optim() will

Re: [R] bivariate vector numerical integration with infinite range

2010-09-21 Thread Hans W Borchers
baptiste auguie baptiste.auguie at googlemail.com writes: Dear list, I'm seeking some advice regarding a particular numerical integration I wish to perform. The integrand f takes two real arguments x and y and returns a vector of constant length N. The range of integration is [0,

Re: [R] bivariate vector numerical integration with infinite range

2010-09-21 Thread Hans W Borchers
baptiste auguie baptiste.auguie at googlemail.com writes: Thanks, adaptIntegrate() seems perfectly suited, I'll just need to figure a transformation rule for the infinite limit. The suggestion of x-1/x does not seem to work here because it also transforms 0 into -infinity. I think exp(pi*

Re: [R] Prime Factorization

2010-09-21 Thread Hans W Borchers
David Winsemius dwinsemius at comcast.net writes: On Sep 21, 2010, at 11:33 AM, Cliff Clive wrote: Hi everyone, I have a very quick question: Is there a ready-made function in R or any R packages to find the prime factorization of an integer? Yes. At least two. The obvious

Re: [R] bivariate vector numerical integration with infinite range

2010-09-21 Thread Hans W Borchers
baptiste auguie baptiste.auguie at googlemail.com writes: Thanks. I am having trouble getting adaptIntegrate to work with a multivalued integrand though, and cannot find a working example. Anyone had better luck with it? The function to be integrated needs a vector as input: f -

Re: [R] Prime Factorization

2010-09-21 Thread Hans W Borchers
David Winsemius dwinsemius at comcast.net writes: A further citation that answers the question I raised (and inaccurately predicted no value) regarding prime.sieve : http://finzi.psych.upenn.edu/R/Rhelp02/archive/49773.html This was found with Barons search facility set for rhelp

Re: [R] optimized value worse than starting Value

2010-09-09 Thread Hans W Borchers
Barry Rowlingson b.rowlingson at lancaster.ac.uk writes: On Wed, Sep 8, 2010 at 1:35 PM, Michael Bernsteiner dethlef1 at hotmail.com wrote: Dear all, I'm optimizing a relatively simple function. Using optimize the optimized parameter value is worse than the starting. why? I would

Re: [R] calculating area between plot lines

2010-09-07 Thread Hans W Borchers
A. Marcia BARBOSA anamarciabarbosa at gmail.com writes: Hi everyone. I have these data: probClass-seq(0,0.9,0.1) prob1-c(0.0070,0.0911,0.1973,0.2949,0.3936,0.5030,0.5985,0.6869,0.7820,0.8822) prob2-c(0.0066,0.0791,0.2358,0.3478,0.3714,0.3860,0.6667,0.6400,0.7000,1.) # which I'm

Re: [R] Regarding naive baysian classifier in R

2010-08-31 Thread Hans W Borchers
Vishnampettai akron_aadhithya at yahoo.com writes: Hi, I have a small doubt regarding naive Bayes. I am able to classify the data's properly but i am just stuck up with how to get the probability values for naive bayes. In the case of SVM we have attr function that helps in displaying the

Re: [R] LSAP and the Hungarian algorithm [was: R-help]

2010-08-27 Thread Hans W Borchers
venkatesh bandaru venkatesh.bandaru at gmail.com writes: If you had searched for it, you would have easily found the 'clue' package: In package clue solve_LSAP() enables the user to solve the linear sum assignment problem (LSAP) using an efficient C implementation of the Hungarian

Re: [R] LSAP and the Hungarian algorithm [was: R-help]

2010-08-27 Thread Hans W Borchers
Ravi Varadhan rvaradhan at jhmi.edu writes: However, The clue package has the solve_LSAP() function (as pointed out by Hans Werner) that solves the linear sum assignment problem, but it uses the Hungarian algorithm that was requested by you. The lp.assign() function in lpSolve package (as

Re: [R] Puzzle

2010-08-26 Thread Hans W Borchers
Ben Holt BHolt at bio.ku.dk writes: I have data similar to this: Location Surveyor Result A1 83 A2 76 A3 45 B1 71 B4 67 C2 23 C5 12 D3 34 E4

Re: [R] Problem to compute a function with very large numbers

2010-08-20 Thread Hans W Borchers
Nan Zhao nzhao at student.ethz.ch writes: Thank you Dennis for your explanations! The results you found are the same as mine. with first an infinity result, followed by NaN. It seems that, when the number becomes too small, R must round it up to 0. Hence I was wondering if there might a

Re: [R] Games

2010-08-13 Thread Hans W Borchers
Silvano silvano at uel.br writes: Hi, I want to build the table of a football league with 11 teams. All play together. So will 55 games. Since there are an odd number of teams in each round a team will not play. The easy solution is moving around a table with one team pausing.

Re: [R] Why use numFold in evaluate_Weka_classifier of RWeka

2010-08-10 Thread Hans W Borchers
s0300851 s0300851 at tp.edu.tw writes: Hi everyone, I have a question about using RWeka package, we know that instruction make_Weka_classifier that can help us to build a model,and evaluate_Weka_classifier instruction can help us to evaluate the performance of the model using on new

Re: [R] Optimization problem with nonlinear constraint

2010-07-28 Thread Hans W Borchers
Uli Kleinwechter u.kleinwechter at uni-hohenheim.de writes: Dear Ravi, As I've already written to you, the problem indeed is to find a solution to the transcendental equation y = x * T^(x-1), given y and T and the optimization problem below only a workaround. I don't think

Re: [R] Hurst Exponent Estimation

2010-07-22 Thread Hans W Borchers
Lorenzo Isella lorenzo.isella at gmail.com writes: Dear All, I am a novice when it comes to time-series analysis and at the moment I am actually interested in calculating the Hurst exponent of a time series. Some time ago I tested some of the classical chaotic time series (such as the

Re: [R] Hurst Exponent Estimation

2010-07-22 Thread Hans W Borchers
, Hans W Borchers wrote: Lorenzo Isellalorenzo.isellaat gmail.com writes: Dear All, I am a novice when it comes to time-series analysis and at the moment I am actually interested in calculating the Hurst exponent of a time series. Some time ago I tested some of the classical chaotic

Re: [R] Double Integration

2010-07-03 Thread Hans W Borchers
Bogaso Christofer bogaso.christofer at gmail.com writes: Hi Ravi, your suggestion helped me as well a lot. If I look into that function, I see this function is calling another function : .Call(doCubature, as.integer(fDim), body(f.check), as.double(lowerLimit),

Re: [R] Double Integration

2010-07-01 Thread Hans W Borchers
Sarah Sanchez sarah_sanchez09 at yahoo.com writes: Dear R helpers I am working on the Bi-variate Normal distribution probabilities. I need to double integrate the following function (actually simplified form of bivariate normal distribution) f(x, y) = exp [ - 0.549451 * (x^2 + y^2 - 0.6

Re: [R] double integral

2010-06-18 Thread Hans W Borchers
suman dhara suman.dhara89 at gmail.com writes: Sir, I want to calculate double integral in R. Is there any function to do this? If your domain of integration is a hypercube, try packages 'cubature' or 'R2cuba'. Otherwise, you have to uncover more information about your specific problem

Re: [R] logarithmic integrals in R?

2010-05-30 Thread Hans W. Borchers
. I am certain I 've seen them in AandS' handbook (where else?), but sure cannot remember in which chapter or page. Which logarithmic integrals do you really need, and on what range? Regards, Hans Werner On Sat, May 29, 2010 at 01:15:29PM +, Hans W. Borchers wrote: Oliver Kullmann

Re: [R] logarithmic integrals in R?

2010-05-29 Thread Hans W. Borchers
Oliver Kullmann O.Kullmann at swansea.ac.uk writes: Hello, I couldn't find information on whether the logarithmic integrals Li_m(x) = integral_0^x log(t)^(-m) dt for x = 0 are available in R? I saw your request only this weekend. The first logarithmic integral can be computed using

Re: [R] best polynomial approximation

2010-05-18 Thread Hans W Borchers
I guess you may be looking for the Remez algorithm. AFAIK there is no implementation in one of the R packages. You can find FORTRAN code in the Collected Algorithms of the ACM (no. 604) which probably could be called from R. There appears to exist a discrete, equi-distant(?) version as function

Re: [R] Find a rectangle of maximal area

2010-03-22 Thread Hans W Borchers
Hans W Borchers hwborchers at googlemail.com writes: For an application in image processing -- using R for statistical purposes -- I need to solve the following task: Given n (e.g. n = 100 or 200) points in the unit square, more or less randomly distributed. Find a rectangle of maximal

Re: [R] Find a rectangle of maximal area

2010-03-22 Thread Hans W Borchers
Barry Rowlingson b.rowlingson at lancaster.ac.uk writes: On Mon, Mar 22, 2010 at 4:28 PM, Hans W Borchers hwborchers at googlemail.com wrote: Still I believe that a clever approach might be possible avoiding the need to call a commercial solver. I am getting this hope from one of Jon

[R] Find a rectangle of maximal area

2010-03-21 Thread Hans W Borchers
For an application in image processing -- using R for statistical purposes -- I need to solve the following task: Given n (e.g. n = 100 or 200) points in the unit square, more or less randomly distributed. Find a rectangle of maximal area within the square that does not contain any of these

Re: [R] Sorting

2010-02-06 Thread Hans W Borchers
David Neu david at davidneu.com writes: Hi, I have a list of vectors (of varying lengths). I'd like to sort this list by applying a function to each pair of vectors in the list and returning information to sorting routine that let's it know which one is larger. To solve problems like

Re: [R] Sorting

2010-02-06 Thread Hans W Borchers
David Neu david at davidneu.com writes: David Neu david at davidneu.com writes: Hi, I have a list of vectors (of varying lengths).  I'd like to sort this list by applying a function to each pair of vectors in the list and returning information to sorting routine that let's it know which

Re: [R] Solving an optimization problem: selecting an quot; optimalquot; subset

2010-02-02 Thread Hans W Borchers
Erwin Kalvelagen-2 wrote: Hans W Borchers hwborchers at googlemail.com writes: # Prepare inputs for MILP solver obj - c(rep(0, n), 0, 1, 1, 0) typ - c(rep(B, n), B, C, C, B) mat - matrix(c(s, -z, -1, 1, 0,# a = a_p + a_m rep(0, n), 1, 0, 0

Re: [R] Solving an optimization problem: selecting an quot;optimalquot; subset

2010-02-01 Thread Hans W Borchers
Dimitri Shvorob dimitri.shvorob at gmail.com writes: Given vector of numbers x, I wish to select an n-subset with sum closest fixed value s. Can anyone advise me how to approach this, in R? I have considered Rcplex package, which handles integer/binary linear/quadratic optimization

Re: [R] Solving an optimization problem: selecting an optimal subset

2010-01-31 Thread Hans W Borchers
Dimitri Shvorob wrote: Same request to Hans: I am afraid I need a little more spoon-feeding following I sent a GAMS script modeling this problem to the NEOS solvers Thanks a lot! If you have access to CPLEX (I mean the commercial program, not Rcplex which is just an interface to

Re: [R] Solving an optimization problem: selecting an quot;optimalquot; subset

2010-01-30 Thread Hans W Borchers
Dimitri Shvorob dimitri.shvorob at gmail.com writes: Is it a subset of a vector containing 100 elements, or 1ths? I need to pick 2-40 elements out of a 50-200-element-long vector. A random number of elements that should be chosen, or the best 10 values which sums up to a defined

Re: [R] Solving an optimization problem: selecting an quot;optimalquot; subset

2010-01-30 Thread Hans W Borchers
Dimitri Shvorob dimitri.shvorob at gmail.com writes: This is a subset sum problem and has been discussed here in December Thanks a lot! Will investigate. Can you settle for an approximate solution? Absolutely. You can use the script from the thread subset sum problem to find

Re: [R] problem with the precision of numbers

2010-01-25 Thread Hans W Borchers
kayj kjaja27 at yahoo.com writes: Hi All, thank you all for your help. I have tried Bill's script and it works! so I am trying to think what was the problem and it looks like it i sthe precision. so you defined a function of the precision and evaluates at precision=500. Bill, I was

Re: [R] problem with the precision of numbers

2010-01-19 Thread Hans W Borchers
Ted.Harding at manchester.ac.uk writes: [...] I suspect this is an invented computation -- the 3456 strikes me as unlikely (it reminds me of my habitual illustrative use of set.seed(54321)). There is a definite problem with the development given by kayj. When k=2000 and i=k, the

Re: [R] optimization problem

2010-01-17 Thread Hans W. Borchers
Ravi Varadhan rvaradhan at jhmi.edu writes: Interesting! Now, if I change the cost matrix, D, in the LSAP formulation slightly such that it is quadratic, it finds the best solution to your example: Dear Ravi, I thought your solution is ingenious, but after the discussion with Erwin

Re: [R] optimization problem

2010-01-17 Thread Hans W. Borchers
Ravi Varadhan rvaradhan at jhmi.edu writes: Dear Hans, I agree with your comments. My intuition was that the quadratic form would be better behaved than the radical form (less nonlinear!?). So, I was hoping to see a change in behavior when the cost function was altered from a radical

Re: [R] how to count the total number of (INCLUDING overl apping) occurrences of a substring within a string ?

2009-12-20 Thread Hans W Borchers
Gabor Grothendieck ggrothendieck at gmail.com writes: Use a zero lookaround expression. It will not consume its match. See ?regexp gregexpr(a(?=a), aaa, perl = TRUE) [[1]] [1] 1 2 attr(,match.length) [1] 1 1 I wonder how you would count the number of occurrences of, for example,

Re: [R] how to count the total number of (INCLUDING overl apping) occurrences of a substring within a string ?

2009-12-20 Thread Hans W Borchers
Hans Werner Correction: I meant the '\G' metacharacter in Perl, not a modifier. On Sun, Dec 20, 2009 at 7:22 AM, Hans W Borchers hwborchers at googlemail.com wrote: Gabor Grothendieck ggrothendieck at gmail.com writes: [Sorry; Gmane forces me to delete more quoted text

Re: [R] Subset sum problem.

2009-12-09 Thread Hans W Borchers
Geert Janssens janssens-geert at telenet.be writes: Hi, I'm quite new to the R-project. I was suggested to look into it because I am trying to solve the Subset sum problem, which basically is: Given a set of integers and an integer s, does any non-empty subset sum to s? (See

Re: [R] Subset sum problem.

2009-12-09 Thread Hans W Borchers
Geert Janssens janssens-geert at telenet.be writes: On Wednesday 9 December 2009, Hans W Borchers wrote: Geert Janssens janssens-geert at telenet.be writes: [ ... ] Has anybody tackled this issue before in R ? If so, I would be very grateful if you could share your solution with me

Re: [R] Modula Generators

2009-12-08 Thread Hans W Borchers
Sam K upperhalfplane at yahoo.co.uk writes: Hi all, Is there function on R for calculating Modula generators? For example for primes above 100, e.g 157, i want to know which number generates the group under multiplication mod 157. i.e i want to find an element whose order is 156. The

Re: [R] Solve linear program without objective function

2009-12-04 Thread Hans W Borchers
Andreas Wittmann andreas_wittmann at gmx.de writes: Dear R-users, i try to solve to following linear programm in R 0 * x_1 + 2/3 * x_2 + 1/3 * x_3 + 1/3 * x_4 = 0.3 x_1 + x_2 + x_3 + x_4 = 1 x_1, x_2, x_3, x_4 0, x_1, x_2, x_3, x_4 1 as you can see i have no objective function

Re: [R] Modular inverses

2009-12-01 Thread Hans W Borchers
SJ Robson-Davis sr6827 at bristol.ac.uk writes: I want to find the inverse of an integer k mod p (prime.) Is there a function that can do this for me? I know i could simply write (k^(p-2)) %% p, but i need to do this for large primes (above 100) and this gives the warning message:

Re: [R] optim or nlminb for minimization, which to believe?

2009-11-29 Thread Hans W Borchers
Your function named 'gradient' is not the correct gradient. Take as an example the following point x0, very near to the true minimum, x0 - c(-0.2517964, 0.4898680, -0.2517962, 0.4898681, 0.7500995) then you get gradient(x0) [1] -0.0372110470 0.0001816991 -0.0372102284

Re: [R] Implementation of the Shuffled Complex Evolution (SCE-UA) Algorithm

2009-11-11 Thread Hans W Borchers
Simon Seibert simon.seibert at mytum.de writes: Good evening list, I'm looking for an R implementation of the Shuffled Complex Evolution” (SCE-UA) algorithm after Duan et al. (1993). Does anybody know if there is an extension/ package existing that contains it? Thanks very much for your

Re: [R] help using R's linprog for LP

2009-10-24 Thread Hans W Borchers
Medha Atre medha.atre at gmail.com writes: Hi, I found the reason. By default it puts a condition for x = 0. Is there a way to get rid of this condition? The constraints x = 0 are used in most linear programming realizations. Some bounds from below are needed. The trick to circumvent the

Re: [R] Ubuntu, Revolutions, R

2009-10-05 Thread Hans W. Borchers
I updated to Ubuntu 9.10 Beta yesterday, and yes I do see the same message and I am a bit irritated. I don't want to read these 'marketing' lines any time I start up R. I simply deleted the lines from /etc/R/Rprofile.site for now, but I am still wondering who put that in. Is there any deeper

Re: [R] Semi continous variable- define bounds using lpsolve

2009-09-22 Thread Hans W. Borchers
I played around a bit with the original 'lp-solve' program --- i.e., not the R package but the program to be downloaded from Sourceforge ---, at least version 5.5.0.15 through its IDE. I was not even able to reproduce the example on semi-continuous variables in the reference documentation at

Re: [R] Semi continous variable- define bounds using lpsolve

2009-09-22 Thread Hans W. Borchers
But of course, it is always possible to emulate a semi-continuous variable by introducing a binary variable and use some big-M trick. That is, with a new binary variable b we add the following two conditions: x3 - 3.6 * b = 0 and x3 - 10 * b = 0 # Big-M trick, here M = 10 (If

Re: [R] Quadratic Constraints

2009-09-20 Thread Hans W. Borchers
The package lpSolve (that I have recommended before) supports so-called 'semi-continuous variables', that is Semi-continuous variables are variables that must take a value between their their minimum and maximum or zero. So these variables are treated the same as regular variables,

Re: [R] computing the radius of an arc

2009-07-28 Thread Hans W Borchers
Nair, Murlidharan T mnair at iusb.edu writes: Alex Brenning, the developer of the RSAGA package told me that and I quote the RSAGA package (which uses functions from the free geographical information system [GIS] SAGA GIS) has a curvature function that is designed to calculate the

Re: [R] re ading jpeg images?

2009-07-28 Thread Hans W. Borchers
I found the 'biOps' package for Image and data analysis quite helpful. (I did some astronomical investigations with it --- counting galaxies in a Hubble picture---and I do recommend this package.) Under Windows you have to unpack the 'libjpeg' and 'libtiff' libraries beforehand somewhere in

Re: [R] computing the radius of an arc

2009-07-24 Thread Hans W Borchers
Nair, Murlidharan T mnair at iusb.edu writes: Hi!! I am interesting in computing the radius of an arc that best approximates a curve. Is there an R function that I can use to draw an arc? Nothing useful came up when I searched help.search. Does anyone have any suggestion to do this?

[R] R interface to PSwarm

2009-07-23 Thread Hans W Borchers
Dear list: Being a bit unsatisfied with global optimization approaches in R such as SANN (in 'optim') or DEoptim, I looked for alternatives on the Web, such as PSwarm http://www.norg.uminho.pt/aivaz/pswarm/ or PIKAIA. There is an R interface to PSwarm (version 1.4) on its home page which I was

Re: [R] Automatic differentiation in R

2009-07-23 Thread Hans W. Borchers
Having given a lecture on Numerical Derivatives just a short time ago, I would like to mention the following: Many functions, especially in engineering, are not available as formulas built simply from arithmetical operators and elementary functions. They are provided as intricate procedures, as

[R] r-project.org address blacklisted by anti-spam software

2009-07-07 Thread Hans W Borchers
Dear List: An e-mail mentioning the r-project.org address and sent to a friend at a German university was considered spam by the local spam filter. Its reasoning: the URL r-project.org is blacklisted at uribl.swinog.ch resp. at antispam.imp.ch. I checked the list

Re: [R] R- NLP on R but ....

2009-07-06 Thread Hans W Borchers
Rahul Varshney itsrahulvarshney at gmail.com writes: I'll appreciate the help on the following problem: I solved many Nonlinear programming problems with nonlinear constraintsRdonlp is working well but i am unable to get INTEGER data with nonlinear constraints in Rdonlp. Is it

Re: [R] Equivalent to Matlab's Ans

2009-06-30 Thread Hans W. Borchers
There is a discussion on this topic under the heading A shorter version of .Last.value on July 7, 2008, see for example http://www.nabble.com/A-shorter-version-of-%22.Last.value%22--to18322831.html#a18322831 --Hans Werner Stephane-18 wrote: Hi everyone, I was just wondering if there is an

Re: [R] Optimization and Linear Programming in R

2009-06-26 Thread Hans W Borchers
Chris.Wilcox at csiro.au writes: Dear List, [...] We are looking for a solver that can deal with this nonlinear integer programming problem. We looked at a number of packages on the CRAN Task View: Optimization and Mathematical Programming, however, we have not been able to locate

Re: [R] newtons method

2009-05-14 Thread Hans W. Borchers
Dear Ravi: Thanks for pointing out the homotopy methods. Coming from Mathematics I was always considering SINGULAR for such a task which is also providing results when the solution set is not isolated points, but an algebraic variety. For single points, homotopy methods appear to be an

Re: [R] MDL - Fayyad Irani Discretization

2009-04-27 Thread Hans W. Borchers
I have asked for this some years ago and did not got a positive hint. To my knowledge this has not changed since then. Brian Ripley proposed to use the 'rpart' algorithm for discretization. I think I applied the RELIEF-F method at that time and also wrote a simplified Fayyad-Irani function on

Re: [R] Constrined dependent optimization.

2009-04-09 Thread Hans W. Borchers
and will not speed up the optimization procedure considerably. Regards, Hans Werner rkevinburton wrote: It has been suggested that I do a cluster analysis. Wouldn't this bet mepart way there? Thank you for your response I am looking up the book now. Kevin Hans W. Borchers hwborch

Re: [R] Constrined dependent optimization.

2009-04-04 Thread Hans W. Borchers
Just in case you are still interested in theoretical aspects: In combinatorial optimization, the problem you describe is known as the Quadratic (Sum) Assignment Problem (QAP or QSAP) and is well known to arise in facility and warehouse layouts. The task itself is considered hard, comparable to

Re: [R] Constrined dependent optimization.

2009-03-30 Thread Hans W. Borchers
Image you want to minimize the following linear function f - function(x) sum( c(1:50, 50:1) * x / (50*51) ) on the set of all permutations of the numbers 1,..., 100. I wonder how will you do that with lpSolve? I would simply order the coefficients and then sort the numbers 1,...,100

Re: [R] asking advice for Integer Programming packages

2009-03-29 Thread Hans W. Borchers
Looking into the 'Optimization' task view you will also identify the Rsymphony package. A short overview of these packages and a benchmark based on a test suite can be found in http://www.rmetrics.org/Meielisalp2008/Presentations/Theussl2.pdf. Different problems may have quite different

Re: [R] Rdonlp2 -Query

2009-03-07 Thread Hans W. Borchers
Leo Guelman leo.guelman at gmail.com writes: Hi, Did anyone used this package? Could you please share your thought on it? What do you, exactly, mean with share your thought on it? It has its pros and cons, as always. Sure Rdonlp2 has been used, and it has been requested and discussed

Re: [R] Self-Organizing Map analysis

2009-03-04 Thread Hans W. Borchers
glaporta glaporta at freeweb.org writes: Dear list, I read the SOM package manual but I don't understand how to perform (for example) 1) the SOM analysis on Iris data 2) with a visualization similar to that of figure 7 in http://www.cis.hut.fi/projects/somtoolbox/package/papers/techrep.pdf

Re: [R] R - MATLAB apply like function

2009-03-04 Thread Hans W. Borchers
ARDIA David david.ardia at unifr.ch writes: Dear all, I very often use the R function apply, for speedup purposes. I am now also using MATLAB, and would like to use the same kind of function. I have already asked MATLAB people, and the answer is : vectorize... but of course, this is not

Re: [R] Subset Regression Package

2009-02-17 Thread Hans W. Borchers
Alex Roy alexroy2008 at gmail.com writes: Dear all , Is there any subset regression (subset selection regression) package in R other than leaps? Lars and Lasso are other 'subset selection' methods, see the corresponding packages 'lars' and 'lasso2' and its description in

Re: [R] Outlier Detection for timeseries

2009-02-14 Thread Hans W. Borchers
Pele drdionc at yahoo.com writes: Hello R users, Can someone tell if there is a package in R that can do outlier detection that give outputs simiilar to what I got from SAS below. Many thanks in advance for any help! I guess you are talking about the OUTLIER procedure in SAS that

Re: [R] general inverse solver?

2009-02-10 Thread Hans W. Borchers
system. Hans W. Borchers That being said, without taking away from yacas there is work going on to interface R to a second CAS. On Tue, Feb 10, 2009 at 2:33 AM, Hans W. Borchers hwborchers at googlemail.com wrote: I know that Ryacas is promoted here whenever requests about symbolic

Re: [R] general inverse solver?

2009-02-10 Thread Hans W. Borchers
Gabor Grothendieck ggrothendieck at gmail.com writes: I am not sure what the point of all this is. The point of all this is: 1. Some users of R here on the list apparently would like to have more powerful CAS functionalities than Yacas can provide. 2. Many of the solution hints to Ryacas

Re: [R] general inverse solver?

2009-02-09 Thread Hans W. Borchers
integrated Maxima into its numerical environment in a way that is really useful for numerical and symbolic computations. I could imagine that in a similar way Maxima can be integrated into R bringing the full power of computer algebra to the R community. Hans W. Borchers ABB Corporate Research

Re: [R] Event sequence analysis

2009-02-03 Thread Hans W. Borchers
, Hans Werner Borchers Le 2 févr. 09 à 13:31, Hans W. Borchers a écrit : Dear R help, I am analyzing sequences of events described by time and a unique event tag. And I am searching for recurring patterns where patterns have to show up in a certain time window, e.g. 5 or 10 minutes

[R] Event sequence analysis

2009-02-02 Thread Hans W. Borchers
Dear R help, I am analyzing sequences of events described by time and a unique event tag. And I am searching for recurring patterns where patterns have to show up in a certain time window, e.g. 5 or 10 minutes. Of course, inbetween these events other events may occur. I have applied basket

Re: [R] How to get solution of following polynomial?

2009-01-12 Thread Hans W. Borchers
RON70 ron_michael70 at yahoo.com writes: Hi Ravi, Thanks for this reply. However I could not understand meaning of vectorizing the function. Can you please be little bit elaborate on that? Secondly the package polynomial is not available in CRAN it seems. What is the alternate package?

Re: [R] Execution of a .sce file through R

2008-12-31 Thread Hans W. Borchers
ykank spicyankit4u at gmail.com writes: Dear R Users Does some has any idea about how to execute a scilab file(.sce file) from the Terminal in R. Any kind of guidance would be highly welcomed and appreciated. Have you looked into your Scilab directory? In mine -- under Windows XP

Re: [R] getting ISO week

2008-12-11 Thread Hans W. Borchers
Gabor Grothendieck ggrothendieck at gmail.com writes: format(d, %U) and format(d, %W) give week numbers using different conventions. See ?strptime Gabor, the results of format(aDate, W) appear to be incorrect anyway, see: format(as.Date(2008-01-01), %W) #- 00 There is never a

Re: [R] getting ISO week

2008-12-11 Thread Hans W. Borchers
Gabor Grothendieck ggrothendieck at gmail.com writes: According to the definition in ?strptime (which is not the same as the ISO definition): format(x, %W) returns Week of the year as decimal number (00–53) using Monday as the first day of week (and typically with the first Monday of

Re: [R] Pre-model Variable Reduction

2008-12-09 Thread Hans W. Borchers
Harsh singhalblr at gmail.com writes: Hello All, I am trying to carry out variable reduction. I do not have information about the dependent variable, and have only the X variables as it were. ... I looked for other R packages that allow me to do variable reduction without considering a

Re: [R] How to solve following equation?

2008-12-02 Thread Hans W. Borchers
I need to solve a equation like this : a = b/(1+x) + c/(1+x)^2 + d/(1+x)^3 where a,b,c,d are known constant. Is there any R-way to do that? Multiplying this expression with (1+x)^3 leads to a polynomial equation. I would certainly recommend the 'PolynomF' package here: #

Re: [R] optimization problem

2008-12-01 Thread Hans W. Borchers
Why not use one of the global optimizers in R, for instance 'DEoptim', and then apply optim() to find the last six decimals? I am relatively sure that the Differential Evolution operator has a better chance to come near a global optimum than a loop over optim(), though 'DEoptim' may be a bit slow

Re: [R] Examples of advanced data visualization

2008-11-29 Thread Hans W. Borchers
Tom Backer Johnsen backer at psych.uib.no writes: [...] The question is interesting, but what I have a somewhat negative reaction to is the next passage: Please answer to my e-mail address. In case enough interesting material comes up, I will enter a summary here. It is nice that

[R] Examples of advanced data visualization

2008-11-28 Thread Hans W. Borchers
Dear R-help, I am looking for ideas and presentations of new and advanced data visualization methods. As an example of what I am searching for, the 'Many Eyes' pages at http://manyeyes.alphaworks.ibm.com/manyeyes/ may provide a good paradigm. I would be interested even if it will not be

Re: [R] Which data structure in R can mimic hash in PERL best?

2008-11-15 Thread Hans W. Borchers
Dear all, Which data structure in R can mimic hash in PERL? I'd like to set up a lookup table, which could be accomplished by HASH if using PERL. Which of the data structures in R is the most efficient for lookup table? Thanks for your help. Best regards, Leon The regular answer to

Re: [R] strsplit (regex)

2008-11-11 Thread Hans W. Borchers
stephen sefick ssefick at gmail.com writes: #how do I break these up into first two letters (RM), number, and then the last part #is there an easily accessible regex tutorial on the internet? For regular expressions, the perl man pages at http://perldoc.perl.org/ perlre.html are quite good

Re: [R] ttda and text-mining

2008-11-10 Thread Hans W. Borchers
LE PAPE Gilles lepape.gilles at neuf.fr writes: The ttda package was devoted to text-mining. It seems to be no more available. Did the name change ? Are other packages devoted to text-mining? The 'ttda' package is deprecated as is noted on the Task View Natural Language Processing. The old

Re: [R] How to do knn regression?

2008-09-28 Thread Hans W. Borchers
This is a summary of discussions between Shengqiao Li and me, entered here as a reference for future requests on knn regression or missing value imputation based on a nearest neighbor approach. There several functions that can be used for 'nearest neighbor' classification such as knn, knn1 (in

Re: [R] How to do knn regression?

2008-09-19 Thread Hans W. Borchers
Shengqiao Li shli at stat.wvu.edu writes: Hello, I want to do regression or missing value imputation by knn. I searched r-help mailing list. This question was asked in 2005. ksmooth and loess were recommended. But my case is different. I have many predictors (p20) and I really want try

Re: [R] How to find a shift between two curves or data sets

2008-09-19 Thread Hans W. Borchers
Dear Sébastien, identifying similarity in curves or time series is one of the main tasks in the quite recent field of 'Functional Data analysis' (FDA). See the 2005 book by Silverman from Springer Verlag or the corresponding Web page at url{http://www.psych.mcgill.ca/misc/fda/}. The 'fda'

Re: [R] Convex optimization in R?

2008-09-11 Thread Hans W. Borchers
Hesen Peng hesen.peng at gmail.com writes: Hi my R buddies, I'm trying to solve a specific group of convex optimization in R. The admissible region is the inside and surface of a multi-dimensional eclipse area and the goal function is the sum of absolution values of the variables. Could

Re: [R] Exponential smoothing?

2008-08-19 Thread Hans W. Borchers
Öhagen Patrik Patrik.Ohagen at mpa.se writes: Dear List, I have used all my resources (i.e. help.search) and I still havn't been able to figure out if there is an Exponential Smoothing command in R. A few weeks ago the book Forecasting with Exponential Smoothing by Hyndman et al. has

Re: [R] nonlinear constrained optimization

2008-08-19 Thread Hans W. Borchers
Paul Smith phhs80 at gmail.com writes: Up to my best knowledge, R cannot deal with optimization problems with nonlinear constraints, unless one uses the penalty method. Outside R, Ipopt and Algencan can solve problems like yours, but one needs to program in AMPL and/or C/Fortran. Paul

Re: [R] exporting adaBoost model

2008-08-18 Thread Hans W. Borchers
One way to port these kinds of models between applications is the Predictive Model Markup Language (PMML). The R package 'PMML' supports linear regression, rpart, SVM, and others, not adaBoost. On the other side, not even the Python machine learning library Orange does have an import function

Re: [R] FastICA

2008-08-12 Thread Hans W. Borchers
Maura E Monville maura.monville at gmail.com writes: Is the FastICA R implementation as good as the MatLab Implementation ? I would appreciate talking to someone who has used FastICA for R. The fastICA packages for Matlab and R (and there is even a version for Python) have a common origin at

Re: [R] FastICA

2008-08-12 Thread Hans W. Borchers
Prof Brian Ripley ripley at stats.ox.ac.uk writes: On Tue, 12 Aug 2008, someone with no signature wrote: Maura E Monville maura.monville at gmail.com writes: Is the FastICA R implementation as good as the MatLab Implementation ? I would appreciate talking to someone who has used

Re: [R] FastICA

2008-08-12 Thread Hans W. Borchers
I may not have been as wrong as Prof. Ripley suggested when I wrote The fastICA packages for Matlab and R (...) have a common origin at the Helsinki University of Technology. Please consider the following lines from the 'fastICA' help page (?fastICA): FastICA algorithm Description:

Re: [R] re cursive root finding

2008-08-08 Thread Hans W. Borchers
As your curve is defined by its points, I don't see any reason to artificially apply functions such as 'uniroot' or 'optim' (being a real overkill in this situation). First smooth the curve with splines, Savitsky-Golay, or Whittacker smoothing, etc., then loop through the sequence of points

<    1   2   3   >