similar to: Inaccuracy in svd() with R ubuntu package

[I thought I'd submitted this bug report some time ago, but it's never showed up on the bug tracking system, so I'm submitting again.] qr.solve() gives incorrect results when dealing with complex matrices or with qr objects that have been computed with LAPACK=TRUE, whenever the b argument has more than one column. This bug flows from qr.coef(), which has a similar problem. I believe

Dear R Development Team, I have encountered the following difficulty in using the function polyroot under either NT4.0 (R version 1.2.1) or linux (R version 0.90.1). In the provided example, the non-zero root of c(0,0,0,1) depends on the results of the previous call of polyroot. R : Copyright 2001, The R Development Core Team Version 1.2.1 (2001-01-15) R is free software and comes with

Hi, I'm responding to a post about finding roots of a cubic or quartic equation non-iteratively. One obviously could create functions using the explicit algebraic solutions. One post on the subject noted that the square-roots in those solutions also require iteration, and one post claimed iterative solutions are more accurate than the explicit solutions. This post, however, is about

Hi I am having difficulty with eigen() on R-devel_2006-01-05.tar.gz Specifically, in R-2.2.0 I get expected behaviour: > eigen(matrix(1:100,10,10),FALSE,TRUE)$values [1] 5.208398e+02+0.000000e+00i -1.583980e+01+0.000000e+00i [3] -4.805412e-15+0.000000e+00i 1.347691e-15+4.487511e-15i [5] 1.347691e-15-4.487511e-15i -4.269863e-16+0.000000e+00i [7] 1.364748e-16+0.000000e+00i

Thanks for telling me that you could not get my message, I hope this work better... so my question was: I built a population matrix to which I applied the fonction eigen in order to find the main parameters about my population. I know that the first eigen value correspond to lambda or exponential growth rate of my population. My problem is that I want to have the 95% confidence interval of the

I'm doing the time-series data analysis using FFT. After fft(), I used spectrum() to calculate the power spectrum, which is used to define the main frequency. Now I want to get the phase information according to this frequency. For example, input 12 time-point data,

El problema del módulo es que pierde el signo. En tu caso sale igual porque has invertido el signo del coeficiente en el polinomio (en realidad se me pasó a a mí advertir que el término independiente debe ir con signo negativo): .> polyroot(z=c(0.5,0,0,0,0,1)) [1] 0.7042902+0.5116968i -0.2690149+0.8279428i -0.2690149-0.8279428i [4] 0.7042902-0.5116968i -0.8705506+0.0000000i .> .>

Hi, with small matrices eigen works as expected: > eigen(cbind(c(1,4),c(4,7)), only.values = TRUE) $values [1] 9 -1 $vectors NULL > eigen(cbind(c(1,4),c(4,7))) $values [1] 9 -1 $vectors [,1] [,2] [1,] 0.4472136 -0.8944272 [2,] 0.8944272 0.4472136 > eigen(cbind(c(1,-1),c(1,-1))) $values [1] -3.25177e-17+1.570092e-16i -3.25177e-17-1.570092e-16i $vectors

Hello, I have 3 columns : a, b and a*b I would like to find the pair (a,b) so that a*b is the minimum but not from the points I measured but from the fit of the curve (I have more points that the ones given below but I fit only on this part because I know that the minimum a*b is in this interval). I thought doing it this way : - to fit a*b=f(a) abfit<-lm(ab ~ poly(a,8,raw=T)) - to use the

Hi, All: Is there a simple way to detect complex conjugates in the roots returned by 'polyroot'? The obvious comparison of each root with the complex conjugate of the next sometimes produces roundoff error, and I don't know how to bound its magnitude: (tst <- polyroot(c(1, -.6, .4))) tst[-1]-Conj(tst[-2]) [1] 3.108624e-15+2.22045e-16i

Hi all, Would anybody care to suggest a third party SATA-RAID card that works out of the box with CentOS 6, without having to jump through hoops to make it work? The card should preferably be able to connect ten harddrives, but I guess three four-port cards should work as well. There's no need for anything fancy really, as long as I can create a single big software-raid on it at CentOS

Hi, is it save to use ext2 / ext3 hashed b-trees feature with linux kernel >= 2.4.25 ? thanx, Gregor -- echo '16i[q]sa[ln0=aln100%Pln100/snlbx]sbA0D3F204445524F42snlbxq'|dc

When upstream released 5.1, everybody wanted to test a new kernel parameter that could adjust the system clock rate at boot time to something else than the standard 1000Hz clock rate. A lot of testings has been done (thanks to Akemi Yagi for her great work) and you can see the results here : http://bugs.centos.org/view.php?id=2189 As you can read at the bottom of the comments, it seems it was a

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi! I'm new to wine, but when I start stronghold (MouseGrab enabled), i can only scroll up and left the screen, but not right and down... any ideas? please email me directly! THX! Mike - -- < "Enjoyment predicts efficiency" -- Eric Steven Raymond > <echo

Hi all! I tried to get WarCraft 3 to run with latest CVS and lastest WineX CVS, but it didn't start at all... I followed the descriptions in the FAQ, so I got this NoCD patch, but something seems to be still wrong... who has got it to work? are you using a stand-alone wine or with an existing/original Windows? Once, somehow, it already started but only played the sound, but nothing was on

The CentOS wiki has a page regarding the nomachine nx client setup ... But i think that a page describing the server part, including ssh key authentication, can be interesting for people wanting to use nx/freenx. Because by default, it will modify sshd_config to allow password authentication (even if configured manually to not do so ...) Range : can you create a new page ? or do we have to edit

Hi, i have two twofold e2image related questions: 1) The man page mentions "cp (1)"'s --sparse=always option. I wonder if there are sparse aware tools which a) display the "real" amount of disk space occupied by a sparse file b) compress sparse files (other than compressing very thightly several GiB of zeros) 2) "E2fsck, mke2fs, etc.

Greetings, We''re trying out a new JBOD here. Multipath (mpxio) is not working, and we could use some feedback and/or troubleshooting advice. The OS is oi151a7, running on an existing server with a 54TB pool of internal drives. I believe the server hardware is not relevant to the JBOD issue, although the internal drives do appear to the OS with multipath device names (despite the fact

Hi, I'm trying to build a bridge using a 486DX2/66 with two 10 Mbit Ethernet NICs. The machine (tsushima) is running Debian stable and kernel 2.4.26. Previously I used it as a router, so I know the hardware (NICs) is working. The NICs are as follows: eth0: WD80x3 at 0x280, 00 00 C0 0A 2C 2F WD8003-old, IRQ 10, shared memory at 0xd0000-0xd1fff. eth16i.c: v0.35 01-Jul-1999 Mika Kuoppala

If I have matrics as follows: > a <- c(1,1,0,0) > b <- c(4,4,0,0) > c <- c(3,5,5,6) How can I use R code to solve the equation ax^2+bx+c=0. thanks! yuying shi [[alternative HTML version deleted]]