15,146 reputation
12986
bio website quantdec.com
location Northeastern US
age 14
visits member for 2 years, 7 months
seen 2 days ago

Consultant (environmental and spatial stats a specialty), expert witness, and teacher. I can be reached through (outdated but still valid) links posted on my web site.

Twitter: @WilliamAHuber // ASA-P website: http://amstatphilly.org/


Why waste time learning, when ignorance is instantaneous?

--T(iger) Hobbes.

For any complex problem there is a simple solution. And it's always wrong.

--[Mis?]attributed to H.L. Mencken by Dava Sobel, Longitude.


Jan
4
comment How to get exact roots of this polynomial?
@b.gates And the next two steps are to let $x\to z/2$ to clear out powers of $2$ and then to take the big factor, $p(z)=1+3 z-3 z^2-4 z^3+z^4+z^5$ and symmetrize it via $p(z+1/z)z^5$: the primitive eleventh roots of unity pop right out.
Jan
2
comment Turn list of edges into a polygon function
@Daniel The example output is an intersection of half-planes, which can describe only a convex figure.
Jan
2
comment Turn list of edges into a polygon function
The type of example you give will work only for convex polygons. Is that a fair assumption to make in your application? If so, may we also assume the vertices have already been sorted in the order they appear around the polygon's boundary, and that the sorting follows a conventional orientation (such as keeping the interior of the polygon always to the left)? A solution for non-convex polygons can be obtained but would require more work (equivalent to triangulating them). Is it possible you only need some procedure to solve the point-in-polygon problem?
Jan
2
reviewed Leave Open Domain Coloring
Jan
2
reviewed Close Problem with SphericalPlot3D plotting
Dec
28
comment How to represent a list as a cycle
@Mr Very nice! I understand that your second method selects the lexicographically earliest representative of the class of all rotations of a list and then replaces any list by its representative. Because that is far superior in timing to the (quadratic) pattern-matching solution I posted, I'm deleting that solution.
Dec
28
comment How to visualize 3D fit
+1. It is useful at this stage to plot the residuals against the fit: that will make your comments more precise. BTW, fit["FitResiduals"] is a built-in way to obtain the residuals.
Dec
28
comment How to visualize 3D fit
A standard--and very effective--way to show the fit is to plot the residuals (equal to actual - estimated) against $x$ and $y$. For a good fit they should cluster evenly and randomly around the $xy$ plane. Because there's no question about where that plane lies, you don't even need the third dimension. For instance, many people map out the residuals in 2D using scaled and/or colored point symbols to represent their sizes and signs.
Dec
26
revised Define Log so that negative reals evaluate on lower edge of branch
added 154 characters in body
Dec
26
comment Using Transpose with a list as the second argument
+1 I am thinking the illustrations might work even better with some color coding. E.g., colors[1, 1] = Darker[Red]; colors[1, 2] = Darker[Green]; colors[1, 3] = Darker[Blue]; colors[2, 1] = Lighter[Red]; colors[2, 2] = Lighter[Green]; colors[2, 3] = Lighter[Blue]; f[x_] := Replace[MatrixForm[x], Subscript[a, i_, j_, k___] -> Style[Subscript[a, i, j, k], colors[i, j]], -1]; f[A], etc.
Dec
26
comment How do I expand a sum?
You are requesting two changes: first to Expand the products in the summation and second to Distribute the action of Sum over the addition. Consulting the help pages for Expand and Distribute will answer your question.
Dec
26
answered Define Log so that negative reals evaluate on lower edge of branch
Dec
24
revised How do I solve this equation?
Removed a potentially deceptive term from the title (this is not a functional equation).
Dec
21
comment Can Depth be used as an equivalent for MatrixQ?
@jVincent I do not know why you are addressing that last comment to me: I have not written anything here about lack of knowledge, intuition, FullForm, f[x_], or anything else you mention, nor have I ever suggested that the status of this question should be resolved based on our guesses about the O.P.'s state of mind. I have explained why I think this question, based on its merits, needs improvement, without resorting to any speculation about the knowledge or background of the person who asked it. (And--unlike 7 other people--I have not downvoted it.)
Dec
21
reviewed Leave Closed Can Depth be used as an equivalent for MatrixQ?
Dec
21
comment Creating an animation illustrating the time-evolution of a pre-computed orbit
+1 The tube is a great idea for giving a 3D visual cue.
Dec
21
comment Can Depth be used as an equivalent for MatrixQ?
I appreciate your comment, @jVincent. I would vote to reopen if the O.P. would at least explain what they think each expression is supposed to be doing. Without that, I maintain there is no question here except as each reader individually chooses to imagine it. Your personal construction of the question, as reflected in your answer, is reasonable, but whether it is a unique interpretation or the intended interpretation is difficult to determine.
Dec
21
comment An recurrence equation that can be solved in version 7 but not in version 8
It is interesting that the same problem, with $5$ changed to $4$, is solved in less than $0.1$ seconds with MMA 8, and with $5$ changed to $7$, the solution also takes less than $0.1$ seconds (although it's trivially expressed in terms of DifferenceRoot, which is just another way of giving up). With $5$ changed to $6$, a useful solution is obtained in $6$ seconds.
Dec
21
comment Can Depth be used as an equivalent for MatrixQ?
I confess to being mystified by this question. First, what is a "phrase," exactly? Second, if we take these to be Mathematica expressions, it is obvious they do completely different things, leaving me wondering what you might mean by "interchangeable." I am voting to close this as a non-question.
Dec
21
awarded  Nice Answer