Skip to content

Commit

Permalink
removed some redundent scripts
Browse files Browse the repository at this point in the history
  • Loading branch information
jeffhsu3 committed Jul 15, 2015
1 parent 0413a1d commit 33ce010
Show file tree
Hide file tree
Showing 6 changed files with 14 additions and 226 deletions.
17 changes: 13 additions & 4 deletions genda/formats/PED.py
Original file line number Diff line number Diff line change
@@ -1,31 +1,40 @@
"""
"""
import Genotype
import pandas as pd
import numpy as np

class PED(Genotype.Genotype):

def __init__(self, PED, MAP, encoder = None):
self.PED = pd.read_table(PED, delim_whitespace = True, header = None)
self.MAP = pd.read_table(MAP, delim_whitespace = True, header = None)
self.PED = pd.read_table(PED,
delim_whitespace = True, header = None)
self.MAP = pd.read_table(MAP,
delim_whitespace = True, header = None)
if encoder != None:
self.geno = self.genotype(self.PED, self.MAP, encoder)
else:
self.geno = None


def genotype(self, PED, MAP, encoder):
geno = PED.ix[:, 6:]
snps = []
for s in MAP.ix[:,1]:
snps.extend([s]*2)
geno.columns = snps
geno = geno.apply(self.apply_encoder, encoder = encoder)
geno = pd.DataFrame([[geno.ix[r,c*2] + geno.ix[r,(c*2)+1]\
for c in range(geno.shape[1]/2)] for r in range(geno.shape[0])])
geno = pd.DataFrame([[geno.ix[r,c*2] +\
geno.ix[r,(c*2)+1]\
for c in range(geno.shape[1]/2)] for r\
in range(geno.shape[0])])

geno = geno.transpose()
geno.index = MAP.ix[:,1]
geno.columns = PED.ix[:,1]
return geno


def apply_encoder(self, PED, encoder):
snp = encoder[PED.name]
encoding_dict = {snp[0] : 0, snp[2] : 1, '0': np.nan}
Expand Down
2 changes: 1 addition & 1 deletion genda/stats/chipSeqnorm.py
Original file line number Diff line number Diff line change
Expand Up @@ -17,7 +17,7 @@ def window_norm(y, window_size):
"""

# Handles the edge cases
pass
raise NotImplementedError


def savitzky_golay(y, window_size, order, deriv=0):
Expand Down
63 changes: 0 additions & 63 deletions scripts/getSeq.py

This file was deleted.

38 changes: 0 additions & 38 deletions scripts/getSequences.py

This file was deleted.

37 changes: 0 additions & 37 deletions scripts/tconsView.py

This file was deleted.

83 changes: 0 additions & 83 deletions scripts/vcf_protin.py

This file was deleted.

0 comments on commit 33ce010

Please sign in to comment.