ankurankan on dev
Updates doc copyright; Also fix… (compare)
ankurankan on dev
Adds CI tests documentation to … Replace the independence symbol… (compare)
def bic(train,test,name,folder,resultlist,address):
array=['Activity']
trainstart=time()
bic=BicScore(train)
hc=HillClimbSearch(train, scoring_method=bic)
best_model=hc.estimate()
edges=best_model.edges()
model=BayesianModel(edges)
model.fit(train,estimator=BayesianEstimator, prior_type="BDeu")
trainend=time()-trainstart
print "nodes", model.nodes()
print "edges", model.edges()
flag=0
if(set(model.nodes())-set(array) ==set(model.nodes())):
flag=1
elif(set(model.nodes())-set(array) == set(test.columns)):
teststart=time()
print "before"
draw(model.edges(),name,"bic",folder)
result=model.predict(test).values.ravel()
print result
testend=time()-teststart
pred=list(result)
print pred
else:
indicator=list(set(test.columns)-set(model.nodes()))
print indicator
testchange=test.copy()
for f in range(len(indicator)):
print f
del testchange[indicator[f]]
print testchange
teststart=time()
result=model.predict(testchange).values.ravel()
testend=time()-teststart
pred=list(result)
model_data = XMLBIFWriter(model)
model_data.write_xmlbif(address+name+'_bic.bif')
if(flag==0):
fscore,accuracy,precision,recall=calscore(resultlist,pred)
draw(model.edges(),name,"bic",folder)
WriteData(address,name+".xlsx",name,pred)
else:
fscore=accuracy=precision=recall=trainend=testend=0
print "fscore :",fscore
print "accuracy :",accuracy
print "precision :",precision
print "recall :",recall
plt.scatter(pred , resultlist)
return fscore,accuracy,precision,recall,trainend,testend
secondtestindicexls1=pd.ExcelFile("D:\OrdonezB_Sensors.xls")
secondtestindicedf1=secondtestindicexls1.parse("OSB")
secondtestindicexls2=pd.ExcelFile("D:\OrdonezA_Sensors.xls")
secondtestindicedf2=secondtestindicexls2.parse("OSA")
secondtestindicexls3=pd.ExcelFile("D:\label.xlsx")
secondtestindicedf3=secondtestindicexls3.parse("Sheet1")
bic(secondtestindicedf1,secondtestindicedf2,"res11","D:\result of run python\",secondtestindicedf3,"D:\result of run python\")
Hello there, I'm trying to work out the ProbModelXML file format by generating a file from a working model I created using the method said on the notebook tutorial (https://github.com/pgmpy/pgmpy_notebook/blob/master/notebooks/8.%20Reading%20and%20Writing%20from%20pgmpy%20file%20formats.ipynb)
model_data = get_probmodel_data(valve_model);
writer = ProbModelXMLWriter(model_data=model_data);
print(writer);
and it give me this error
Traceback (most recent call last):
File "/Users/-/Desktop/Valve.py", line 49, in <module>
model_data = get_probmodel_data(valve_model);
File "/Users/-/anaconda3/lib/python3.5/site-packages/pgmpy/readwrite/ProbModelXML.py", line 266, in get_probmodel_data
model_data['probnet']['edges'][str(edge)] = model.edge[edge[0]][edge[1]]
AttributeError: 'BayesianModel' object has no attribute 'edge'
I tried to fix it by modify line 266 in ProbModelXML.py to
model_data['probnet']['edges'][str(edge)] = model.edges[edge[0]][edge[1]]
and now I'm receiving this error:
Traceback (most recent call last):
File "/Users/-/Desktop/Valve.py", line 49, in <module>
model_data = get_probmodel_data(valve_model);
File "/Users/-/anaconda3/lib/python3.5/site-packages/pgmpy/readwrite/ProbModelXML.py", line 266, in get_probmodel_data
model_data['probnet']['edges'][str(edge)] = model.edges[edge[0]][edge[1]]
File "/Users/-/anaconda3/lib/python3.5/site-packages/networkx/classes/reportviews.py", line 930, in getitem
return self._adjdict[u][v]
KeyError: 'n'
It seems like there's an error associated with the indexing (model.edges[edge[0]][edge[1]]), but I'm not sure how to fix it. Would appreciate if anyone could help.
+-------------------------------------+-----------+
| project_type(Fast Track Onboarding) | 0.0299222 |
+-------------------------------------+-----------+
| project_type(Innovation) | 0.0113704 |
+-------------------------------------+-----------+
| project_type(governance) | 0.0388989 |
+-------------------------------------+-----------+
| project_type(innovation) | 0.032316 |
+-------------------------------------+-----------+
| project_type(other) | 0.831837 |
+-------------------------------------+-----------+
| project_type(performance) | 0.0359066 |
+-------------------------------------+-----------+
| project_type(productivity) | 0.0197487 |
+-------------------------------------+-----------+
In this example I want to get Fast Track Onboarding, Innovation...