def fit_mammals():
model = LMMBasicAnimal(ylist = ['Body'],plist=['Brain'])
dataFile = "data/bzdata.txt"
treeFile = "data/bz.trees.newick"
# loadFromFile(self,dataFile,treeFile,taxa_column='taxa',sep='TAB'):
model.loadFromFile(dataFile,treeFile,taxa_column='Name',sep='SPACE')
model.fitGLS()
# model.fitLMMBasic()
print "--------------------"
model.fitGLSLambda()
def fit_test():
print "Testing GLS"
model = LMMBasicAnimal(ylist = ['length'],plist=['size'])
dataFile = "data/test.data"
treeFile = "data/test.newick"
# loadFromFile(self,dataFile,treeFile,taxa_column='taxa',sep='TAB'):
model.loadFromFile(dataFile,treeFile,sep='SPACE')
model.fitGLS()
# model.fitLMMBasic()
print "--------------------"
