comp.soft-sys.math.mathematica Google Group

ML estimators for Box-Cox Extreme Value distribution

2009-11-23T12:11:03Z

Dear all,
I have a set of data, and I want to fit a distribution on this data. I
am particularily interested in the Box-Cox General Extreme Value
distribution, for which I have the CDF and PDF:
---
BoxCoxGEVDistribution /:
CDF[BoxCoxGEVDistribution[\[Mu ]_, \[Sigma]_, \[Xi]_, \[Phi]_], x_] :=
1 + (((Exp[-(1 + \[Xi] ( (x - \

Re: Mathematica scoping / contexts / variable localization

2009-11-23T12:10:53Z

It's called SubValues. One advantage is that it cleanly separates parameters
from variables. Another is that you can write expressions such as:
f[a, b]'[x]
You can also use it conveniently as a postfix operator:
x//f[a, b]
There would be many operations in Mathematica that would be easier to use if
they were defined with SubValues. For example the various

Re: Setting global InputAutoReplacements

2009-11-23T12:10:38Z

John, much thanks for the quick reply. Allow me to elaborate, for the
benefits of whoever might be interested in setting up replacement
aliases (I think they are very handy!)
First, back up the existing Core.nb.
Then, open the options inspector, select "Core.nb" from the first
drop-down box, and under Notebook Options -> File Options select

how to read in a number in hex and then get it in binary?

2009-11-23T12:10:27Z

Hi, newbie question I'm afraid, but I couldn't work out how to do this
from the help. I'm using Mathematica 5.2.
I've got a textfile containing a long list of numbers in hex. I can
read them in one by one, using e.g.
p=Read["hello.txt",Expression]
Then I would like to turn each number into a fixed-length binary

Re: compositions with the restricted set of integers

2009-11-23T12:00:04Z

Hi, Michael.
Your solution is indeed very memory - hungry due to the mapping of
Permutations, as you mentioned. The total number of permutations can be
easily deduced from the list of multiplicities of elements in a given
partition: n!/(n1!n2!...nk!), where n1, ..., nk are multiplicities of
elements, and n is the length of the partition: n=n1+...+nk. The

Re: Labelling a plot with maximum

2009-11-23T11:59:53Z

Hi Shalin,
here is an example. I set ContinuousAction -> False due to the delay
cause by NMaximize:
Manipulate[Column[{
Plot[f[x, \[CapitalDelta]], {x, -1.22, 1.22},
PlotRange -> {{-0.5, 0.5}, {0, 1.5}},
Epilog -> {Red,
Point[{x /. #[[2]], peak = #[[1]]}] & @
NMaximize[f[x, \[CapitalDelta]], x]}]

How to get the output of FindFit from MathLink in C program?

2009-11-23T11:59:43Z

Hi, all,
I want to write a C program that calls MathLink to use "FindFit" in
Mathematica. The output of "FindFit" is a list of rules. How can I get
those packets? Thanks for help!
Best regards,
DONG Li

Re: A little C program calling MathLink

2009-11-23T11:59:32Z

¶Àí Ð´µÀ:
Hi, all,
I found what the problem is. The output of "Print" is Null. I do not
quite understand why it is this. Where is the symbol output going?
Best regards,
DONG Li

Re: More Efficient Method

2009-11-23T11:59:11Z

If you want something less complex than Leonid's code, here are two
simple speedups for the IntervalMemberQ approach. First, it seems
to run a little faster with integer limits than with real limits.
datalist = Table[Random[Integer,11*^5],{m = 1*^6}];
targetlist = Table[{n,n+20},{n,100,m,50}];
Block[{realreject = Interval@@({.5,-.5}+#&)/@targe tlist},

Re: Setting global InputAutoReplacements

2009-11-23T11:59:21Z

Format->Edit Stylesheet...
Paste the cell I included. For a single notebook, you're done.
Click Install in the toolbar. Choose a filename and install it. That name will
now show up in the Stylesheets menu. Choose that stylesheet and it works.
Or, choose the name of an existing stylesheet (like Default.nb). This will,

Re: bug in Mathematica 7?

2009-11-23T11:58:49Z

It looks you're right. You may consider reporting this at
supp...@wolfram.com.
Cheers -- Sjoerd

Not all points plot on my graph...

2009-11-23T11:59:00Z

Why is it when I write this...
Clear[f, x]
f[x_] = (x + 2)^(2/3)
Plot[f[x], {x, -3, 0}, PlotRange -> {{-3, 0}, {-4, 4}}]
f[-3]
.. I get no values show in the graph for x<-2 even though f[-3]
evaluates to (-1)^(2/3).
(-1)^(2/3) evaluates as f[-3]=1 so why wouldn't the point (-3,1) plot
on my graph?

Re: Labelling a plot with maximum

2009-11-23T11:58:38Z

Manipulate[
max = NMaximize[f[x, \[CapitalDelta]], x];
Plot[f[x, \[CapitalDelta]], {x, -1.22, 1.22},
PlotRange -> {{-0.5, 0.5}, {0, 1.5}},
Epilog -> {Text[max, Scaled[{0, 1}], {-1, 1}], PointSize[0.025],
Point[{max[[2, 1, 2]], 0}]}], {\[CapitalDelta], 0.125, 0.5}]
As far as I know Export does not export to Quicktime. Try an .avi

Re: Color depth - Wikipedia, the free encyclopedia

2009-11-23T11:57:58Z

This looks like a troll or I am missing something here? For someone
who has been around for a while I cannot imagine that you haven't been
working with colors before.
Manipulate[
Graphics[{RGBColor[r, g, b], Disk[]}], {r, 0, 1, 1/256}, {g, 0, 1,
1/256}, {b, 0, 1, 1/256}]
Cheers -- Sjoerd
.google.com/profiles/Roger.Bag ula

how to read in a number in hex and then get it in binary?

2009-11-23T11:57:28Z