Antag att man i ett anfall av ostyrbar nyfikenhet får för sig att man ska väga alla godisbitar i sin fredagsgodispåse, för att man t ex har en ny våg, och för att man bara måste få veta. Man väger bitarna och skriver upp resultatet, bit för bit, i en tabell. Och så sitter man där med en massa siffror.
Då tänker man förstås att man skulle vilja veta lite mer om egenskapen hos mängden godis. T ex vill man kanske veta vad medelgodisen väger, och vilken spridning det är på godisarna, dvs vilken standardavvikelse man har i godispåsen. Och den absolut brännande frågan är förstås är godismassor normalfördelade, eller kommer de från någon annan fördelning? Exponential? Weibull? Rektangelfördelning?
Det funderade jag på idag. Att hitta fördelningen för en datamängd. Frågan kom nämligen upp på jobbet. Man kan ju göra en hypotestest och kolla om man har en normalfördelning t ex. Men sådär lite ingenjörsmässigt göra en anpassning, fit, till data. Jag har ingen funktion för det. Hade inte.
Nu använde jag bara godisexemplet för att anknyta till en vardagstillämpning, istället för abstrakta datapunkter från god knows where. En vardagstillämpningar som mina barn väl känner till, eftersom deras fredagsgodis får väga 400 gram, med en avvikelse uppåt på 10 %. Som i ett svagt ögonblick förvandlades till 10% + 6g. Deras fredagsgodis väger alltså 446 gram. Exakt. De vet vad alla godissorter väger. De plockar och anpassar. Vår fredagsgodisshopping är en orgie i matematik, eller snarare i statistik.
Hursomhelst hittade jag funktionen allfitdist på matlab central. Jag testade den på enkelt vis.
for i = 1:1000
E(i) = random(‘exp’,2);
N(i) = random(‘norm’,2,1);
end
[DE,PDE]=allfitdist(E,’PDF’);
[DN,PDN]=allfitdist(N,’PDF’);

För exponentialfördelningen, E, fick jag tillbaka förslaget Inverse Gaussian. Inte helt nöjd med det.
För normalfördelningen fick jag förslaget Normalfördelning. Med rätt medelvärde och standardavvikelse. Bra.
I badrumsskåpet hittade jag de sista resterna av ett badskum från Lush. Frågan är nu om jag ska bada skumbad, eller om jag ska undersöka allfitdist lite noggrannare?

Och det får väl illustrera hur matematiken kommer in i vardagen på ett sätt som man knappt lägger märke till.