1  # -*- coding: utf-8 -*-
   2  """
   3  Created on Sun May 23 00:36:02 2021
   4  
   5  @author: PC
   6  """
   7  
   8  import pandas as pd
   9  import os
  10  import scipy.stats as stats
  11  
  12  os.chdir(r'C:\Users\PC\Pytho')
  13  
  14  df=pd.read_csv('NISPUF17.csv')
  15  df=df[['P_NUMVRC', 'HAD_CPOX']].dropna()
  16  df=df[df['HAD_CPOX']< 3]
  17  
  18  corr, pval=stats.pearsonr(df["HAD_CPOX"],df["P_NUMVRC"])
  19      
  20  corr1 = df.corr()['HAD_CPOX']['P_NUMVRC']
  21  
  22  
  23  if round(corr,10) == round(corr1,10):
  24      print(corr1)