예제 #1
0
    def __init__(self, name, a, d, mu_r, mu_v, delta_r, delta_v):
        self.name = name
        self.a = a
        self.d = d

        self.mu_r = mu_r
        self.mu_v = mu_v

        self.delta_r = delta_r
        self.delta_v = delta_v

        self.log_pi_r = self._get_log_mix_wts(delta_r)
        self.log_pi_v = self._get_log_mix_wts(delta_v)

        self._log_bin_norm_const = [u.log_bin_coeff(self.d[tp], self.a[tp]) for tp in arange(len(self.a))]
예제 #2
0
파일: data.py 프로젝트: 16NWallace/phylowgs
	def __init__(self, name, id, a, d, mu_r=0, mu_v=0):
		self.name = name # SSM name, blank for CNV
		self.id = id
		self.a = a
		self.d = d
		self.mu_r = mu_r # 1-p_error
		self.mu_v = mu_v
		self._log_bin_norm_const = [u.log_bin_coeff(self.d[tp], self.a[tp]) for tp in arange(len(self.a))]

		## all variables below are used by cnv related computations
		self.nr = 0
		self.nv = 0
		self.node = None # this is the node where the datum resides
		self.cnv = [] # for SSM, this is [(cnv,cp,cm)]
		
		self.tssb = None # this is just a pointer to tssb (tree object), gets initialized in evolve.py
예제 #3
0
파일: data.py 프로젝트: levvim/gdocker
	def __init__(self, name, id, a, d, mu_r=0, mu_v=0):
		self.name = name # SSM name, blank for CNV
		self.id = id
		self.a = a
		self.d = d
		self.mu_r = mu_r # 1-p_error
		self.mu_v = mu_v
		self._log_bin_norm_const = [u.log_bin_coeff(self.d[tp], self.a[tp]) for tp in arange(len(self.a))]

		## all variables below are used by cnv related computations
		self.nr = 0
		self.nv = 0
		self.node = None # this is the node where the datum resides
		self.cnv = [] # for SSM, this is [(cnv,cp,cm)]
		
		self.tssb = None # this is just a pointer to tssb (tree object), gets initialized in evolve.py
예제 #4
0
파일: data.py 프로젝트: Yeung678/phylosub
	def __init__(self, name,id,a, d, mu_r, mu_v, delta_r, delta_v):
		self.name = name; self.id=id;
		self.a=a
		self.d=d
			
		self.mu_r=mu_r
		self.mu_v=mu_v
		
		self.delta_r = delta_r
		self.delta_v = delta_v
		
		self.log_pi_r = self._get_log_mix_wts(delta_r)
		self.log_pi_v = self._get_log_mix_wts(delta_v)
		
		self._log_bin_norm_const = [u.log_bin_coeff(self.d[tp], self.a[tp]) for tp in arange(len(self.a))]		
		
		self.pairs = [(mu_r,mu_v,log_pi_r+log_pi_v) for mu_r, log_pi_r in zip(self.mu_r, self.log_pi_r) for mu_v, log_pi_v in zip(self.mu_v, self.log_pi_v)]
예제 #5
0
파일: data.py 프로젝트: svembu/phylosub
    def __init__(self, name, a, d, mu_r, mu_v, delta_r, delta_v):
        self.name = name
        self.a = a
        self.d = d

        self.mu_r = mu_r
        self.mu_v = mu_v

        self.delta_r = delta_r
        self.delta_v = delta_v

        self.log_pi_r = self._get_log_mix_wts(delta_r)
        self.log_pi_v = self._get_log_mix_wts(delta_v)

        self._log_bin_norm_const = [
            u.log_bin_coeff(self.d[tp], self.a[tp])
            for tp in arange(len(self.a))
        ]
예제 #6
0
파일: data.py 프로젝트: hrk2109/phylosub
    def __init__(self, name, id, a, d, mu_r, mu_v, delta_r, delta_v):
        self.name = name
        self.id = id
        self.a = a
        self.d = d

        self.mu_r = mu_r
        self.mu_v = mu_v

        self.delta_r = delta_r
        self.delta_v = delta_v

        self.log_pi_r = self._get_log_mix_wts(delta_r)
        self.log_pi_v = self._get_log_mix_wts(delta_v)

        self._log_bin_norm_const = [
            u.log_bin_coeff(self.d[tp], self.a[tp])
            for tp in arange(len(self.a))
        ]

        self.pairs = [(mu_r, mu_v, log_pi_r + log_pi_v)
                      for mu_r, log_pi_r in zip(self.mu_r, self.log_pi_r)
                      for mu_v, log_pi_v in zip(self.mu_v, self.log_pi_v)]