Seminaariin 30.1.2007 osallistuneille

[vastaus aiempaan viestiin]

Palaan lämpimältään eilisessä seminaarissa käsiteltyyn
kahvikuppiesimerkkiin.

En yllätyksekseni saanut "paremman", kahden eksponentiaalisen komponen-
tin mallin parametreja estimoitua "oikein". Syynä oli, että olin poista-
nut vahingossa ensimmäisen malliyrityksen (alla riveillä 17-26) alku-
arvon T=25 (rivillä 20).
Kun se puuttui, "paremman" mallin estimoinnissa käytettiin
alkuarvoina ensimmäisen mallin tuloksia, jolloin uusi komponentti ei
pääse esille estimoinnin juuttuessa ensimmäisen mallin tuottamaan
"kuoppaan".

Aika ei muutenkaan (erilaisten sivujuonteiden johdosta :) riittänyt
kohteen täysimääräiseen käsittelyyn. Olen valmis jatkamaan aiheesta
vielä uudelleen keväämmällä, mikäli sitä toivotaan.

Ote toimituskentästä, jossa riveillä 34-39 "oikea" tulos:

  4 *This data set is from an experiment made in a Finnish secondary school
  5 *(Pohjoispuiston lukio, Hyvinkää) by Rolf Dahlström and his pupils.
  6 *
  7 *Temperature T (in C) of a coffee cup as a function of time t (in min)
  8 *
  9 *DATA CUP:(t,T) 0,69.5 1,65.0 2,62.5 3,60.0 4,58.0 5,56.0 6,54.5 7,
 10 *53.0 8,51.5 9,50.0 10,48.5 11,47.5 12,46.5 13,45.5 14,44.5 15,43.5 16,
 11 *43.0 17,42.0 18,41.0 19,40.5 20,39.5 21,39.0 22,38.5 23,38.0 24,37.5 25,
 12 *36.5 END
 13 *
 14 *GPLOT CUP,t,T
 15 *XSCALE=-5(5)30 YSCALE=25(5)75
 16 *
 17 *MODEL M1
 18 *T=T0+a1*exp(-b1*t)
 19 *
 20 *T=25
 21 *ESTIMATE CUP,M1,CUR+1 / RESULTS=0 METHOD=M
 22 *Estimated parameters of model M1:
 23 *T0=32.004 (0.65557)
 24 *a1=36.078 (0.553967)
 25 *b1=0.0774714 (0.00322618)
 26 *n=26 rss=4.540517 R^2=0.99788 nf=24
 27 *
 28 *
 29 *MODEL M2
 30 *T=T0+a1*exp(-b1*t)+a2*exp(-b2*t)
 31 *
 32 *ESTIMATE CUP,M2,CUR+1 / RESULTS=0 METHOD=M
 33 *Estimated parameters of model M2:
 34 *T0=29.2965 (0.664366)
 35 *a1=3.4019 (0.439223)
 36 *b1=0.875556 (0.222983)
 37 *a2=36.7547 (0.347846)
 38 *b2=0.0632299 (0.00282834)
 39 *n=26 rss=0.641183 R^2=0.99970 nf=963
 40 *

• Re: Seminaariin 30.1.2007 osallistuneille (Petri Palmu 31.1.2007 22:22)

Luettelo keskustelunaiheista

Survo-keskustelupalstan (2001-2013) viestit arkistoitiin aika ajoin sukrolla, joka automaattisesti rakensi viesteistä (yli 1600 kpl) HTML-muotoisen sivukokonaisuuden. Vuoden 2013 alusta Survo-keskustelua on jatkettu entistäkin aktiivisemmin osoitteessa forum.survo.fi. Tervetuloa mukaan!