# -*- coding: utf-8 -*-
"""
Created on Sun May 23 00:36:02 2021

@author: PC
"""

import pandas as pd
import os
import scipy.stats as stats

os.chdir(r'C:\Users\PC\Pytho')

df=pd.read_csv('NISPUF17.csv')
df=df[['P_NUMVRC', 'HAD_CPOX']].dropna()
df=df[df['HAD_CPOX']< 3]

corr, pval=stats.pearsonr(df["HAD_CPOX"],df["P_NUMVRC"])
    
corr1 = df.corr()['HAD_CPOX']['P_NUMVRC']


if round(corr,10) == round(corr1,10):
    print(corr1)