comp.soft-sys.math.mathematica Google Group

ML estimators for Box-Cox Extreme Value distribution

petri.totter...@hanken.fi (Petri Tötterman) — Mon, 11 Nov 2009 12:11:03 UT

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

djmp...@comcast.net (David Park) — Mon, 11 Nov 2009 12:10:53 UT

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

dnqu...@gmail.com (Leo Alekseyev) — Mon, 11 Nov 2009 12:10:38 UT

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?

hanrahan...@yahoo.co.uk — Mon, 11 Nov 2009 12:10:27 UT

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

lsh...@gmail.com (Leonid Shifrin) — Mon, 11 Nov 2009 12:00:04 UT

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

d...@metrohm.com (dh) — Mon, 11 Nov 2009 11:59:53 UT

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?

dongli2...@gmail.com (董理) — Mon, 11 Nov 2009 11:59:43 UT

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

dongli2...@gmail.com (董理) — Mon, 11 Nov 2009 11:59:32 UT

¶Àí Ð´µÀ:
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

koop...@sfu.ca (Ray Koopman) — Mon, 11 Nov 2009 11:59:11 UT

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

jfu...@wolfram.com (John Fultz) — Mon, 11 Nov 2009 11:59:21 UT

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?

sjoerd.c.devr...@gmail.com (Sjoerd C. de Vries) — Mon, 11 Nov 2009 11:58:49 UT

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

Not all points plot on my graph...

davidfrick2...@yahoo.com (davef) — Mon, 11 Nov 2009 11:59:00 UT

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

sjoerd.c.devr...@gmail.com (Sjoerd C. de Vries) — Mon, 11 Nov 2009 11:58:38 UT

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

sjoerd.c.devr...@gmail.com (Sjoerd C. de Vries) — Mon, 11 Nov 2009 11:57:58 UT

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?

hanrahan...@yahoo.co.uk — Mon, 11 Nov 2009 11:57:28 UT