function [mu_hat sigma_hat] = estimBrownian(data)

mu_hat = sum(log(data))/length(data);

sigma_hat = sqrt(sum((log(data) - mu_hat).^2)/length(data));