similar to: Inverting a square matrix using solve() with LAPACK=TRUE (PR#13765)

Full_Name: Ravi Varadhan Version: 2.8.1 OS: Windows Submission from: (NULL) (162.129.251.19) Inverting a matrix with solve(), but using LAPACK=TRUE, gives erroneous results: Here is an example: hilbert <- function(n) { i <- 1:n; 1 / outer(i - 1, i, "+") } h5 <- hilbert(5) hinv1 <- solve(qr(h5)) hinv2 <- solve(qr(h5, LAPACK=TRUE)) all.equal(hinv1, hinv2) #

Refiling this. The actual fix was slightly more complicated. Will soon be committed to R-Patched (aka 2.9.1 beta). -p rvaradhan at jhmi.edu wrote: > Full_Name: Ravi Varadhan > Version: 2.8.1 > OS: Windows > Submission from: (NULL) (162.129.251.19) >=20 >=20 > Inverting a matrix with solve(), but using LAPACK=3DTRUE, gives erroneo= us > results: Thanks, but there seems

The bug I submitted yesterday (It's not entered in the bug data base, so I have no ID for it) included a suggested fix that is not correct. It worked for the examples I gave because there was no pivoting in fact, or only pivot permutations that were idempotent. A correction that works in general on the examples I gave makes these two changes in qr.coef(): ## coef[qr$pivot, ]

Hello. I am trying to invert a matrix, and I am finding that I can get different answers depending on whether I set LAPACK true or false using "qr". I had understood that LAPACK is, in general more robust and faster than LINPACK, so I am confused as to why I am getting what seems to be invalid answers. The matrix is ostensibly the Hessian for a function I am optimizing. I want to get

Dear all, I am trying to get an estimate of the intercept for a linear model. In this case, I know the slope of the model, can anyone tell me how to constrain the formula in lm() so that it only estimates the intercept not the slope? Many thanks in advance, Sincerely, Liqing Zhang Dept. of Eco. Evol. Biol. Univ. of CA, Irvine email: lzhang at uci.edu >From VM Mon Apr 30 08:18:45 2001

Hi everyone, I am looking to use R as a MATLAB replacement for linear algebra. I've done a fairly good job for finding replacements for most of the functions I'm interested in, I John Fox wrote a program for implementing the reduced row echelon form of a matrix (by doing the Gauss-Jordan elimination). I modified it a bit: rref <- function(A,

Hi everyone, I've noticed a problem with Debian wheezy + samba 3.6.6 configured with acl_xattr in my configuration. The following test sequence causes Windows Explorer to report incorrectly ordered permission entries: 1) Map a share as with "admin" user credentials to a drive letter on a Windows client 2) Create a folder at the root of the share "rootfolder" 3) Create a

Hi, I have successfully setup a tinc network between five hosts (in switch mode). Two of the hosts have static and known IP addresses (S1 and S2). Other hosts (H3-H5) connect one (or both) of them. The traffic flows nicely between all hosts. The initial edges (ConnectTo configuration directives) in my test network are: S1<->S2 H3 -> S1 and S2 H4 -> S1 H5 -> S2 As far as I have

Quantitative Risk Analysis with R May 1-4, 2018 Fort Collins, Colorado, USA Join us this spring for our QRA with R training. Our class will focus on applied risk modeling methods using the R statistical language and will cover the core principles of QRA and Monte Carlo simulation modeling. Both Bayesian and frequentist methods will be discussed. This class is very popular with a variety of

>>>>> Suharto Anggono Suharto Anggono via R-devel <r-devel at r-project.org> >>>>> on Fri, 11 Aug 2017 17:11:06 +0000 writes: >>>>> Suharto Anggono Suharto Anggono via R-devel <r-devel at r-project.org> >>>>> on Fri, 11 Aug 2017 17:11:06 +0000 writes: > See

On 26/12/2017 9:40 AM, Dirk Eddelbuettel wrote: > > On 26 December 2017 at 22:14, Sun Yijiang wrote: > | Thanks for the solution. Now I know the work-arounds, but still don't > | quite get it. Why does R_DEFAULT_PACKAGES has anything to do with > | library(methods)? > > Because it governs which packages are loaded by default. And while R also > loads

Examples similar to pretty(c(-1,1)*1e300, n = 1e9, min.n = 1) with smaller 'n': pretty(c(-1,1)*1e304, n = 1e5, min.n = 1) pretty(c(-1,1)*1e306, n = 1e3, min.n = 1) A report on 'pretty' when working with integers, similar to what led to change of 'seq' fuzz, is https://bugs.r-project.org/bugzilla3/show_bug.cgi?id=15137 -------------------------------------------- On Tue,

Hi Dirk, Thanks for the solution. Now I know the work-arounds, but still don't quite get it. Why does R_DEFAULT_PACKAGES has anything to do with library(methods)? If library(h5) works, it should just work, not depend on an environment variable. Rscript is not consistent with R, that's my confusion. Steve 2017-12-26 20:46 GMT+08:00 Dirk Eddelbuettel <edd at debian.org>: >

See https://stat.ethz.ch/pipermail/r-devel/2017-August/074746.html for the origin of the example here. That pretty(c(-1,1)*1e300, n = 1e9, min.n = 1) gave 20 intervals, far from 1e9, but pretty(c(-1,1)*1e300, n = 1e6, min.n = 1) gave 1000000 intervals (on a machine), made me trace through the code to function 'R_pretty' in https://svn.r-project.org/R/trunk/src/appl/pretty.c . *lo is

Consider this script (with h5 installed): $ cat test.R library(h5) name <- tempfile() f <- h5file(name) file.remove(name) $ Rscript test.R Error in initialize(value, ...) : cannot use object of class "character" in new(): class "H5File" does not extend that class Calls: h5file -> new -> initialize -> initialize Execution halted $ /usr/lib64/R/bin/R --slave

On 26 December 2017 at 15:24, Sun Yijiang wrote: | After looking into C source code, I found that Rscript by default fills | environment variable R_DEFAULT_PACKAGES with | "datasets,utils,grDevices,graphics,stats", and it somehow fails some | package like h5. | | The problem here is, not setting R_DEFAULT_PACKAGES is equivalent to | setting it to a magic value, it's really

Hello All, I just signed up to this list, so I apologize if this question has been asked before. I am trying to load an hdf5 file into R and running into some problems. Here are the steps I took to configure my environment: * R 2.10.0 (x64) on Mac OS X 10.6 * hdf5 1.8.3 installed via macports * hdf5_1.6.9.tar.gz from CRAN I suspect the problem I am having relates to

*Hi everyone, I am a new subscriber, so I hope I'm posting in the right way. I need to load a hdf5 file, so I installed the hdf5 libraries ad the hdf5 R package; in downloaded a small sample file from the hdf5 website and I could correctly load it. But when I try to load my file: * require(hdf5) rm(list=ls()) ls() hdf5load("my_file.h5", verbosity=3) *I get this error message: *

On 8/8/20 12:21 AM, Rowland penny via samba wrote: > On 07/08/2020 22:44, Simon Matthews via samba wrote: >> >>>> >>>>> >>>> >>>> This is where your problems start, you do not have enough lines, I >>>> would expect something like this: >>>> >>>> idmap config * : backend = tdb >>>>

On 08/08/2020 19:43, Simon Matthews wrote: > On 8/8/20 12:21 AM, Rowland penny via samba wrote: >> On 07/08/2020 22:44, Simon Matthews via samba wrote: >>> >>>>> >>>>>> >>>>> >>>>> This is where your problems start, you do not have enough lines, I >>>>> would expect something like this: