Re: [R] eigen() may fail for some symmetric matrices, affects mvrnorm()

2005-05-01 Thread Prof Brian Ripley
I do not begin to understand this. mvrnorm() uses EISPACK=TRUE and always has done, explicitly since there was such an argument to eigen() (since 15 Jan 2003: the argument was introduced in R 1.7.0 in April 2003). Just take a look at the code. The claimed `patch' (which is not a patch at all but

Re: [R] legend(): how to put variable in subscript?

2005-05-01 Thread Peter Dalgaard
Aleksey Naumov [EMAIL PROTECTED] writes: Dear List, I would like to plot a simple legend with two math expressions, e.g. plot(0) legend(1, 0.5, expression(sigma[i], sigma[j])) The difficulty is that i and j should be variables rather than strings i and j. In other words I'd like to

Re: [R] User-defined random variable

2005-05-01 Thread Matthias Kohl
Achim Zeileis wrote: On Sat, 30 Apr 2005, Peter Dalgaard wrote: Paul Smith [EMAIL PROTECTED] writes: Dear All I would like to know whether it is possible with R to define a discrete random variable different from the ones already defined inside R and generate random numbers from that

[R] Problems using html and dvips of Hmisc

2005-05-01 Thread Marco Chiarandini
Dear all, I am experiencing the two following problems trying to produce a comfortable output using the package Hmisc. If I use: html(latex(D,file=here.tex),file=here.html) where D is a data.frame the function html does not write in here.html but only in a temporary file somewhere else in my

Re: [R] eigen() may fail for some symmetric matrices, affects mvrnorm()

2005-05-01 Thread Jouni Kerman
Hi all, Professor Ripley is right. This problem does not affect newer versions of mvrnorm. We had an older version of mvrnorm which did not set EISPACK=TRUE as it should. Sorry for the confusion. Thank you for the correction. Jouni Kerman On May 1, 2005, at 6:08 AM, Prof Brian Ripley [EMAIL

[R] International Association for Statistical Computing Conference (session on Diagnostics)

2005-05-01 Thread Stephen Nyangoma
Dear Everyone, I am organizing a session on Diagnostics (Regression etc.) for the IASC on Computational Statistics and Data Analysis which will take place at the Amathus Beach Hotel in Limassol Cyprus, October 28-31, 2005. Those of you who would like to attend and present a paper may submit

[R] Roots of quadratic system.

2005-05-01 Thread John Janmaat
Hello, I have a system of quadratic equations (results of a Hamiltonian optimization) which I need to find the roots for. Is there a package and/or function which will find the roots for a quadratic system? Note that I am not opimizing, but rather solving the first order conditions which come

[R] simulate zero-truncated Poisson distribution

2005-05-01 Thread Glazko, Galina
Dear All I would like to know whether it is possible with R to generate random numbers from zero-truncated Poisson distribution. Thanks in advance, Galina __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] Programa de email e contatos

2005-05-01 Thread Programas de email
Convidamos você a experimentar o Magic Seven, o melhor programa de envio de emails para listas de contatos. Com o Magic Seven você pode enviar emails personalizados e com recursos multimídia, além de ser o programa mais rápido e confiável do mercado. O Magic Seven está no mercado há mais de 5

RE: [R] simulate zero-truncated Poisson distribution

2005-05-01 Thread Ted Harding
On 01-May-05 Glazko, Galina wrote: Dear All I would like to know whether it is possible with R to generate random numbers from zero-truncated Poisson distribution. Thanks in advance, Galina Maybe someone has written an efficient function for this. If not, you anyway have at least

[R] opimization problem

2005-05-01 Thread Gottfried Gruber
hi, i want to execute the following opimization problem: max r*w s.t.: w*z=1 # sum of w is 1 r, w are [nx1] vectors, z is a [nx1] vector consisting of 1 so far so good, works fine with lp the problem arises with the additional restriction w' * V * w where V is a [nxn] matrix how can i

Re: [R] simulate zero-truncated Poisson distribution

2005-05-01 Thread Peter Dalgaard
(Ted Harding) [EMAIL PROTECTED] writes: On 01-May-05 Glazko, Galina wrote: Dear All I would like to know whether it is possible with R to generate random numbers from zero-truncated Poisson distribution. Thanks in advance, Galina . (B) This has a deeper theoretical

[R] Partial F-test comparing full and reduced regression models

2005-05-01 Thread Jim Milks
Dear all: I have a regression model that has collinearity problems (between three regressor variables). I need a F-test that will allow me to compare between full (with all variables) and partial models (minus 1= variables). The general F-test formula I'm using is: F = {[SS(full model) -

Re: [R] Partial F-test comparing full and reduced regression models

2005-05-01 Thread Peter Dalgaard
Jim Milks [EMAIL PROTECTED] writes: Dear all: I have a regression model that has collinearity problems (between three regressor variables). I need a F-test that will allow me to compare between full (with all variables) and partial models (minus 1= variables). The general F-test formula

Re: [R] simulate zero-truncated Poisson distribution

2005-05-01 Thread Spencer Graves
Hi, Peter: Your solution is simple, elegant, very general yet sufficiently subtle that I (and I suspect many others) might not have thought of it until you mentioned it. It is worth remembering. Thanks. spencer graves Peter Dalgaard wrote: (Ted Harding) [EMAIL PROTECTED] writes: On 01-May-05

Re: [R] simulate zero-truncated Poisson distribution

2005-05-01 Thread Ted Harding
On 01-May-05 Peter Dalgaard wrote: (Ted Harding) [EMAIL PROTECTED] writes: Something like the following code will do the job vectorially: n-1000 # desired size of sample T-3.5# pre-truncation mean of Poisson U-runif(n) # the uniform sample t = -log(1 - U*(1 -

Re: [R] Roots of quadratic system.

2005-05-01 Thread Kjetil Brinchmann Halvorsen
John Janmaat wrote: Hello, I have a system of quadratic equations (results of a Hamiltonian optimization) which I need to find the roots for. Is there a package and/or function which will find the roots for a quadratic system? Certainly you cxould use solve, see ?solve Alternatively you could

RE: [R] Roots of quadratic system.

2005-05-01 Thread Bill.Venables
Are you looking for a unique solution or families of solutions? Can't you turn a root-finding problem for a system of equations with a unique solution into an optimisation problem, anyway? E.g. You want to solve f1(x) = g1 f2(x) = g2 ... Why not optimise L(x) = (f1(x) - g1)^2 + (f2(x) -

Re: [R] how to replace text...

2005-05-01 Thread Sean Davis
Jonathan, Just for future reference, you can send replies to the list--everyone benefits then. I'm not sure how pricelts works, but you can type ?pricelts to see. It is possible that pricelts takes a vector of stock names. You need to look at the Value part of the help to see what pricelts

[R] eigenvalues of a circulant matrix

2005-05-01 Thread Globe Trotter
Hi, It is my understanding that the eigenvectors of a circulant matrix are given as follows: 1,omega,omega^2,,omega^{p-1} where the matrix has dimension given by p x p and omega is one of p complex roots of unity. (See Bellman for an excellent discussion on this). The matrix created by

[R] help : Bisect( )

2005-05-01 Thread mingan
I am using the package of Icens but for the function Bisect I do not know what ndir is it just says the direction of pvec, but what kind of values should I use ? thanks for reply pls give me an example.. Bisect(tA, pvec, ndir, Meps, tolbis=1e-07) Arguments: tA: The transpose of the