Christiaan Swanepoel christiaanjs

## coalescent-softsort.ipynb

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                christiaanjs
                / coalescent-softsort.ipynb
            
            
              Created
              October 12, 2021 04:14
            
              
                Relaxing the coalescent with SoftSort
              
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## product-bijector.ipynb

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                christiaanjs
                / product-bijector.ipynb
            
            
              Last active
              July 24, 2020 04:22
            
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## tf-likelihood-test.ipynb

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                christiaanjs
                / tf-likelihood-test.ipynb
            
            
              Last active
              May 14, 2020 09:33
            
              
                Making Tensorflow likelihood gradients work in function mode
              
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## phylo-model.py
model = tfd.JointDistributionNamed(dict(
        frequencies=tfd.Dirichlet(concentration=[4,4,4,4]),
        kappa=tfd.LogNormal(loc=0, scale=1),
        pop_size=tfd.LogNormal(loc=0, scale=1),
        tree=lambda pop_size: treeflow.coalescent.ConstantCoalescent(pop_size=pop_size, sampling_times=sampling_times),
        sequences=lambda tree, kappa, frequencies: treeflow.sequences.PhyloCTMC(tree=tree, q=HKY(frequencies=frequencies, kappa=kappa))
    ))

## quantile-test.ipynb

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                christiaanjs
                / quantile-test.ipynb
            
            
              Created
              January 16, 2020 03:37
            
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## bijectors.ipynb

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                christiaanjs
                / bijectors.ipynb
            
            
              Created
              September 18, 2019 23:47
            
              
                Truncated Lognormal distribution in Tensorflow Probability
              
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## skygrid-test.ipynb

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                christiaanjs
                / skygrid-test.ipynb
            
            
              Last active
              July 27, 2019 00:55
            
              
                Skygrid + Beagle test
              
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## xml-sequence.py
from xml.etree import ElementTree

data_xml_filename = '311spp_5genes.xml'
xml_root = ElementTree.parse(data_xml_filename)
flat_dict = [{ 'alignment': alignment.attrib['id'], 'taxon': sequence.find('taxon').attrib['idref'], 'sequence': sequence.find('taxon').tail.strip() } for alignment in xml_root.findall('alignment') for sequence in alignment.findall('sequence') ]
taxa_dict = { taxon: { value['alignment']: value['sequence'] for value in values } for taxon, values in groupby(sorted(flat_dict, key=key_func), key=key_func) }
print(taxa_dict)

## operator-jacobian.ipynb

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                christiaanjs
                / operator-jacobian.ipynb
            
            
              Created
              December 14, 2018 04:04
            
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## distribution-test.ipynb

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                christiaanjs
                / distribution-test.ipynb
            
            
              Created
              November 22, 2018 19:53
            
              
                Bootstrap multivariate two-sample test for distribution
              
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
	model = tfd.JointDistributionNamed(dict(
	frequencies=tfd.Dirichlet(concentration=[4,4,4,4]),
	kappa=tfd.LogNormal(loc=0, scale=1),
	pop_size=tfd.LogNormal(loc=0, scale=1),
	tree=lambda pop_size: treeflow.coalescent.ConstantCoalescent(pop_size=pop_size, sampling_times=sampling_times),
	sequences=lambda tree, kappa, frequencies: treeflow.sequences.PhyloCTMC(tree=tree, q=HKY(frequencies=frequencies, kappa=kappa))
	))
	from xml.etree import ElementTree

	data_xml_filename = '311spp_5genes.xml'
	xml_root = ElementTree.parse(data_xml_filename)
	flat_dict = [{ 'alignment': alignment.attrib['id'], 'taxon': sequence.find('taxon').attrib['idref'], 'sequence': sequence.find('taxon').tail.strip() } for alignment in xml_root.findall('alignment') for sequence in alignment.findall('sequence') ]
	taxa_dict = { taxon: { value['alignment']: value['sequence'] for value in values } for taxon, values in groupby(sorted(flat_dict, key=key_func), key=key_func) }
	print(taxa_dict)