import sys sys.argv[1:2] = [str(x) for x in sys.argv[1:2]] output_folder = sys.argv[1] threshold_dir = (int)(sys.argv[2]) num_blocks = (int)(sys.argv[3]) compute_mse_current = (int)(sys.argv[4]) cycles = (int)(sys.argv[5]) max_active_set = (int)(sys.argv[6]) ####load data file = "ENEL_2014/Enel_dataset.npz" results = Result(file, "lasso") XTrain, YTrain, XTest, YTest = results.extract_train_test() enel_dict = results.extract_dict() Coord, Coord_turb, power_curve = results.extract_coords() angles_coord_turb, _ = compute_angle(Coord, Coord_turb) ##transformation of data X = np.concatenate((XTrain, XTest), axis=0) enel_transf = Enel_powerCurveTransformation() X_angle, _, _ = enel_transf.compute_angle_matrix(X) output_dict = dict.fromkeys(np.arange(0, 49), np.array([[]], dtype="int64")) k_levels = np.arange(0, 12).reshape([12, 1]) for key in np.arange(0, 49): current_values = np.arange(key * 12, key * 12 + 12).reshape([12, 1]) output_dict[key] = np.concatenate((current_values, k_levels), axis=1)
import numpy as np from LASSOModel import Shooting, LASSOEstimator from Lasso_utils import compute_weightedLASSO, compute_lasso from utility import get_current_data, assign_weights, assign_weights_ordered, \ get_beta_div_zeros, print_features_active import sys sys.argv[1:] = [str(x) for x in sys.argv[1:]] file_name = sys.argv[1] ext = ".npz" file = "ENEL_2014/"+file_name+ext results = Result(file, "lasso") dict_ = results.extract_dict() XTrain, YTrain, XVal, YVal = results.extract_train_val() score = "mean_squared_error" if score=="r2_score": score_f = r2_score scoring = "r2" else: score_f = mean_squared_error scoring = "mean_squared_error" verbose = True ###compute ranking
import sys sys.argv[1:] = [str(x) for x in sys.argv[1:]] file_name = sys.argv[1] file_name_dict = sys.argv[2] ext = ".npz" file = "ENEL_2014/"+file_name+ext file_dict = "ENEL_2014/"+file_name_dict+ext results_dict = Result(file_dict, "lasso") results = Result(file, "lasso") #dict_ = results.extract_dict() dict_ = results_dict.extract_dict() XTrain, YTrain, XVal, YVal = results.extract_train_val() new_loss, beta = compute_lasso(XTrain, YTrain, XVal, YVal,score = "mean_squared_error") print("loss lineare", new_loss) score = "mean_squared_error" if score=="r2_score": score_f = r2_score scoring = "r2" else: score_f = mean_squared_error scoring = "mean_squared_error" verbose = True
file_cross_val = folder+file_name+ext results_cross_val = Result(file_cross_val, "lasso") ##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) ##ranking verbose = True dict_ = results_cross_val.extract_dict() weights_data = results_cross_val.extract_weights() index_mse = len(weights_data)-1 weights_data = weights_data[index_mse] final_weights = np.zeros(len(dict_.keys())) keys_ = np.array((list)(dict_.keys())).astype("int64") for key in keys_: final_weights[key] += np.sum(weights_data[dict_.get(key).astype("int64")]) ordered_final_weights = np.argsort(final_weights)[::-1] if verbose: print("-------------")
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) if filename=="Enel_cross_val_blocks_direction_single_plus_verso.npz": print("trasformazione direzione verso") enel_transf = Enel_directionVersoPowerCurveTransformation()