Tietääkö joku millä softalla voi sovittaa vaikkapa normaalijakauma datapisteisiin, jossa näkyy vaikkapa noin normaalijakautunut huippu? Matlab ei käy, koska ei ole rahaa kaupallisiin softiin, pitäisi olla Linuxilla pyörivä softa. Ohjelman täytyy myös osata laskea sovituksen parametrien virhearviot.
Tiedän vastauksen itse, mutta luulen että tällä palstalla pyörii niin vitun epäpätevää porukkaa, ettei ole tosikaan.
Tiedän vastauksen itse, mutta luulen että tällä palstalla pyörii niin vitun epäpätevää porukkaa, ettei ole tosikaan.
Vastaukset 12
Vanhimmat
No, itse olen laskenut matematiikka noin 20 vuotta, enkä osaa tuohon vastata.
Tunnustetaan, olen käyttänyt Exceliä laskujen tekoon, se antaa vääriä vastauksia, mutta ei ole ihan ammattiylpeyttä korjata asiaa.
Vitun pellet!!!
Itselläni tohtorin tutkinto matematiikassa, mutta ongelma on minulle ylivoimainen. Ehdottaisin että vaikka Banachin kiintopistelausetta käytät, ja siinä se spiraali menee sisäänpäin, ja sovituksen saa samalla. L2-normia ehdotan mitaksi. Älä vedä data logaritmin läpi, ja polynomisovituksella sitten, koska en osaa sanoa toimiiko.
Laske otoksen keskiarvo ja hajonta ja ota niitä vastaava normaalijakaantuma.
Voi perkele!
Eihän sitä noin tehdä. Vaaditaan epälineaarinen sovitus, se on monen muuttujan optimointiongelma.
Varmaan Sagella onnistuu.
Käviskö tämmöinen: Laske ensin otoskeskiarvo ja otosvarianssi. Sitten käytä khii^2 testiä tarkistamaan sovituksen virhe. Sitten Sagella sovitus: http://ask.sagemath.org/question/1164/gauss-distribution-fit .
Ei käy. Te epäonnistuitte kaikki, kuten arvasinkin. Homma käy näin yksinkertaisesta gnuplotilla:
esim:
f1(x) = 1/((sqrt(2*pi)*rho))*exp(-((x-my)**2)/(2*rho**2))+c
rho = 1; my = 0.5; c = 2
fit f1(x) 'foo' using 1:2 via rho, my, c
missä 'foo' on data (x,y-tiedosto), johon normaalijakauma sovitetaan. Alkuarvaukset on tuossa annettu myös.
Tulokseksi tulee esim.:
Final set of parameters Asymptotic Standard Error
======================= ==========================
rho = 0.537129 +/- 0.007918 (1.474%)
my = 0.00391748 +/- 0.009035 (230.6%)
c = 0.277775 +/- 0.00601 (2.164%)
No, tuloksia voi laskea monella ohjelmalla. Minusta Sage on ihan yhtä pätevä laskuihin kuin Gnuplot. Olisi kiva kuulla perustelut sille, miksi Sage ei käy tuohon tarkoitukseen.
Taisit itsekin epäonnistua. :) Jotta f1(x) = 1/((sqrt(2*pi)*rho))*exp(-((x-my)**2)/(2*rho**2)) c olisi normaalijakauma, olisi oltava c=0. Eli ilmeisesti kyseessä on joku henkilö, joka ei tajua edes normaalijakauman määritelmää. Tällöin on turha tulla palstalle avautumaan siitä, miksi kukaan ei löydä softaa, jos kerran vääriä tuloksia pitää tehdä.
Ja muuten, Gnuplot sisältyy Sageen.
R
Tuo on ilmainen ja kohtuullisen helppokäyttöinen. Hyvin saa piirreltyä vaikka minkälaisia kuvioita.