Mercurial > repos > vipints > rdiff

annotate rDiff/src/variance/estimate_variance.m @ 2:233c30f91d66

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
updated python based GFF parsing module which will handle GTF/GFF/GFF3 file types
author vipints <vipin@cbio.mskcc.org>
date Tue, 08 Oct 2013 07:15:44 -0400
parents 0f80a5141704
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
1 function [VARIANCE]=estimate_variance(MEANS,VARS)
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
2 %estimates the variance using the locfit package
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
3
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
4 IDX_ZERO=and(VARS>=0,MEANS>0);
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
5 VARIANCE =locfit(log(MEANS(IDX_ZERO)),VARS(IDX_ZERO),'family','gamma','acri','cp','deg',2,'maxk',1000);
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
6
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
7 return