sc_readtsvfile.m

function [X, genelist, celllist] = sc_readtsvfile(filename, genecolnum)
%Read TSV/TXT file
if nargin < 2
    genecolnum = 1;
end
if nargin < 1
    [filename, pathname] = uigetfile( ...
        {'*.csv;*.tsv;*.tab;*.txt', 'Expression Matrix Files (*.csv, *.tsv, *.tab, *.txt)'; ...
        '*.*', 'All Files (*.*)'}, ...
        'Pick a Exprssion Matrix file');
    if isequal(filename, 0), X = [];
        genelist = [];
        return;
    end
    filename = [pathname, filename];
end

if exist(filename, 'file') ~= 2
    error('FileNotFound');
end
%     try
if nargout > 2
    [X, genelist, celllist] = i_read_exprmat(filename, genecolnum);
    if size(X, 2) ~= length(celllist)
        celllist = strrep(celllist, '"', '');
        celllist(strlength(celllist) == 0) = [];
    end
    if size(X, 2) ~= length(celllist)
        warning('LENGTH(BARCODELIST) is not equal to SIZE(X,2)')
        if size(X, 2) - length(celllist) == -1
            % celllist(1:5)
            celllist = celllist(2:end);
            % warning('Correction was made.');
        end
    end
else
    [X, genelist] = i_read_exprmat(filename, genecolnum);
end
%     catch
%         fprintf('An Expression Matrix file looks like:\n');
%         fprintf('===================\n');
%         fprintf('Gene,c1,c2,c3,c4,c5\n');
%         fprintf('HES4,45,50,15,19,50\n');
%         fprintf('ISG15,168,279,312,425,180\n');
%         fprintf('AGRN,2,3,4,9,5\n');
%         fprintf('C1orf159,0,0,0,1,1\n');
%         fprintf('===================\n');
%         error('Error hapens when reading Expression Matrix file.');
%     end

end

function [X, genelist, sampleid] = i_read_exprmat(filename, genecolnum, verbose)
% Validate input args
narginchk(1, Inf);

% Get Filename
if ~ischar(filename) && ~(isstring(filename) && isscalar(filename))
    error('FileNameMustBeString');
end
filename = char(filename);

% Make sure file exists
if exist(filename, 'file') ~= 2
    error('FileNotFound');
end


if nargin < 2, genecolnum = 1; end
if nargin < 3, verbose = true; end
if verbose, fprintf('Reading %s ...... ', filename); end
warning('off', 'MATLAB:table:ModifiedAndSavedVarnames');

T = readtable(filename, 'filetype', 'text', 'HeaderLines', 1);
% T=readtable(filename,'filetype','text');

warning('on', 'MATLAB:table:ModifiedAndSavedVarnames');
X = table2array(T(:, (1 + genecolnum):end));
genelist = string(table2array(T(:, 1:genecolnum)));
if nargout > 2
    % sampleid=string(T.Properties.VariableNames(1+genecolnum:end)');
    fid = fopen(filename);
    sampleid = string(strsplit(fgetl(fid), {'\t', ','})');
    fclose(fid);
end
if verbose, fprintf('done.\n'); end
end