Displaying 3 results from an estimated 3 matches for "det2".
Did you mean:
dat2
1999 Apr 20
1
eigenvalue calculations
...anyone else is interested, I redid the determinant
calculations in
Version 0.64.0 Patched (unreleased snapshot) (April 19, 1999)
using the method from Stephan Steinhaus's script (det0), the method
based on calculating the eigenvalues alone (det1) and the method based
on the QR decomposition (det2). The results from the three methods
are consistent except that the eigenvalue methods both produce a
non-trivial imaginary component.
I enclose the timing results below. Notice that the determinant
becomes extremely large as the matrix size increases and that the
error introduced by the complex...
1999 Apr 20
2
eigenvalue/eigenvector calculations
...r products are noticeably
different.
R : Copyright 1999, The R Development Core Team
Version 0.64.0 (April 8, 1999)
> x <- array( runif( 10000 ), c(100, 100))
> determinant <- function(x) prod(eigen(x)$values)
> det1 <- function(x) prod(eigen(x, only.values = TRUE)$values)
> det2 <- function(x) prod(diag(qr(x)$qr)) * (-1)^(ncol(x) - 1)
> system.time(print(determinant(x)))
[1] -27320631-8.050401e-10i
[1] 0.06 0.04 0.00 0.00 0.00
> system.time(print(det1(x)))
[1] 36072865748-5.155196e-06i
[1] 0.06 0.00 0.00 0.00 0.00
> system.time(print(det2(x)))
[1] -9.854994e+25...
2023 Feb 01
1
Detpack package
I did use "detpack" ie not with a capital
detpack::chi2testuniform(vals,0.05)
gives this:
Error: 'chi2testuniform' is not an exported object from 'namespace:detpack'
??
Thanks Nick
On Wed, 1 Feb 2023 at 16:29, Eric Berger <ericjberger at gmail.com> wrote:
> Detpack or detpack?
>
> What happens when you try detpack::chi2testuniform(...) ?
>
>
>