Skip to main content

I am working on airport data and I am new to python. I have a code in STATA and I need help to convert it into python [closed]

I am working on airport data and I am new to python. I have a code in STATA and I need help to convert it into python. The code is as follows:

gen FlightID= origin+ dest+ op_unique_carrier+ string(op_carrier_fl_num)

egen N_Flights = count( FlightID ), by ( FlightID )

 

keep fl_date op_unique_carrier tail_num origin dest crs_dep_time dep_time dep_delay taxi_out wheels_off wheels_on taxi_in crs_arr_time arr_time FlightID N_Flights

drop taxi_in crs_arr_time //crs_dep_time

drop if dep_time==.

drop if wheels_off==.

keep if origin=="EWR"

drop if wheels_off<dep_time //revisit to bring back those flights around midnight

sort fl_date dep_time

gen fl_id=_n

//gen Departed = 1

//replace Departed= Departed[_n-1]+1 if _n>1

sort fl_date wheels_off

//gen Tookoff = 0

//replace Tookoff= Tookoff[_n-1]+1 if _n>1

//gen DepartedQ=Departed-Tookoff

gen date=date(fl_date, "YMD")

format date %td

gen dow=dow( date)+1

gen QuJump=0

replace QuJump=1 if dep_time>dep_time[_n+1] & date==date[_n+1]

gen QuJumpMin=(dep_time-dep_time[_n+1])+(wheels_off[_n+1]-wheels_off)

summarize QuJumpMin if QuJumpMin>0

gen id = _n

gen obs=2

expand obs

drop obs

sort id

gen time_stamp=dep_time

replace time_stamp=wheels_off if id==id[_n-1]

gen stamp="dep_time"

replace stamp="wheels_off" if id==id[_n-1]

sort fl_date time_stamp

gen Qid = _n

gen QLength=1

replace QLength=cond(stamp=="dep_time" & _n>=2,QLength[_n-1]+1,cond(stamp=="wheels_off" & _n>=2,QLength[_n-1]-1,QLength))

gen Tookoff = 0

replace Tookoff=cond(stamp=="wheels_off" & _n>=2,Tookoff[_n-1]+1,cond(stamp=="dep_time" & _n>=2,Tookoff[_n-1],Tookoff))

sort fl_id stamp

gen QJumpSpots=QLength-(Tookoff[_n+1]-Tookoff) if stamp=="dep_time" & fl_id==fl_id[_n+1]


source https://stackoverflow.com/questions/72382902/i-am-working-on-airport-data-and-i-am-new-to-python-i-have-a-code-in-stata-and

Labels

Labels:

Comments

Post a Comment

Popular posts from this blog

How to show number of registered users in Laravel based on usertype?

i'm trying to display data from the database in the admin dashboard i used this: <?php use Illuminate\Support\Facades\DB; $users = DB::table('users')->count(); echo $users; ?> and i have successfully get the correct data from the database but what if i want to display a specific data for example in this user table there is "usertype" that specify if the user is normal user or admin i want to user the same code above but to display a specific usertype i tried this: <?php use Illuminate\Support\Facades\DB; $users = DB::table('users')->count()->WHERE usertype =admin; echo $users; ?> but it didn't work, what am i doing wrong? source https://stackoverflow.com/questions/68199726/how-to-show-number-of-registered-users-in-laravel-based-on-usertype
Post a Comment
Read more

Why is my reports service not connecting?

I am trying to pull some data from a Postgres database using Node.js and node-postures but I can't figure out why my service isn't connecting. my routes/index.js file: const express = require('express'); const router = express.Router(); const ordersCountController = require('../controllers/ordersCountController'); const ordersController = require('../controllers/ordersController'); const weeklyReportsController = require('../controllers/weeklyReportsController'); router.get('/orders_count', ordersCountController); router.get('/orders', ordersController); router.get('/weekly_reports', weeklyReportsController); module.exports = router; My controllers/weeklyReportsController.js file: const weeklyReportsService = require('../services/weeklyReportsService'); const weeklyReportsController = async (req, res) => { try { const data = await weeklyReportsService; res.json({data}) console...
Post a Comment
Read more

ValueError: X has 10 features, but LinearRegression is expecting 1 features as input

So, I am trying to predict the model but its throwing error like it has 10 features but it expacts only 1. So I am confused can anyone help me with it? more importantly its not working for me when my friend runs it. It works perfectly fine dose anyone know the reason about it? cv = KFold(n_splits = 10) all_loss = [] for i in range(9): # 1st for loop over polynomial orders poly_order = i X_train = make_polynomial(x, poly_order) loss_at_order = [] # initiate a set to collect loss for CV for train_index, test_index in cv.split(X_train): print('TRAIN:', train_index, 'TEST:', test_index) X_train_cv, X_test_cv = X_train[train_index], X_test[test_index] t_train_cv, t_test_cv = t[train_index], t[test_index] reg.fit(X_train_cv, t_train_cv) loss_at_order.append(np.mean((t_test_cv - reg.predict(X_test_cv))**2)) # collect loss at fold all_loss.append(np.mean(loss_at_order)) # collect loss at order plt.plot(np.log(al...
Post a Comment
Read more