Re: Survo-ristikoiden lukumäärien mallintamisesta

[vastaus aiempaan viestiin]

Joonas Kauppisen kokeilu innosti minua laskemaan luvun S(2,10),
jolle ratkaisuohjelmani antoi noin puolen tunnin työn jälkeen
arvoksi 843476. Kun sen logaritmi log(843476)=13.645286727688
pannaan aikaisemman datan KOE2B jatkoksi, saadaan tulokset:

*DATA KOE2B,a,b,v,m
m11 11.111111  RRR.RRRRR
v X  LX        RES
a 3  2.890372    0.27858
* 4  4.127134   -0.03119
* 5  5.627621   -0.07723
* 6  7.044033   -0.20734
* 7  8.649274   -0.14863
* 8 10.264896   -0.07954
* 9 11.948512    0.05755
b10 13.645287    0.20780
*
*MODEL M2B
*LX=a+b*X+c*X^2
*
*ESTIMATE KOE2B,M2B,CUR+1
*Estimated parameters of model M2B:
*a=-0.776581 (0.13844)
*b=1.10691 (0.0461171)
*c=0.0338165 (0.00350622)
*n=8 rss=0.010327 R^2=0.99990 nf=11
*Correlations:
*                 a      b      c
* a           1.000 -0.976  0.937
* b          -0.976  1.000 -0.988
* c           0.937 -0.988  1.000
*
*
*Piirretään lineaarinen sovite punaisena viivana ja kvadraattinen mustana:
*........................................
*OUTFILE=A MODE=800,600
*XSCALE=0:_,3(1)10,11:_ YSCALE=0(2)16
*GPLOT KOE2B,X,LX   / TREND=[RED],0
*........................................
*INFILE=A  MODE=800,600
*a=-0.776581
*b=1.10691
*c=0.0338165
*XSCALE=0:_,3(1)10,11:_ YSCALE=0(2)16
*GPLOT Y(X)=a+b*X+c*X^2
*........................................
*

Viimeistä kuvaa katsomalla huomaa, että kvadraattinen malli
toimii lineaarista paremmin.
Edelleen on kuitenkin syytä olla varovainen johtopäätöksiä
tehdessä.

- Seppo

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!