Mercurial > repos > vipints > rdiff

view rDiff/src/variance/estimate_variance.m @ 3:29a698dc5c7e default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge multiple heads.
author Dave Bouvier <dave@bx.psu.edu>
date Mon, 27 Jan 2014 14:15:36 -0500
parents 0f80a5141704
children
line wrap: on
line source

function [VARIANCE]=estimate_variance(MEANS,VARS) 
%estimates the variance using the locfit package

IDX_ZERO=and(VARS>=0,MEANS>0);
VARIANCE =locfit(log(MEANS(IDX_ZERO)),VARS(IDX_ZERO),'family','gamma','acri','cp','deg',2,'maxk',1000);

return