Displaying 3 results from an estimated 3 matches for "det1".
Did you mean:
dat1
1999 Apr 20
1
eigenvalue calculations
...at the CRAN sites) that bug has been fixed.
In case 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 in...
1999 Apr 20
2
eigenvalue/eigenvector calculations
...'t think it is just a matter
of the ordering of the results. Their 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...
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(...) ?
>
>
>