search for: det

...ates? Thanks a lot for your help. Maomao ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ > p<-optim(c(0,0,0), f, method ="BFGS", hessian =T, y=y,X=X,W=W) There were 31 warnings (use warnings() to see them) > warnings() Warning messages: 1: In log(det(I_N - pd * wd - po * wo - pw * ww)) : NaNs produced 2: In log(det(I_N - pd * wd - po * wo - pw * ww)) : NaNs produced 3: In log(det(I_N - pd * wd - po * wo - pw * ww)) : NaNs produced 4: In log(det(I_N - pd * wd - po * wo - pw * ww)) : NaNs produced 5: In log(det(I_N - pd * wd - po * wo - pw *...

Hello, I am trying to write an R package for doing analysis of speaker recognition systems. The big thing in this line of research is a DET-plot (detection Error Trade-off, a ROC plot with qnorm() warped axes). My approach has been to make a class "det" and also a function "det()" that will prepare the data into the right class. Now loading the library doesn't like me re-defining the determinant function &quo...

Hello, Currently, determinant(A) calculates the determinant of 'A' by factorizing A=LU and computing prod(diag(U)) [or the logarithm of the absolute value]. The factorization is done by LAPACK routine DGETRF, which gives a status code INFO, documented [1] as follows: *> INFO is INTEGER *>...

I just detected, that det() is not working on complex matrices any more, due to the fix to the bug reports noted above. I am not happy with this, as determinants are perfectly usable on complex matrices. AFAIUI the bugs resulted from less than optimal behaviour of qr() in certain cases. IMHO this is due to t...

I have been running into difficulties with dispatching on an S4 class defined in the SparseM package, when the method calls are inside a function passed as the f= argument to optimize() in functions in the spdep package. The S4 methods are typically defined as: setMethod("det","matrix.csr", function(x, ...) det(chol(x))^2) that is within setMethod() rather than by name before the setMethod(). When called from within functions passed as the f= argument to optimize, the S3 generics for det() and chol() get picked up, not the S4 generics for the S4 Spars...

Given the following data, I am plotting log.det ~ norm.beta, where the points depend on a parameter, lambda (but there is no functional form). I want to find the (x,y) positions along this curve corresponding to two special values of lambda lambda.HKB <- 0.004275357 lambda.LW <- 0.03229531 and draw reference lines at ~ -45 degrees (or...

...st, I have recently written an R package that solves a linear least squares problem, and computes the multivariate normal density function. The bulk of the code is written in C, with interfacing code to the BLAS and Lapack libraries. The motivation here is speed. I ran into a problem computing the determinant of a symmetric matrix in packed storage. Apparently, there are no explicit routines for this as part of Lapack. While there IS an explicit routine for this in Linpack, I did not want to use the older library. Also, right before I needed the determinant of the matrix A, I had used the Lapack...

Full_Name: Krzysztof Podgorski Version: R version 2.4.1 (2006-12-18) OS: Windows XP Submission from: (NULL) (130.235.3.79) The function ''det'' works improperly for a singular matrix and returns a non-zero value even if ''solve'' reports singularity. The matrix is very simple as shown below. A <- diag(rep(c(64,8), c(8,8))) A[9:16,1] <- 8 A[1,9:16] <- 8 det(A) #[1] -196608 solve(A) #Error in solve.defaul...

Full_Name: Armin Roehrl Version: 1.4.0(2001-12-19) OS: Linux; Suse 7.3 Submission from: (NULL) (217.84.18.204) Hello, On a given matrix det returns different results whether I use the method "qr" or "eivenvalues". The matrix: > m2 [,1] [,2] [,3] [1,] 9.822616e+09 3.841723e+09 79790.09 [2,] 3.841723e+09 1.502536e+09 31251.82 [3,] 7.979009e+04 3.125182e+04 64156419.36 > det(m2...

Hello, When I use det() and qr() on complex matrices the result is in some cases indeterministic. The documentation speaks of numeric matrices (and not of complex matrices) but det() uses qr() which should handle complex matrices correctly. I've also tried using only qr() with similar results. det() returns a valu...

...e(1:100, m^2), m, m) C[[i]] <- matrix(sample(1:100, m^2), m, m) } } result <- list(m=m, A=A, B=B, C=C) class(result) <- "foo" result } I'd like to define an as.data.frame function for class "foo" object that applies either det() or tr() to the matrices in the case m>1. The function below works, but I'd rather define the argument FUN as c(det, tr) and use something like match.arg(FUN) to determine which was supplied. Is there some easy way to do this? # matrix trace tr <- function(M) sum(diag(M)) # why can&...

Hello, Is there a way to compare the relative relevance of fixed and random effects in mixed models? I have in mind measures of effect size in ANOVAs, and would like to obtain similar information with mixed models. Are there information criteria that allow to compare the relevance of each of the effects in a mixed model to the overall fit? Thank you, Bruno

Hi, Which R function I should use to obtain determinant of a matrix with real(and complex) numbers? Kalin --------------------------------- Never miss a thing. Make Yahoo your homepage. [[alternative HTML version deleted]]

...I should report to the SAMBA team? or is this a known bug? (I've not seen this specific problem referenced) This problem is driving us crazy! Help! -David *-Attachment--- log.smbd [2005/04/05 04:39:22, 2] auth/auth.c:check_ntlm_password(305) check_ntlm_password: authentication for user [det-liv-cd1] -> [det-liv-cd1] -> [det-liv-cd1] succeeded [2005/04/05 04:39:22, 1] smbd/service.c:make_connection_snum(642) mpc-cdirect1 (10.1.5.84) connect to service test initially as user det-liv-cd1 (uid=1001, gid=100) (pid 28067) [2005/04/05 04:39:22, 2] smbd/dosmode.c:unix_mode(60) unix_...

Greetings; I've posted the following to R's bug tracking system (at https://bugs.r-project.org/bugzilla/show_bug.cgi?id=17210 ) and Martin Maechler requested that I post to this list as well. If I start R from the command line with --vanilla, then repeatedly execute the following line: det(matrix(nrow=10,ncol=10)) ... I eventually get a crash, with error: *** Error in `/usr/lib/R/bin/exec/R': malloc(): memory corruption: 0x0000000002399400 *** The specific address in memory that is referenced varies. The number of times I need to execute the above line before getting a crash a...

-----Original Message----- From: R-SIG-Debian [mailto:r-sig-debian-bounces at r-project.org] On Behalf Of Rolf Turner Sent: Thursday, 19 January 2017 10:11 AM To: Ian Erickson Cc: r-sig-debian at r-project.org Subject: Re: [R-sig-Debian] [FORGED] Taking determinant of a matrix of NAs results in intermittent memory corruption >On 19/01/17 11:54, Ian Erickson wrote: >> Greetings; I've posted the following to R's bug tracking system (at https://bugs.r-project.org/bugzilla/show_bug.cgi?id=17210 ) and Martin Maechler requested that I post...

Hi, I have the following hypothetical optimization problem: -det(t(x-A%*%x1)%*%(x-A%*%x1)) where A,x,x1 are matrices. A coefficients and x and x1 are variable matrices or vectors. I tried to apply optim and nlm functions but I kept receive the following massage: Error in A%*%x1 : non-conformable arguments. The massage appears even the -det() can be calculated an...

...problem. Suppose I have a matrix of size nxn say, generated as follows   z<-matrix(rnorm(n*n,0,1),nrow=n)   I want to write a function such that for i in 1:n, I will remove the rows and columns corresponding to i (so, will be left with n-1*n-1 submatrix in each cases). Now I need the sum of the determinant of each of this submatrices. As an example, if n=3, it means I will have det(1strow and 1stcolum removed) + det(2ndrow and 2ndcolum removed) + det(3rdrow and 3rdcolum removed).   Any help will be appreciated. Thanks   John [[alternative HTML version deleted]]

I'm confused by this: > identical(function() {}, function() {}) [1] FALSE Yet, after loading the Matrix package (which redefines det), the following is checked (in library.checkConflicts): > identical(get("det", baseenv()), get("det", asNamespace("Matrix")), ignore.environment=T) [1] TRUE I've looked at the code in identical.c and for closures it seems to compare the FORMALS and the BOD...

...a square matrix of dimension k by k with alpha on the diagonal and beta everywhee else. This symmetric matrix is called symmetric compound matrix and has the form a( I + cJ), where I is the k by k identity matrix J is the k by k matrix of all ones a = alpha - beta c = beta/a I need to evaluate the determinant of this matrix. Is there any algebric formula for that? thank you for your help Stefano [[alternative HTML version deleted]]