Last updated on 2021-04-20 23:47:40 CEST.
| Flavor | Version | Tinstall | Tcheck | Ttotal | Status | Flags |
|---|---|---|---|---|---|---|
| r-devel-linux-x86_64-debian-clang | 0.3-0 | 4.33 | 167.94 | 172.27 | OK | |
| r-devel-linux-x86_64-debian-gcc | 0.3-0 | 3.97 | 127.19 | 131.16 | OK | |
| r-devel-linux-x86_64-fedora-clang | 0.3-0 | 209.19 | NOTE | |||
| r-devel-linux-x86_64-fedora-gcc | 0.3-0 | 200.11 | NOTE | |||
| r-devel-windows-ix86+x86_64 | 0.3-0 | 16.00 | 168.00 | 184.00 | OK | |
| r-devel-windows-x86_64-gcc10-UCRT | 0.3-0 | NOTE | ||||
| r-patched-linux-x86_64 | 0.3-0 | 6.45 | 159.18 | 165.63 | WARN | |
| r-patched-solaris-x86 | 0.3-0 | 229.50 | NOTE | |||
| r-release-linux-x86_64 | 0.3-0 | 4.33 | 159.51 | 163.84 | OK | |
| r-release-macos-x86_64 | 0.3-0 | NOTE | ||||
| r-release-windows-ix86+x86_64 | 0.3-0 | 11.00 | 217.00 | 228.00 | OK | |
| r-oldrel-macos-x86_64 | 0.3-0 | ERROR | ||||
| r-oldrel-windows-ix86+x86_64 | 0.3-0 | 8.00 | 189.00 | 197.00 | ERROR |
Version: 0.3-0
Check: dependencies in R code
Result: NOTE
Namespace in Imports field not imported from: ‘graphics’
All declared Imports should be used.
Flavors: r-devel-linux-x86_64-fedora-clang, r-devel-linux-x86_64-fedora-gcc, r-devel-windows-x86_64-gcc10-UCRT, r-patched-solaris-x86, r-release-macos-x86_64
Version: 0.3-0
Check: whether package can be installed
Result: WARN
Found the following significant warnings:
Rd warning: /tmp/RtmpQxImo6/R.INSTALL1b6d64b658e1/DPQmpfr/man/hyperGeomQ.Rd:9: file link ‘chooseZ’ in package ‘gmp’ does not exist and so has been treated as a topic
Rd warning: /tmp/RtmpQxImo6/R.INSTALL1b6d64b658e1/DPQmpfr/man/hyperGeomQ.Rd:35: file link ‘chooseZ’ in package ‘gmp’ does not exist and so has been treated as a topic
Flavor: r-patched-linux-x86_64
Version: 0.3-0
Check: whether package can be installed
Result: ERROR
Installation failed.
Flavor: r-oldrel-macos-x86_64
Version: 0.3-0
Check: tests
Result: ERROR
Running 'beta-Ding94.R' [107s]
Running 'dpqnorm-tst.R' [3s]
Running 'hypergeom.R' [3s]
Running 't-nonc_tst.R' [6s]
Running the tests in 'tests/dpqnorm-tst.R' failed.
Complete output:
> require(DPQmpfr)
Loading required package: DPQmpfr
> require(Rmpfr)
Loading required package: Rmpfr
Loading required package: gmp
Attaching package: 'gmp'
The following objects are masked from 'package:base':
%*%, apply, crossprod, matrix, tcrossprod
C code of R package 'Rmpfr': GMP using 64 bits per limb
Attaching package: 'Rmpfr'
The following object is masked from 'package:gmp':
outer
The following objects are masked from 'package:stats':
dbinom, dgamma, dnbinom, dnorm, dpois, pnorm
The following objects are masked from 'package:base':
cbind, pmax, pmin, rbind
>
> ## very accurate set of *large* abscissa values
> x <- mpfr(2,1024)^ seq(3, 100, by=1/2)
> xN <- asNumeric(x)
>
> ## Looking at pnormU_S53() = upper bound to pnorm(), for large x,
> ## the term log(4 / (3 + sqrt(1 + 8/x^2))) early underflows to 0:
> eps_S53 <- function(x, method = c("direct", "log1pExct", "deltaFull",
+ "delta1", "delta2", "delta3")) {
+ method <- match.arg(method)
+ delta <- function(x) { xm2 <- x^-2; (.5*(1 - sqrt(1 + 8*xm2)) + xm2) / (1 + xm2) }
+ switch(method,
+ "direct" = log(4 / (3 + sqrt(1 + 8*x^-2))),
+ "log1pExct" = log1p(delta(x)),
+ "deltaFull" = delta(x),
+ "delta1" = -1/(x^2+1) , ## 1 term for delta(x)
+ "delta2" = -(1 - (2/x)^2) / (x^2+1) , ## 2 terms for delta(x)
+ "delta3" = {tx2 <- (2/x)^2 ; -(1 - tx2*(1 - tx2)) / (x^2+1) }, ## 3 terms
+ stop("invalid 'method': ", method))
+ }
>
>
> (S53meth <- eval(formals(eps_S53)$method))
[1] "direct" "log1pExct" "deltaFull" "delta1" "delta2" "delta3"
> ## "direct" "log1pExct" "deltaFull" "delta1" "delta2" "delta3"
>
> epsAllS53 <- function(x, meth = eval(formals(eps_S53)$method)) {
+ stopifnot(is.character(meth), is.numeric(x) || is(x,"mpfr"))
+ Rl <- lapply(meth, function(m) eps_S53(x, method=m))
+ new("mpfrMatrix", unlist(Rl),
+ Dim = c(length(x), length(meth)),
+ Dimnames = list(NULL, meth))
+ }
>
> EA <- epsAllS53(x)
> asNumeric(EA) # looks +- ok
direct log1pExct deltaFull delta1 delta2
[1,] -1.505120e-02 -1.458470e-02 -1.447885e-02 -1.538462e-02 -1.442308e-02
[2,] -7.664653e-03 -7.545372e-03 -7.516977e-03 -7.751938e-03 -7.509690e-03
[3,] -3.868707e-03 -3.838543e-03 -3.831185e-03 -3.891051e-03 -3.830253e-03
[4,] -1.943664e-03 -1.936080e-03 -1.934207e-03 -1.949318e-03 -1.934089e-03
[5,] -9.741879e-04 -9.722861e-04 -9.718136e-04 -9.756098e-04 -9.717988e-04
[6,] -4.876864e-04 -4.872103e-04 -4.870916e-04 -4.880429e-04 -4.870897e-04
[7,] -2.439918e-04 -2.438726e-04 -2.438429e-04 -2.440810e-04 -2.438427e-04
[8,] -1.220331e-04 -1.220033e-04 -1.219958e-04 -1.220554e-04 -1.219958e-04
[9,] -6.102585e-05 -6.101840e-05 -6.101653e-05 -6.103143e-05 -6.101653e-05
[10,] -3.051525e-05 -3.051339e-05 -3.051292e-05 -3.051665e-05 -3.051292e-05
[11,] -1.525821e-05 -1.525774e-05 -1.525762e-05 -1.525856e-05 -1.525762e-05
[12,] -7.629249e-06 -7.629133e-06 -7.629104e-06 -7.629336e-06 -7.629103e-06
[13,] -3.814661e-06 -3.814632e-06 -3.814625e-06 -3.814683e-06 -3.814625e-06
[14,] -1.907340e-06 -1.907332e-06 -1.907330e-06 -1.907345e-06 -1.907330e-06
[15,] -9.536720e-07 -9.536702e-07 -9.536698e-07 -9.536734e-07 -9.536698e-07
[16,] -4.768366e-07 -4.768361e-07 -4.768360e-07 -4.768369e-07 -4.768360e-07
[17,] -2.384184e-07 -2.384183e-07 -2.384183e-07 -2.384185e-07 -2.384183e-07
[18,] -1.192093e-07 -1.192092e-07 -1.192092e-07 -1.192093e-07 -1.192092e-07
[19,] -5.960464e-08 -5.960463e-08 -5.960463e-08 -5.960464e-08 -5.960463e-08
[20,] -2.980232e-08 -2.980232e-08 -2.980232e-08 -2.980232e-08 -2.980232e-08
[21,] -1.490116e-08 -1.490116e-08 -1.490116e-08 -1.490116e-08 -1.490116e-08
[22,] -7.450580e-09 -7.450580e-09 -7.450580e-09 -7.450581e-09 -7.450580e-09
[23,] -3.725290e-09 -3.725290e-09 -3.725290e-09 -3.725290e-09 -3.725290e-09
[24,] -1.862645e-09 -1.862645e-09 -1.862645e-09 -1.862645e-09 -1.862645e-09
[25,] -9.313226e-10 -9.313226e-10 -9.313226e-10 -9.313226e-10 -9.313226e-10
[26,] -4.656613e-10 -4.656613e-10 -4.656613e-10 -4.656613e-10 -4.656613e-10
[27,] -2.328306e-10 -2.328306e-10 -2.328306e-10 -2.328306e-10 -2.328306e-10
[28,] -1.164153e-10 -1.164153e-10 -1.164153e-10 -1.164153e-10 -1.164153e-10
[29,] -5.820766e-11 -5.820766e-11 -5.820766e-11 -5.820766e-11 -5.820766e-11
[30,] -2.910383e-11 -2.910383e-11 -2.910383e-11 -2.910383e-11 -2.910383e-11
[31,] -1.455192e-11 -1.455192e-11 -1.455192e-11 -1.455192e-11 -1.455192e-11
[32,] -7.275958e-12 -7.275958e-12 -7.275958e-12 -7.275958e-12 -7.275958e-12
[33,] -3.637979e-12 -3.637979e-12 -3.637979e-12 -3.637979e-12 -3.637979e-12
[34,] -1.818989e-12 -1.818989e-12 -1.818989e-12 -1.818989e-12 -1.818989e-12
[35,] -9.094947e-13 -9.094947e-13 -9.094947e-13 -9.094947e-13 -9.094947e-13
[36,] -4.547474e-13 -4.547474e-13 -4.547474e-13 -4.547474e-13 -4.547474e-13
[37,] -2.273737e-13 -2.273737e-13 -2.273737e-13 -2.273737e-13 -2.273737e-13
[38,] -1.136868e-13 -1.136868e-13 -1.136868e-13 -1.136868e-13 -1.136868e-13
[39,] -5.684342e-14 -5.684342e-14 -5.684342e-14 -5.684342e-14 -5.684342e-14
[40,] -2.842171e-14 -2.842171e-14 -2.842171e-14 -2.842171e-14 -2.842171e-14
[41,] -1.421085e-14 -1.421085e-14 -1.421085e-14 -1.421085e-14 -1.421085e-14
[42,] -7.105427e-15 -7.105427e-15 -7.105427e-15 -7.105427e-15 -7.105427e-15
[43,] -3.552714e-15 -3.552714e-15 -3.552714e-15 -3.552714e-15 -3.552714e-15
[44,] -1.776357e-15 -1.776357e-15 -1.776357e-15 -1.776357e-15 -1.776357e-15
[45,] -8.881784e-16 -8.881784e-16 -8.881784e-16 -8.881784e-16 -8.881784e-16
[46,] -4.440892e-16 -4.440892e-16 -4.440892e-16 -4.440892e-16 -4.440892e-16
[47,] -2.220446e-16 -2.220446e-16 -2.220446e-16 -2.220446e-16 -2.220446e-16
[48,] -1.110223e-16 -1.110223e-16 -1.110223e-16 -1.110223e-16 -1.110223e-16
[49,] -5.551115e-17 -5.551115e-17 -5.551115e-17 -5.551115e-17 -5.551115e-17
[50,] -2.775558e-17 -2.775558e-17 -2.775558e-17 -2.775558e-17 -2.775558e-17
[51,] -1.387779e-17 -1.387779e-17 -1.387779e-17 -1.387779e-17 -1.387779e-17
[52,] -6.938894e-18 -6.938894e-18 -6.938894e-18 -6.938894e-18 -6.938894e-18
[53,] -3.469447e-18 -3.469447e-18 -3.469447e-18 -3.469447e-18 -3.469447e-18
[54,] -1.734723e-18 -1.734723e-18 -1.734723e-18 -1.734723e-18 -1.734723e-18
[55,] -8.673617e-19 -8.673617e-19 -8.673617e-19 -8.673617e-19 -8.673617e-19
[56,] -4.336809e-19 -4.336809e-19 -4.336809e-19 -4.336809e-19 -4.336809e-19
[57,] -2.168404e-19 -2.168404e-19 -2.168404e-19 -2.168404e-19 -2.168404e-19
[58,] -1.084202e-19 -1.084202e-19 -1.084202e-19 -1.084202e-19 -1.084202e-19
[59,] -5.421011e-20 -5.421011e-20 -5.421011e-20 -5.421011e-20 -5.421011e-20
[60,] -2.710505e-20 -2.710505e-20 -2.710505e-20 -2.710505e-20 -2.710505e-20
[61,] -1.355253e-20 -1.355253e-20 -1.355253e-20 -1.355253e-20 -1.355253e-20
[62,] -6.776264e-21 -6.776264e-21 -6.776264e-21 -6.776264e-21 -6.776264e-21
[63,] -3.388132e-21 -3.388132e-21 -3.388132e-21 -3.388132e-21 -3.388132e-21
[64,] -1.694066e-21 -1.694066e-21 -1.694066e-21 -1.694066e-21 -1.694066e-21
[65,] -8.470329e-22 -8.470329e-22 -8.470329e-22 -8.470329e-22 -8.470329e-22
[66,] -4.235165e-22 -4.235165e-22 -4.235165e-22 -4.235165e-22 -4.235165e-22
[67,] -2.117582e-22 -2.117582e-22 -2.117582e-22 -2.117582e-22 -2.117582e-22
[68,] -1.058791e-22 -1.058791e-22 -1.058791e-22 -1.058791e-22 -1.058791e-22
[69,] -5.293956e-23 -5.293956e-23 -5.293956e-23 -5.293956e-23 -5.293956e-23
[70,] -2.646978e-23 -2.646978e-23 -2.646978e-23 -2.646978e-23 -2.646978e-23
[71,] -1.323489e-23 -1.323489e-23 -1.323489e-23 -1.323489e-23 -1.323489e-23
[72,] -6.617445e-24 -6.617445e-24 -6.617445e-24 -6.617445e-24 -6.617445e-24
[73,] -3.308722e-24 -3.308722e-24 -3.308722e-24 -3.308722e-24 -3.308722e-24
[74,] -1.654361e-24 -1.654361e-24 -1.654361e-24 -1.654361e-24 -1.654361e-24
[75,] -8.271806e-25 -8.271806e-25 -8.271806e-25 -8.271806e-25 -8.271806e-25
[76,] -4.135903e-25 -4.135903e-25 -4.135903e-25 -4.135903e-25 -4.135903e-25
[77,] -2.067952e-25 -2.067952e-25 -2.067952e-25 -2.067952e-25 -2.067952e-25
[78,] -1.033976e-25 -1.033976e-25 -1.033976e-25 -1.033976e-25 -1.033976e-25
[79,] -5.169879e-26 -5.169879e-26 -5.169879e-26 -5.169879e-26 -5.169879e-26
[80,] -2.584939e-26 -2.584939e-26 -2.584939e-26 -2.584939e-26 -2.584939e-26
[81,] -1.292470e-26 -1.292470e-26 -1.292470e-26 -1.292470e-26 -1.292470e-26
[82,] -6.462349e-27 -6.462349e-27 -6.462349e-27 -6.462349e-27 -6.462349e-27
[83,] -3.231174e-27 -3.231174e-27 -3.231174e-27 -3.231174e-27 -3.231174e-27
[84,] -1.615587e-27 -1.615587e-27 -1.615587e-27 -1.615587e-27 -1.615587e-27
[85,] -8.077936e-28 -8.077936e-28 -8.077936e-28 -8.077936e-28 -8.077936e-28
[86,] -4.038968e-28 -4.038968e-28 -4.038968e-28 -4.038968e-28 -4.038968e-28
[87,] -2.019484e-28 -2.019484e-28 -2.019484e-28 -2.019484e-28 -2.019484e-28
[88,] -1.009742e-28 -1.009742e-28 -1.009742e-28 -1.009742e-28 -1.009742e-28
[89,] -5.048710e-29 -5.048710e-29 -5.048710e-29 -5.048710e-29 -5.048710e-29
[90,] -2.524355e-29 -2.524355e-29 -2.524355e-29 -2.524355e-29 -2.524355e-29
[91,] -1.262177e-29 -1.262177e-29 -1.262177e-29 -1.262177e-29 -1.262177e-29
[92,] -6.310887e-30 -6.310887e-30 -6.310887e-30 -6.310887e-30 -6.310887e-30
[93,] -3.155444e-30 -3.155444e-30 -3.155444e-30 -3.155444e-30 -3.155444e-30
[94,] -1.577722e-30 -1.577722e-30 -1.577722e-30 -1.577722e-30 -1.577722e-30
[95,] -7.888609e-31 -7.888609e-31 -7.888609e-31 -7.888609e-31 -7.888609e-31
[96,] -3.944305e-31 -3.944305e-31 -3.944305e-31 -3.944305e-31 -3.944305e-31
[97,] -1.972152e-31 -1.972152e-31 -1.972152e-31 -1.972152e-31 -1.972152e-31
[98,] -9.860761e-32 -9.860761e-32 -9.860761e-32 -9.860761e-32 -9.860761e-32
[99,] -4.930381e-32 -4.930381e-32 -4.930381e-32 -4.930381e-32 -4.930381e-32
[100,] -2.465190e-32 -2.465190e-32 -2.465190e-32 -2.465190e-32 -2.465190e-32
[101,] -1.232595e-32 -1.232595e-32 -1.232595e-32 -1.232595e-32 -1.232595e-32
[102,] -6.162976e-33 -6.162976e-33 -6.162976e-33 -6.162976e-33 -6.162976e-33
[103,] -3.081488e-33 -3.081488e-33 -3.081488e-33 -3.081488e-33 -3.081488e-33
[104,] -1.540744e-33 -1.540744e-33 -1.540744e-33 -1.540744e-33 -1.540744e-33
[105,] -7.703720e-34 -7.703720e-34 -7.703720e-34 -7.703720e-34 -7.703720e-34
[106,] -3.851860e-34 -3.851860e-34 -3.851860e-34 -3.851860e-34 -3.851860e-34
[107,] -1.925930e-34 -1.925930e-34 -1.925930e-34 -1.925930e-34 -1.925930e-34
[108,] -9.629650e-35 -9.629650e-35 -9.629650e-35 -9.629650e-35 -9.629650e-35
[109,] -4.814825e-35 -4.814825e-35 -4.814825e-35 -4.814825e-35 -4.814825e-35
[110,] -2.407412e-35 -2.407412e-35 -2.407412e-35 -2.407412e-35 -2.407412e-35
[111,] -1.203706e-35 -1.203706e-35 -1.203706e-35 -1.203706e-35 -1.203706e-35
[112,] -6.018531e-36 -6.018531e-36 -6.018531e-36 -6.018531e-36 -6.018531e-36
[113,] -3.009266e-36 -3.009266e-36 -3.009266e-36 -3.009266e-36 -3.009266e-36
[114,] -1.504633e-36 -1.504633e-36 -1.504633e-36 -1.504633e-36 -1.504633e-36
[115,] -7.523164e-37 -7.523164e-37 -7.523164e-37 -7.523164e-37 -7.523164e-37
[116,] -3.761582e-37 -3.761582e-37 -3.761582e-37 -3.761582e-37 -3.761582e-37
[117,] -1.880791e-37 -1.880791e-37 -1.880791e-37 -1.880791e-37 -1.880791e-37
[118,] -9.403955e-38 -9.403955e-38 -9.403955e-38 -9.403955e-38 -9.403955e-38
[119,] -4.701977e-38 -4.701977e-38 -4.701977e-38 -4.701977e-38 -4.701977e-38
[120,] -2.350989e-38 -2.350989e-38 -2.350989e-38 -2.350989e-38 -2.350989e-38
[121,] -1.175494e-38 -1.175494e-38 -1.175494e-38 -1.175494e-38 -1.175494e-38
[122,] -5.877472e-39 -5.877472e-39 -5.877472e-39 -5.877472e-39 -5.877472e-39
[123,] -2.938736e-39 -2.938736e-39 -2.938736e-39 -2.938736e-39 -2.938736e-39
[124,] -1.469368e-39 -1.469368e-39 -1.469368e-39 -1.469368e-39 -1.469368e-39
[125,] -7.346840e-40 -7.346840e-40 -7.346840e-40 -7.346840e-40 -7.346840e-40
[126,] -3.673420e-40 -3.673420e-40 -3.673420e-40 -3.673420e-40 -3.673420e-40
[127,] -1.836710e-40 -1.836710e-40 -1.836710e-40 -1.836710e-40 -1.836710e-40
[128,] -9.183550e-41 -9.183550e-41 -9.183550e-41 -9.183550e-41 -9.183550e-41
[129,] -4.591775e-41 -4.591775e-41 -4.591775e-41 -4.591775e-41 -4.591775e-41
[130,] -2.295887e-41 -2.295887e-41 -2.295887e-41 -2.295887e-41 -2.295887e-41
[131,] -1.147944e-41 -1.147944e-41 -1.147944e-41 -1.147944e-41 -1.147944e-41
[132,] -5.739719e-42 -5.739719e-42 -5.739719e-42 -5.739719e-42 -5.739719e-42
[133,] -2.869859e-42 -2.869859e-42 -2.869859e-42 -2.869859e-42 -2.869859e-42
[134,] -1.434930e-42 -1.434930e-42 -1.434930e-42 -1.434930e-42 -1.434930e-42
[135,] -7.174648e-43 -7.174648e-43 -7.174648e-43 -7.174648e-43 -7.174648e-43
[136,] -3.587324e-43 -3.587324e-43 -3.587324e-43 -3.587324e-43 -3.587324e-43
[137,] -1.793662e-43 -1.793662e-43 -1.793662e-43 -1.793662e-43 -1.793662e-43
[138,] -8.968310e-44 -8.968310e-44 -8.968310e-44 -8.968310e-44 -8.968310e-44
[139,] -4.484155e-44 -4.484155e-44 -4.484155e-44 -4.484155e-44 -4.484155e-44
[140,] -2.242078e-44 -2.242078e-44 -2.242078e-44 -2.242078e-44 -2.242078e-44
[141,] -1.121039e-44 -1.121039e-44 -1.121039e-44 -1.121039e-44 -1.121039e-44
[142,] -5.605194e-45 -5.605194e-45 -5.605194e-45 -5.605194e-45 -5.605194e-45
[143,] -2.802597e-45 -2.802597e-45 -2.802597e-45 -2.802597e-45 -2.802597e-45
[144,] -1.401298e-45 -1.401298e-45 -1.401298e-45 -1.401298e-45 -1.401298e-45
[145,] -7.006492e-46 -7.006492e-46 -7.006492e-46 -7.006492e-46 -7.006492e-46
[146,] -3.503246e-46 -3.503246e-46 -3.503246e-46 -3.503246e-46 -3.503246e-46
[147,] -1.751623e-46 -1.751623e-46 -1.751623e-46 -1.751623e-46 -1.751623e-46
[148,] -8.758115e-47 -8.758115e-47 -8.758115e-47 -8.758115e-47 -8.758115e-47
[149,] -4.379058e-47 -4.379058e-47 -4.379058e-47 -4.379058e-47 -4.379058e-47
[150,] -2.189529e-47 -2.189529e-47 -2.189529e-47 -2.189529e-47 -2.189529e-47
[151,] -1.094764e-47 -1.094764e-47 -1.094764e-47 -1.094764e-47 -1.094764e-47
[152,] -5.473822e-48 -5.473822e-48 -5.473822e-48 -5.473822e-48 -5.473822e-48
[153,] -2.736911e-48 -2.736911e-48 -2.736911e-48 -2.736911e-48 -2.736911e-48
[154,] -1.368456e-48 -1.368456e-48 -1.368456e-48 -1.368456e-48 -1.368456e-48
[155,] -6.842278e-49 -6.842278e-49 -6.842278e-49 -6.842278e-49 -6.842278e-49
[156,] -3.421139e-49 -3.421139e-49 -3.421139e-49 -3.421139e-49 -3.421139e-49
[157,] -1.710569e-49 -1.710569e-49 -1.710569e-49 -1.710569e-49 -1.710569e-49
[158,] -8.552847e-50 -8.552847e-50 -8.552847e-50 -8.552847e-50 -8.552847e-50
[159,] -4.276424e-50 -4.276424e-50 -4.276424e-50 -4.276424e-50 -4.276424e-50
[160,] -2.138212e-50 -2.138212e-50 -2.138212e-50 -2.138212e-50 -2.138212e-50
[161,] -1.069106e-50 -1.069106e-50 -1.069106e-50 -1.069106e-50 -1.069106e-50
[162,] -5.345529e-51 -5.345529e-51 -5.345529e-51 -5.345529e-51 -5.345529e-51
[163,] -2.672765e-51 -2.672765e-51 -2.672765e-51 -2.672765e-51 -2.672765e-51
[164,] -1.336382e-51 -1.336382e-51 -1.336382e-51 -1.336382e-51 -1.336382e-51
[165,] -6.681912e-52 -6.681912e-52 -6.681912e-52 -6.681912e-52 -6.681912e-52
[166,] -3.340956e-52 -3.340956e-52 -3.340956e-52 -3.340956e-52 -3.340956e-52
[167,] -1.670478e-52 -1.670478e-52 -1.670478e-52 -1.670478e-52 -1.670478e-52
[168,] -8.352390e-53 -8.352390e-53 -8.352390e-53 -8.352390e-53 -8.352390e-53
[169,] -4.176195e-53 -4.176195e-53 -4.176195e-53 -4.176195e-53 -4.176195e-53
[170,] -2.088097e-53 -2.088097e-53 -2.088097e-53 -2.088097e-53 -2.088097e-53
[171,] -1.044049e-53 -1.044049e-53 -1.044049e-53 -1.044049e-53 -1.044049e-53
[172,] -5.220244e-54 -5.220244e-54 -5.220244e-54 -5.220244e-54 -5.220244e-54
[173,] -2.610122e-54 -2.610122e-54 -2.610122e-54 -2.610122e-54 -2.610122e-54
[174,] -1.305061e-54 -1.305061e-54 -1.305061e-54 -1.305061e-54 -1.305061e-54
[175,] -6.525304e-55 -6.525304e-55 -6.525304e-55 -6.525304e-55 -6.525304e-55
[176,] -3.262652e-55 -3.262652e-55 -3.262652e-55 -3.262652e-55 -3.262652e-55
[177,] -1.631326e-55 -1.631326e-55 -1.631326e-55 -1.631326e-55 -1.631326e-55
[178,] -8.156631e-56 -8.156631e-56 -8.156631e-56 -8.156631e-56 -8.156631e-56
[179,] -4.078315e-56 -4.078315e-56 -4.078315e-56 -4.078315e-56 -4.078315e-56
[180,] -2.039158e-56 -2.039158e-56 -2.039158e-56 -2.039158e-56 -2.039158e-56
[181,] -1.019579e-56 -1.019579e-56 -1.019579e-56 -1.019579e-56 -1.019579e-56
[182,] -5.097894e-57 -5.097894e-57 -5.097894e-57 -5.097894e-57 -5.097894e-57
[183,] -2.548947e-57 -2.548947e-57 -2.548947e-57 -2.548947e-57 -2.548947e-57
[184,] -1.274474e-57 -1.274474e-57 -1.274474e-57 -1.274474e-57 -1.274474e-57
[185,] -6.372368e-58 -6.372368e-58 -6.372368e-58 -6.372368e-58 -6.372368e-58
[186,] -3.186184e-58 -3.186184e-58 -3.186184e-58 -3.186184e-58 -3.186184e-58
[187,] -1.593092e-58 -1.593092e-58 -1.593092e-58 -1.593092e-58 -1.593092e-58
[188,] -7.965460e-59 -7.965460e-59 -7.965460e-59 -7.965460e-59 -7.965460e-59
[189,] -3.982730e-59 -3.982730e-59 -3.982730e-59 -3.982730e-59 -3.982730e-59
[190,] -1.991365e-59 -1.991365e-59 -1.991365e-59 -1.991365e-59 -1.991365e-59
[191,] -9.956824e-60 -9.956824e-60 -9.956824e-60 -9.956824e-60 -9.956824e-60
[192,] -4.978412e-60 -4.978412e-60 -4.978412e-60 -4.978412e-60 -4.978412e-60
[193,] -2.489206e-60 -2.489206e-60 -2.489206e-60 -2.489206e-60 -2.489206e-60
[194,] -1.244603e-60 -1.244603e-60 -1.244603e-60 -1.244603e-60 -1.244603e-60
[195,] -6.223015e-61 -6.223015e-61 -6.223015e-61 -6.223015e-61 -6.223015e-61
delta3
[1,] -1.448317e-02
[2,] -7.517260e-03
[3,] -3.831203e-03
[4,] -1.934208e-03
[5,] -9.718137e-04
[6,] -4.870916e-04
[7,] -2.438429e-04
[8,] -1.219958e-04
[9,] -6.101653e-05
[10,] -3.051292e-05
[11,] -1.525762e-05
[12,] -7.629104e-06
[13,] -3.814625e-06
[14,] -1.907330e-06
[15,] -9.536698e-07
[16,] -4.768360e-07
[17,] -2.384183e-07
[18,] -1.192092e-07
[19,] -5.960463e-08
[20,] -2.980232e-08
[21,] -1.490116e-08
[22,] -7.450580e-09
[23,] -3.725290e-09
[24,] -1.862645e-09
[25,] -9.313226e-10
[26,] -4.656613e-10
[27,] -2.328306e-10
[28,] -1.164153e-10
[29,] -5.820766e-11
[30,] -2.910383e-11
[31,] -1.455192e-11
[32,] -7.275958e-12
[33,] -3.637979e-12
[34,] -1.818989e-12
[35,] -9.094947e-13
[36,] -4.547474e-13
[37,] -2.273737e-13
[38,] -1.136868e-13
[39,] -5.684342e-14
[40,] -2.842171e-14
[41,] -1.421085e-14
[42,] -7.105427e-15
[43,] -3.552714e-15
[44,] -1.776357e-15
[45,] -8.881784e-16
[46,] -4.440892e-16
[47,] -2.220446e-16
[48,] -1.110223e-16
[49,] -5.551115e-17
[50,] -2.775558e-17
[51,] -1.387779e-17
[52,] -6.938894e-18
[53,] -3.469447e-18
[54,] -1.734723e-18
[55,] -8.673617e-19
[56,] -4.336809e-19
[57,] -2.168404e-19
[58,] -1.084202e-19
[59,] -5.421011e-20
[60,] -2.710505e-20
[61,] -1.355253e-20
[62,] -6.776264e-21
[63,] -3.388132e-21
[64,] -1.694066e-21
[65,] -8.470329e-22
[66,] -4.235165e-22
[67,] -2.117582e-22
[68,] -1.058791e-22
[69,] -5.293956e-23
[70,] -2.646978e-23
[71,] -1.323489e-23
[72,] -6.617445e-24
[73,] -3.308722e-24
[74,] -1.654361e-24
[75,] -8.271806e-25
[76,] -4.135903e-25
[77,] -2.067952e-25
[78,] -1.033976e-25
[79,] -5.169879e-26
[80,] -2.584939e-26
[81,] -1.292470e-26
[82,] -6.462349e-27
[83,] -3.231174e-27
[84,] -1.615587e-27
[85,] -8.077936e-28
[86,] -4.038968e-28
[87,] -2.019484e-28
[88,] -1.009742e-28
[89,] -5.048710e-29
[90,] -2.524355e-29
[91,] -1.262177e-29
[92,] -6.310887e-30
[93,] -3.155444e-30
[94,] -1.577722e-30
[95,] -7.888609e-31
[96,] -3.944305e-31
[97,] -1.972152e-31
[98,] -9.860761e-32
[99,] -4.930381e-32
[100,] -2.465190e-32
[101,] -1.232595e-32
[102,] -6.162976e-33
[103,] -3.081488e-33
[104,] -1.540744e-33
[105,] -7.703720e-34
[106,] -3.851860e-34
[107,] -1.925930e-34
[108,] -9.629650e-35
[109,] -4.814825e-35
[110,] -2.407412e-35
[111,] -1.203706e-35
[112,] -6.018531e-36
[113,] -3.009266e-36
[114,] -1.504633e-36
[115,] -7.523164e-37
[116,] -3.761582e-37
[117,] -1.880791e-37
[118,] -9.403955e-38
[119,] -4.701977e-38
[120,] -2.350989e-38
[121,] -1.175494e-38
[122,] -5.877472e-39
[123,] -2.938736e-39
[124,] -1.469368e-39
[125,] -7.346840e-40
[126,] -3.673420e-40
[127,] -1.836710e-40
[128,] -9.183550e-41
[129,] -4.591775e-41
[130,] -2.295887e-41
[131,] -1.147944e-41
[132,] -5.739719e-42
[133,] -2.869859e-42
[134,] -1.434930e-42
[135,] -7.174648e-43
[136,] -3.587324e-43
[137,] -1.793662e-43
[138,] -8.968310e-44
[139,] -4.484155e-44
[140,] -2.242078e-44
[141,] -1.121039e-44
[142,] -5.605194e-45
[143,] -2.802597e-45
[144,] -1.401298e-45
[145,] -7.006492e-46
[146,] -3.503246e-46
[147,] -1.751623e-46
[148,] -8.758115e-47
[149,] -4.379058e-47
[150,] -2.189529e-47
[151,] -1.094764e-47
[152,] -5.473822e-48
[153,] -2.736911e-48
[154,] -1.368456e-48
[155,] -6.842278e-49
[156,] -3.421139e-49
[157,] -1.710569e-49
[158,] -8.552847e-50
[159,] -4.276424e-50
[160,] -2.138212e-50
[161,] -1.069106e-50
[162,] -5.345529e-51
[163,] -2.672765e-51
[164,] -1.336382e-51
[165,] -6.681912e-52
[166,] -3.340956e-52
[167,] -1.670478e-52
[168,] -8.352390e-53
[169,] -4.176195e-53
[170,] -2.088097e-53
[171,] -1.044049e-53
[172,] -5.220244e-54
[173,] -2.610122e-54
[174,] -1.305061e-54
[175,] -6.525304e-55
[176,] -3.262652e-55
[177,] -1.631326e-55
[178,] -8.156631e-56
[179,] -4.078315e-56
[180,] -2.039158e-56
[181,] -1.019579e-56
[182,] -5.097894e-57
[183,] -2.548947e-57
[184,] -1.274474e-57
[185,] -6.372368e-58
[186,] -3.186184e-58
[187,] -1.593092e-58
[188,] -7.965460e-59
[189,] -3.982730e-59
[190,] -1.991365e-59
[191,] -9.956824e-60
[192,] -4.978412e-60
[193,] -2.489206e-60
[194,] -1.244603e-60
[195,] -6.223015e-61
>
> ## testing cbind() i.e., the Rmpfr-method a bit here !
> chkDn <- function(r) stopifnot(identical(
+ dimnames(r),
+ list(NULL, c("x", "log2x", eval(formals(eps_S53)$method)))))
> chkDn2 <- function(r) stopifnot(identical(
+ colnames(r)[-(1:2)], eval(formals(eps_S53)$method)))
> chkDn ( r1 <- cbind(x, log2x = log2(x), EA))
> chkDn2( r2 <- selectMethod("cbind","mpfr")(x, log2x = log2(x), EA) )
> ## are the same if names are *made* to match:
> stopifnot(identical(r1, local({t <- r2; dimnames(t) <- dimnames(r1); t})))
> chkDn2( r3 <- selectMethod("cbind","mpfr")(x, log2(x), epsAllS53(x)) )
> chkDn2( r4 <- cbind(x, log2(x), epsAllS53(x)) )
> stopifnot(identical(r1, local({t <- r4; dimnames(t) <- dimnames(r1); t})))
>
> U_S53mat <- function(x) cbind(x, log2x = log2(x), epsAllS53(x))
> chkDn(Rx <- U_S53mat(x))
> dimnames(Rx) # .. all empty, i.e., list(NULL, NULL) ... even though epsAll..() gives nice colnames
[[1]]
NULL
[[2]]
[1] "x" "log2x" "direct" "log1pExct" "deltaFull" "delta1"
[7] "delta2" "delta3"
> ##
> if(FALSE) ## no longer needed !
+ ## A version with correct colnames :
+ U_S53mat <- function(x) {
+ eA <- epsAllS53(x)
+ r <- cbind(x, log2(x), eA)
+ colnames(r) <- c("x", "log2.x", colnames(eA))
+ r
+ }
>
> stopifnot(Rx[,-(1:2)] < 0)
>
> op <- options(width = 100)
> print(asNumeric(Rx), digits=5) ; options(op)
x log2x direct log1pExct deltaFull delta1 delta2 delta3
[1,] 8.0000e+00 3.0 -1.5051e-02 -1.4585e-02 -1.4479e-02 -1.5385e-02 -1.4423e-02 -1.4483e-02
[2,] 1.1314e+01 3.5 -7.6647e-03 -7.5454e-03 -7.5170e-03 -7.7519e-03 -7.5097e-03 -7.5173e-03
[3,] 1.6000e+01 4.0 -3.8687e-03 -3.8385e-03 -3.8312e-03 -3.8911e-03 -3.8303e-03 -3.8312e-03
[4,] 2.2627e+01 4.5 -1.9437e-03 -1.9361e-03 -1.9342e-03 -1.9493e-03 -1.9341e-03 -1.9342e-03
[5,] 3.2000e+01 5.0 -9.7419e-04 -9.7229e-04 -9.7181e-04 -9.7561e-04 -9.7180e-04 -9.7181e-04
[6,] 4.5255e+01 5.5 -4.8769e-04 -4.8721e-04 -4.8709e-04 -4.8804e-04 -4.8709e-04 -4.8709e-04
[7,] 6.4000e+01 6.0 -2.4399e-04 -2.4387e-04 -2.4384e-04 -2.4408e-04 -2.4384e-04 -2.4384e-04
[8,] 9.0510e+01 6.5 -1.2203e-04 -1.2200e-04 -1.2200e-04 -1.2206e-04 -1.2200e-04 -1.2200e-04
[9,] 1.2800e+02 7.0 -6.1026e-05 -6.1018e-05 -6.1017e-05 -6.1031e-05 -6.1017e-05 -6.1017e-05
[10,] 1.8102e+02 7.5 -3.0515e-05 -3.0513e-05 -3.0513e-05 -3.0517e-05 -3.0513e-05 -3.0513e-05
[11,] 2.5600e+02 8.0 -1.5258e-05 -1.5258e-05 -1.5258e-05 -1.5259e-05 -1.5258e-05 -1.5258e-05
[12,] 3.6204e+02 8.5 -7.6292e-06 -7.6291e-06 -7.6291e-06 -7.6293e-06 -7.6291e-06 -7.6291e-06
[13,] 5.1200e+02 9.0 -3.8147e-06 -3.8146e-06 -3.8146e-06 -3.8147e-06 -3.8146e-06 -3.8146e-06
[14,] 7.2408e+02 9.5 -1.9073e-06 -1.9073e-06 -1.9073e-06 -1.9073e-06 -1.9073e-06 -1.9073e-06
[15,] 1.0240e+03 10.0 -9.5367e-07 -9.5367e-07 -9.5367e-07 -9.5367e-07 -9.5367e-07 -9.5367e-07
[16,] 1.4482e+03 10.5 -4.7684e-07 -4.7684e-07 -4.7684e-07 -4.7684e-07 -4.7684e-07 -4.7684e-07
[17,] 2.0480e+03 11.0 -2.3842e-07 -2.3842e-07 -2.3842e-07 -2.3842e-07 -2.3842e-07 -2.3842e-07
[18,] 2.8963e+03 11.5 -1.1921e-07 -1.1921e-07 -1.1921e-07 -1.1921e-07 -1.1921e-07 -1.1921e-07
[19,] 4.0960e+03 12.0 -5.9605e-08 -5.9605e-08 -5.9605e-08 -5.9605e-08 -5.9605e-08 -5.9605e-08
[20,] 5.7926e+03 12.5 -2.9802e-08 -2.9802e-08 -2.9802e-08 -2.9802e-08 -2.9802e-08 -2.9802e-08
[21,] 8.1920e+03 13.0 -1.4901e-08 -1.4901e-08 -1.4901e-08 -1.4901e-08 -1.4901e-08 -1.4901e-08
[22,] 1.1585e+04 13.5 -7.4506e-09 -7.4506e-09 -7.4506e-09 -7.4506e-09 -7.4506e-09 -7.4506e-09
[23,] 1.6384e+04 14.0 -3.7253e-09 -3.7253e-09 -3.7253e-09 -3.7253e-09 -3.7253e-09 -3.7253e-09
[24,] 2.3170e+04 14.5 -1.8626e-09 -1.8626e-09 -1.8626e-09 -1.8626e-09 -1.8626e-09 -1.8626e-09
[25,] 3.2768e+04 15.0 -9.3132e-10 -9.3132e-10 -9.3132e-10 -9.3132e-10 -9.3132e-10 -9.3132e-10
[26,] 4.6341e+04 15.5 -4.6566e-10 -4.6566e-10 -4.6566e-10 -4.6566e-10 -4.6566e-10 -4.6566e-10
[27,] 6.5536e+04 16.0 -2.3283e-10 -2.3283e-10 -2.3283e-10 -2.3283e-10 -2.3283e-10 -2.3283e-10
[28,] 9.2682e+04 16.5 -1.1642e-10 -1.1642e-10 -1.1642e-10 -1.1642e-10 -1.1642e-10 -1.1642e-10
[29,] 1.3107e+05 17.0 -5.8208e-11 -5.8208e-11 -5.8208e-11 -5.8208e-11 -5.8208e-11 -5.8208e-11
[30,] 1.8536e+05 17.5 -2.9104e-11 -2.9104e-11 -2.9104e-11 -2.9104e-11 -2.9104e-11 -2.9104e-11
[31,] 2.6214e+05 18.0 -1.4552e-11 -1.4552e-11 -1.4552e-11 -1.4552e-11 -1.4552e-11 -1.4552e-11
[32,] 3.7073e+05 18.5 -7.2760e-12 -7.2760e-12 -7.2760e-12 -7.2760e-12 -7.2760e-12 -7.2760e-12
[33,] 5.2429e+05 19.0 -3.6380e-12 -3.6380e-12 -3.6380e-12 -3.6380e-12 -3.6380e-12 -3.6380e-12
[34,] 7.4146e+05 19.5 -1.8190e-12 -1.8190e-12 -1.8190e-12 -1.8190e-12 -1.8190e-12 -1.8190e-12
[35,] 1.0486e+06 20.0 -9.0949e-13 -9.0949e-13 -9.0949e-13 -9.0949e-13 -9.0949e-13 -9.0949e-13
[36,] 1.4829e+06 20.5 -4.5475e-13 -4.5475e-13 -4.5475e-13 -4.5475e-13 -4.5475e-13 -4.5475e-13
[37,] 2.0972e+06 21.0 -2.2737e-13 -2.2737e-13 -2.2737e-13 -2.2737e-13 -2.2737e-13 -2.2737e-13
[38,] 2.9658e+06 21.5 -1.1369e-13 -1.1369e-13 -1.1369e-13 -1.1369e-13 -1.1369e-13 -1.1369e-13
[39,] 4.1943e+06 22.0 -5.6843e-14 -5.6843e-14 -5.6843e-14 -5.6843e-14 -5.6843e-14 -5.6843e-14
[40,] 5.9316e+06 22.5 -2.8422e-14 -2.8422e-14 -2.8422e-14 -2.8422e-14 -2.8422e-14 -2.8422e-14
[41,] 8.3886e+06 23.0 -1.4211e-14 -1.4211e-14 -1.4211e-14 -1.4211e-14 -1.4211e-14 -1.4211e-14
[42,] 1.1863e+07 23.5 -7.1054e-15 -7.1054e-15 -7.1054e-15 -7.1054e-15 -7.1054e-15 -7.1054e-15
[43,] 1.6777e+07 24.0 -3.5527e-15 -3.5527e-15 -3.5527e-15 -3.5527e-15 -3.5527e-15 -3.5527e-15
[44,] 2.3727e+07 24.5 -1.7764e-15 -1.7764e-15 -1.7764e-15 -1.7764e-15 -1.7764e-15 -1.7764e-15
[45,] 3.3554e+07 25.0 -8.8818e-16 -8.8818e-16 -8.8818e-16 -8.8818e-16 -8.8818e-16 -8.8818e-16
[46,] 4.7453e+07 25.5 -4.4409e-16 -4.4409e-16 -4.4409e-16 -4.4409e-16 -4.4409e-16 -4.4409e-16
[47,] 6.7109e+07 26.0 -2.2204e-16 -2.2204e-16 -2.2204e-16 -2.2204e-16 -2.2204e-16 -2.2204e-16
[48,] 9.4906e+07 26.5 -1.1102e-16 -1.1102e-16 -1.1102e-16 -1.1102e-16 -1.1102e-16 -1.1102e-16
[49,] 1.3422e+08 27.0 -5.5511e-17 -5.5511e-17 -5.5511e-17 -5.5511e-17 -5.5511e-17 -5.5511e-17
[50,] 1.8981e+08 27.5 -2.7756e-17 -2.7756e-17 -2.7756e-17 -2.7756e-17 -2.7756e-17 -2.7756e-17
[51,] 2.6844e+08 28.0 -1.3878e-17 -1.3878e-17 -1.3878e-17 -1.3878e-17 -1.3878e-17 -1.3878e-17
[52,] 3.7963e+08 28.5 -6.9389e-18 -6.9389e-18 -6.9389e-18 -6.9389e-18 -6.9389e-18 -6.9389e-18
[53,] 5.3687e+08 29.0 -3.4694e-18 -3.4694e-18 -3.4694e-18 -3.4694e-18 -3.4694e-18 -3.4694e-18
[54,] 7.5925e+08 29.5 -1.7347e-18 -1.7347e-18 -1.7347e-18 -1.7347e-18 -1.7347e-18 -1.7347e-18
[55,] 1.0737e+09 30.0 -8.6736e-19 -8.6736e-19 -8.6736e-19 -8.6736e-19 -8.6736e-19 -8.6736e-19
[56,] 1.5185e+09 30.5 -4.3368e-19 -4.3368e-19 -4.3368e-19 -4.3368e-19 -4.3368e-19 -4.3368e-19
[57,] 2.1475e+09 31.0 -2.1684e-19 -2.1684e-19 -2.1684e-19 -2.1684e-19 -2.1684e-19 -2.1684e-19
[58,] 3.0370e+09 31.5 -1.0842e-19 -1.0842e-19 -1.0842e-19 -1.0842e-19 -1.0842e-19 -1.0842e-19
[59,] 4.2950e+09 32.0 -5.4210e-20 -5.4210e-20 -5.4210e-20 -5.4210e-20 -5.4210e-20 -5.4210e-20
[60,] 6.0740e+09 32.5 -2.7105e-20 -2.7105e-20 -2.7105e-20 -2.7105e-20 -2.7105e-20 -2.7105e-20
[61,] 8.5899e+09 33.0 -1.3553e-20 -1.3553e-20 -1.3553e-20 -1.3553e-20 -1.3553e-20 -1.3553e-20
[62,] 1.2148e+10 33.5 -6.7763e-21 -6.7763e-21 -6.7763e-21 -6.7763e-21 -6.7763e-21 -6.7763e-21
[63,] 1.7180e+10 34.0 -3.3881e-21 -3.3881e-21 -3.3881e-21 -3.3881e-21 -3.3881e-21 -3.3881e-21
[64,] 2.4296e+10 34.5 -1.6941e-21 -1.6941e-21 -1.6941e-21 -1.6941e-21 -1.6941e-21 -1.6941e-21
[65,] 3.4360e+10 35.0 -8.4703e-22 -8.4703e-22 -8.4703e-22 -8.4703e-22 -8.4703e-22 -8.4703e-22
[66,] 4.8592e+10 35.5 -4.2352e-22 -4.2352e-22 -4.2352e-22 -4.2352e-22 -4.2352e-22 -4.2352e-22
[67,] 6.8719e+10 36.0 -2.1176e-22 -2.1176e-22 -2.1176e-22 -2.1176e-22 -2.1176e-22 -2.1176e-22
[68,] 9.7184e+10 36.5 -1.0588e-22 -1.0588e-22 -1.0588e-22 -1.0588e-22 -1.0588e-22 -1.0588e-22
[69,] 1.3744e+11 37.0 -5.2940e-23 -5.2940e-23 -5.2940e-23 -5.2940e-23 -5.2940e-23 -5.2940e-23
[70,] 1.9437e+11 37.5 -2.6470e-23 -2.6470e-23 -2.6470e-23 -2.6470e-23 -2.6470e-23 -2.6470e-23
[71,] 2.7488e+11 38.0 -1.3235e-23 -1.3235e-23 -1.3235e-23 -1.3235e-23 -1.3235e-23 -1.3235e-23
[72,] 3.8874e+11 38.5 -6.6174e-24 -6.6174e-24 -6.6174e-24 -6.6174e-24 -6.6174e-24 -6.6174e-24
[73,] 5.4976e+11 39.0 -3.3087e-24 -3.3087e-24 -3.3087e-24 -3.3087e-24 -3.3087e-24 -3.3087e-24
[74,] 7.7747e+11 39.5 -1.6544e-24 -1.6544e-24 -1.6544e-24 -1.6544e-24 -1.6544e-24 -1.6544e-24
[75,] 1.0995e+12 40.0 -8.2718e-25 -8.2718e-25 -8.2718e-25 -8.2718e-25 -8.2718e-25 -8.2718e-25
[76,] 1.5549e+12 40.5 -4.1359e-25 -4.1359e-25 -4.1359e-25 -4.1359e-25 -4.1359e-25 -4.1359e-25
[77,] 2.1990e+12 41.0 -2.0680e-25 -2.0680e-25 -2.0680e-25 -2.0680e-25 -2.0680e-25 -2.0680e-25
[78,] 3.1099e+12 41.5 -1.0340e-25 -1.0340e-25 -1.0340e-25 -1.0340e-25 -1.0340e-25 -1.0340e-25
[79,] 4.3980e+12 42.0 -5.1699e-26 -5.1699e-26 -5.1699e-26 -5.1699e-26 -5.1699e-26 -5.1699e-26
[80,] 6.2198e+12 42.5 -2.5849e-26 -2.5849e-26 -2.5849e-26 -2.5849e-26 -2.5849e-26 -2.5849e-26
[81,] 8.7961e+12 43.0 -1.2925e-26 -1.2925e-26 -1.2925e-26 -1.2925e-26 -1.2925e-26 -1.2925e-26
[82,] 1.2440e+13 43.5 -6.4623e-27 -6.4623e-27 -6.4623e-27 -6.4623e-27 -6.4623e-27 -6.4623e-27
[83,] 1.7592e+13 44.0 -3.2312e-27 -3.2312e-27 -3.2312e-27 -3.2312e-27 -3.2312e-27 -3.2312e-27
[84,] 2.4879e+13 44.5 -1.6156e-27 -1.6156e-27 -1.6156e-27 -1.6156e-27 -1.6156e-27 -1.6156e-27
[85,] 3.5184e+13 45.0 -8.0779e-28 -8.0779e-28 -8.0779e-28 -8.0779e-28 -8.0779e-28 -8.0779e-28
[86,] 4.9758e+13 45.5 -4.0390e-28 -4.0390e-28 -4.0390e-28 -4.0390e-28 -4.0390e-28 -4.0390e-28
[87,] 7.0369e+13 46.0 -2.0195e-28 -2.0195e-28 -2.0195e-28 -2.0195e-28 -2.0195e-28 -2.0195e-28
[88,] 9.9516e+13 46.5 -1.0097e-28 -1.0097e-28 -1.0097e-28 -1.0097e-28 -1.0097e-28 -1.0097e-28
[89,] 1.4074e+14 47.0 -5.0487e-29 -5.0487e-29 -5.0487e-29 -5.0487e-29 -5.0487e-29 -5.0487e-29
[90,] 1.9903e+14 47.5 -2.5244e-29 -2.5244e-29 -2.5244e-29 -2.5244e-29 -2.5244e-29 -2.5244e-29
[91,] 2.8147e+14 48.0 -1.2622e-29 -1.2622e-29 -1.2622e-29 -1.2622e-29 -1.2622e-29 -1.2622e-29
[92,] 3.9807e+14 48.5 -6.3109e-30 -6.3109e-30 -6.3109e-30 -6.3109e-30 -6.3109e-30 -6.3109e-30
[93,] 5.6295e+14 49.0 -3.1554e-30 -3.1554e-30 -3.1554e-30 -3.1554e-30 -3.1554e-30 -3.1554e-30
[94,] 7.9613e+14 49.5 -1.5777e-30 -1.5777e-30 -1.5777e-30 -1.5777e-30 -1.5777e-30 -1.5777e-30
[95,] 1.1259e+15 50.0 -7.8886e-31 -7.8886e-31 -7.8886e-31 -7.8886e-31 -7.8886e-31 -7.8886e-31
[96,] 1.5923e+15 50.5 -3.9443e-31 -3.9443e-31 -3.9443e-31 -3.9443e-31 -3.9443e-31 -3.9443e-31
[97,] 2.2518e+15 51.0 -1.9722e-31 -1.9722e-31 -1.9722e-31 -1.9722e-31 -1.9722e-31 -1.9722e-31
[98,] 3.1845e+15 51.5 -9.8608e-32 -9.8608e-32 -9.8608e-32 -9.8608e-32 -9.8608e-32 -9.8608e-32
[99,] 4.5036e+15 52.0 -4.9304e-32 -4.9304e-32 -4.9304e-32 -4.9304e-32 -4.9304e-32 -4.9304e-32
[100,] 6.3691e+15 52.5 -2.4652e-32 -2.4652e-32 -2.4652e-32 -2.4652e-32 -2.4652e-32 -2.4652e-32
[101,] 9.0072e+15 53.0 -1.2326e-32 -1.2326e-32 -1.2326e-32 -1.2326e-32 -1.2326e-32 -1.2326e-32
[102,] 1.2738e+16 53.5 -6.1630e-33 -6.1630e-33 -6.1630e-33 -6.1630e-33 -6.1630e-33 -6.1630e-33
[103,] 1.8014e+16 54.0 -3.0815e-33 -3.0815e-33 -3.0815e-33 -3.0815e-33 -3.0815e-33 -3.0815e-33
[104,] 2.5476e+16 54.5 -1.5407e-33 -1.5407e-33 -1.5407e-33 -1.5407e-33 -1.5407e-33 -1.5407e-33
[105,] 3.6029e+16 55.0 -7.7037e-34 -7.7037e-34 -7.7037e-34 -7.7037e-34 -7.7037e-34 -7.7037e-34
[106,] 5.0952e+16 55.5 -3.8519e-34 -3.8519e-34 -3.8519e-34 -3.8519e-34 -3.8519e-34 -3.8519e-34
[107,] 7.2058e+16 56.0 -1.9259e-34 -1.9259e-34 -1.9259e-34 -1.9259e-34 -1.9259e-34 -1.9259e-34
[108,] 1.0190e+17 56.5 -9.6296e-35 -9.6296e-35 -9.6296e-35 -9.6296e-35 -9.6296e-35 -9.6296e-35
[109,] 1.4412e+17 57.0 -4.8148e-35 -4.8148e-35 -4.8148e-35 -4.8148e-35 -4.8148e-35 -4.8148e-35
[110,] 2.0381e+17 57.5 -2.4074e-35 -2.4074e-35 -2.4074e-35 -2.4074e-35 -2.4074e-35 -2.4074e-35
[111,] 2.8823e+17 58.0 -1.2037e-35 -1.2037e-35 -1.2037e-35 -1.2037e-35 -1.2037e-35 -1.2037e-35
[112,] 4.0762e+17 58.5 -6.0185e-36 -6.0185e-36 -6.0185e-36 -6.0185e-36 -6.0185e-36 -6.0185e-36
[113,] 5.7646e+17 59.0 -3.0093e-36 -3.0093e-36 -3.0093e-36 -3.0093e-36 -3.0093e-36 -3.0093e-36
[114,] 8.1524e+17 59.5 -1.5046e-36 -1.5046e-36 -1.5046e-36 -1.5046e-36 -1.5046e-36 -1.5046e-36
[115,] 1.1529e+18 60.0 -7.5232e-37 -7.5232e-37 -7.5232e-37 -7.5232e-37 -7.5232e-37 -7.5232e-37
[116,] 1.6305e+18 60.5 -3.7616e-37 -3.7616e-37 -3.7616e-37 -3.7616e-37 -3.7616e-37 -3.7616e-37
[117,] 2.3058e+18 61.0 -1.8808e-37 -1.8808e-37 -1.8808e-37 -1.8808e-37 -1.8808e-37 -1.8808e-37
[118,] 3.2610e+18 61.5 -9.4040e-38 -9.4040e-38 -9.4040e-38 -9.4040e-38 -9.4040e-38 -9.4040e-38
[119,] 4.6117e+18 62.0 -4.7020e-38 -4.7020e-38 -4.7020e-38 -4.7020e-38 -4.7020e-38 -4.7020e-38
[120,] 6.5219e+18 62.5 -2.3510e-38 -2.3510e-38 -2.3510e-38 -2.3510e-38 -2.3510e-38 -2.3510e-38
[121,] 9.2234e+18 63.0 -1.1755e-38 -1.1755e-38 -1.1755e-38 -1.1755e-38 -1.1755e-38 -1.1755e-38
[122,] 1.3044e+19 63.5 -5.8775e-39 -5.8775e-39 -5.8775e-39 -5.8775e-39 -5.8775e-39 -5.8775e-39
[123,] 1.8447e+19 64.0 -2.9387e-39 -2.9387e-39 -2.9387e-39 -2.9387e-39 -2.9387e-39 -2.9387e-39
[124,] 2.6088e+19 64.5 -1.4694e-39 -1.4694e-39 -1.4694e-39 -1.4694e-39 -1.4694e-39 -1.4694e-39
[125,] 3.6893e+19 65.0 -7.3468e-40 -7.3468e-40 -7.3468e-40 -7.3468e-40 -7.3468e-40 -7.3468e-40
[126,] 5.2175e+19 65.5 -3.6734e-40 -3.6734e-40 -3.6734e-40 -3.6734e-40 -3.6734e-40 -3.6734e-40
[127,] 7.3787e+19 66.0 -1.8367e-40 -1.8367e-40 -1.8367e-40 -1.8367e-40 -1.8367e-40 -1.8367e-40
[128,] 1.0435e+20 66.5 -9.1835e-41 -9.1835e-41 -9.1835e-41 -9.1835e-41 -9.1835e-41 -9.1835e-41
[129,] 1.4757e+20 67.0 -4.5918e-41 -4.5918e-41 -4.5918e-41 -4.5918e-41 -4.5918e-41 -4.5918e-41
[130,] 2.0870e+20 67.5 -2.2959e-41 -2.2959e-41 -2.2959e-41 -2.2959e-41 -2.2959e-41 -2.2959e-41
[131,] 2.9515e+20 68.0 -1.1479e-41 -1.1479e-41 -1.1479e-41 -1.1479e-41 -1.1479e-41 -1.1479e-41
[132,] 4.1740e+20 68.5 -5.7397e-42 -5.7397e-42 -5.7397e-42 -5.7397e-42 -5.7397e-42 -5.7397e-42
[133,] 5.9030e+20 69.0 -2.8699e-42 -2.8699e-42 -2.8699e-42 -2.8699e-42 -2.8699e-42 -2.8699e-42
[134,] 8.3480e+20 69.5 -1.4349e-42 -1.4349e-42 -1.4349e-42 -1.4349e-42 -1.4349e-42 -1.4349e-42
[135,] 1.1806e+21 70.0 -7.1746e-43 -7.1746e-43 -7.1746e-43 -7.1746e-43 -7.1746e-43 -7.1746e-43
[136,] 1.6696e+21 70.5 -3.5873e-43 -3.5873e-43 -3.5873e-43 -3.5873e-43 -3.5873e-43 -3.5873e-43
[137,] 2.3612e+21 71.0 -1.7937e-43 -1.7937e-43 -1.7937e-43 -1.7937e-43 -1.7937e-43 -1.7937e-43
[138,] 3.3392e+21 71.5 -8.9683e-44 -8.9683e-44 -8.9683e-44 -8.9683e-44 -8.9683e-44 -8.9683e-44
[139,] 4.7224e+21 72.0 -4.4842e-44 -4.4842e-44 -4.4842e-44 -4.4842e-44 -4.4842e-44 -4.4842e-44
[140,] 6.6784e+21 72.5 -2.2421e-44 -2.2421e-44 -2.2421e-44 -2.2421e-44 -2.2421e-44 -2.2421e-44
[141,] 9.4447e+21 73.0 -1.1210e-44 -1.1210e-44 -1.1210e-44 -1.1210e-44 -1.1210e-44 -1.1210e-44
[142,] 1.3357e+22 73.5 -5.6052e-45 -5.6052e-45 -5.6052e-45 -5.6052e-45 -5.6052e-45 -5.6052e-45
[143,] 1.8889e+22 74.0 -2.8026e-45 -2.8026e-45 -2.8026e-45 -2.8026e-45 -2.8026e-45 -2.8026e-45
[144,] 2.6714e+22 74.5 -1.4013e-45 -1.4013e-45 -1.4013e-45 -1.4013e-45 -1.4013e-45 -1.4013e-45
[145,] 3.7779e+22 75.0 -7.0065e-46 -7.0065e-46 -7.0065e-46 -7.0065e-46 -7.0065e-46 -7.0065e-46
[146,] 5.3427e+22 75.5 -3.5032e-46 -3.5032e-46 -3.5032e-46 -3.5032e-46 -3.5032e-46 -3.5032e-46
[147,] 7.5558e+22 76.0 -1.7516e-46 -1.7516e-46 -1.7516e-46 -1.7516e-46 -1.7516e-46 -1.7516e-46
[148,] 1.0685e+23 76.5 -8.7581e-47 -8.7581e-47 -8.7581e-47 -8.7581e-47 -8.7581e-47 -8.7581e-47
[149,] 1.5112e+23 77.0 -4.3791e-47 -4.3791e-47 -4.3791e-47 -4.3791e-47 -4.3791e-47 -4.3791e-47
[150,] 2.1371e+23 77.5 -2.1895e-47 -2.1895e-47 -2.1895e-47 -2.1895e-47 -2.1895e-47 -2.1895e-47
[151,] 3.0223e+23 78.0 -1.0948e-47 -1.0948e-47 -1.0948e-47 -1.0948e-47 -1.0948e-47 -1.0948e-47
[152,] 4.2742e+23 78.5 -5.4738e-48 -5.4738e-48 -5.4738e-48 -5.4738e-48 -5.4738e-48 -5.4738e-48
[153,] 6.0446e+23 79.0 -2.7369e-48 -2.7369e-48 -2.7369e-48 -2.7369e-48 -2.7369e-48 -2.7369e-48
[154,] 8.5484e+23 79.5 -1.3685e-48 -1.3685e-48 -1.3685e-48 -1.3685e-48 -1.3685e-48 -1.3685e-48
[155,] 1.2089e+24 80.0 -6.8423e-49 -6.8423e-49 -6.8423e-49 -6.8423e-49 -6.8423e-49 -6.8423e-49
[156,] 1.7097e+24 80.5 -3.4211e-49 -3.4211e-49 -3.4211e-49 -3.4211e-49 -3.4211e-49 -3.4211e-49
[157,] 2.4179e+24 81.0 -1.7106e-49 -1.7106e-49 -1.7106e-49 -1.7106e-49 -1.7106e-49 -1.7106e-49
[158,] 3.4194e+24 81.5 -8.5528e-50 -8.5528e-50 -8.5528e-50 -8.5528e-50 -8.5528e-50 -8.5528e-50
[159,] 4.8357e+24 82.0 -4.2764e-50 -4.2764e-50 -4.2764e-50 -4.2764e-50 -4.2764e-50 -4.2764e-50
[160,] 6.8387e+24 82.5 -2.1382e-50 -2.1382e-50 -2.1382e-50 -2.1382e-50 -2.1382e-50 -2.1382e-50
[161,] 9.6714e+24 83.0 -1.0691e-50 -1.0691e-50 -1.0691e-50 -1.0691e-50 -1.0691e-50 -1.0691e-50
[162,] 1.3677e+25 83.5 -5.3455e-51 -5.3455e-51 -5.3455e-51 -5.3455e-51 -5.3455e-51 -5.3455e-51
[163,] 1.9343e+25 84.0 -2.6728e-51 -2.6728e-51 -2.6728e-51 -2.6728e-51 -2.6728e-51 -2.6728e-51
[164,] 2.7355e+25 84.5 -1.3364e-51 -1.3364e-51 -1.3364e-51 -1.3364e-51 -1.3364e-51 -1.3364e-51
[165,] 3.8686e+25 85.0 -6.6819e-52 -6.6819e-52 -6.6819e-52 -6.6819e-52 -6.6819e-52 -6.6819e-52
[166,] 5.4710e+25 85.5 -3.3410e-52 -3.3410e-52 -3.3410e-52 -3.3410e-52 -3.3410e-52 -3.3410e-52
[167,] 7.7371e+25 86.0 -1.6705e-52 -1.6705e-52 -1.6705e-52 -1.6705e-52 -1.6705e-52 -1.6705e-52
[168,] 1.0942e+26 86.5 -8.3524e-53 -8.3524e-53 -8.3524e-53 -8.3524e-53 -8.3524e-53 -8.3524e-53
[169,] 1.5474e+26 87.0 -4.1762e-53 -4.1762e-53 -4.1762e-53 -4.1762e-53 -4.1762e-53 -4.1762e-53
[170,] 2.1884e+26 87.5 -2.0881e-53 -2.0881e-53 -2.0881e-53 -2.0881e-53 -2.0881e-53 -2.0881e-53
[171,] 3.0949e+26 88.0 -1.0440e-53 -1.0440e-53 -1.0440e-53 -1.0440e-53 -1.0440e-53 -1.0440e-53
[172,] 4.3768e+26 88.5 -5.2202e-54 -5.2202e-54 -5.2202e-54 -5.2202e-54 -5.2202e-54 -5.2202e-54
[173,] 6.1897e+26 89.0 -2.6101e-54 -2.6101e-54 -2.6101e-54 -2.6101e-54 -2.6101e-54 -2.6101e-54
[174,] 8.7536e+26 89.5 -1.3051e-54 -1.3051e-54 -1.3051e-54 -1.3051e-54 -1.3051e-54 -1.3051e-54
[175,] 1.2379e+27 90.0 -6.5253e-55 -6.5253e-55 -6.5253e-55 -6.5253e-55 -6.5253e-55 -6.5253e-55
[176,] 1.7507e+27 90.5 -3.2627e-55 -3.2627e-55 -3.2627e-55 -3.2627e-55 -3.2627e-55 -3.2627e-55
[177,] 2.4759e+27 91.0 -1.6313e-55 -1.6313e-55 -1.6313e-55 -1.6313e-55 -1.6313e-55 -1.6313e-55
[178,] 3.5014e+27 91.5 -8.1566e-56 -8.1566e-56 -8.1566e-56 -8.1566e-56 -8.1566e-56 -8.1566e-56
[179,] 4.9518e+27 92.0 -4.0783e-56 -4.0783e-56 -4.0783e-56 -4.0783e-56 -4.0783e-56 -4.0783e-56
[180,] 7.0028e+27 92.5 -2.0392e-56 -2.0392e-56 -2.0392e-56 -2.0392e-56 -2.0392e-56 -2.0392e-56
[181,] 9.9035e+27 93.0 -1.0196e-56 -1.0196e-56 -1.0196e-56 -1.0196e-56 -1.0196e-56 -1.0196e-56
[182,] 1.4006e+28 93.5 -5.0979e-57 -5.0979e-57 -5.0979e-57 -5.0979e-57 -5.0979e-57 -5.0979e-57
[183,] 1.9807e+28 94.0 -2.5489e-57 -2.5489e-57 -2.5489e-57 -2.5489e-57 -2.5489e-57 -2.5489e-57
[184,] 2.8011e+28 94.5 -1.2745e-57 -1.2745e-57 -1.2745e-57 -1.2745e-57 -1.2745e-57 -1.2745e-57
[185,] 3.9614e+28 95.0 -6.3724e-58 -6.3724e-58 -6.3724e-58 -6.3724e-58 -6.3724e-58 -6.3724e-58
[186,] 5.6023e+28 95.5 -3.1862e-58 -3.1862e-58 -3.1862e-58 -3.1862e-58 -3.1862e-58 -3.1862e-58
[187,] 7.9228e+28 96.0 -1.5931e-58 -1.5931e-58 -1.5931e-58 -1.5931e-58 -1.5931e-58 -1.5931e-58
[188,] 1.1205e+29 96.5 -7.9655e-59 -7.9655e-59 -7.9655e-59 -7.9655e-59 -7.9655e-59 -7.9655e-59
[189,] 1.5846e+29 97.0 -3.9827e-59 -3.9827e-59 -3.9827e-59 -3.9827e-59 -3.9827e-59 -3.9827e-59
[190,] 2.2409e+29 97.5 -1.9914e-59 -1.9914e-59 -1.9914e-59 -1.9914e-59 -1.9914e-59 -1.9914e-59
[191,] 3.1691e+29 98.0 -9.9568e-60 -9.9568e-60 -9.9568e-60 -9.9568e-60 -9.9568e-60 -9.9568e-60
[192,] 4.4818e+29 98.5 -4.9784e-60 -4.9784e-60 -4.9784e-60 -4.9784e-60 -4.9784e-60 -4.9784e-60
[193,] 6.3383e+29 99.0 -2.4892e-60 -2.4892e-60 -2.4892e-60 -2.4892e-60 -2.4892e-60 -2.4892e-60
[194,] 8.9636e+29 99.5 -1.2446e-60 -1.2446e-60 -1.2446e-60 -1.2446e-60 -1.2446e-60 -1.2446e-60
[195,] 1.2677e+30 100.0 -6.2230e-61 -6.2230e-61 -6.2230e-61 -6.2230e-61 -6.2230e-61 -6.2230e-61
>
> matplot(Rx[,"x"], - Rx[,-(1:2)], type = "b", log="xy")
Error in xy.coords(x, y, xlabel, ylabel, log = log) :
'list' object cannot be coerced to type 'double'
Calls: matplot -> xy.coords
Execution halted
Flavor: r-oldrel-windows-ix86+x86_64