-
Notifications
You must be signed in to change notification settings - Fork 3
/
abundance.m
74 lines (57 loc) · 2.03 KB
/
abundance.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
hist(score)
hist(score,[100])
hist(score(score(:,1)>0),[100])
hist(score(log(score(:,1)>0)),[100])
hist(log(score(score(:,1)>0)),[100])
hist(mw(mw(:,1)<500),[40])
hist(mw(mw(:,1)<500),[50])
hist(mw(mw(:,1)<500),[200])
hist(mw(mw(:,1)<500),[100])
xlabel('Molecular Weight')
ylabel('# of proteins')
hist(log10(abd),[100])
trad=xlsread('L:\Plasma\BlooPlasmaProt_TradMethod_Anna_3CombinedFiLes.xlsx')
hist(direct(:,11),[40])
xlabel('Molecular Weight (Da.)')
ylabel('# Peptides')
direct=xlsread('L:\Plasma\BlooPlasmaProt_DirectMethod_Anna_3CombinedFiLes.xlsx')
hist((direct(direct(:,4)>0,4)),[105])
xlabel('Molecular Weight (Da.)')
ylabel('# Peptides')
bins = linspace(0,120,25)
%bins = unique([trad(:,4);direct(:,4)]);
y1 = hist(trad(:,4), bins);
y2 = hist(direct(:,4), bins);
bar(bins, [y1;y2]');
xlabel('# Unique Peptides')
ylabel('# Proteins')
[h p v]=ttest(trad(:,4),direct(:,4))
trad=xlsread('L:\Plasma\BlooPlasmaProt_TradMethod_Anna_3CombinedFiLesPep.xlsx')
hist(trad(:,11),[40])
xlabel('Molecular Weight (Da.)')
ylabel('# Peptides')
%% synaptos
trad=xlsread('L:\Plasma\SynaptospmalProteinComplexes_TraditionalMethod_16hrsDigestion_02_2CombFiLesPep.xlsx')
hist(trad(:,11),[40])
xlabel('Molecular Weight (Da.)')
ylabel('# Peptides')
direct=xlsread('L:\Plasma\SynaptospmalProteinComplexes_DirectMethod_16hrsDigestion_2CombFiLesPep.xlsx')
hist(direct(:,11),[40])
xlabel('Molecular Weight (Da.)')
ylabel('# Peptides')
%% comp
trad=xlsread('L:\Plasma\SynaptospmalProteinComplexes_TraditionalMethod_16hrsDigestion_02_2CombFiLes.xlsx')
hist(trad(:,11),[40])
xlabel('Molecular Weight (Da.)')
ylabel('# Peptides')
direct=xlsread('L:\Plasma\SynaptospmalProteinComplexes_DirectMethod_16hrsDigestion_2CombFiLes.xlsx')
hist(direct(:,11),[40])
xlabel('Molecular Weight (Da.)')
ylabel('# Peptides')
bins = linspace(1,50,50)
%bins = unique([trad(:,4);direct(:,4)]);
y1 = hist(trad(:,4), bins);
y2 = hist(direct(:,4), bins);
bar(bins, [y1;y2]');
xlabel('# Unique Peptides')
ylabel('# Proteins')