Re: Survo-ristikoiden lukumäärien mallintamisesta

[vastaus aiempaan viestiin]

Kirjoittaja: Seppo Mustonen
Sähköposti:    -
Päiväys: 17.10.2007 19:13

Olen tehnyt samantapaisia kokeiluja kuin Joonas, mutta
mielestäni esim. S(2,n)-lukujen kasvu on jopa eksponentiaalista
nopeampaa eli niiden logaritmien kasvu ei ole lineaarista vaan
saattaa olla kvadraattista tai jotain sen kaltaista.

Tämä näkyy mm. Joonaksen tarkastelun toisesta kuvasta

GPLOT KOE2,M,LX / XSCALE=0:_,1(1)10 YSCALE=-1:_,0(1)13 POINT=0,3 TREND=0
                  HEADER=2_x_M_-ristikoiden_lkm_(exp(LX))

josta (varsinkin kun jättää kaksi ensimmäistä havaintoa pois)
ilmenee selvä käyräviivainen kasvu.

Kasvu lienee (logaritmeilla) aika lähellä kvadraattista, mutta
saattaa olla jopa sitäkin voimakkaampaa, mikä ilmenee seuraavasta
jatkotarkastelusta:

*DATA KOE2B,a,b,v,m
m1 11.111111  RRR.RRRRR
vX  LX        RES
a3  2.890372    0.02947
*4  4.127134   -0.06087
*5  5.627621    0.03662
*6  7.044033   -0.02586
*7  8.649274    0.02459
*8 10.264896    0.00952
b9 11.948512   -0.01344
*
*MODEL M2B
*LX=a+b*X+c*X^2
*
*ESTIMATE KOE2B,M2B,CUR+1
*Estimated parameters of model M2B:
*a=-0.66504 (0.159366)
*b=1.06147 (0.0571271)
*c=0.0379477 (0.00471177)
*n=7 rss=0.007459 R^2=0.99988 nf=10
*Correlations:
*                 a      b      c
* a           1.000 -0.980  0.946
* b          -0.980  1.000 -0.990
* c           0.946 -0.990  1.000
*
*........................................
*OUTFILE=A MODE=800,800
*XSCALE=0:_,3(1)10,11:_ YSCALE=0(2)16
*GPLOT KOE2B,X,LX
*........................................
*INFILE=A  MODE=800,800
*a=-0.66504
*b=1.06147
*c=0.0379477
*XSCALE=0:_,3(1)10,11:_ YSCALE=0(2)16
*GPLOT Y(X)=a+b*X+c*X^2
*

Havaintojen vähyys vie kuitenkin pohjaa pois miltään varmoilta
päätelmiltä.

-Seppo

Vastaukset:

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!

Etusivu  |  Keskustelu
Copyright © Survo Systems 2001-2013. All rights reserved.
Updated 2013-06-15.