diff --git a/functions/@mmo/checkworkspace.m b/functions/@mmo/checkworkspace.m
index 96392b53f..1e6cde6d4 100644
--- a/functions/@mmo/checkworkspace.m
+++ b/functions/@mmo/checkworkspace.m
@@ -29,7 +29,7 @@
             end
         end
                         
-        % if subfunction, must be a copie
+        % if subfunction, must be a copy
         if ~isempty(obj.workspace) && strcmpi(stack(end).file, obj.workspace(end).file) && ...
                 ~strcmpi(stack(end).name, obj.workspace(end).name)
             % We are within a subfunction. The MMO must have
diff --git a/functions/@mmo/subsasgn_old.m b/functions/@mmo/subsasgn_old.m
index c9326beba..64a31c95d 100644
--- a/functions/@mmo/subsasgn_old.m
+++ b/functions/@mmo/subsasgn_old.m
@@ -40,7 +40,7 @@
 if ~isempty(stack)
     % check if we are in a different workspace
     if ~isequal(stack, obj.workspace)
-        % if subfunction, must be a copie
+        % if subfunction, must be a copy
         if ~isempty(obj.workspace) && strcmpi(stack(end).file, obj.workspace(end).file) && ...
                 ~strcmpi(stack(end).name, obj.workspace(end).name)
             % We are within a subfunction. The MMO must have
diff --git a/functions/adminfunc/eeg_checkset.m b/functions/adminfunc/eeg_checkset.m
index abfedbbf0..190ca1b6f 100644
--- a/functions/adminfunc/eeg_checkset.m
+++ b/functions/adminfunc/eeg_checkset.m
@@ -165,7 +165,7 @@
 % 01-26-02 chandeg events and trial condition format -ad
 % 01-27-02 debug when trial condition is empty -ad
 % 02-15-02 remove icawinv recompute for pop_epoch -ad & ja
-% 02-16-02 remove last modification and test icawinv separatelly -ad
+% 02-16-02 remove last modification and test icawinv separately -ad
 % 02-16-02 empty event and epoch check -ad
 % 03-07-02 add the eeglab options -ad
 % 03-07-02 corrected typos and rate/point calculation -ad & ja
diff --git a/functions/adminfunc/gethelpvar.m b/functions/adminfunc/gethelpvar.m
index 325b1e5d7..99bb0096e 100644
--- a/functions/adminfunc/gethelpvar.m
+++ b/functions/adminfunc/gethelpvar.m
@@ -10,7 +10,7 @@
 %
 % Outputs:
 %   vartext     - variable text
-%   varnames    - name of the varialbe returned;
+%   varnames    - name of the variable returned;
 %
 % Notes: see help2html() for file format
 %
diff --git a/functions/miscfunc/corrimage.m b/functions/miscfunc/corrimage.m
index c7bdeeede..d750fdd4f 100755
--- a/functions/miscfunc/corrimage.m
+++ b/functions/miscfunc/corrimage.m
@@ -33,7 +33,7 @@
 %                percentile of sortvar values (and associated data) before
 %                computing statistics. Default is [0 100].
 %   'align'    - [float] same as 'align' parameter of erpimage(). This 
-%                parameter is used to contrain the 'times' parameter so
+%                parameter is used to constrain the 'times' parameter so
 %                correlation with data trials containing 0-values (as a
 %                result of data alignment) are avoided: computing these
 %                correlations would produce spurious significant results.
diff --git a/functions/miscfunc/eegplotgold.m b/functions/miscfunc/eegplotgold.m
index 4c2a60763..4610add6e 100644
--- a/functions/miscfunc/eegplotgold.m
+++ b/functions/miscfunc/eegplotgold.m
@@ -12,7 +12,7 @@
 %   yscaling   - initial y scaling factor (0 - default is 300)
 %   range      - how many seconds to display in window (0 -> 10)
 %
-% Note: this version of eegplotgold() reguires that your data matrix 
+% Note: this version of eegplotgold() requires that your data matrix 
 %       be defined as a global variable before running this routine. 
 %
 % Example:  >> global dataname
diff --git a/functions/miscfunc/makehtml.m b/functions/miscfunc/makehtml.m
index e825abdbf..acbfbaab5 100644
--- a/functions/miscfunc/makehtml.m
+++ b/functions/miscfunc/makehtml.m
@@ -133,7 +133,7 @@
 end
 
 % ------------------------------------------- 
-% scrips which generate a web page for eeglab
+% scripts which generate a web page for eeglab
 % ------------------------------------------- 
 STYLEHEADER = '<BR><a name="%s"></a><H2>%s</H2>\n'; 
 OPENWIN = [ '<script language="JavaScript"><!--' 10 'function openhelp(fnc){' 10 ... 
diff --git a/functions/popfunc/eeg_eegrej.m b/functions/popfunc/eeg_eegrej.m
index ddacfec42..5a0a226c6 100644
--- a/functions/popfunc/eeg_eegrej.m
+++ b/functions/popfunc/eeg_eegrej.m
@@ -188,7 +188,7 @@
     if 100*differs/length(EEG.event) > 50
         db = dbstack;
         if length(db) > 1 
-            fprintf(['bug 1971 warning for scritps older than 2017: see https://sccn.ucsd.edu/wiki/EEGLAB_bug1971\n' ]);
+            fprintf(['bug 1971 warning for scripts older than 2017: see https://sccn.ucsd.edu/wiki/EEGLAB_bug1971\n' ]);
         end
     end
     
@@ -196,7 +196,7 @@
     if ~isempty(event2)
         otherlatencies = [event2.latency];
         if ~isequal(alllats, otherlatencies)
-            warning([ 'Minor discrepency when recomputing event latency.' 10 '(this will not affect the accuracy of analyses).' 10 'Try to reproduce the problem and send us your dataset' ]);
+            warning([ 'Minor discrepancy when recomputing event latency.' 10 '(this will not affect the accuracy of analyses).' 10 'Try to reproduce the problem and send us your dataset' ]);
         end
     end
 end
diff --git a/functions/popfunc/eeg_laplac.m b/functions/popfunc/eeg_laplac.m
index 412e3b750..684828d8e 100644
--- a/functions/popfunc/eeg_laplac.m
+++ b/functions/popfunc/eeg_laplac.m
@@ -103,7 +103,7 @@
 %---------------------------------------------------------------------
 % internal functions g2 gm_1
 function valor=g2(r1,r2,r3,r4)
-%g entrega la sumatoria hasta 7 con m=3 de la ec. (2) de la tecnical note del
+%g entrega la sumatoria hasta 7 con m=3 de la ec. (2) de la technical note del
 %Electrical Geodesics, salvo que esta esta bien
 
 x=cos(r1).*cos(r2).*cos(r3).*cos(r4) + ...
@@ -134,7 +134,7 @@
     valor = 1/4/pi*sum(aux2.*P);
  %-----------------------------------------------------   
 function valor=gm_1(r1,r2,r3,r4)
-%g entrega la sumatoria hasta 7 con m=3 de la ec. (2) de la tecnical note del
+%g entrega la sumatoria hasta 7 con m=3 de la ec. (2) de la technical note del
 %Electrical Geodesics
 
 x=cos(r1).*cos(r2).*cos(r3).*cos(r4) + ...
diff --git a/functions/popfunc/pop_biosig.m b/functions/popfunc/pop_biosig.m
index cfbd95750..5d9fe14b2 100644
--- a/functions/popfunc/pop_biosig.m
+++ b/functions/popfunc/pop_biosig.m
@@ -18,7 +18,7 @@
 %  'overflow'    - ['on'|'off'] overflow detection. Default is 'off'
 %  'uncalibrated' - ['on'|'off'] import uncalibrated data. Default is 'off'
 %  'blockepoch'  - ['on'|'off'] force importing continuous data. Default is 'on'
-%  'bdfeventmode' - [interger] see bdf2biosig_events function help. Default is 4.
+%  'bdfeventmode' - [integer] see bdf2biosig_events function help. Default is 4.
 %  'ref'         - [integer] channel index or index(s) for the reference.
 %                  Reference channels are not removed from the data,
 %                  allowing easy re-referencing. If more than one
@@ -147,7 +147,7 @@
                  { 'style' 'checkbox' 'String' 'Import as memory mapped file (use if out of memory error)' 'value' option_memmapdata 'tag' 'memorymapped' } ...
                  { 'style' 'checkbox' 'String' 'Enable overflow detection' 'value' 0 'tag' 'overflow' } ...
                  { 'style' 'checkbox' 'String' 'Import uncalibrated data' 'value' 0 'tag' 'uncalibrated' } ...
-                 { 'style' 'checkbox' 'string' 'Import EDF+ anotations (try also mexSLOAD below)' 'value' 1 'enable' 'on' 'tag' 'importannot' } ...
+                 { 'style' 'checkbox' 'string' 'Import EDF+ annotations (try also mexSLOAD below)' 'value' 1 'enable' 'on' 'tag' 'importannot' } ...
                  { 'style' 'checkbox' 'string' 'Import using alternative BIOSIG method (mexSLOAD)' 'value' 0 'callback' checkmex 'tag' 'importmex' } };
     geom = { [3 1] [3 1] [3 1] [3 1] [1] [1] [1] [1] [1] [1] [1] };
 
diff --git a/functions/popfunc/pop_chanevent.m b/functions/popfunc/pop_chanevent.m
index 13ab5b24f..9ed3de82c 100644
--- a/functions/popfunc/pop_chanevent.m
+++ b/functions/popfunc/pop_chanevent.m
@@ -51,7 +51,7 @@
 %                    at each time point (returning 1 if the data channel value 
 %                    is larger than 3, and 0 otherwise). You may also use 
 %                    any function (Ex: 'myfunction(X)'). If an equal character
-%                    is not present, this function preprend 'X=' to your 
+%                    is not present, this function prepend 'X=' to your 
 %                    expression before evaluating it. Otherwise it just 
 %                    evaluate the expression. For example, one may use
 %                    'if X(1)>100,X(1)=0;end; X(find(X>100))=X(find(X>100)-1);'
diff --git a/functions/popfunc/pop_fusechanrej.m b/functions/popfunc/pop_fusechanrej.m
index ae711ff5c..45deaa268 100644
--- a/functions/popfunc/pop_fusechanrej.m
+++ b/functions/popfunc/pop_fusechanrej.m
@@ -80,17 +80,17 @@
     if length(dats{index}) == 1
         TMPALLEEG = EEG(dats{index}(1));
         if ~isempty(TMPALLEEG(1).session)
-            fprintf('Skipping selecting common channels accross datasets for subject %s session %s (only 1 dataset)\n', TMPALLEEG(1).subject, num2str(TMPALLEEG(1).session));
+            fprintf('Skipping selecting common channels across datasets for subject %s session %s (only 1 dataset)\n', TMPALLEEG(1).subject, num2str(TMPALLEEG(1).session));
         else
-            fprintf('Skipping selecting common channels accross datasets for subject %s (only 1 dataset)\n', TMPALLEEG(1).subject);
+            fprintf('Skipping selecting common channels across datasets for subject %s (only 1 dataset)\n', TMPALLEEG(1).subject);
         end
     else
         TMPALLEEG = EEG(dats{index});
         commonChans = myintersect(TMPALLEEG.chanlocs);
         if ~isempty(TMPALLEEG(1).session)
-            fprintf('Selecting common channels accross datasets for subject %s session %s ***************\n', TMPALLEEG(1).subject, num2str(TMPALLEEG(1).session));
+            fprintf('Selecting common channels across datasets for subject %s session %s ***************\n', TMPALLEEG(1).subject, num2str(TMPALLEEG(1).session));
         else
-            fprintf('Selecting common channels accross datasets for subject %s ***************\n', TMPALLEEG(1).subject);
+            fprintf('Selecting common channels across datasets for subject %s ***************\n', TMPALLEEG(1).subject);
         end
         for iSet = 1:length(TMPALLEEG)
             if length(commonChans) ~= TMPALLEEG(iSet).nbchan
diff --git a/functions/popfunc/pop_newtimef.m b/functions/popfunc/pop_newtimef.m
index 9d4b00dc3..6968cf040 100644
--- a/functions/popfunc/pop_newtimef.m
+++ b/functions/popfunc/pop_newtimef.m
@@ -16,7 +16,7 @@
 %   "Sub-epoch time limits" - [edit box] sub epochs may be extracted (note that
 %              this function aims at plotting data epochs not continuous data).
 %              You may select the new epoch limits in this edit box.
-%   "Use n time points" - [muliple choice list] this is the number of time
+%   "Use n time points" - [multiple choice list] this is the number of time
 %              points to use for the time-frequency decomposition. The more
 %              time points, the longer the time-frequency decomposition
 %              takes to compute.
@@ -25,7 +25,7 @@
 %              of limits, you may also enter a sequence of frequencies. For
 %              example to compute the time-frequency decomposition at all
 %              frequency between 5 and 50 hertz with 1 Hz increment, enter "1:50"
-%   "Use limits, padding n" - [muliple choice list] "using limits" means
+%   "Use limits, padding n" - [multiple choice list] "using limits" means
 %              to use the upper and lower limits in "Frequency limits" with
 %              a specific padding ratio (padratio argument of newtimef).
 %              The last option "use actual frequencies" forces newtimef to
@@ -35,7 +35,7 @@
 %              frequencies. Note that this is only relevant if you specify
 %              frequency limits (in case you specify actual frequencies,
 %              this parameter is ignored).
-%   "Use divisive baseline" - [muliple choice list] there are two types of
+%   "Use divisive baseline" - [multiple choice list] there are two types of
 %              baseline correction, additive (the baseline is subtracted)
 %              or divisive (the data is divided by the baseline values).
 %              The choice is yours. There is also the option to perform 
diff --git a/functions/popfunc/pop_rejchan.m b/functions/popfunc/pop_rejchan.m
index 0dee47db5..f5b076083 100644
--- a/functions/popfunc/pop_rejchan.m
+++ b/functions/popfunc/pop_rejchan.m
@@ -178,7 +178,7 @@
     indFiltered = find(mean(measure) < -20);
     if ~isempty(indFiltered) && indFiltered(1) > 11, 
         measure = measure(:,1:indFiltered(1)-10); 
-        myfprintf('Removing spectrum data below -20dB (most likelly filtered out)\n'); 
+        myfprintf('Removing spectrum data below -20dB (most likely filtered out)\n'); 
     end
     meanSpec = mean(measure);
     stdSpec  = std( measure);
diff --git a/functions/popfunc/pop_reref.m b/functions/popfunc/pop_reref.m
index 3cfc4cbaa..0c400da67 100644
--- a/functions/popfunc/pop_reref.m
+++ b/functions/popfunc/pop_reref.m
@@ -414,7 +414,7 @@
             newLabelPos = strmatch(oldLabel, { newchanlocs.labels }, 'exact');
             
             if length(newLabelPos) > 1
-                warning('More than one match for specified reference channel; First one selected. This may cause eratic behavior. If the 2 channels are identical, delete one of them.');
+                warning('More than one match for specified reference channel; First one selected. This may cause erratic behavior. If the 2 channels are identical, delete one of them.');
             end
             if ~isempty( newLabelPos )
                 icachansind(i) = newLabelPos(1);
diff --git a/functions/sigprocfunc/mri3dplot.m b/functions/sigprocfunc/mri3dplot.m
index 3258679db..6eda0c8b0 100644
--- a/functions/sigprocfunc/mri3dplot.m
+++ b/functions/sigprocfunc/mri3dplot.m
@@ -36,7 +36,7 @@
 %                 Default is 'on'.
 %   'mixfact'   - [float] factor for mixing the background image with the
 %                 array3d information. Default is 0.5.
-%   'mixmode'   - ['add'|'overwrite'] 'add' will allow for trasnparency
+%   'mixmode'   - ['add'|'overwrite'] 'add' will allow for transparency
 %                 while 'overwrite' will preserve the original MRI image 
 %                 and overwrite the pixel showind density changes.
 %
@@ -406,7 +406,7 @@
     maxdens = max(prob3d(:));
     ncolors = size(g.cmap,1);
     
-    prob3d    = round((prob3d-g.cmin)/(g.cmax - g.cmin)*(ncolors-1))+1; % project desnity image into the color space: [1:ncolors]
+    prob3d    = round((prob3d-g.cmin)/(g.cmax - g.cmin)*(ncolors-1))+1; % project density image into the color space: [1:ncolors]
     prob3d( find(prob3d > ncolors) ) = ncolors;
     prob3d( find(prob3d < 1))        = 1; % added by Makoto
     newprob3d = zeros(size(prob3d,1), size(prob3d,2), size(prob3d,3), 3);
diff --git a/functions/sigprocfunc/plotsphere.m b/functions/sigprocfunc/plotsphere.m
index b2c1fb4da..ec29ff1a3 100644
--- a/functions/sigprocfunc/plotsphere.m
+++ b/functions/sigprocfunc/plotsphere.m
@@ -20,7 +20,7 @@
 %   'colormap' - [real] sphere colormap { default: jet }
 %
 % Output:
-%   handle    - sphere graphic handle(s). If projected sphere are ploted
+%   handle    - sphere graphic handle(s). If projected sphere are plotted
 %               handles of plotted projected spheres are also returned.
 %
 % Example:
diff --git a/functions/sigprocfunc/readlocs.m b/functions/sigprocfunc/readlocs.m
index be8bb5b7d..7d8f96821 100644
--- a/functions/sigprocfunc/readlocs.m
+++ b/functions/sigprocfunc/readlocs.m
@@ -36,7 +36,7 @@
 %                  coordinates, import either in EEGLAB format (nose pointing toward +X). 
 %                  This may not always be possible since EEGLAB might not be able to 
 %                  determine the nose direction for scanned electrode files. 'native' import
-%                  original carthesian coordinates (user can then specify the position of
+%                  original cartesian coordinates (user can then specify the position of
 %                  the nose when calling the topoplot() function; in EEGLAB the position
 %                  of the nose is stored in the EEG.chaninfo structure). {default 'eeglab'}
 %   'format'    -  [cell array] Format of a 'custom' channel location file (see above).
diff --git a/functions/sigprocfunc/readneurodat.m b/functions/sigprocfunc/readneurodat.m
index 0197d576b..9f60bc180 100644
--- a/functions/sigprocfunc/readneurodat.m
+++ b/functions/sigprocfunc/readneurodat.m
@@ -9,7 +9,7 @@
 % Outputs:
 %   CHANLOCS       - [structure] EEGLAB channel location data structure. See
 %                    help readlocs()
-%   labels         - [cell arrya] channel labels
+%   labels         - [cell array] channel labels
 %   theta          - [float array]array containing 3-D theta angle electrode
 %                    position (in degree)
 %   phi            - [float array]array containing 3-D phi angle electrode
diff --git a/functions/sigprocfunc/spectopo.m b/functions/sigprocfunc/spectopo.m
index 53e79e691..e17bed853 100644
--- a/functions/sigprocfunc/spectopo.m
+++ b/functions/sigprocfunc/spectopo.m
@@ -43,7 +43,7 @@
 %   'blckhn'   = [integer] Parameter to scale the windows length when Blackman-Harris 
 %                windows is used in 'wintype' {default: 2}
 %   'reref'    = ['averef'|'off'] convert data to average reference {default: 'off'}
-%   'mapnorm'  = [float vector] If 'data' contain the activity of an independant 
+%   'mapnorm'  = [float vector] If 'data' contain the activity of an independent 
 %                component, this parameter should contain its scalp map. In this case
 %                the spectrum amplitude will be scaled to component RMS scalp power.
 %                Useful for comparing component strengths {default: none}
@@ -72,7 +72,7 @@
 %   'icamode'  = ['normal'|'sub'] in 'sub' mode, instead of computing the spectra of
 %                individual ICA components, the function computes the spectrum of
 %                the data minus their contributions {default: 'normal'}
-%   'icamaps'  = [integer array] force plotting of selected ICA compoment maps 
+%   'icamaps'  = [integer array] force plotting of selected ICA component maps 
 %                {default: [] = the 'nicamaps' largest contributing components}.
 %   'icawinv'  = [float array] inverse component weight or mixing matrix. Normally,
 %                this is computed by inverting the ICA unmixing matrix 'weights' (above).
@@ -422,7 +422,7 @@
         if strcmp(g.memory, 'high') && strcmp(g.icamode, 'normal')
             myfprintf(g.verbose, 'Computing component spectra: ')
             [compeegspecdB freqs] = spectcomp( newweights*data(:,:), frames, srate, epoch_subset, g);
-        else % in case out of memory error, multiply conmponent sequencially
+        else % in case out of memory error, multiply conmponent sequentially
             if strcmp(g.icamode, 'sub') && ~isempty(g.plotchan) && g.plotchan == 0
                 % scan all electrodes
                 myfprintf(g.verbose, 'Computing component spectra at each channel: ')
@@ -695,7 +695,7 @@
 			[realpos(index) allaxuse] = closestplot( freqnormpos, allaxcoords, allaxuse );
 		end
 	
-		% put the channel plot a liitle bit higher
+		% put the channel plot a little bit higher
 		tmppos = get(headax(realpos(1)), 'position');
 		tmppos(2) = tmppos(2)+0.04;
 		set(headax(realpos(1)), 'position', tmppos);
@@ -918,7 +918,7 @@ function plotcolbar(g)
     
     if ~usepwelch
         myfprintf(g.verbose, '\nSignal processing toolbox (SPT) absent: spectrum computed using the pwelch()\n');
-        myfprintf(g.verbose, 'function from Octave which is suposedly 100%% compatible with the Matlab pwelch function\n');
+        myfprintf(g.verbose, 'function from Octave which is supposedly 100%% compatible with the Matlab pwelch function\n');
     end
     myfprintf(g.verbose,' (window length %d; fft length: %d; overlap %d):\n', winlength, fftlength, g.overlap);	
         
diff --git a/functions/sigprocfunc/transformcoords.m b/functions/sigprocfunc/transformcoords.m
index 8cf685923..48c0bbb74 100644
--- a/functions/sigprocfunc/transformcoords.m
+++ b/functions/sigprocfunc/transformcoords.m
@@ -100,7 +100,7 @@
         % pitch roll yaw rotation
         % -----------------------
         % pitch (x-axis); roll = y axis rotation; yaw = z axis
-        % see http://bishopw.loni.ucla.edu/AIR5/homogenous.html
+        % see http://bishopw.loni.ucla.edu/AIR5/homogeneous.html
         cp = cos(pitch); sp = sin(pitch);
         cr = cos(roll);  sr = sin(roll);
         cy = cos(yaw);   sy = sin(yaw);
diff --git a/functions/statistics/statcond.m b/functions/statistics/statcond.m
index 999faa8d1..0ae3d63eb 100644
--- a/functions/statistics/statcond.m
+++ b/functions/statistics/statcond.m
@@ -49,7 +49,7 @@
 %   'variance' = ['homegenous'|'inhomogenous'] this option is exclusively
 %                for parametric statistics using unpaired t-test. It allows
 %                to compute a more accurate value for the degree of freedom
-%                using the formula for inhomogenous variance (see
+%                using the formula for inhomogeneous variance (see
 %                ttest2_cell function). Default is 'inhomegenous'.
 %   'surrog'   = surrogate data array (see output).
 %   'stats'    = F- or T-value array (see output).
@@ -559,11 +559,11 @@
 
 % compute t-test
 % -------------------
-function [t, df] = ttest_cell_select( res, paired, homogenous)
+function [t, df] = ttest_cell_select( res, paired, homogeneous)
     if strcmpi(paired,'on')
         [t, df] = ttest_cell( res{1}, res{2});
     else
-        [t, df] = ttest2_cell( res{1}, res{2}, homogenous);
+        [t, df] = ttest2_cell( res{1}, res{2}, homogeneous);
     end
 
 % function to compute the number of dimensions
diff --git a/functions/statistics/ttest2_cell.m b/functions/statistics/ttest2_cell.m
index e22fc7433..443687ccd 100644
--- a/functions/statistics/ttest2_cell.m
+++ b/functions/statistics/ttest2_cell.m
@@ -10,8 +10,8 @@
 %   a,b       = data consisting of UNPAIRED arrays to be compared. The last 
 %               dimension of the data array is used to compute the t-test.
 %   'inhomogenous' = use computation for the degree of freedom using 
-%                    inhomogenous variance. By default the computation of
-%                    the degree of freedom is done with homogenous
+%                    inhomogeneous variance. By default the computation of
+%                    the degree of freedom is done with homogeneous
 %                    variances.
 %
 % Outputs:
@@ -36,7 +36,7 @@
 %
 % Author: Arnaud Delorme, SCCN/INC/UCSD, La Jolla, 2005
 %         (thank you to G. Rousselet for providing the formula for 
-%         inhomogenous variances).
+%         inhomogeneous variances).
 %
 % Reference:
 %   Schaum's outlines in statistics (3rd edition). 1999. Mc Graw-Hill.
@@ -76,19 +76,19 @@
         return;
     end
     
-    homogenous = 'homogenous';
+    homogeneous = 'homogeneous';
     if nargin > 1 && ischar(b)
-        homogenous = b;
+        homogeneous = b;
     end
     if nargin > 2 && ischar(c)
-        homogenous = c;
+        homogeneous = c;
     end
     if iscell(a), 
         b = a{2}; 
         a = a{1}; 
     end
-    if ~strcmpi(homogenous, 'inhomogenous') && ~strcmpi(homogenous, 'homogenous')
-        error('Value for homogenous parameter can only be ''homogenous'' or ''inhomogenous''');
+    if ~strcmpi(homogeneous, 'inhomogenous') && ~strcmpi(homogeneous, 'homogeneous')
+        error('Value for homogeneous parameter can only be ''homogeneous'' or ''inhomogenous''');
     end
 
     nd    = myndims(a);
@@ -97,8 +97,8 @@
     meana = mymean(a, nd);
     meanb = mymean(b, nd);
     
-    if strcmpi(homogenous, 'inhomogenous')
-        % inhomogenous variance from Howel, 2009, "Statistical Methods for Psychology"
+    if strcmpi(homogeneous, 'inhomogenous')
+        % inhomogeneous variance from Howel, 2009, "Statistical Methods for Psychology"
         % thank you to G. Rousselet for providing these formulas
         m  = meana - meanb;
         s1 = var(a,0,nd) ./ na;
diff --git a/functions/studyfunc/std_dipplot.m b/functions/studyfunc/std_dipplot.m
index 0452e8111..6b8962d3c 100644
--- a/functions/studyfunc/std_dipplot.m
+++ b/functions/studyfunc/std_dipplot.m
@@ -21,7 +21,7 @@
 %                {default: 'all'}.
 %   'mode'     - ['together'|'apart'|'multicolor'] Display all requested cluster on one 
 %                figure ('together') or separate figures ('apart'). 
-%                'together'-> plot all 'clusters' individuall in one multi-pane figure (without the gui).
+%                'together'-> plot all 'clusters' individually in one multi-pane figure (without the gui).
 %                'apart'   -> plot each cluster in a separate figure. 
 %                'multicolor' -> plot all clusters in one figure, 
 %                Note that this parameter has no effect if the 'comps' option (above) is used.
diff --git a/functions/studyfunc/std_plotdmat.m b/functions/studyfunc/std_plotdmat.m
index 39f06bd03..a6830cacb 100644
--- a/functions/studyfunc/std_plotdmat.m
+++ b/functions/studyfunc/std_plotdmat.m
@@ -243,7 +243,7 @@ function callback_popup_subject(src,eventdata,handles)
 check = find(sum(isnan(tmpdmat),2));
 tmpdmat(check,:) = [];
 
-% Checking checbox to sort/unsort
+% Checking checkbox to sort/unsort
 dmatsortindex = 1:size(tmpdmat, 1);
 if get(handles.popup_sort, 'Value') ~= 1
     [tmp,dmatsortindex] = sortrows(tmpdmat,get(handles.popup_sort, 'Value')-1);
diff --git a/functions/studyfunc/std_propplot.m b/functions/studyfunc/std_propplot.m
index 5fbefad08..2775d6787 100644
--- a/functions/studyfunc/std_propplot.m
+++ b/functions/studyfunc/std_propplot.m
@@ -2,7 +2,7 @@
 %                  properties for a STUDY set. 
 %                  Displays mean cluster scalp map, ERP, ERSP; 
 %                  dipole model, spectrum, and ITC in one figure 
-%                  per cluster. Only meaasures computed during 
+%                  per cluster. Only measures computed during 
 %                  pre-clustering (by pop_preclust() or std_preclust()) 
 %                  are plotted. Called by pop_clustedit().
 %                  Leaves the plotted grand mean cluster measures 
diff --git a/functions/studyfunc/std_rmdat.m b/functions/studyfunc/std_rmdat.m
index 8456d3686..372baadab 100644
--- a/functions/studyfunc/std_rmdat.m
+++ b/functions/studyfunc/std_rmdat.m
@@ -1,4 +1,4 @@
-% std_rmdat() - remvoe datasets from STUDY
+% std_rmdat() - remove datasets from STUDY
 %
 % Usage: 
 %   >> [STUDY,ALLEEG] = std_rmdat(STUDY, ALLEEG, 'key', val);  
diff --git a/functions/timefreqfunc/bootstat.m b/functions/timefreqfunc/bootstat.m
index 5d03b5ec6..1bce3aaf4 100644
--- a/functions/timefreqfunc/bootstat.m
+++ b/functions/timefreqfunc/bootstat.m
@@ -448,7 +448,7 @@
     end
 
 % this shuffling preserve the number of -1 and 1
-% for cloumns and rows (assuming matrix size is multiple of 2
+% for columns and rows (assuming matrix size is multiple of 2
 % -----------------------------------------------------------
 function array = supershuffle(array, dim)
     if size(array, 1) == 1 || size(array,2) == 1
diff --git a/functions/timefreqfunc/correctfit.m b/functions/timefreqfunc/correctfit.m
index ce78bb926..288acc666 100644
--- a/functions/timefreqfunc/correctfit.m
+++ b/functions/timefreqfunc/correctfit.m
@@ -10,7 +10,7 @@
 %
 % Optional inputs:
 %    'allpval'   - [float array] collection of p values drawn from random
-%                  distributions (theoritically uniform).
+%                  distributions (theoretically uniform).
 %    'gamparams' - [phat pci zerofreq] parameter for gamma function fitting.
 %                  zerofreq is the frequency of occurrence of p=0.
 %    'zeromode'  - ['on'|'off'] enable estimation of frequency of pval=0
diff --git a/functions/timefreqfunc/newtimef.m b/functions/timefreqfunc/newtimef.m
index 109e52d06..612502b14 100644
--- a/functions/timefreqfunc/newtimef.m
+++ b/functions/timefreqfunc/newtimef.m
@@ -111,7 +111,7 @@
 %       'wletmethod' = ['dftfilt'|'dftfilt2'|'dftfilt3'] Wavelet type to use.
 %                     'dftfilt2' -> Morlet-variant wavelets, or Hanning DFT.
 %                     'dftfilt3' -> Morlet wavelets.  See the timefreq() function 
-%                     for more detials {default: 'dftfilt3'}
+%                     for more details {default: 'dftfilt3'}
 %       'cycleinc'    ['linear'|'log'] mode of cycles increase when [min max] cycles 
 %                     are provided in 'cycle' parameter. Applies only to 
 %                     'wletmethod','dftfilt'  {default: 'linear'}
@@ -366,7 +366,7 @@
 % 02-28-00 added NOTE on formula derivation below -sm
 % 03-16-00 added axcopy() feature -sm && tpj
 % 04-16-00 added multiple marktimes loop -sm
-% 04-20-00 fixed ITC cbar limits when spcified in input -sm
+% 04-20-00 fixed ITC cbar limits when specified in input -sm
 % 07-29-00 changed frequencies displayed msg -sm
 % 10-12-00 fixed bug in freqs when cycles>0 -sm
 % 02-07-01 fixed inconsistency in BASE_BOOT use -sm
diff --git a/functions/timefreqfunc/timef.m b/functions/timefreqfunc/timef.m
index c11ad0529..a61fabf8d 100644
--- a/functions/timefreqfunc/timef.m
+++ b/functions/timefreqfunc/timef.m
@@ -181,7 +181,7 @@
 % 02-28-00 added NOTE on formula derivation below -sm
 % 03-16-00 added axcopy() feature -sm & tpj
 % 04-16-00 added multiple marktimes loop -sm
-% 04-20-00 fixed ITC cbar limits when spcified in input -sm
+% 04-20-00 fixed ITC cbar limits when specified in input -sm
 % 07-29-00 changed frequencies displayed msg -sm
 % 10-12-00 fixed bug in freqs when cycles>0 -sm
 % 02-07-01 fixed inconsistency in BASE_BOOT use -sm
diff --git a/sample_data/test_data/testbva.vhdr b/sample_data/test_data/testbva.vhdr
index d5b7f129e..ba613ab51 100644
--- a/sample_data/test_data/testbva.vhdr
+++ b/sample_data/test_data/testbva.vhdr
@@ -20,7 +20,7 @@ BinaryFormat=IEEE_FLOAT_32
 [Channel Infos]
 ; Each entry: Ch<Channel number>=<Name>,<Reference channel name>,
 ; <Resolution in microvolts>,<Future extensions..
-; Fields are delimited by commas, some fields might be omited (empty).
+; Fields are delimited by commas, some fields might be omitted (empty).
 ; Commas in channel names are coded as "\1".
 Ch1=Fp1,,
 Ch2=Fp2,,
diff --git a/sample_data/test_data/testbva.vmrk b/sample_data/test_data/testbva.vmrk
index f1e2193e8..9c15b7e7a 100644
--- a/sample_data/test_data/testbva.vmrk
+++ b/sample_data/test_data/testbva.vmrk
@@ -9,7 +9,7 @@ DataFile=testbva.dat
 ; Each entry: Mk<Marker number>=<Type>,<Description>,<Position in data points>,
 ; <Size in data points>, <Channel number (0 = marker is related to all channels)>,
 ; <Date (YYYYMMDDhhmmssuuuuuu)>
-; Fields are delimited by commas, some fields might be omited (empty).
+; Fields are delimited by commas, some fields might be omitted (empty).
 ; Commas in type or description text are coded as "\1".
 Mk1=New Segment,,1,1,0,20030924105038119829
 Mk2=Stimulus,S  4,108,1,0