-
Notifications
You must be signed in to change notification settings - Fork 2
/
TODO-MM-nolonger.html
132 lines (117 loc) · 4.51 KB
/
TODO-MM-nolonger.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
<!doctype html public "-//w3c//dtd html 4.0 transitional//en">
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<title>No longer on MM's TODO list for R</title>
</head>
<h1 align=center>Things no longer on M.M's TODO list for R</h1>
<p>
<i> Here are things that earlier where on my
<a href="TODO-MM.html"> TODO list</a>
with the reasons why they are there no longer ...
</i>
<br>
Martin Maechler
<tt><a href="mailto:maechler@r-project.org">maechler@r-project.org</a></tt>
$Date: 2004/05/15 13:52:29 $ UTC
</p>
<hr>
<!-- ----------------------------------- -->
<H2>Things done by myself (at least partly) </H2>
<!---------------------------------------------------------------------------->
<TABLE Border width=100%>
<TR>
<TH width=38.2%> Item <TH width=61.8%> How and when done
</TR> <!------------------------------------------>
<TR>
<TD> Running Medians for library(modreg) :
We need a fast robust smoother; <tt>lowess()</tt> is <em>not</em> robust. </TD>
<TD>done for R 1.7 "unstable", in Feb.2003.</TD>
</TR>
<TR>
<TD> <A HREF="http://www.statistik.uni-dortmund.de/leute/ligges/R/scatterplot3d/">
scatterplot3d</A> from Uwe Ligges): A simple 3D scatterplot with somewhat
<i>orthogonal functionality</i> to persp()</TD>
<TD> The package has has been well available from CRAN, is mature and
documented, and we have submitted a paper (as of end 2002).
In the longer term, we will probably see real 3D graphics much more widely
available from R. </TD>
</TR>
<TR>
<TD>Implement C.Loader's dbinom(), dpois() improvements</TD>
<TD>done for R 1.2 "unstable", as of Sept.29, 2000, together with more,
using CL's bd0() function</TD>
</TR> <!------------------------------------------>
<TR>
<TD> Replace <tt>qnorm()</tt>'s algorithm with AS 241 </TD>
<TD> done for R 1.1.0 (June 15, 2000) </TD>
</TR> <!------------------------------------------>
<TR>
<TD> <tt><b>cor</b>(.)</tt> should get a new <tt>method</tt> argument, with
a `match.arg' default of <tt>c("pearson", "kendall", "spearman")</tt>
completely in line with the <tt><b>cor.test</b>(.)</tt> function in the
<tt>ctest</tt> package --- need to move C code for Kendall from ctest to
base (and coordinate all this with Kurt!).
</TD>
<TD> done for R 1.8.0 (Oct 2003) </TD>
</TR> <!------------------------------------------>
<TR>
<TD> Implement the final Newton step in <tt><b>qgamma</b>()</tt>'s C
code, suggested Morten Welinder, in bug report PR #2214
</TD>
<TD> done for R 1.9.0 (early 2004) </TD>
</TR> <!------------------------------------------>
<TR>
<TD> </TD>
<TD> </TD>
</TR> <!------------------------------------------>
</TABLE>
<H2>Things done by others </H2>
<!---------------------------------------------------------------------------->
<TABLE Border width=100%>
<TR>
<TH width=38.2%> Item <TH width=61.8%> How and when done
</TR> <!------------------------------------------>
<TR>
<TD> Numerical Integration at least as good as Quadpack (see also
<A HREF="http://sourceware.cygnus.com/gsl/">GSL -- The GNU Scientific Library
</A>)</TD>
<TD> R function integrate() since R 1.3.0 (2001-06-22);<br>
C API to same code, for R 1.5.0 (2002-04-26) </TD>
</TR> <!------------------------------------------>
<TR>
<TD> </TD>
<TD> </TD>
</TR> <!------------------------------------------>
</TABLE>
<H2>Things no longer ``interesting'' to me </H2>
<!---------------------------------------------------------------------------->
<TABLE Border width=100%>
<TR>
<TH width=38.2%> Topic <TH width=61.8%> Reasoning
</TR> <!------------------------------------------>
<TR>
<TD>
An option which would allow to use the old
(quite less precise, but a tad faster) version of <tt> pnorm( )</tt>.
<em>Not</em> a new argument to qnorm(), but rather a global
new entry in <tt>options(..)</tt>, something like
<pre>
options(prob.precision = 1e-10)
options(dpq.precision = 1e-10) </pre>
(where 1e-10 is somewhat arbitrary)
which would use faster less precise algorithms when the required
precision was set high (i.e., 1e-4 instead of 1e-10).
</TD>
<TD>
Got no feedback on the proposition;
The (considerable) drawback would have been that a global option
could change numerical results (not just the printing of them),
which should only ever happen if really deemed important.
</TD>
</TR> <!------------------------------------------>
<TR>
<TD> </TD>
<TD> </TD>
</TR> <!------------------------------------------>
</TABLE>