diff options
author | elavington <elavington@hotmail.com> | 2017-07-27 15:33:17 -0400 |
---|---|---|
committer | GitHub <noreply@github.com> | 2017-07-27 15:33:17 -0400 |
commit | 228f8f203eac1b5881d890e266ac10d46bb1b024 (patch) | |
tree | 6ee415846c8889e94480676c991581b92fa44de6 /cgi-bin | |
parent | 3e5406bbd3e05ac59c1cc954de7bf8187baf9a39 (diff) | |
download | cressdna-backup.tar.gz cressdna-backup.zip |
Add files via uploadbackup
Diffstat (limited to 'cgi-bin')
-rw-r--r-- | cgi-bin/SVM_linear_aa_clf.pkl | bin | 0 -> 187597 bytes | |||
-rw-r--r-- | cgi-bin/UniqRepsGemys_6089_StSCALER.pkl | bin | 0 -> 980 bytes | |||
-rw-r--r-- | cgi-bin/classifier.py | 52 |
3 files changed, 52 insertions, 0 deletions
diff --git a/cgi-bin/SVM_linear_aa_clf.pkl b/cgi-bin/SVM_linear_aa_clf.pkl Binary files differnew file mode 100644 index 0000000..1afce0a --- /dev/null +++ b/cgi-bin/SVM_linear_aa_clf.pkl diff --git a/cgi-bin/UniqRepsGemys_6089_StSCALER.pkl b/cgi-bin/UniqRepsGemys_6089_StSCALER.pkl Binary files differnew file mode 100644 index 0000000..3a098bd --- /dev/null +++ b/cgi-bin/UniqRepsGemys_6089_StSCALER.pkl diff --git a/cgi-bin/classifier.py b/cgi-bin/classifier.py new file mode 100644 index 0000000..ec2b634 --- /dev/null +++ b/cgi-bin/classifier.py @@ -0,0 +1,52 @@ +#!/usr/bin/python + +#import packages to be used +from sklearn.svm import SVC +from sklearn.feature_extraction.text import CountVectorizer +from sklearn.preprocessing import StandardScaler +from sklearn.externals import joblib +import cgi, cgitb + +cgitb.enable() +form=cgi.FieldStorage() +if form.getvalue('fasta'): + alignment = form.getvalue('fasta') + alignment=[alignment] + name=form.getvalue('seqname') + size=len(alignment[0]) +else: + alignment = ["MPSKKSGPQPHKRWVFTLNNPSEEEKNKIRELPISLFDYFVCGEEGLEEGRTAHLQGFANFAKKQTFNKVKWYFGARCHIEKAKGTDQQNKEYCSKEGHILIECGAPRNQGKRSDLSTAYFDYQQSGPPGMVLLNCCPSCRSSLSEDYYFAILEDCWRTINGGTRRPI"] + name='demo' + size=len(alignment[0]) + +html = open("./www.html/CRESSresults.html") +page=html.read() + + +AAs=['a','c','d','e','f','g','h','i','k','l','m','n','p','q','r','s','t','v','w','y'] +clf=joblib.load("./cgi-bin/SVM_linear_aa_clf.pkl") +StSc=joblib.load("./cgi-bin/UniqRepsGemys_6089_StSCALER.pkl") +cv=CountVectorizer(analyzer='char',ngram_range=(1,1),vocabulary=AAs) + + +#initialize text data vectorizer + +dataVect=cv.transform(alignment) + +#Scale the data to the training set +X=StSc.transform(dataVect.astype("float64")) + +#make predictions for the original dataset +results=",".join([name,clf.predict(X)[0]]) +results=",".join([results,str(size)]) +#for i in results: + #print(i[0],"\t",i[1]) + +output = page.format(prediction=results) +"""f=open('test.html','w') +f.write(output) +f.close()""" +print (output) + + +quit()
\ No newline at end of file |