search for: qra

Displaying 12 results from an estimated 12 matches for "qra".

Did you mean: qa
2009 Jun 17
1
Inverting a square matrix using solve() with LAPACK=TRUE (PR#13762)
...# They are not equal Here is a function that I wrote to correct this problem: solve.lapack <- function(A, LAPACK=TRUE, tol=1.e-07) { # A function to invert a matrix using "LAPACK" or "LINPACK" if (nrow(A) != ncol(A)) stop("Matrix muxt be square") qrA <- qr(A, LAPACK=LAPACK, tol=tol) if (LAPACK) { apply(diag(1, ncol(A)), 2, function(x) solve(qrA, x)) } else solve(qrA) } hinv3 <- solve.lapack(h5) all.equal(hinv1, hinv3) # Now, they are equal
2017 Nov 06
0
QRA with R Course 12-4 to 12-7-17 - Ft. Collins, CO
Quantitative Risk Analysis with R Fort Collins, Colorado, USA December 4-7, 2017 This course will cover the core principles of quantitative risk analysis and Monte Carlo simulation modeling. The focus of the course is on how to conduct accurate and effective quantitative risk analyses, including best practices of risk modeling, selecting the appropriate distribution, using data and expert
2016 Dec 20
0
EpiX Analytics-QRA Consultant Position
Position: Quantitative Risk Analysis Consultant APPLICATIONS RECEIVED UNTIL JANUARY 9TH, 2017 Location: Fort Collins, CO (physical presence required) Position type: Full time We are a Colorado-based specialist risk and decision modeling consulting and research firm, with diverse experience and skills. We are seeking to fill a full-time position as a Quantitative Risk Analysis Consultant, to
2009 Jun 18
1
Inverting a square... (PR#13762)
...unction that I wrote to correct this problem: >=20 > solve.lapack <- function(A, LAPACK=3DTRUE, tol=3D1.e-07) { > # A function to invert a matrix using "LAPACK" or "LINPACK" > if (nrow(A) !=3D ncol(A)) stop("Matrix muxt be square") > qrA <- qr(A, LAPACK=3DLAPACK, tol=3Dtol) > if (LAPACK) { > apply(diag(1, ncol(A)), 2, function(x) solve(qrA, x))=20 > } else solve(qrA) > } >=20 > hinv3 <- solve.lapack(h5) > all.equal(hinv1, hinv3) # Now, they are equal >=20 > _____________...
2009 Jun 18
0
Inverting a square matrix using solve() with LAPACK=TRUE (PR#13765)
...e is a function that I wrote to correct this problem: > > solve.lapack <- function(A, LAPACK=TRUE, tol=1.e-07) { > # A function to invert a matrix using "LAPACK" or "LINPACK" > if (nrow(A) != ncol(A)) stop("Matrix muxt be square") > qrA <- qr(A, LAPACK=LAPACK, tol=tol) > if (LAPACK) { > apply(diag(1, ncol(A)), 2, function(x) solve(qrA, x)) > } else solve(qrA) > } > > hinv3 <- solve.lapack(h5) > all.equal(hinv1, hinv3) # Now, they are equal > > _______________________...
2007 Sep 01
1
row echelon form
...getting the lower triangular part to be zero, along with the first non-zero entry in each non-zero row being a one. Here is what I wrote: ref <- function(A) { ## Implements gaussian elimination using the QR factorization. ## Note: ref form is not unique. ## written by S. K. Hyde ## qrA <- qr(A) r <- qrA$rank ##computed rank of the matrix R <- qr.R(qrA) if (r < nrow(R)) R[-(1:r),] <- 0 #zero out rows that should be zero (according to the rank). #Get ones as the first nonzero entry in each row and return result. diag(c(1/diag(R)[1:r],rep(1,nrow(R)-r)...
2010 Dec 06
1
use pcls to solve least square fitting with constraints
...ely, it is stuck with the following error: > M <- list() > M$y = Dmat[,1] > M$X = Cmat > M$Ain = as.matrix(Amat) > M$bin = rep(0, dim(Amat)[1]) > M$p=qr.solve(as.matrix(Cmat), Dmat[,1]) > M$w = rep(1, length(M$y)) > M$C = matrix(0,0,0) > p<-pcls(M) Error in t(qr.qty(qra, t(M$X))[(j + 1):k, ]) : error in evaluating the argument 'x' in selecting a method for function 't' After some searches, I still couldn't find any solution, any help and/or advice will be highly appreciated! Baoqiang
2018 Jan 30
0
Quantitative Risk Analysis with R Course 5/1/18 to 5/4/18
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 participants,...
2019 Nov 05
3
newbie question on a dovecot buffer
...used Dovecot before so its a vertical learning curve Thanks -- Best wishes /73 Richard Bown Email : richard at g8jvm.com HTTP : http://www.g8jvm.com ###################################################################### Ham Call: G8JVM . QRV: 50-432 MHz + Microwave Maidenhead QRA: IO82SP38, LAT. 52 39.720' N LONG. 2 28.171 W QRV 6mtrs 200W, 4mtrs 150W, 2mtrs 300W, 70cms 200W, Microwave 1296MHz 110W, 2320MHz 100W, 5760MHz 10W & 10368MHz 5W OS: Linux Mint 19.2 x86_64 Tina, on a Dell Inspiron 3580 laptop #########################################################...
2002 Jun 12
3
How does R compares for speed?
Hi, For those who are interested, I have update my R benchmark (to version 1.5.0) and also to Splus 6. They are available at: http://www.sciviews.org/other/benchmark.htm. A comparison is made between Matlab (5.3 & 6), R 1.5.0, Splus 6 rel 2, O-Matrix 5.1, Octave 2.1.31, Scilab 2.6, Rlab 2.1 and OX 3.00 under Windows 2000 pro. Overall, R is not the fastest package, but it is one of the
2001 Feb 16
1
Sub_scribe and a question
...})[3] cumulate <- cumulate + timing } timing <- cumulate/runs times[4, 1] <- timing cat(c("550x550 cross-product matrix (b =3D a' * a)___________ (sec): ", = timing, "\n")) remove("a", "b") flush.console() # (5) cumulate <- 0; c <- 0; qra <-0 for (i in 1:runs) { a <- rnorm(700*700); dim(a) <- c(700,700) b <- 1:700 timing <- system.time({ qra <- qr(a, tol =3D 1e-10); c <- qr.coef(qra, b) #Rem: a little faster than c <- lsfit(a, b, inter=3DF)$coefficients })[3] cumulate <- cumulate + ti...
2008 Jun 30
4
Rebuild of kernel 2.6.9-67.0.20.EL failure
Hello list. I'm trying to rebuild the 2.6.9.67.0.20.EL kernel, but it fails even without modifications. How did I try it? Created a (non-root) build environment (not a mock ) Installed the kernel.scr.rpm and did a rpmbuild -ba --target=`uname -m` kernel-2.6.spec 2> prep-err.log | tee prep-out.log The build failed at the end: Processing files: kernel-xenU-devel-2.6.9-67.0.20.EL Checking