# yum install postgis2_96postgres=# create database geodb owner iyo; postgres=# \c geodb postgres=# create extension postgis; geodb=> CREATE TABLE geotable ( id SERIAL NOT NULL PRIMARY KEY, name varchar (255) NOT NULL, geom geography(POINT…

読み込んで、列順を変えて書き込む import csvcsv_r = open("./test.csv", "r", encoding="utf-8", errors="", newline="" ) csv_w = open('output.csv', 'w') rf = csv.reader(csv_r, delimiter=",", doublequote=True, lineterminator="\r\n", quotechar='…

mport numpy as np import pandas as pd %matplotlib inline from matplotlib import pyplot as pltx = np.random.normal(50, 10, 1000) dataframe_value = pd.DataFrame(x) dataframe_value.plot(kind="kde", subplots=True, figsize=(10, 6))#集約 std_dat…

# cd /usr/local/bin #curl -O https://raw.githubusercontent.com/spawngrid/kerl/master/kerl # cd kerl # chmod a+x kerl $ kerl list releases $ kerl build 20.3 build20.3 $ mkdir -p ~/bin/erlang/20.3 $ kerl install build20.3 ~/bin/erlang/20.3/ …

$ git clone https://github.com/riywo/anyenv ~/.anyenv$ vim ~/.bash_profile export PATH="$HOME/.anyenv/bin:$PATH" eval "$(anyenv init -)"$ anyenv install pyenv $ pyenv install -l | grep ana $ pyenv install anaconda3-5.1.0 $ pyenv rehash $ p…

%matplotlib inlinefrom scipy.stats import norm import matplotlib.pyplot as plt import numpy as npx = np.arange(-3, 3, 0.001)plt.plot(x, norm.pdf(x)) plt.show()#ファイルに保存 plt.plot(x, norm.pdf(x)) plt.plot(x, norm.pdf(x, 1.0, 0.5)) plt.…

現在価値 n年後に受け取るX円の現在価値、割引率r% X/(1+r%)**n 5年後にもらえる100万円の現在価値、（国債などの利率0.8%） >>> np.pv(0.008, 5, 0, -100) 96.094236264559825将来価値 X円を利率r%でn年運用したn年後の将来価値 X*(1+r%)**n >>> n = np.pv(…

二項分布 n回の連続した試行 コインを10回連続で投げる それぞれの試行は互いに排他的 表か裏がでる 表の確率がpならば裏の確率は1-p 表が0.5ならば裏は1-0.5 それぞれの試行は独立 1回目の結果は２回目に影響しない コインを10回連続で投げた場合の表がでる…

共分散とはxの偏差とyの偏差の積の平均値 一方が増加すると一方も増加するのが正の相関 一方が増加すると一方は減少するのが負の相関 0に近づくほど相関は弱い x=[3, 8, 9, 7, 4, 5, 8, 10, 9, 7] y=[4, 6, 8, 4, 5, 4, 7, 9, 7, 6] xの平均7 yの平均6 xの偏…

#配列 list = np.array([1,2,3,4]) my_list1 = [1,2,3,4] my_array1 = np.array(my_list1) #始点、終点（含まない）、間隔 >>>np.arange(0,10,2) array([0, 2, 4, 6, 8]) ＃ 0~10の区間を15等分 >>> np.linspace(0,10,15) array([ 0. , 0.71428571, 1.428571…