clear

A=xlsread('pdata.xlsx',1);


%time series plot
figure
plot(A(:,1),A(:,2),'ro')

lsline

grid on
xlabel('year')
ylabel('precipitation depth (mm)')

figure
hist(A(:,2))

figure
ksdensity(A(:,2))

figure
probplot(A(:,2))


figure
cdfplot(A(:,2))


% KS test

[h, p]= kstest(zscore(A(:,2)));



%summary statistics

for j=2:7
figure

subplot(2,2,1)
probplot(A(:,j));
grid on

subplot(2,2,2)
hist(A(:,j));
grid on

subplot(2,2,3)
plot(A(:,j))
grid on

subplot(2,2,4)
autocorr(A(:,j))
grid on

end










