Package: varclushi
VarClusHiThis is a Python module to perform variable clustering (varclus) with a hierarchical structure. Varclus is a nice dimension reduction algorithm. Here is a short description:A cluster is chosen for splitting.The chosen cluster is split into two clusters by finding the first two principal components, performing an orthoblique rotation, and assigning each variable to the rotated component with which it has the higher squared correlation.Variables are iteratively reassigned to clusters to maximize the variance accounted for by the cluster components.Indented Audience:Those who are familar with the usage of varclus algorithm in other analytical software like SAS, but always feel distressed when trying to find a RIGHT python module.Pythoners who are new to varclus algorithm. The source code could help you gain a deep understanding of the math behind this algorithm.INSIGHTS & HIGHLIGHTS:(this is a pure theoretical part, ignore this bullet point does not affect the usage of this package) Existing literatures always mention we need principal components (refer step 2-3 above). Actually, implementing this algorithm DOES NOT require principle components to be calulated, correlation matrix and its eigenvectors are enough to get the squared correlation between component and variable (this can be proved by math). If our dataset has millions of observations and hundreds of variables, not using principal components will save time and memory.Python package VarClusHi can produce very similar results, if we use SAS VARCLUS Procedure as a benchmark. This gurantees the correctness of the code.:)ExampleSeedemo.ipynbfor more details.importpandasaspdfromvarclushiimportVarClusHiCreate a VarClusHi object and pass the dataframe (df) to be analyzed as a parameter, you can also specifya feature list (feat_list, default all columns of df)max second eigenvalue (maxeigval2, default 1)max clusters (maxclus, default None)Then call method varclus(), which performs hierachical variable clustering algorithmdemo1_df=pd.read_csv('https://archive.ics.uci.edu/ml/machine-learning-databases/wine-quality/winequality-red.csv',sep=';')demo1_df.drop('quality',axis=1,inplace=True)demo1_vc=VarClusHi(demo1_df,maxeigval2=1,maxclus=None)demo1_vc.varclus()<varclushi.varclushi.VarClusHi at 0x15f96e35e10>Call info, you can get the number of clusters, number of variables in each cluster (N_vars), variance explained by each cluster (Eigval1), etc.demo1_vc.infoClusterN_VarsEigval1Eigval2VarProp0032.1413570.6584130.7137861131.7668850.9009910.5889622221.3712600.6287400.6856303321.5524960.4475040.7762484411.0000000.0000001.000000Call rsquare, you can get the (1 - rsquare) ratio of each variabledemo1_vc.rsquareClusterVariableRS_OwnRS_NCRS_Ratio00fixedacidity0.8822100.2772560.16297610density0.6220700.2461940.50136220pH0.6370760.1943590.45047831freesulfurdioxide0.7777960.0103580.22453041totalsulfurdioxide0.7866600.0422940.22276151residualsugar0.2024280.0454240.83552562sulphates0.6856300.1060220.35165372chlorides0.6856300.0489030.33053483citricacid0.7762480.3982080.37181093volatileacidity0.7762480.0409200.233299104alcohol1.0000000.0820550.000000InstallationRequirements: Python 3.4+Install by pip:pip install varclushiOther Comments:The parameters controlling this algorithm only include second eigenvalues and max number of clusters. I do not develop other functions because it is enough for my use. If you have a need for more flexibility, you can reach out to me viaxuanjing@hotmail.com.Comments for source code will be added once I have time.ThanksThank my former manager ***, I first heard of this method from him. Thank my current manager Mr. Mingsong Li, who gave me enough encouragement and support to complete this project.