similar to: Catching errors from solve() with near-singular matrices

dear R experts: ?I am writing my own little newey-west standard error function, with heteroskedasticity and arbitrary x period autocorrelation corrections. ?including my function in this post here may help others searching for something similar. it is working quite well, except on occasion, it complains that Error in solve.default(crossprod(x.na.omitted, x.na.omitted)) : system is

thank you, achim. I will try chol2inv. sandwich is a very nice package, but let me make some short suggestions. I am not a good econometrician, so I do not know what prewhitening is, and the vignette did not explain it. "?coeftest" did not work after I loaded the library. automatic bandwidth selection can be a good thing, but is not always. as to my own little function, I like the

I maintain the geometry package, which integrates the Qhull C library (http://qhull.org) into R. The Qhull function I hook into requires an open FILE handle as one of its arguments. I had set this file handle to stdout, but now R check NOTEs the presence of stdout, and the CRAN maintainers asked me to get rid of these NOTEs. Including the following defines means the checks are passed on CRAN:

I?m new to R so maybe this issue has been asked before and I still could not read the complete set of past messages sent to the list. I found a weird behabiour that I will explain with a simple example. Lets consider the following block of commands: > x <- diag(c(1,4,10)) > x [,1] [,2] [,3] [1,] 1 0 0 [2,] 0 4 0 [3,] 0 0 10 > invx <- x^-1 > invx

Hi uRsers, when inverting a 2 by 2 matrix using solve, I encountered a error message: solve.default(sigma, tol = 1e-07) : system is computationally singular: reciprocal condition number = 1.7671e-017 and then I test the determinant of this matrix: 6.341393e-06. In my program, I have a condition block that whether a matrix is invertible like this: if(det(sigma)<1e-7) return NULL;

With a matrix such as C I do not have any problem: >library(gmp) > C V1 V2 V3 V4 V5 V6 V7 [1,] 1 0 0 0 1 0 0 [2,] 0 1 0 0 0 1 0 [3,] 0 0 1 0 0 0 1 [4,] 0 0 0 1 0 0 0 [5,] 0 0 0 0 1 0 0 [6,] 0 0 0 0 0 1 0 [7,] 0 0 0 0 0 0 1 > solve.bigz(C) [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] "1" "0" "0"

Dear Arne, thank you very much upfront for your time that you invest in answering my question. I'm using your systemfit R package to analyze if companies risk appetite has a relationship to their performance and how it evolves over time. We collected data for ca. 80 companies, from 6 segments over 15 years that gives us a final data set with 1,312 firm-year observations. We have set up 3

Hi there everyone, I have the following model (this is naturally a simplified version just for showing my problem, in case you're wondering this is a translog cost function with the associated cost share equations): C ~ á + â1 log X + â2 log Y + ã1 log Z + ã2 log XX C1 ~ â1 + â2 log YY + ã1 log ZZ Then I have some restrictions on the coefficients, namely that the sum of â equal 1 and the

Bug #2150: Don't call VBEGetSetDACPaletteFormat when in direct colour or YUV modes. (David Sterratt) CVS tag: xf86-video-vesa-1_2_1 http://xorg.freedesktop.org/releases/individual/driver/xf86-video-vesa-1.2.1.tar.gz: MD5: 1bf4a36af8422245bcdd99c2aa2f8f4f xf86-video-vesa-1.2.1.tar.gz SHA1: 92323676608ff9ec357c8d94f14133794943b542 xf86-video-vesa-1.2.1.tar.gz

Perhaps someone else could confirm this, I have an APC UPS with a built-in network management card with environmental sensor for temp/humidity. The OID index for the sensors appears to be .1 (the last component for each of the OIDs listed below). The original source code has them at .0 drivers/apccmib.h #define APCC_OID_IEM_TEMP ".1.3.6.1.4.1.318.1.1.10.2.3.2.1.4.1" #define

Hello list. I'm contemplating buying new batteries and putting an older APC Smart-UPS 3000XL back into service. Would that work well with NUT? Any drawbacks to consider? Thank you, Tobias -- Tobias Balle-Petersen Kontrapunkt Group / System Administrator

We recently had to replace an old APC SmartUPS 3000 with a newer 3000VA (the 2U one). It has both serial and USB in the back (with an RJ-45 on the end of the serial cable). I can't get the apcsmart driver to connect. I figured it should work the exact same as the old UPS, with /dev/ttyS0 as the port (I'm using the same port on the server that the last one was plugged into). I even

Hi, I have a set of matrix data named “invest” consists of 450 observations (75 countries, 6 years) with 7 variables (set as I, pop, inv, gov, c, life, d; which each is “numeric[450]”). The procedure is modify from code provided by B.E. Hansen at http://www.ssc.wisc.edu/~bhansen/progs/ecnmt_00.html. *Then the variable is being transformed to* y <- lag_v(i,0) cf <-

Hello, I am running a simple plspm for a class project due later today and I am receiving the following error despite following along exactly with Gaston Sanchez's directions in PLS Path Modeling with R: Error in solve.qr(qr(X.blok), Z[, j]) : singular matrix 'a' in 'solve' I would greatly appreciate any help resolving this matter. I got the same error after changing the

Hi there! I´m a noob when it comes to R and I´m using it to run statisc analysis. With the code for ARIMA below I´m getting this error: Error in solve.default(res$hessian * n.used) :Lapack routine dgesv: system is exactly singular The code is: > s.ts <- ts(x[,7], start = 2004, fre=12) > get.best.arima <- function (x.ts, maxord=c(1,1,1,1,1,1)) + { + best.aic <- 1e8 + n <-

Hi R-users: When I run the R script (as the following), I got the error message: Error in solve.default(C[ind, ind]) : Lapack routine dgesv: system is exactly singular. Any help is appreciated. Ying library(sem) R.pw <- matrix(c( 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2137356, 0.2137356, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,

Hello everyone, I'm trying to estimate the parameters of the returns series attached using the GARCH code below, but I get the following error message: Error in solve.default(Hessian) : system is computationally singular: reciprocal condition number = 0 Error in diag(solve(Hessian)) : error in evaluating the argument 'x' in selecting a method for function 'diag' Can

* Is there a way to configure NUT so that it recognizes the huge additional battery pack that has been attached to the 2200XL? The 51minutes runtime is a number consistent with the load we have placed on a stock standalone 2200XL. With the extra large battery pack we've added, the runtime should be more than 2.5 hours. * I see the "Bizzare configuration" graphic on the

Dear R, If I have remembered correctly, a square matrix is singular if and only if its determinant is zero. I am a bit confused by the following code error. Can someone give me a hint? > a <- matrix(c(1e20,1e2,1e3,1e3),2) > det(a) [1] 1e+23 > solve(a) Error in solve.default(a) : system is computationally singular: reciprocal condition number = 1e-17 Thanks in advance! Feng --

Here's a patch (against 2.2.2) which adds battery.temperature to the ietf-mib snmp driver. While looking at ietfmib.h I noticed that the value for IETF_OID_IN_LINES is actually that for upsOutputNumLines.0. Does anyone know whether that's deliberate or a typo? -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: application/x-patch