# 相關分析 cd C:\Users\admin\Desktop import excel 相關分析.xlsx, firstrow clear corr LNY LNX1 LNX2 LNX3 LNX4 LNX5 estpost summarize LNY LNX1 LNX2 LNX3 LNX4 LNX5, detail esttab using 計量經濟學服務中心.rtf, cells("count mean(fmt(3)) p50 sd(fmt(2)) min max") noobs append
import numpy as np import pandas as pd
df=pd.read_excel(r'C:\Users\admin\Desktop\相關分析.xlsx')df.head
df.describe()
df.corr()
corr=df[['LNY','LNX1','LNX2','LNX3','LNX4','LNX5']].corr()corr
df[['LNY','LNX1','LNX2','LNX3','LNX4','LNX5']].corr()變量LNY與LNX1-LNX5進行相關分析結果為:
LNYLNX1LNX2LNX3LNX4LNX5LNY1.0000000.9988400.9983920.9326280.9920380.992266LNX10.9988401.0000000.9953620.9269930.9901630.989255LNX20.9983920.9953621.0000000.9210520.9909270.991783LNX30.9326280.9269930.9210521.0000000.9259650.930723LNX40.9920380.9901630.9909270.9259651.0000000.999477LNX50.9922660.9892550.9917830.9307230.9994771.000000library(readxl)data=read_excel( 'C:/Users/admin/Desktop/data/相關分析.xlsx' )View(data)
cor1=cor(data[,c('LNY','LNX1','LNX2','LNX3','LNX4','LNX5')])cor1◆◆◆◆