% Lag box plot av snodybde og snotetthet for Rybekken og Sylsjodammen
% INPUT: .txt-filer med observasjoner fra Rybekken og Sylsj?dammen med
% aarstall, snodybde, snotetthet og SWE i henholdsvis forste, andre, tredje og fjerde kolonne
clc, clear all, close all

%% Last inn data
ry = load('rybekken.txt');
sy = load('sylsjodammen.txt');

ry_depth = ry(:,2);
ry_density = ry(:,3);
sy_depth = sy(:,2);
sy_density = sy(:,3);

%% Sett sammen vektorene
depth = [ry_depth; sy_depth];
density = [ry_density; sy_density];

%% Lagre antall observasjoner
ry_num = length(ry);
sy_num = length(sy);

%% Opprett grupperingsvariabel
grp = [repmat(1,ry_num,1); repmat(2,sy_num,1)];

%% Lag box plot


boxplot(depth,grp,'Labels',{'Rybekken','Sylsjodammen'})
ylabel('Snodybde (cm)','FontSize',18)
title('Boksplott av snodybde (cm)','FontSize',16)


