I want to check a difference between measured and estimated by MLIT(: Ministry of Land, Infrastructure, Transport and Tourism) in seminar as part of research in university.
Data mean the number of medium and large size category tracks include trailers in Japan.
# -*- coding: utf-8 -*- """ Created on Sat Jan 25 17:08:42 2020 @author: aoki """ import numpy as np import matplotlib.pyplot as plt import pandas as pd #%% jikenkyo = pd.read_excel('nums.xlsx') mlit = pd.read_excel('MLITest.xlsx') #reforme jikenkyo.heisei = jikenkyo.heisei+1988 mlit = np.transpose(mlit) jikenkyo = jikenkyo.loc[:,['heisei','ICE']] jikenkyo = jikenkyo.drop(0) mlit = mlit.drop(['year','Unnamed: 1']) mlitmiddle = (mlit.astype(float)+mlit.astype(float))*10**5 mlitall = (mlit.astype(float))*10**5 plt.plot(jikenkyo.heisei,jikenkyo.ICE,label='measured') plt.plot(mlit.index,mlitmiddle,label='MLIT estimate: over middle') plt.plot(mlit.index,mlitall,label='MLIT estimate: all') plt.legend() plt.ylabel('number of freight vehicles') plt.ylim(0,) plt.xlabel('time /year') plt.show()
A difference is larger than I think.
This is difference of category in each survey or statics maybe.
I'll check it.