Regresjonsanalyse matlab

[Nebuchadnezzar]

Har igjen noen målepunkter som jeg skal vurdere hvor godt passer med
den tilpassede grafen. Tenker det bare er greit å spørre tidlig, da jeg er borte i hele dag, og fristen er i morgo.

Uansett vanligvis når jeg skal bedømme hvor godt samsvar det er mellom målte data og en graf benytter jeg meg av [tex]R[/tex] eller korrelasjons koeffisienten.
Men denne fant jeg ikke automagisk i matlab.

Her er koden min

http://pastebin.com/tqeMDLKC

Den kjørte ikke når jeg kortet den ned =(
Og de mest relevante tallene er vel

Kode: Velg alt

k= 1.0e+03 *[  0.3397  ,  0.5096  ,  0.6589  ,  0.7907  ,  0.8917  ,  1.0191

    1.1227  ,  1.2267  ,  1.3270  ,  1.3488  ,  1.5055]

målt omega =  1.0e+05 * [ 0.0395 ,   0.0888  ,  0.1579  ,  0.2467    0.3553 ,   0.4836  ,  0.6317 ,   0.7994  ,  0.9870  ,  1.1942 , 1.4212 ]

teoretisk omega=  1.0e+05 * [
    [ 0.0468  ,  0.0953  ,  0.1626  ,  0.2469  ,  0.3303  ,  0.4624  ,  0.5948  ,  0.7524  ,  0.9302  ,  0.9726 1.3160 ]

Plottet mitt ser slik ut

http://i.imgur.com/z4Rff.png

Jeg prøvde å kjøre et standardavvikk på måleserien, og fikk da

Kode: Velg alt

p = (sum((omegakvadratteori2-omegakvadrat).^2)/(length(n)-1))^(1/2)
P = 8.2570e+03

(Ja, jeg vet det finnes en [tex]\text{std(x)}[/tex] funsksjon...
Mens når jeg prøvde å finne R,

så fikk jeg bare

Kode: Velg alt

P = corrcoef(omegakvadratteori2)
P = 1

og dette kan jo ikke stemme. Hvilken feil er det jeg gjør?
Jeg klarer ikke få ut noen fornuftige tall som sier at den korrigerte teoretiske verdien samsvarer godt med målepunktene ( eg, lav R eller innenfor ett standardavik)

[drgz]

Hint: prøv å tenk ut hva som er et naturlig resultat av å bruke corrcoef() på en vektor, eventuelt sjekk det som står på 'help corrcoef' / 'doc corrcoef'.

[Nebuchadnezzar]

Fant feilen to min senere, men måtte stikke å være lærerassistent

Selvsagt må jeg ta inn både x og y, da funker det for eksempel å skrive

Kode: Velg alt

R = \corrcoef(omegakvadrat,omegakvadratteori2)
R(2,1)

Dumme meg, sånn går det når en prøver å gjøre matematikk på morgenkvisten.