|
|
|
@ -1,5 +1,6 @@ |
|
|
|
|
import os |
|
|
|
|
import time as tm |
|
|
|
|
from line_profiler import LineProfiler |
|
|
|
|
|
|
|
|
|
import network_graph as ng |
|
|
|
|
import sample_path as sp |
|
|
|
@ -20,22 +21,22 @@ class ParametersEstimator: |
|
|
|
|
self.net_graph.get_ordered_by_indx_parents_values_for_all_nodes()) |
|
|
|
|
|
|
|
|
|
def parameters_estimation(self): |
|
|
|
|
print("Starting computing") |
|
|
|
|
t0 = tm.time() |
|
|
|
|
#print("Starting computing") |
|
|
|
|
#t0 = tm.time() |
|
|
|
|
for indx, trajectory in enumerate(self.sample_path.trajectories): |
|
|
|
|
self.parameters_estimation_single_trajectory(trajectory.get_trajectory()) |
|
|
|
|
#print("Finished Trajectory number", indx) |
|
|
|
|
t1 = tm.time() - t0 |
|
|
|
|
print("Elapsed Time ", t1) |
|
|
|
|
#t1 = tm.time() - t0 |
|
|
|
|
#print("Elapsed Time ", t1) |
|
|
|
|
|
|
|
|
|
def parameters_estimation_single_trajectory(self, trajectory): |
|
|
|
|
#print(type(trajectory[0][0])) |
|
|
|
|
t0 = tm.time() |
|
|
|
|
for indx, row in enumerate(trajectory): |
|
|
|
|
if trajectory[indx][1] == -1: |
|
|
|
|
break |
|
|
|
|
if trajectory[indx + 1][1] != -1: |
|
|
|
|
transition = self.find_transition(trajectory[indx], trajectory[indx + 1]) |
|
|
|
|
which_node = self.net_graph.get_node_by_index(transition[0]) |
|
|
|
|
which_node = transition[0] |
|
|
|
|
# print(which_node) |
|
|
|
|
which_matrix = self.which_matrix_to_update(row, transition[0]) |
|
|
|
|
which_element = transition[1] |
|
|
|
@ -43,17 +44,22 @@ class ParametersEstimator: |
|
|
|
|
|
|
|
|
|
#changed_node = which_node |
|
|
|
|
time = self.compute_time_delta(trajectory[indx], trajectory[indx + 1]) |
|
|
|
|
which_element = transition[1][0] |
|
|
|
|
self.amalgamated_cims_struct.update_state_residence_time_for_matrix(which_node, which_matrix, which_element, |
|
|
|
|
time) |
|
|
|
|
|
|
|
|
|
for node_indx, node in enumerate(self.net_graph.get_nodes()): |
|
|
|
|
#if node != changed_node: |
|
|
|
|
for node_indx in range(0, 3): |
|
|
|
|
if node_indx != transition[0]: |
|
|
|
|
# print(node) |
|
|
|
|
which_node = node |
|
|
|
|
which_node = node_indx |
|
|
|
|
which_matrix = self.which_matrix_to_update(row, node_indx) |
|
|
|
|
which_element = row[node_indx + 1] |
|
|
|
|
# print("State res time element " + str(which_element) + node) |
|
|
|
|
# print("State res time matrix indx" + str(which_matrix)) |
|
|
|
|
self.amalgamated_cims_struct.update_state_residence_time_for_matrix(which_node, which_matrix, which_element, time) |
|
|
|
|
|
|
|
|
|
self.amalgamated_cims_struct.update_state_residence_time_for_matrix(which_node, which_matrix, |
|
|
|
|
which_element, time) |
|
|
|
|
t1 = tm.time() - t0 |
|
|
|
|
print("Elapsed Time ", t1) |
|
|
|
|
|
|
|
|
|
def find_transition(self, current_row, next_row): |
|
|
|
|
for indx in range(1, len(current_row)): |
|
|
|
@ -85,14 +91,18 @@ g1.init_graph() |
|
|
|
|
|
|
|
|
|
pe = ParametersEstimator(s1, g1) |
|
|
|
|
pe.init_amalgamated_cims_struct() |
|
|
|
|
print(pe.amalgamated_cims_struct.get_set_of_cims('X').get_cims_number()) |
|
|
|
|
print(pe.amalgamated_cims_struct.get_set_of_cims('Y').get_cims_number()) |
|
|
|
|
print(pe.amalgamated_cims_struct.get_set_of_cims('Z').get_cims_number()) |
|
|
|
|
print(pe.amalgamated_cims_struct.get_set_of_cims(0).get_cims_number()) |
|
|
|
|
print(pe.amalgamated_cims_struct.get_set_of_cims(1).get_cims_number()) |
|
|
|
|
print(pe.amalgamated_cims_struct.get_set_of_cims(2).get_cims_number()) |
|
|
|
|
#pe.parameters_estimation_single_trajectory(pe.sample_path.trajectories[0].get_trajectory()) |
|
|
|
|
pe.parameters_estimation() |
|
|
|
|
for matrix in pe.amalgamated_cims_struct.get_set_of_cims('Y').actual_cims: |
|
|
|
|
lp = LineProfiler() |
|
|
|
|
lp_wrapper = lp(pe.parameters_estimation_single_trajectory) |
|
|
|
|
lp_wrapper(pe.sample_path.trajectories[0].get_trajectory()) |
|
|
|
|
lp.print_stats() |
|
|
|
|
#pe.parameters_estimation() |
|
|
|
|
"""for matrix in pe.amalgamated_cims_struct.get_set_of_cims(1).actual_cims: |
|
|
|
|
print(matrix.state_residence_times) |
|
|
|
|
print(matrix.state_transition_matrix) |
|
|
|
|
matrix.compute_cim_coefficients() |
|
|
|
|
print(matrix.cim) |
|
|
|
|
print(matrix.cim)""" |
|
|
|
|
|
|
|
|
|