% Matlabkommandoer til ekstraoppgave 2.
% 23.01.2008 Steffen Gr?nneberg (steffeng@math.uio.no)
% For ? l?se denne oppgaven er det klart at man ikke trenger ? bruke
% Matlab! Dette datasettet er jo ganske s? lite. Jeg har alikevel insistet
% p? ? gj?re dette, siden det illustrerer en veldig typisk situasjon i
% anvendt statistikk: man har en  tekstfil med data, har lyst til ? lese
% disse tallene inn for s? ? regne p? dem og visualisere visse aspekter
% ved dem.
%
% Merk at hvis man skulle befinne seg i en situasjon der man heller har
% f.eks en Excel-fil med data, kan Matlab importere dette direkte. Ellers
% kunne man f? Excel til ? eksportere datasettet til f.eks csv-formatet som
% er et "flatt" tekstformat -- dvs uten noe s?rlig struktur.

% Last inn filen:

% For unix:
X = importdata('~/math/stk1100/h08/uke4/foedsler.csv')
% For windows:

% Beregn relativ frekvens:
total = X.data(:,1);
antall.tvillinger = X.data(:,2);
antall.trillinger = X.data(:,3);
frekvens.tvillinger = antall.tvillinger./total;
frekvens.trillinger = antall.trillinger./total;

% Viser de resulterende tallene direkte:
frekvens.tvillinger
frekvens.trillinger

% Viser de resulterende tallene ved ? visualisere dem,
% spesifikt blir plottes de relative frekvensene sortert etter tid:
aar = X.textdata(:,1)

plot(frekvens.tvillinger)
set(gca,'XTickLabel',aar);
title('Relativ frekvens for tvillingf?dsler i forhold til totalt antall f?dsler')
xlabel('?r'), ylabel('Relativ Frekvens')

figure
plot(frekvens.trillinger)
set(gca,'XTickLabel',aar);
title('Relativ frekvens for trillingf?dsler i forhold til totalt antall f?dsler')
xlabel('?r'), ylabel('Relativ Frekvens')

% Vi tar s? og beregner samlet gjennomsnitt for de f?rste syv ?rene.
empirisk_sannsynlighet_tvillinger = sum(antall.tvillinger(1:7))/sum(total(1:7))
empirisk_sannsynlighet_tvillinger = sum(antall.trillinger(1:7))/sum(total(1:7))

X = 

          data: [10x3 double]
      textdata: {10x1 cell}
    rowheaders: {10x1 cell}


ans =

    0.0126
    0.0116
    0.0109
    0.0099
    0.0093
    0.0095
    0.0098
    0.0111
    0.0136
    0.0161


ans =

   1.0e-03 *

    0.1441
    0.1111
    0.1250
    0.1050
    0.0652
    0.0773
    0.1579
    0.3341
    0.3987
    0.4032


aar = 

    ''1951-55''
    ''1956-60''
    ''1961-65''
    ''1966-70''
    ''1971-75''
    ''1976-80''
    ''1981-85''
    ''1986-90''
    ''1991-95''
    ''1996-00''


empirisk_sannsynlighet_tvillinger =

    0.0106


empirisk_sannsynlighet_tvillinger =

   1.1191e-04

Published with MATLAB® 7.5