score = "mean_squared_error"
if score=="r2_score":
score_f = r2_score
scoring = "r2"
else:
score_f = mean_squared_error
scoring = "mean_squared_error"
ext = ".npz"
file_cross_val = file_name+ext
fine_name_weights = file_name+"ranking"+ext
results_cross_val = Result(file_cross_val, "lasso")
results_weighted_lasso = Result(fine_name_weights, "lasso")
mses = results_weighted_lasso.extract_mses()
mses_int = list(map(int, mses))
iter = np.argmin(mses_int)
print ("iter chosen:",iter, "with mse:",mses_int[iter])
print("--------------")
indexes_beta = results_weighted_lasso.extract_beta_div_zeros()[iter]
##get transformed data
XTrain, XTest = results_cross_val.extract_data_transf()
_,YTrain,_, YTest = results_cross_val.extract_train_test()
### centratura dei dati
XTrain, YTrain, X_mean, y_mean, X_std = center_data(XTrain, YTrain, fit_intercept=True, normalize = True)
XTest, YTest = center_test(XTest,YTest,X_mean,y_mean,X_std)
import sys
from Enel_utils import find_nearest_turbine, compute_angle
from ExtractResult import Result
from Transformation import Enel_powerCurveTransformation, EnelWindSpeedTransformation, \
Enel_directionVersoPowerCurveTransformation, Enel_directionPowerCurveTransformation
sys.argv[1:3] = [str(x) for x in sys.argv[1:]]
folder_name = sys.argv[1]
filename = sys.argv[2]
compute_dict = (int)(sys.argv[3])
tot_filename = folder_name+"/"+filename
results = Result(tot_filename, "lasso")
weights_list = results.extract_weights()
mses = results.extract_mses()
XTrain_transf, _ = results.extract_data_transf()
XTrain_, YTrain_, XVal_, YVal_ = results.extract_train_val()
saved_indexes_list = results.get_saved_indexes()
XTrain_ValNoCenter, YTrainVal_noCenter, XVal_noCenter,YVal_noCenter = results.extract_train_val_no_centered()
output_dict = results.extract_dict()
file = "ENEL_2014/Enel_dataset.npz"
results = Result(file, "lasso")
XTrain, YTrain, XTest, YTest = results.extract_train_test()
X_speed,_ = EnelWindSpeedTransformation().transform(XTest)
Coord, Coord_turb, power_curve = results.extract_coords()
angles_coord_turb,verso_turb_point = compute_angle(Coord, Coord_turb)
enel_transf = Enel_powerCurveTransformation()
X_angle,x_verso,_ = enel_transf.compute_angle_matrix(XTest)
