myawesomepkg 0.1.8__py3-none-any.whl

myawesomepkg/TSAPY1/Practical No 4 B.py

@@ -0,0 +1,137 @@
+PRCTICAL 4 B
+
+
+
+def test_stationarity(dataFrame, var):
+  dataFrame['rollMean']=dataFrame[var].rolling(window=12).mean()
+  dataFrame['rollStd']=dataFrame[var].rolling(window=12).std()
+
+  from statsmodels.tsa.stattools import adfuller
+  import seaborn as sns
+  adfTest = adfuller(dataFrame[var], autolag='AIC')
+  stats=pd.Series(adfTest[0:4],index=['Test Statistic','p-value', '#lags used', 'number of observations used'])
+  print(stats)
+
+  for key, value in adfTest[4].items():
+    print('\t%s: %.3f' % (key, value))
+
+  sns.lineplot(data=dataFrame, x=dataFrame.index, y=var)
+  sns.lineplot(data=dataFrame, x=dataFrame.index, y='rollMean')
+  sns.lineplot(data=dataFrame, x=dataFrame.index, y='rollStd')
+
+
+********
+
+
+
+import pandas as pd
+import numpy as np
+
+#Reading the airline-passengers data
+
+data = pd.read_csv('/content/drive/MyDrive/MScDS TSA/AirPassengers.csv', index_col='Month')
+
+#Checking for some values of the data.
+
+data.head()
+
+
+
+************
+
+
+
+air_df=data[['Passengers']]
+air_df.head()
+
+
+
+
+**********
+
+air_df['shift']=air_df.Passengers.shift()
+air_df['shiftDiff']=air_df.Passengers - air_df['shift']
+air_df.head()
+
+
+**********
+
+test_stationarity(air_df.dropna(),'shiftDiff')
+
+
+***********
+
+log_df=air_df[['Passengers']]
+log_df['log']=np.log(log_df['Passengers'])
+log_df.head()
+
+
+************
+
+test_stationarity(log_df,'log')
+
+
+
+sqrt_df=air_df[['Passengers']]
+sqrt_df['sqrt']=np.sqrt(air_df['Passengers'])
+sqrt_df.head()
+
+
+********
+
+test_stationarity(sqrt_df,'sqrt')
+
+
+***********
+
+cbrt_df=air_df[['Passengers']]
+cbrt_df['cbrt']=np.cbrt(air_df['Passengers'])
+cbrt_df.head()
+
+***********
+
+test_stationarity(cbrt_df,'cbrt')
+
+
+
+
+************
+
+
+log_df2=log_df[['Passengers','log']]
+log_df2['log_sqrt']=np.sqrt(log_df['log'])
+log_df2.head()
+
+**********
+
+
+test_stationarity(log_df2,'log_sqrt')
+
+
+
+********
+
+
+log_df2=log_df[['Passengers','log']]
+log_df2['log_sqrt']=np.sqrt(log_df['log'])
+log_df2['logShiftDiff']=log_df2['log_sqrt']-log_df2['log_sqrt'].shift()
+log_df2.head()
+
+
+
+*********
+
+
+test_stationarity(log_df2.dropna(),'logShiftDiff')
+
+*************88
+
+
+
+
+
+
+
+
+
+

myawesomepkg/TSAPY1/Practical No 5.py

@@ -0,0 +1,52 @@
+Aim: Implementing auto correlation and partial auto-correlation on timeseries
+
+
+# ACF plot of time series
+from pandas import read_csv
+from matplotlib import pyplot
+#from statsmodels.graphics.tsaplots import plot_acf
+from pandas.plotting import autocorrelation_plot
+series = read_csv('/content/drive/MyDrive/MScDS TSA/daily-min-temperatures.csv', header=0, index_col=0,parse_dates=True, squeeze=True)
+#plot_acf(series)
+autocorrelation_plot(series)
+pyplot.show()
+
+
+
+
+*********
+
+
+# zoomed-in ACF plot of time series
+from pandas import read_csv
+from matplotlib import pyplot
+from statsmodels.graphics.tsaplots import plot_acf
+series = read_csv('/content/drive/MyDrive/MScDS TSA/daily-min-temperatures.csv', header=0, index_col=0,parse_dates=True, squeeze=True)
+plot_acf(series, lags=50)
+pyplot.show()
+
+
+
+**************
+
+
+
+# PACF plot of time series
+from pandas import read_csv
+from matplotlib import pyplot
+from statsmodels.graphics.tsaplots import plot_pacf
+series = read_csv('/content/drive/MyDrive/MScDS TSA/daily-min-temperatures.csv', header=0, index_col=0,
+parse_dates=True, squeeze=True)
+plot_pacf(series, lags=50)
+pyplot.show()
+
+
+
+***************8
+
+
+
+
+
+
+

myawesomepkg/TSAPY1/Practical No 6.py

@@ -0,0 +1,29 @@
+Aim: Perform autoregression on time series data
+
+
+# create and evaluate a static autoregressive model
+from pandas import read_csv
+from matplotlib import pyplot
+from statsmodels.tsa.ar_model import AutoReg
+from sklearn.metrics import mean_squared_error
+from math import sqrt
+# load dataset
+series = read_csv('/content/drive/MyDrive/MScDS TSA/daily-min-temperatures.csv', header=0, index_col=0,parse_dates=True, squeeze=True)
+# split dataset
+X = series.values
+train, test = X[1:len(X)-7], X[len(X)-7:]
+# train autoregression
+model = AutoReg(train,30)
+model_fit = model.fit()
+print('Lag: %s' % model_fit.ar_lags)
+print('Coefficients: %s' % model_fit.params)
+# make predictions
+predictions = model_fit.predict(start=len(train), end=len(train)+len(test)-1, dynamic=False)
+for i in range(len(predictions)):
+  print('predicted=%f, expected=%f' % (predictions[i], test[i]))
+rmse = sqrt(mean_squared_error(test, predictions))
+print('Test RMSE: %.3f' % rmse)
+# plot results
+pyplot.plot(test)
+pyplot.plot(predictions, color='red')
+pyplot.show()

myawesomepkg/TSAPY1/Practical No 7.py

@@ -0,0 +1,67 @@
+Aim: Forecasting using MA model.
+
+
+# correct forecasts with a model of forecast residual errors
+from pandas import read_csv
+from pandas import DataFrame
+from pandas import concat
+from statsmodels.tsa.ar_model import AutoReg
+from matplotlib import pyplot
+from sklearn.metrics import mean_squared_error
+from math import sqrt
+# load data
+series = read_csv('/content/drive/MyDrive/MScDS TSA/daily-total-female-births-CA.csv', header=0, index_col=0, parse_dates=True,squeeze=True)
+# create lagged dataset
+values = DataFrame(series.values)
+dataframe = concat([values.shift(1), values], axis=1)
+dataframe.columns = ['t', 't+1']
+print(dataframe)
+X = dataframe.values
+
+
+*******
+
+
+# split into train and test sets
+X = dataframe.values
+train_size = int(len(X) * 0.66)
+train, test = X[1:train_size], X[train_size:]
+train_X, train_y = train[:,0], train[:,1]
+test_X, test_y = test[:,0], test[:,1]
+# persistence model on training set
+train_pred = [x for x in train_X]
+# calculate residuals
+train_resid = [train_y[i]-train_pred[i] for i in range(len(train_pred))]
+# model the training set residuals
+model = AutoReg(train_resid,20)
+model_fit = model.fit()
+window = len(model_fit.ar_lags)
+coef = model_fit.params
+# walk forward over time steps in test
+history = train_resid[len(train_resid)-window:]
+history = [history[i] for i in range(len(history))]
+predictions = list()
+for t in range(len(test_y)):
+  # persistence
+  yhat = test_X[t]
+  error = test_y[t] - yhat
+  # predict error
+  length = len(history)
+  lag = [history[i] for i in range(length-window,length)]
+  pred_error = coef[0]
+  for d in range(window):
+    pred_error += coef[d+1] * lag[window-d-1]
+  # correct the prediction
+  yhat = yhat + pred_error
+  predictions.append(yhat)
+  history.append(error)
+  print('predicted=%f, expected=%f' % (yhat, test_y[t]))
+# error
+rmse = sqrt(mean_squared_error(test_y, predictions))
+print('Test RMSE: %.3f' % rmse)
+# plot predicted error
+pyplot.plot(test_y)
+pyplot.plot(predictions, color='red')
+pyplot.show()
+
+

myawesomepkg/TSAPY1/Practical No 8.py

@@ -0,0 +1,108 @@
+Aim: Forecasting using ARIMA model --TEMPERATURE 
+Time Series Forecasting With ARIMA Model in Python for Temperature Prediction.
+
+1) Reading Time Series Data in Python using Pandas library
+
+import pandas as pd
+df=pd.read_csv('/content/drive/MyDrive/MScDS TSA/MaunaLoaDailyTemps.csv',index_col='DATE',parse_dates=True)
+df=df.dropna()
+print('Shape of data',df.shape)
+df.head()
+df
+
+**********
+
+df['AvgTemp'].plot(figsize=(12,5))
+
+
+*******
+
+2) Checking for stationarity of time series model
+
+
+from statsmodels.tsa.stattools import adfuller
+def adf_test(dataset):
+     dftest = adfuller(dataset, autolag = 'AIC')
+     print("1. ADF : ",dftest[0])
+     print("2. P-Value : ", dftest[1])
+     print("3. Num Of Lags : ", dftest[2])
+     print("4. Num Of Observations Used For ADF Regression:",      dftest[3])
+     print("5. Critical Values :")
+     for key, val in dftest[4].items():
+         print("\t",key, ": ", val)
+adf_test(df['AvgTemp'])
+
+
+*************
+
+
+ Auto Arima Function to select order of Auto Regression Model
+
+pip install pmdarima
+
+
+
+from pmdarima import auto_arima
+import warnings
+warnings.filterwarnings("ignore")
+stepwise_fit=auto_arima(df['AvgTemp'],trace=True,suppress_warnings=True)
+stepwise_fit.summary()
+
+
+
+************8
+
+
+Split Your Dataset
+
+print(df.shape)
+train=df.iloc[:-30]
+test=df.iloc[-30:]
+print(train.shape,test.shape)
+
+
+
+
+
+from statsmodels.tsa.arima.model import ARIMA
+model=ARIMA(train['AvgTemp'],order=(1,0,5))
+model=model.fit()
+model.summary()
+
+
+
+
+**************
+
+
+
+Check How Good Your Model Is
+
+
+start=len(train)
+end=len(train)+len(test)-1
+pred=model.predict(start=start,end=end,typ='levels').rename('ARIMA Predictions')
+print(pred)
+pred.index=df.index[start:end+1]
+pred.plot(legend=True)
+test['AvgTemp'].plot(legend=True)
+
+
+
+
+
+***********8
+
+
+Check your Accuracy Metric
+
+
+from sklearn.metrics import mean_squared_error
+from math import sqrt
+test['AvgTemp'].mean()
+rmse=sqrt(mean_squared_error(pred,test['AvgTemp']))
+print(rmse)
+
+
+
+

myawesomepkg/TSAPY1/Print_R.py

@@ -0,0 +1,123 @@
+ 1. Print in R
+r
+Copy
+Edit
+print("Hello Boss!")
+🔹 2. Comments in R
+r
+Copy
+Edit
+# This is a single-line comment
+🔹 3. Variables in R
+r
+Copy
+Edit
+x <- 10
+y <- "Data"
+🔹 4. Concatenate Elements
+r
+Copy
+Edit
+v <- c(1, 2, 3, 4)
+print(v)
+🔹 5. Multiple Variables
+r
+Copy
+Edit
+a <- 5
+b <- 10
+c <- a + b
+print(c)
+🔹 6. Variable Names
+r
+Copy
+Edit
+user_name <- "Boss"
+user_age <- 25
+🔹 7. Data Types
+r
+Copy
+Edit
+num <- 10.5        # Numeric
+str <- "Hello"     # Character
+bool <- TRUE       # Logical
+vec <- c(1, 2, 3)  # Vector
+🔹 8. Strings
+r
+Copy
+Edit
+name <- "R Programming"
+paste("Welcome to", name)
+🔹 9. Boolean
+r
+Copy
+Edit
+is_true <- TRUE
+is_false <- FALSE
+🔹 10. Operators
+r
+Copy
+Edit
+a <- 10
+b <- 3
+a + b   # Addition
+a > b   # Comparison
+a == b  # Equal
+a %% b  # Modulus
+🔹 11. If Else
+r
+Copy
+Edit
+x <- 10
+if (x > 5) {
+  print("Greater than 5")
+} else {
+  print("5 or less")
+}
+🔹 12. List
+r
+Copy
+Edit
+my_list <- list(name="Boss", age=25, scores=c(90, 85))
+print(my_list)
+🔹 13. Matrices
+r
+Copy
+Edit
+matrix_data <- matrix(1:6, nrow=2, ncol=3)
+print(matrix_data)
+🔹 14. Data Frames
+r
+Copy
+Edit
+df <- data.frame(Name=c("A", "B"), Age=c(20, 25))
+print(df)
+🔹 15. Functions
+r
+Copy
+Edit
+add_numbers <- function(x, y) {
+  return(x + y)
+}
+🔹 16. Call a Function
+r
+Copy
+Edit
+result <- add_numbers(5, 3)
+print(result)
+🔹 17. Global Variable
+r
+Copy
+Edit
+x <- 5
+my_func <- function() {
+  x <<- 10  # Modify global x
+}
+my_func()
+print(x)
+🔹 18. Vectors
+r
+Copy
+Edit
+my_vector <- c(1, 2, 3, 4, 5)
+print(my_vector)

myawesomepkg/TSAPY1/R_Graph.py

@@ -0,0 +1,32 @@
+ 1. Line Plot
+r
+Copy
+Edit
+x <- c(1, 2, 3, 4, 5)
+y <- c(2, 4, 6, 8, 10)
+
+plot(x, y, type="l", col="blue", main="Line Plot", xlab="X-axis", ylab="Y-axis")
+🔹 2. Scatter Plot
+r
+Copy
+Edit
+x <- c(1, 2, 3, 4, 5)
+y <- c(5, 3, 6, 2, 7)
+
+plot(x, y, main="Scatter Plot", xlab="X", ylab="Y", col="red", pch=19)
+🔹 3. Pie Chart
+r
+Copy
+Edit
+slices <- c(10, 20, 30, 40)
+labels <- c("A", "B", "C", "D")
+
+pie(slices, labels=labels, main="Pie Chart")
+🔹 4. Bar Chart
+r
+Copy
+Edit
+values <- c(5, 10, 15, 20)
+names <- c("A", "B", "C", "D")
+
+barplot(values, names.arg=names, col="green", main="Bar Chart", ylab="Values")

myawesomepkg/TSAPY1/Working_Ggplot.py

@@ -0,0 +1,53 @@
+✅ Step 1: Install & Load ggplot2
+r
+Copy
+Edit
+install.packages("ggplot2")  # Run once
+library(ggplot2)
+✅ Step 2: Sample Data
+r
+Copy
+Edit
+data <- data.frame(
+  category = rep(c("A", "B", "C"), each=4),
+  subcat = rep(c("X", "Y"), times=6),
+  value = c(4, 7, 6, 9, 5, 3, 8, 4, 7, 5, 6, 2)
+)
+✅ Step 3: Basic ggplot
+r
+Copy
+Edit
+ggplot(data, aes(x=subcat, y=value)) +
+  geom_bar(stat="identity", fill="steelblue") +
+  ggtitle("Basic Bar Chart")
+✅ Step 4: Facets
+r
+Copy
+Edit
+ggplot(data, aes(x=subcat, y=value)) +
+  geom_bar(stat="identity", fill="tomato") +
+  facet_wrap(~ category) +
+  ggtitle("Faceted by Category")
+✅ Step 5: Geometric Objects
+r
+Copy
+Edit
+ggplot(data, aes(x=subcat, y=value, fill=category)) +
+  geom_bar(stat="identity", position="dodge") +   # Bar chart
+  geom_point(aes(color=category), size=3, shape=21) +  # Add points
+  ggtitle("Geometric Objects: Bars + Points")
+✅ Step 6: Position Adjustment
+r
+Copy
+Edit
+ggplot(data, aes(x=subcat, y=value, fill=category)) +
+  geom_bar(stat="identity", position=position_dodge(width=0.7)) +
+  ggtitle("Position: Dodge for Side-by-Side Bars")
+✅ Step 7: Coordinate System (Flip Axis)
+r
+Copy
+Edit
+ggplot(data, aes(x=subcat, y=value, fill=category)) +
+  geom_bar(stat="identity") +
+  coord_flip() +
+  ggtitle("Flipped Coordinates")