search for: dgesv

Dear R-devel, I am trying to use Fortran DGESV subroutine into C. Here it is the relevant part of the C file I am currently writing #include<stdio.h> #include<R.h> #include<Rmath.h> #include<math.h> void F77_NAME(DGESV)( int*, int*, double*, int*, int*, double*, int*, int*); void solve( int *p, double *A, double *Ainv...

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 <- length(x.ts) + for (p in 0:maxord[1]) for (d in 0:maxord[2]) for (q in 0:maxord[3]) + for (P in 0:maxord[4]...

...I am trying to do Log_rank Survival analysis, I have included tables and str command, is it a factor/integer problem? If so how do I correct this, as all my attempt to recode the data have failed. > survdiff(Surv(f2$days.alive , f2$censored)~group, data=f2) Error in drop(.Call("La_dgesv", a, as.matrix(b), tol, PACKAGE = "base")) : Lapack routine dgesv: system is exactly singular > head(f2) group days.alive censored 1 PRI_CAS_5_NODU 1826 1 2 PRI_CAS_5_NODU 1488 0 3 PRI_CAS_5_NODU 1826 1 4 PRI_CAS_5_NODU...

....names(Dataset) <- as.character(Dataset$cntry) Dataset$cntry <- NULL cl <- cmeans(Dataset, 4, 20, verbose=TRUE, method="cmeans") resultindexes <- fclustIndex(cl, Dataset, index= "all") resultindexes I am receiving the following error message: ERROR: Lapack routine dgesv I have tried running the example for fclustIndex and it works fine. The code above is based on the example from the manual (below). I also checked that I have the latest version installed. I have also had no trouble running similar commands from other packages which leads me to believe...

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, 0.09792217, 0.09792217, -0.2361812, 1, 0, 0,...

...silly question. But I really couldn't understand why cov() and ccov() don't work for my calculation of covariance matrix. a <- matrix(1:8, 2, 4) a [,1] [,2] [,3] [,4] [1,] 1 3 5 7 [2,] 2 4 6 8 > ccov(a) Error in solve.default(cov, ...) : Lapack routine dgesv: system is exactly singular I also tried colume bind, but it still doesn't work. > a1<- a[,1] > a2 <- a[,2] > a3<- a[,3] > a4 <- a[,4] > acomb <- cbind(a1,a2,a3,a4) > acomb a1 a2 a3 a4 [1,] 1 3 5 7 [2,] 2 4 6 8 > ccov(acomb) Error in solve.defaul...

...tAttrib(A, R_DimSymbol), INTSXP)); Bdims = INTEGER(coerceVector(getAttrib(B, R_DimSymbol), INTSXP)); int nrA, ncB; double *Ap, *Bp; nrA = Adims[0]; ncB = Bdims[1]; Ap = REAL(A); Bp = REAL(B); /* int info, *ipiv = (int *) R_alloc(nrA, sizeof(int)); */ /* F77_CALL(dgesv)(&nrA, &ncB, Ap, &nrA, ipiv, Bp, &nrA, &info); */ mysolve2(Ap, &nrA, Bp, &ncB, &tol); UNPROTECT(2); return B; } void mysolve2(double *A, int *nrA, double *B, int *ncB, double *tolin) { int info; int *ipiv = (int *) R_alloc((int) nrA, sizeof(int));...

Dear R-experts, I am trying to fit my model but I couldn't because of this error. here is the error "Error in solve.default(dial(m) -A) : Lapack routine dgesv: system is exactly singular" thank you all. sonam

Hi, I'm trying to use the solve() function in R to invert a matrix. I get the following error, "Lapack routine dgesv: system is exactly singular" However, My matrix doesn't appear to be singular. [,1] [,2] [,3] [,4] [1,] 0.99252358 0.93715047 0.7540535 0.4579895 [2,] 0.01607797 0.09616267 0.2452471 0.3088614 [3,] 0.09772828 0.58451468 1.4907090 1.8773815 [4,] -0.01000000...

...default. Can this precision be set to values above 22? The problem occurs especially when trying to perform 2SLS with the 'systemfit' package. There appears always an error message like the following from the inverting routine: solve(tcross) Error in solve.default(tcross) : Lapack routine dgesv: system is exactly singular Or is there another source of error? We would like to embed R-routines in a non-commercial web application for which we need to employ 'systemfit' together with individual, user-submitted data. So the problem needs a general solution and not a special one for th...

...anted to assess autocorrelation characteristics and tried to run a Durbin-Watson using: library(lmtest) dwtest(formula=Index1stdiff~Comps1stdiff,alternative=c("greater")) I am getting the following error: Error in solve.default(crossprod(X), tol = tol) : Lapack routine dgesv: system is exactly singular Can anyone assess why my attempt crashes? I tried this earlier on what I thought was similar data and was returned an answer for "d" and a "p-value". Thanks in advance for your assistance. Rick

When trying to fit a generalized linear mixed model using glmmPQL: > fit0 <- glmmPQL(ifelse(response=="A",1,0)~gender,data=set1, random=~1|subject,family=binomial) iteration 1 Error in solve.default(pdMatrix(a, fact = TRUE)) : Lapack routine dgesv: system is exactly singular Could this be occuring because the paired responses for each subject are always the same? If this were the case, what would be the best way to handle this situation? I replaced the response variable with fake data that are not always the same for each pair and then I d...

Dear R group, I have to solve a hessian matrix 40*40, called M, in order to obtain the standart deviations of estimators. When I use the function solve(M), I have the following error message: "Error in solve.default(M) : Lapack routine dgesv: system is exactly singular" Do you know an alternative approach which could succeed? I have found some information about the function "ginv" (library MASS), has someone already used it in such a situation? Thank you very much for your responses. Eve MATHIEU.

Dear All, While inverting a matrix the following error appears on my console: Error in solve.default(my_matrix) : Lapack routine dgesv: system is exactly singular With this respect, I have been replacing the solve() function with ginv(): the Moore-Penrose generalized inverse of a matrix. These are the questions I would like to ask you: 1. Would you also replace solve() with ginv() in these scenarios and using R? 2. Or is t...

...e error. mydata <- read.table("C:/dataForR/radiology/wordFrequencies.csv", header=TRUE, sep=",") > dim(mydata) [1] 982 924 mydataN<-mydata[,-923] > dim(mydataN) [1] 982 923 cor(mydataN) factanal(mydataN, factors=3) Error in solve.default(cv) : Lapack routine dgesv: system is exactly singular Thank you for looking at this! -- View this message in context: http://www.nabble.com/factanal-error-tp23147009p23147009.html Sent from the R help mailing list archive at Nabble.com.

...ot;,m=1.7,rate.par=NULL,weights=1) I would like to calculate the value of all the fuzzy validity measures using the flcustIndex function. However it returned the following error: > fclustIndex(ra.flcust6,ra,index="all") Error in solve.default(scatter[, , i]) : Lapack routine dgesv: system is exactly singular Please could anyone explain what this means and what I have done wrong? thanks Caroline Wallis [[alternative HTML version deleted]]

...7,] 2 5 [18,] 2 5 [19,] 2 5 [20,] 2 5 [21,] 2 6 [22,] 2 6 [23,] 2 6 [24,] 2 6 > Y = A + B + rnorm(a*b*n) > fr = data.frame(Y=Y,A=as.factor(A),B=as.factor(B)) > afit=aov(Y ~ A / B - 1,fr) > X=model.matrix(afit) > solve(t(X) %*% X) Error in solve.default(t(X) %*% X) : Lapack routine dgesv: system is exactly singular

...} z=which(mat!=0) mat[z]=0 hitt=SVAR(var , estmethod = "scoring", Amat = mat, Bmat = NULL,max.iter = 100, maxls = 1000, conv.crit = 1.0e-8) hit[j]=hitt$LR$statistic } } There will be an error comment,e.g. : Error in solve.default(BinvA) : Lapack routine dgesv: system is exactly singular" But I still want to continue to compute the SVAR.. How can I do it? Thanks for answering it.. Regards, Arif _________________________________________________________________ Facebook. k-basics.aspx?ocid=PID23461::T:WLMTAGL:ON:WL:en-id:SI_SB_2:092...

Hey Everyone, I'm running R 2.4.0 on Debian etch 4.0, and I'm trying to call some LAPACK functions from the C code in my package. Actually, to be honest I'm not really having trouble using commands such as La_dgesv from within my C code, but I do get warning when compiling the package saying: ***.c: In function '***': ***.c:37: warning: implicit declaration of function 'La_dgesv' ***.c:37: warning: assignment makes pointer from integer without a cast I tried using: #include <Rmodules/Rla...

...em package, I also have an error message:   supplyreg <- tsls(tscn ~ tsyn + tsqn + tsksn + tsucp, ~ tsupp1 + tsupp2 + tsupp3 + tsyus1 + tsupf + tsd1 + tsyud1 + tspo + tssn + tsupp + tsus + tsus1 + tsyn + tsqn + tsksn) summary(supplyreg)   Error in solve.default(crossprod(Z))"Lapack routine dgesv: system is exactly singular"   Do the two error messages mean the same think?   Thank you in advance for your help. [[alternative HTML version deleted]]