Re: [R] eigenvalues of a circulant matrix

2005-05-01 Thread Prof Brian Ripley
On Sun, 1 May 2005, someone who didn't give his name wrote: 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 exc

[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

[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 th

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

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) - g2)^2

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 go

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(

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-0

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

[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) - SS(

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

[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 include

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

[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 a

[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 d

[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 fr

[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 th

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 PR

[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

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 user

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 wo