similar to: Determinant function (PR#9715)

Dear R users, even if this question is not related to an issue about R, probably some of you will be able to help me. I have 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 =

r-devel list, 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

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

Dear R list members, I have a vector of Cholesky parameterization of a matrix let say A. I would like to compute the determinant and inverse of the original matrix A from the vector of cholesky parameters , would you suggest an R function to do the task. I have tried hard but unable to find anything like that. Please direct me any package or please share R code that can do the task. Thanks and

Dear R-users, I computed determinant of a square matrix "var.r" using the SVD output: detr _ 1 d _ svd(var.r)$d for (i in 1:length(d)) { detr _ detr*d[i] } print(detr) 30.20886 BUT when I tried : det(var.r) I got : -30.20886 Is this because SVD output will only give absolute of the eigenvalues ?, If this is the case how can I get the original eigenvalues? Thanks, Agus

Dear all, I just noticed that the 0.9 update for FrF2 did not work out for Mac OS due to an error in an example that ran without error on all other platforms. I do not find any reason for this. In the past, umlauts or tab characters have sometimes been an issue, but I didn't find any of these. The function definition is FrF2(nruns = NULL, nfactors = NULL, factor.names = if

Full_Name: Mike Meredith Version: 2.3.0 OS: WinXP SP2 Submission from: (NULL) (210.195.228.29) Using dt() with a non-centrality parameter and near-zero values for 'x' results in erratic output. Try this: tst <- c(1e-12, 1e-13, 1e-14, 1e-15, 1e-16, 1e-17, 0) dt(tst,16,1) I get: 0.2381019 0.2385462 0.2296557 0.1851817 0.6288373 3.8163916 (!!) 0.2382217 The 0.238 values are okay,

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

Hi, I have data with multiple sub-second entries: 2011/03/15 09:32:15.035619,-0.403103,1.09664,48.6,126.92,117.32 2011/03/15 09:32:15.069331,-0.39851,1.09874,48.6,126.92,117.32 2011/03/15 09:32:15.289135,-0.402463,1.10084,48.59,126.92,117.32 2011/03/15 09:32:15.296110,-0.450244,1.10063,48.59,126.92,117.32 2011/03/15 09:32:15.451358,-0.438813,1.10273,48.59,126.93,117.32 2011/03/15

Dear R Users, I am having trouble with lmer. I am looking at recombinant versus non recombinant individuals. In the response variable recombinant individuals are coded as 1's and non-recombinant as 0's. I built a model with 2 fixed factors and 1 random effect. Sex (males/females) is the first fixed effect and sexual genotype (XY, YY, WX and WY) the second one. Sexual Genotype is

Dear R users, Simple question. Can anyone help with the code that would allow me to view only the variables who's correlation output is >0.8? This is the code I'm using to date >cor(data, method="spearman") Kind regards Krys -------------------------------------------- _________________________________________________________________ Save time by

Hi R-List, I have a question regarding R-language formats, I think. I am producing a series of graphs (using plot, barplot, barchart, and bwplot, using either text or mtext to place values on the graphs) and tables for a Francophone country. In fact, I have already done so. However, while they are pleased with the results they've requested I convert all of my decimal points into the French

Is there a function in R to calculate the generalized determinant of a singular matrix? - similar to the ginv() used to compute the generalized inverse. I can't seem to find any R related posts at all. Thanks in advance, Sean O'Riordain Trinity College Dublin -- View this message in context: http://r.789695.n4.nabble.com/Moore-Penrose-Generalized-determinant-tp4471629p4471629.html Sent

Hi, What Lapack routine(s) should I use to calculate the determinant of a symmetric matrix of real numbers? I would appreciate any guidance on this issue. Thanks, Kosuke --------------------------------------------------------- Kosuke Imai Office: Corwin Hall 041 Assistant Professor Phone: 609-258-6601 (Direct) Department of Politics Fax: 609-258-1110 (Department)

On 19 January 2017 at 02:49, Klint Gore wrote: | -----Original Message----- | From: Dirk Eddelbuettel [mailto:dirk.eddelbuettel at gmail.com] On Behalf Of Dirk Eddelbuettel | > Sorry, what part of '14.04' is current? | > | > Ubuntu is at 16.10. And release 16.04, which as a LTS replaces the LTS 14.04 you use, also passes. Can you upgrade? | | That version of the openblas package

Hi, What determines the umask of sieve_pipe_bin_dir scripts ? The results from my script are always being set to 0600. My script is simple and shown below, even if I adjust the right line to add " && chmod 644", the actual resulting file still remains at 0600 ?!? #!/bin/bash # Usage: imapsieve_copy <email> <spam|ham> MSG_USER="$1" MSG_TYPE="$2"

> On 27/05/2020 13:28 Laura Smith <n5d9xq3ti233xiyif2vp at protonmail.ch> wrote: > > > Hi, > What determines the umask of sieve_pipe_bin_dir scripts ? > The results from my script are always being set to 0600. > My script is simple and shown below, even if I adjust the right line to add " && chmod 644", the actual resulting file still remains at 0600

-----Original Message----- From: Dirk Eddelbuettel [mailto:dirk.eddelbuettel at gmail.com] On Behalf Of Dirk Eddelbuettel Sent: Thursday, 19 January 2017 11:21 AM To: Klint Gore Cc: r-sig-debian at r-project.org Subject: Re: [R-sig-Debian] Taking determinant of a matrix of NAs results in intermittent memory corruption >So this converges towards 'old versions bad, new versions fine' ?

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

what do you mean by "more optimal instructions" ? -omer On Wed, Nov 16, 2011 at 1:28 AM, Joe Abbey <jabbey at arxan.com> wrote: > I've done a first pass over the past 6 months of changes and some notable > things stood out: > > * The ARM backend has reworked Set Jump Long Jump EH Lowering. > * The ARM backend includes improved support for Cortex-M > *