function [result] = MonteCarlo(data, indices, title)
    data = data(indices);    
    data = data(~isnan(data));                
    mu = mean(data);
    s = std(data);
    % 1000 random trials of same length as data
    % to get an idea of the error in the predicted mean.
    means = mean(normrnd(mu, s, length(data), 1000)); 
    sigma = std(means);
    %fprintf('%-20s mean %8.4f , std %8.4f , (s %8.4f)\n', title, mu, sigma, std(data));
    result = [mu sigma];
end