|
71 | 71 | CM = ColorModel(beadfile, blankfile, channels, colorfiles, colorpairfiles); |
72 | 72 | CM = set_ERF_channel_name(CM, 'FITC-A'); % Name the channel we'll use for ERF units |
73 | 73 | CM=set_dequantization(CM, 1); % important at low levels |
74 | | -CM=set_bead_plot(CM, 2); % 2 = show beads for all channels, even though only FITC will be used |
75 | | -CM=set_bead_min(CM, 1); % Don't consider beads less than this amount |
| 74 | +TASBEConfig.set('beads.rangeMin', 1); % Don't consider beads less than this (log10) amount |
76 | 75 | % Things we'll talk about in the next section... |
77 | 76 | CM=set_translation_plot(CM, true); |
78 | 77 | CM=set_translation_channel_min(CM,[2,2,2]); |
|
88 | 87 | % The FITC channel is the one that we will key off |
89 | 88 | % The identification is pretty terrible, though! |
90 | 89 | % Two things are going wrong here: |
91 | | -% 1) We're getting smearing of peaks from autofluorescence: bead_min should be raised |
| 90 | +% 1) We're getting smearing of peaks from autofluorescence: beads.rangeMin should be raised |
92 | 91 | % 2) Automatic threshold detection is not finding the right value, because there |
93 | 92 | % is not a distinct enough "valley" in the FITC graph |
94 | 93 | % Notice also that you got a warning that: "Warning: Bead calibration probably incorrect" |
95 | 94 | % When peaks are mis-detected, this typically leads to a bad fit against the |
96 | 95 | % expected sequence gaps, giving warning of failures. |
97 | 96 |
|
98 | | -CM=set_bead_min(CM, 2); % Don't consider beads less than this amount |
99 | | -CM=set_bead_peak_threshold(CM, 200); % override default peak threshold |
| 97 | +TASBEConfig.set('beads.rangeMin', 2); % Don't consider beads less than this (log10) amount |
| 98 | +TASBEConfig.set('beads.peakThreshold', 200); % override default peak threshold |
100 | 99 | CM = resolve(CM); |
101 | 100 |
|
102 | 101 | % We get a new warning: "Warning: Only one bead peak found, assuming brightest" |
|
0 commit comments