HyPhy message board - constrain variables in HyPhy

	Welcome, Guest. Please Login

Home

Help

HyPhy message board › HYPHY Package › HyPhy feedback › constrain variables in HyPhy

(Moderators: Sergei, Simon)

‹ Previous Topic | Next Topic ›

Pages: 1

Send Topic

constrain variables in HyPhy (Read 2094 times)

student_of_evolution

YaBB Newbies

Offline

Curious HyPhy user

Posts: 1

constrain variables in HyPhy
Sep 9^th, 2013 at 1:37pm

Hi,

I am trying to devise a LRT (or perhaps use the AIC difference) where the null model is neutral evolution and the alternate model is selection.

For example Yang's M1a vs M2a does this (thus allows use of chisquare with df = 2)
M1a (neutral):
parameters: p0, w_0 < 1, w_1 = 1
M2a (selection):
parameters: p0, p1, w_0 < 1, w_1 = 1, w_2 > 1

I was thinking about setting the neutral model (using the general discrete bi variate distribution) to be:
beta_s <= alpha_s for all s
And alternate model is:
beta_s > alpha_s for at least one s

Can you constrain these variables in hyphy?

IP Logged

Sergei

YaBB Administrator

Offline

Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male

Re: constrain variables in HyPhy
Reply #1 - Sep 9^th, 2013 at 2:22pm

Hi there,

The types of constraints you need can be implemented through ratios, e.g.

Code:

global R_0;

beta_0 := alpha_0 * R_0;

R_0 :< 1;

...

Alternatively, you can use the ratio parameterization directly (e.g. see Multimedia File Viewing and Clickable Links are available for Registered Members only!! You need to

Associate Professor
Division of Infectious Diseases
Division of Biomedical Informatics
School of Medicine
University of California San Diego

WWW

IP Logged

Pages: 1

Send Topic

‹ Previous Topic | Next Topic ›

« Home

‹ Board

Top of this page