statslibx 0.1.7__py3-none-any.whl → 0.1.8__py3-none-any.whl

statslibx/inferential.py

@@ -1,32 +1,39 @@
 from dataclasses import dataclass
 import numpy as np
 import pandas as pd
-from typing import Optional, Union, Literal, List, Dict, Any
+import polars as pl
+from typing import Optional, Union, Literal, List, Dict, Any, Tuple
 from datetime import datetime
 from scipy import stats
 import os
 
 class InferentialStats:
     """
-    Clase para estadística inferencial (pruebas de hipótesis, intervalos de confianza, etc.)
+    Class for inferential statistics (hypothesis tests, confidence intervals, etc.)
     """
     
     def __init__(self, data: Union[pd.DataFrame, np.ndarray],
-                backend: Literal['pandas', 'polars'] = 'pandas'):
+                backend: Literal['pandas', 'polars'] = 'pandas',
+                lang: Literal['es-ES', 'en-US'] = 'es-ES'):
         """
-        Inicializar con DataFrame o array numpy
+        Initialize DataFrame
 
         Parameters:
         -----------
         data : DataFrame o ndarray
-            Datos a analizar
+            Data to analyze
         backend : str
-            'pandas' o 'polars' para procesamiento
+            'pandas' or 'polars' for processing
         """
 
         if isinstance(data, str) and os.path.exists(data):
                 data = InferentialStats.from_file(data).data
 
+        if isinstance(data, pl.DataFrame):
+            raise TypeError(
+                "Polars aún no soportado. Use pandas.DataFrame."
+            )
+
         if isinstance(data, np.ndarray):
             if data.ndim == 1:
                 data = pd.DataFrame({'var': data})
@@ -36,8 +43,9 @@ class InferentialStats:
         self.data = data
         self.backend = backend
         self._numeric_cols = data.select_dtypes(include=[np.number]).columns.tolist()
+        self.lang = lang
 
-    @staticmethod
+    @classmethod
     def from_file(path: str):
         """
         Carga automática de archivos y devuelve instancia de Intelligence.
@@ -76,14 +84,14 @@ class InferentialStats:
     def confidence_interval(self, column: str, confidence: float = 0.95,
                             statistic: Literal['mean', 'median', 'proportion'] = 'mean') -> tuple:
         """
-        Intervalo de confianza para diferentes estadísticos
+        Confidence interval for different statistics
         
         Parameters:
         -----------
         column : str
-            Columna a analizar
+            Column to analyze
         confidence : float
-            Nivel de confianza (default 0.95 = 95%)
+            Confidence level (default 0.95 = 95%)
         statistic : str
             'mean', 'median' o 'proportion'
         
@@ -128,9 +136,10 @@ class InferentialStats:
     
     def t_test_1sample(self, column: str, popmean: float = None, 
                         popmedian: float = None,
-                        alternative: Literal['two-sided', 'less', 'greater'] = 'two-sided') -> 'TestResult':
+                        alternative: Literal['two-sided', 'less', 'greater'] = 'two-sided',
+                        alpha: float = 0.05) -> 'TestResult':
         """
-        Prueba t de una muestra (para media o mediana)
+        One sample t test (for mean or median)
         
         Parameters:
         -----------
@@ -149,7 +158,7 @@ class InferentialStats:
         
         if popmean is not None:
             statistic, pvalue = stats.ttest_1samp(data, popmean, alternative=alternative)
-            
+
             return TestResult(
                 test_name='T-Test de Una Muestra (Media)',
                 statistic=statistic,
@@ -160,13 +169,14 @@ class InferentialStats:
                     'sample_mean': data.mean(), 
                     'n': len(data),
                     'df': len(data) - 1
-                }
+                },
+                alpha=alpha
             )
         
         elif popmedian is not None:
             # Wilcoxon signed-rank test para mediana
             statistic, pvalue = stats.wilcoxon(data - popmedian, alternative=alternative)
-            
+
             return TestResult(
                 test_name='Wilcoxon Signed-Rank Test (Mediana)',
                 statistic=statistic,
@@ -184,7 +194,7 @@ class InferentialStats:
     
     def t_test_2sample(self, column1: str, column2: str,
                         equal_var: bool = True,
-                        alternative: Literal['two-sided', 'less', 'greater'] = 'two-sided') -> 'TestResult':
+                        alternative: Literal['two-sided', 'less', 'greater'] = 'two-sided', alpha: float = 0.05) -> 'TestResult':
         """
         Prueba t de dos muestras independientes
         
@@ -214,11 +224,12 @@ class InferentialStats:
                 'std1': data1.std(), 'std2': data2.std(),
                 'n1': len(data1), 'n2': len(data2),
                 'equal_var': equal_var
-            }
+            },
+            alpha=alpha
         )
     
     def t_test_paired(self, column1: str, column2: str,
-                        alternative: Literal['two-sided', 'less', 'greater'] = 'two-sided') -> 'TestResult':
+                        alternative: Literal['two-sided', 'less', 'greater'] = 'two-sided', alpha: float = 0.05) -> 'TestResult':
         """
         Prueba t pareada
 
@@ -241,11 +252,12 @@ class InferentialStats:
             statistic=statistic,
             pvalue=pvalue,
             alternative=alternative,
-            params={'mean_diff': (data1 - data2).mean(), 'n': len(data1)}
+            params={'mean_diff': (data1 - data2).mean(), 'n': len(data1)},
+            alpha=alpha
         )
     
     def mann_whitney_test(self, column1: str, column2: str,
-                            alternative: Literal['two-sided', 'less', 'greater'] = 'two-sided') -> 'TestResult':
+                            alternative: Literal['two-sided', 'less', 'greater'] = 'two-sided', alpha: float = 0.05) -> 'TestResult':
         """
         Prueba de Mann-Whitney U (alternativa no paramétrica al t-test)
         
@@ -273,10 +285,12 @@ class InferentialStats:
                 'median2': data2.median(),
                 'n1': len(data1),
                 'n2': len(data2)
-            }
+            },
+            alpha=alpha
         )
     
-    def chi_square_test(self, column1: str, column2: str) -> 'TestResult':
+    def chi_square_test(self, column1: str, column2: str,
+                        alpha: float = 0.05) -> 'TestResult':
         """
         Prueba Chi-cuadrado de independencia
         
@@ -295,10 +309,12 @@ class InferentialStats:
             statistic=chi2,
             pvalue=pvalue,
             alternative='two-sided',
-            params={'dof': dof, 'contingency_table': contingency_table}
+            params={'dof': dof, 'contingency_table': contingency_table},
+            alpha=alpha
         )
     
-    def anova_oneway(self, column: str, groups: str) -> 'TestResult':
+    def anova_oneway(self, column: str, groups: str,
+                        alpha: float = 0.05) -> 'TestResult':
         """
         ANOVA de un factor
         
@@ -310,11 +326,16 @@ class InferentialStats:
             Variable de agrupación (categórica)
         """
         from scipy import stats
+        clean_data = self.data[[column, groups]].dropna()
         
-        groups_data = [group[column].values for name, group in self.data.groupby(groups)]
+        groups_data = [group[column].values
+                            for _, group in clean_data.groupby(groups)
+                            if len(group) > 1 and group[column].var() > 0
+                    ]
+
         statistic, pvalue = stats.f_oneway(*groups_data)
         
-        return TestResult(
+        return TestResult(  
             test_name='ANOVA de Un Factor',
             statistic=statistic,
             pvalue=pvalue,
@@ -322,10 +343,12 @@ class InferentialStats:
             params={
                 'groups': len(groups_data),
                 'n_total': sum(len(g) for g in groups_data)
-            }
+            },
+            alpha=alpha
         )
     
-    def kruskal_wallis_test(self, column: str, groups: str) -> 'TestResult':
+    def kruskal_wallis_test(self, column: str, groups: str,
+                            alpha: float = 0.05) -> 'TestResult':
         """
         Prueba de Kruskal-Wallis (ANOVA no paramétrico)
         
@@ -337,8 +360,13 @@ class InferentialStats:
             Variable de agrupación (categórica)
         """
         from scipy import stats
+
+        clean_data = self.data[[column, groups]].dropna()
         
-        groups_data = [group[column].values for name, group in self.data.groupby(groups)]
+        groups_data = [group[column].values
+                            for _, group in clean_data.groupby(groups)
+                            if len(group) > 1 and group[column].var() > 0
+                    ]
         statistic, pvalue = stats.kruskal(*groups_data)
         
         return TestResult(
@@ -349,12 +377,14 @@ class InferentialStats:
             params={
                 'groups': len(groups_data),
                 'n_total': sum(len(g) for g in groups_data)
-            }
+            },
+            alpha=alpha
         )
     
     def normality_test(self, column: str, 
                         method: Literal['shapiro', 'ks', 'anderson', 'jarque_bera', 'all'] = 'shapiro',
-                        test_statistic: Literal['mean', 'median', 'mode'] = 'mean') -> Union['TestResult', dict]:
+                        test_statistic: Literal['mean', 'median', 'mode'] = 'mean',
+                        alpha: float = 0.05) -> Union['TestResult', dict]:
         """
         Prueba de normalidad con múltiples métodos y estadísticos
         
@@ -396,6 +426,9 @@ class InferentialStats:
             scale = np.std(data, ddof=1)
         else:
             raise ValueError(f"test_statistic '{test_statistic}' no reconocido")
+
+        critical_values = None
+        significance_levels = None
         
         if method == 'all':
             results = {}
@@ -423,13 +456,13 @@ class InferentialStats:
             
             # Anderson-Darling
             anderson_result = stats.anderson(data, dist='norm')
-            results['anderson_darling'] = {
-                'test_name': f'Anderson-Darling ({test_statistic})',
-                'statistic': anderson_result.statistic,
-                'critical_values': anderson_result.critical_values,
-                'significance_levels': anderson_result.significance_level,
-                'params': {'n': n, 'test_statistic': test_statistic, 'loc': loc, 'scale': scale}
-            }
+            results['anderson_darling'] = TestResult(
+                test_name=f'Anderson-Darling ({test_statistic})',
+                statistic=anderson_result.statistic,
+                critical_values=anderson_result.critical_values,
+                significance_levels=anderson_result.significance_level,
+                params={'n': n, 'test_statistic': test_statistic, 'loc': loc, 'scale': scale}
+            )
             
             # Jarque-Bera
             stat_jb, p_jb = stats.jarque_bera(data)
@@ -462,14 +495,12 @@ class InferentialStats:
         
         elif method == 'anderson':
             anderson_result = stats.anderson(data, dist='norm')
-            return {
-                'test_name': f'Anderson-Darling ({test_statistic})',
-                'statistic': anderson_result.statistic,
-                'critical_values': anderson_result.critical_values,
-                'significance_levels': anderson_result.significance_level,
-                'params': {'n': n, 'test_statistic': test_statistic, 'loc': loc, 'scale': scale},
-                'interpretation': self._interpret_anderson(anderson_result)
-            }
+            test_name = f'Anderson-Darling ({test_statistic})'
+            pvalue = None
+            statistic = anderson_result.statistic
+            critical_values = anderson_result.critical_values
+            significance_levels = anderson_result.significance_level
+            params = {'n': n, 'test_statistic': test_statistic, 'loc': loc, 'scale': scale}
         
         elif method == 'jarque_bera':
             statistic, pvalue = stats.jarque_bera(data)
@@ -489,25 +520,19 @@ class InferentialStats:
             statistic=statistic,
             pvalue=pvalue,
             alternative='two-sided',
-            params=params
+            params=params,
+            critical_values=critical_values,
+            significance_levels=significance_levels,
+            alpha=alpha
         )
-    
-    def _interpret_anderson(self, anderson_result):
-        """Interpreta resultados de Anderson-Darling"""
-        interpretations = []
-        for i, (crit_val, sig_level) in enumerate(zip(anderson_result.critical_values, 
-                                                    anderson_result.significance_level)):
-            if anderson_result.statistic < crit_val:
-                interpretations.append(f"No se rechaza normalidad al {sig_level}% de significancia")
-            else:
-                interpretations.append(f"Se RECHAZA normalidad al {sig_level}% de significancia")
-        return interpretations
 
     def hypothesis_test(
             self,
             method: Literal["mean", "difference_mean", "proportion", "variance"] = "mean",
             column1: str = None,
             column2: str = None,
+            pop_mean: float = None,
+            pop_proportion: Union[float, Tuple[float, float]] = 0.5,
             alpha: float = 0.05,
             homoscedasticity: Literal["levene", "bartlett", "var_test"] = "levene") -> Dict[str, Any]:
             
@@ -522,11 +547,14 @@ class InferentialStats:
             Columnas numéricas a comparar
         alpha : float
             Nivel de significancia (default 0.05)
+        pop_mean : float
+            Media poblacional
+        pop_proportion : float
+            Proporción poblacional (default 0.5)
         homoscedasticity : str
             Método de homocedasticidad
             'levene', 'bartlett' o 'var_test' 
         """
-
         data = self.data
 
         if column1 is None:
@@ -547,7 +575,7 @@ class InferentialStats:
         # --- MAIN HYPOTHESIS TESTS ---
         if method == "mean":
             # One-sample t-test
-            t_stat, p_value = stats.ttest_1samp(x, popmean=np.mean(x))
+            t_stat, p_value = stats.ttest_1samp(x, popmean=pop_mean)
             test_name = "One-sample t-test"
 
         elif method == "difference_mean":
@@ -558,13 +586,46 @@ class InferentialStats:
 
         elif method == "proportion":
             # Proportion test (z-test)
-            p_hat = np.mean(x)
+
+            x = np.asarray(x)
+
+            # Caso 1: datos ya binarios
+            unique_vals = np.unique(x)
+            if set(unique_vals).issubset({0, 1}):
+
+                if pop_proportion is None:
+                    raise ValueError("Debe especificarse pop_proportion")
+
+                pop_p = pop_proportion
+
+            # Caso 2: datos continuos → binarizar
+            else:
+                if not isinstance(pop_proportion, tuple):
+                    raise ValueError(
+                        "Para datos continuos, pop_proportion debe ser (p0, binizar_value)"
+                    )
+
+                pop_p, binizar_value = pop_proportion
+                x = (x > binizar_value).astype(int)
+
+            if not (0 < pop_p < 1):
+                raise ValueError("pop_proportion debe estar entre 0 y 1")
+
             n = len(x)
-            z_stat = (p_hat - 0.5) / np.sqrt(0.5 * 0.5 / n)
+            p_hat = np.mean(x)
+
+            if n * pop_p < 5 or n * (1 - pop_p) < 5:
+                raise ValueError(
+                    "Condiciones del Z-test no cumplidas: np0 y n(1-p0) deben ser ≥ 5"
+                )
+
+            z_stat = (p_hat - pop_p) / np.sqrt(pop_p * (1 - pop_p) / n)
             p_value = 2 * (1 - stats.norm.cdf(abs(z_stat)))
+
             t_stat = z_stat
             test_name = "Proportion Z-test"
 
+
         elif method == "variance":
             # Classic F-test
             var_x = np.var(x, ddof=1)
@@ -577,15 +638,19 @@ class InferentialStats:
             t_stat = F
             test_name = "Variance F-test"
 
-        return {
-            "test": test_name,
-            "statistic": t_stat,
-            "p_value": p_value,
-            "alpha": alpha,
-            "reject_H0": p_value < alpha,
-            "homoscedasticity_test": homo_result
-        }
-
+        if p_value < alpha:
+            self.interpretation = "Se RECHAZA la hipótesis nula"
+        else:
+            self.interpretation = ("Se RECHAZA la hipotesis alternativa")
+        return TestResult(
+            test_name=test_name,
+            statistic=t_stat,
+            pvalue=p_value,
+            alternative='two-sided',
+            alpha=alpha,
+            homo_result=homo_result
+        )
+    
     def _homoscedasticity_test(
         self,
         x,
@@ -617,8 +682,8 @@ class InferentialStats:
     
     def variance_test(self, column1: str, column2: str,
                     method: Literal['levene', 'bartlett', 'var_test'] = 'levene',
-                    center: Literal['mean', 'median', 'trimmed'] = 'median'
-                    ) -> 'TestResult':
+                    center: Literal['mean', 'median', 'trimmed'] = 'median',
+                    alpha: float = 0.05) -> 'TestResult':
         """
         Prueba de igualdad de varianzas entre dos columnas.
 
@@ -690,23 +755,38 @@ class InferentialStats:
             statistic=statistic,
             pvalue=pvalue,
             alternative='two-sided',
-            params=params
+            params=params,
+            alpha=alpha
         )
 
     
     def help(self):
         """
-        Muestra ayuda completa de la clase InferentialStats
+        Muestra ayuda completa de la clase DescriptiveStats
+
+        Parametros / Parameters:
+        ------------------------
+        lang: str
+            Idioma Usuario: Codigo de Idioma (es-Es) o "Español"
+            User Language: Languaje Code (en-Us) or "English"
         """
-        help_text = """
+
+        if self.lang in ["en-US", "English", "english"]:
+            self.lang = "en-US"
+        else:
+            self.lang = "es-ES"
+        help_text = " "
+        match self.lang:
+            case "es-ES":
+                help_text = """
 ╔════════════════════════════════════════════════════════════════════════════╗
 ║                   🔬 CLASE InferentialStats - AYUDA COMPLETA               ║
 ╚════════════════════════════════════════════════════════════════════════════╝
 
 📝 DESCRIPCIÓN:
-   Clase para estadística inferencial: pruebas de hipótesis, intervalos de
-   confianza y pruebas de normalidad. Permite realizar inferencias sobre
-   poblaciones a partir de muestras de datos.
+    Clase para estadística inferencial: pruebas de hipótesis, intervalos de
+    confianza y pruebas de normalidad. Permite realizar inferencias sobre
+    poblaciones a partir de muestras de datos.
 
 ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
 
@@ -716,291 +796,595 @@ class InferentialStats:
 │ 1. 📊 INTERVALOS DE CONFIANZA                                              │
 └────────────────────────────────────────────────────────────────────────────┘
 
-  • .confidence_interval(column, confidence=0.95, statistic='mean')
+    • .confidence_interval(column, confidence=0.95, statistic='mean')
     
     Calcula intervalos de confianza para diferentes estadísticos
     
-    Parámetros:
-      column      : Columna a analizar (str)
-      confidence  : Nivel de confianza (float, default 0.95 = 95%)
-      statistic   : 'mean', 'median' o 'proportion'
-    
-    Retorna: (lower_bound, upper_bound, point_estimate)
+        Parámetros:
+            column      : Columna a analizar (str)
+            confidence  : Nivel de confianza (float, default 0.95 = 95%)
+            statistic   : 'mean', 'median' o 'proportion'
+        
+        Retorna: (lower_bound, upper_bound, point_estimate)
 
 ┌────────────────────────────────────────────────────────────────────────────┐
 │ 2. 🧪 PRUEBAS DE HIPÓTESIS - UNA MUESTRA                                   │
 └────────────────────────────────────────────────────────────────────────────┘
 
-  • .t_test_1sample(column, popmean=None, popmedian=None, 
-                   alternative='two-sided')
-    
-    Prueba t de una muestra (o Wilcoxon para mediana)
+    • .t_test_1sample(column, popmean=None, popmedian=None, 
+                        alternative='two-sided')
     
-    Parámetros:
-      column      : Columna a analizar
-      popmean     : Media poblacional hipotética (para t-test)
-      popmedian   : Mediana poblacional hipotética (para Wilcoxon)
-      alternative : 'two-sided', 'less', 'greater'
+        Prueba t de una muestra (o Wilcoxon para mediana)
+        
+        Parámetros:
+            column      : Columna a analizar
+            popmean     : Media poblacional hipotética (para t-test)
+            popmedian   : Mediana poblacional hipotética (para Wilcoxon)
+            alternative : 'two-sided', 'less', 'greater'
 
 ┌────────────────────────────────────────────────────────────────────────────┐
 │ 3. 🧪 PRUEBAS DE HIPÓTESIS - DOS MUESTRAS                                  │
 └────────────────────────────────────────────────────────────────────────────┘
 
-  🔹 Pruebas Paramétricas:
-  
-  • .t_test_2sample(column1, column2, equal_var=True, 
-                   alternative='two-sided')
-    Prueba t de dos muestras independientes
-  
-  • .t_test_paired(column1, column2, alternative='two-sided')
-    Prueba t pareada (muestras dependientes)
-
-  🔹 Pruebas No Paramétricas:
-  
-  • .mann_whitney_test(column1, column2, alternative='two-sided')
-    Alternativa no paramétrica al t-test de dos muestras
-
-  🔹 Pruebas Extras:
-  • .hypothesis_test(method='mean', column1=None, column2=None, 
-                   alpha=0.05, homoscedasticity='levene')
-  • .variance_test(column1, column2, method='levene', center='median')
+    🔹 Pruebas Paramétricas:
+
+    • .t_test_2sample(column1, column2, equal_var=True, 
+                        alternative='two-sided')
+        Prueba t de dos muestras independientes
+
+    • .t_test_paired(column1, column2, alternative='two-sided')
+        Prueba t pareada (muestras dependientes)
+
+    🔹 Pruebas No Paramétricas:
+
+    • .mann_whitney_test(column1, column2, alternative='two-sided')
+        Alternativa no paramétrica al t-test de dos muestras
+
+    🔹 Pruebas Extras:
+    • .hypothesis_test(method='mean', column1=None, column2=None, 
+                        alpha=0.05, homoscedasticity='levene')
+    • .variance_test(column1, column2, method='levene', center='median')
     
 
 ┌────────────────────────────────────────────────────────────────────────────┐
 │ 4. 🧪 PRUEBAS PARA MÚLTIPLES GRUPOS                                        │
 └────────────────────────────────────────────────────────────────────────────┘
 
-  🔹 Pruebas Paramétricas:
-  
-  • .anova_oneway(column, groups)
-    ANOVA de un factor para comparar múltiples grupos
-  
-  🔹 Pruebas No Paramétricas:
-  
-  • .kruskal_wallis_test(column, groups)
-    Alternativa no paramétrica a ANOVA
+    🔹 Pruebas Paramétricas:
+
+    • .anova_oneway(column, groups)
+        ANOVA de un factor para comparar múltiples grupos
+
+    🔹 Pruebas No Paramétricas:
+
+    • .kruskal_wallis_test(column, groups)
+        Alternativa no paramétrica a ANOVA
 
 ┌────────────────────────────────────────────────────────────────────────────┐
 │ 5. 🧪 PRUEBAS PARA VARIABLES CATEGÓRICAS                                   │
 └────────────────────────────────────────────────────────────────────────────┘
 
-  • .chi_square_test(column1, column2)
-    Prueba Chi-cuadrado de independencia entre variables categóricas
+    • .chi_square_test(column1, column2)
+        Prueba Chi-cuadrado de independencia entre variables categóricas
 
 ┌────────────────────────────────────────────────────────────────────────────┐
 │ 6. 📈 PRUEBAS DE NORMALIDAD                                                │
 └────────────────────────────────────────────────────────────────────────────┘
 
-  • .normality_test(column, method='shapiro', test_statistic='mean')
+    • .normality_test(column, method='shapiro', test_statistic='mean')
     
-    Prueba si los datos siguen una distribución normal
+        Prueba si los datos siguen una distribución normal
     
-    Métodos disponibles:
-      'shapiro'      : Shapiro-Wilk (mejor para n ≤ 5000)
-      'ks'           : Kolmogorov-Smirnov
-      'anderson'     : Anderson-Darling
-      'jarque_bera'  : Jarque-Bera (basado en asimetría y curtosis)
-      'all'          : Ejecuta todos los tests
+        Métodos disponibles:
+            'shapiro'      : Shapiro-Wilk (mejor para n ≤ 5000)
+            'ks'           : Kolmogorov-Smirnov
+            'anderson'     : Anderson-Darling
+            'jarque_bera'  : Jarque-Bera (basado en asimetría y curtosis)
+            'all'          : Ejecuta todos los tests
     
-    test_statistic: 'mean', 'median' o 'mode' para centrar la distribución
+        test_statistic: 'mean', 'median' o 'mode' para centrar la distribución
 
 ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
 
 💡 EJEMPLOS DE USO:
 
-  ┌─ Ejemplo 1: Intervalos de Confianza ────────────────────────────────────┐
-  │ from inferential import InferentialStats                                │
-  │ import pandas as pd                                                      │
-  │                                                                          │
-  │ df = pd.read_csv('datos.csv')                                           │
-  │ inf_stats = InferentialStats(df)                                        │
-  │                                                                          │
-  │ # IC para la media (95%)                                                 │
-  │ lower, upper, mean = inf_stats.confidence_interval(                     │
-  │     'salario',                                                           │
-  │     confidence=0.95,                                                    │
-  │     statistic='mean'                                                    │
-  │ )                                                                        │
-  │ print(f"IC 95%: [{lower:.2f}, {upper:.2f}]")                            │
-  │                                                                          │
-  │ # IC para la mediana (bootstrap)                                         │
-  │ lower, upper, median = inf_stats.confidence_interval(                   │
-  │     'edad',                                                              │
-  │     confidence=0.99,                                                    │
-  │     statistic='median'                                                  │
-  │ )                                                                        │
-  └──────────────────────────────────────────────────────────────────────────┘
-
-  ┌─ Ejemplo 2: Prueba t de Una Muestra ────────────────────────────────────┐
-  │ # H0: μ = 50000 (la media salarial es 50000)                            │
-  │ # H1: μ ≠ 50000 (prueba bilateral)                                      │
-  │                                                                          │
-  │ resultado = inf_stats.t_test_1sample(                                   │
-  │     column='salario',                                                   │
-  │     popmean=50000,                                                      │
-  │     alternative='two-sided'                                             │
-  │ )                                                                        │
-  │                                                                          │
-  │ print(resultado)                                                         │
-  │ # Muestra: estadístico t, valor p, interpretación                       │
-  │                                                                          │
-  │ # Prueba unilateral                                                      │
-  │ resultado = inf_stats.t_test_1sample(                                   │
-  │     column='salario',                                                   │
-  │     popmean=50000,                                                      │
-  │     alternative='greater'  # H1: μ > 50000                              │
-  │ )                                                                        │
-  └──────────────────────────────────────────────────────────────────────────┘
-
-  ┌─ Ejemplo 3: Comparación de Dos Grupos ──────────────────────────────────┐
-  │ # Prueba t independiente                                                 │
-  │ resultado = inf_stats.t_test_2sample(                                   │
-  │     column1='salario_hombres',                                          │
-  │     column2='salario_mujeres',                                          │
-  │     equal_var=True,                                                     │
-  │     alternative='two-sided'                                             │
-  │ )                                                                        │
-  │ print(resultado)                                                         │
-  │                                                                          │
-  │ # Prueba Mann-Whitney (no paramétrica)                                   │
-  │ resultado = inf_stats.mann_whitney_test(                                │
-  │     column1='salario_grupo_a',                                          │
-  │     column2='salario_grupo_b',                                          │
-  │     alternative='two-sided'                                             │
-  │ )                                                                        │
-  │                                                                          │
-  │ # Prueba t pareada (mediciones antes/después)                            │
-  │ resultado = inf_stats.t_test_paired(                                    │
-  │     column1='peso_antes',                                               │
-  │     column2='peso_despues',                                             │
-  │     alternative='two-sided'                                             │
-  │ )                                                                        │
-  └──────────────────────────────────────────────────────────────────────────┘
-
-  ┌─ Ejemplo 4: ANOVA y Kruskal-Wallis ─────────────────────────────────────┐
-  │ # ANOVA para comparar múltiples grupos                                   │
-  │ resultado = inf_stats.anova_oneway(                                     │
-  │     column='rendimiento',                                               │
-  │     groups='departamento'                                               │
-  │ )                                                                        │
-  │ print(resultado)                                                         │
-  │                                                                          │
-  │ # Kruskal-Wallis (alternativa no paramétrica)                            │
-  │ resultado = inf_stats.kruskal_wallis_test(                              │
-  │     column='satisfaccion',                                              │
-  │     groups='categoria'                                                  │
-  │ )                                                                        │
-  └──────────────────────────────────────────────────────────────────────────┘
-
-  ┌─ Ejemplo 5: Chi-Cuadrado ───────────────────────────────────────────────┐
-  │ # Probar independencia entre variables categóricas                       │
-  │ resultado = inf_stats.chi_square_test(                                  │
-  │     column1='genero',                                                   │
-  │     column2='preferencia_producto'                                      │
-  │ )                                                                        │
-  │ print(resultado)                                                         │
-  │                                                                          │
-  │ # El resultado incluye la tabla de contingencia                          │
-  └──────────────────────────────────────────────────────────────────────────┘
-
-  ┌─ Ejemplo 6: Pruebas de Normalidad ──────────────────────────────────────┐
-  │ # Shapiro-Wilk (recomendado para n ≤ 5000)                              │
-  │ resultado = inf_stats.normality_test(                                   │
-  │     column='edad',                                                      │
-  │     method='shapiro',                                                   │
-  │     test_statistic='mean'                                               │
-  │ )                                                                        │
-  │ print(resultado)                                                         │
-  │                                                                          │
-  │ # Kolmogorov-Smirnov                                                     │
-  │ resultado = inf_stats.normality_test(                                   │
-  │     column='salario',                                                   │
-  │     method='ks'                                                         │
-  │ )                                                                        │
-  │                                                                          │
-  │ # Ejecutar todos los tests                                               │
-  │ resultados = inf_stats.normality_test(                                  │
-  │     column='ingresos',                                                  │
-  │     method='all',                                                       │
-  │     test_statistic='median'                                             │
-  │ )                                                                        │
-  │                                                                          │
-  │ # Acceder a cada test                                                    │
-  │ print(resultados['shapiro'])                                            │
-  │ print(resultados['kolmogorov_smirnov'])                                 │
-  │ print(resultados['anderson_darling'])                                   │
-  │ print(resultados['jarque_bera'])                                        │
-  └──────────────────────────────────────────────────────────────────────────┘
+    ┌─ Ejemplo 1: Intervalos de Confianza ────────────────────────────────────┐
+    │ from inferential import InferentialStats                                │
+    │ import pandas as pd                                                     │
+    │                                                                         │
+    │ df = pd.read_csv('datos.csv')                                           │
+    │ inf_stats = InferentialStats(df)                                        │
+    │                                                                         │
+    │ # IC para la media (95%)                                                │
+    │ lower, upper, mean = inf_stats.confidence_interval(                     │
+    │     'salario',                                                          │
+    │     confidence=0.95,                                                    │
+    │     statistic='mean'                                                    │
+    │ )                                                                       │
+    │ print(f"IC 95%: [{lower:.2f}, {upper:.2f}]")                            │
+    │                                                                         │
+    │ # IC para la mediana (bootstrap)                                        │
+    │ lower, upper, median = inf_stats.confidence_interval(                   │
+    │     'edad',                                                             │
+    │     confidence=0.99,                                                    │
+    │     statistic='median'                                                  │
+    │ )                                                                       │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Ejemplo 2: Prueba t de Una Muestra ────────────────────────────────────┐
+    │ # H0: μ = 50000 (la media salarial es 50000)                            │
+    │ # H1: μ ≠ 50000 (prueba bilateral)                                      │
+    │                                                                         │
+    │ resultado = inf_stats.t_test_1sample(                                   │
+    │     column='salario',                                                   │
+    │     popmean=50000,                                                      │
+    │     alternative='two-sided'                                             │
+    │ )                                                                       │
+    │                                                                         │
+    │ print(resultado)                                                        │
+    │ # Muestra: estadístico t, valor p, interpretación                       │
+    │                                                                         │
+    │ # Prueba unilateral                                                     │
+    │ resultado = inf_stats.t_test_1sample(                                   │
+    │     column='salario',                                                   │
+    │     popmean=50000,                                                      │
+    │     alternative='greater'  # H1: μ > 50000                              │
+    │ )                                                                       │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Ejemplo 3: Comparación de Dos Grupos ──────────────────────────────────┐
+    │ # Prueba t independiente                                                │
+    │ resultado = inf_stats.t_test_2sample(                                   │
+    │     column1='salario_hombres',                                          │
+    │     column2='salario_mujeres',                                          │
+    │     equal_var=True,                                                     │
+    │     alternative='two-sided'                                             │
+    │ )                                                                       │
+    │ print(resultado)                                                        │
+    │                                                                         │
+    │ # Prueba Mann-Whitney (no paramétrica)                                  │
+    │ resultado = inf_stats.mann_whitney_test(                                │
+    │     column1='salario_grupo_a',                                          │
+    │     column2='salario_grupo_b',                                          │
+    │     alternative='two-sided'                                             │
+    │ )                                                                       │
+    │                                                                         │
+    │ # Prueba t pareada (mediciones antes/después)                           │
+    │ resultado = inf_stats.t_test_paired(                                    │
+    │     column1='peso_antes',                                               │
+    │     column2='peso_despues',                                             │
+    │     alternative='two-sided'                                             │
+    │ )                                                                       │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Ejemplo 4: ANOVA y Kruskal-Wallis ─────────────────────────────────────┐
+    │ # ANOVA para comparar múltiples grupos                                  │
+    │ resultado = inf_stats.anova_oneway(                                     │
+    │     column='rendimiento',                                               │
+    │     groups='departamento'                                               │
+    │ )                                                                       │
+    │ print(resultado)                                                        │
+    │                                                                         │
+    │ # Kruskal-Wallis (alternativa no paramétrica)                           │
+    │ resultado = inf_stats.kruskal_wallis_test(                              │
+    │     column='satisfaccion',                                              │
+    │     groups='categoria'                                                  │
+    │ )                                                                       │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Ejemplo 5: Chi-Cuadrado ───────────────────────────────────────────────┐
+    │ # Probar independencia entre variables categóricas                      │
+    │ resultado = inf_stats.chi_square_test(                                  │
+    │     column1='genero',                                                   │
+    │     column2='preferencia_producto'                                      │
+    │ )                                                                       │
+    │ print(resultado)                                                        │
+    │                                                                         │
+    │ # El resultado incluye la tabla de contingencia                         │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Ejemplo 6: Pruebas de Normalidad ──────────────────────────────────────┐
+    │ # Shapiro-Wilk (recomendado para n ≤ 5000)                              │
+    │ resultado = inf_stats.normality_test(                                   │
+    │     column='edad',                                                      │
+    │     method='shapiro',                                                   │
+    │     test_statistic='mean'                                               │
+    │ )                                                                       │
+    │ print(resultado)                                                        │
+    │                                                                         │
+    │ # Kolmogorov-Smirnov                                                    │
+    │ resultado = inf_stats.normality_test(                                   │
+    │     column='salario',                                                   │
+    │     method='ks'                                                         │
+    │ )                                                                       │
+    │                                                                         │
+    │ # Ejecutar todos los tests                                              │
+    │ resultados = inf_stats.normality_test(                                  │
+    │     column='ingresos',                                                  │
+    │     method='all',                                                       │
+    │     test_statistic='median'                                             │
+    │ )                                                                       │
+    │                                                                         │
+    │ # Acceder a cada test                                                   │
+    │ print(resultados['shapiro'])                                            │
+    │ print(resultados['kolmogorov_smirnov'])                                 │
+    │ print(resultados['anderson_darling'])                                   │
+    │ print(resultados['jarque_bera'])                                        │
+    └─────────────────────────────────────────────────────────────────────────┘
 
 ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
 
 📊 GUÍA DE SELECCIÓN DE PRUEBAS:
 
-  ┌─ Comparar Una Muestra vs Valor de Referencia ───────────────────────────┐
-  │ Datos normales        → t_test_1sample (con popmean)                    │
-  │ Datos no normales     → t_test_1sample (con popmedian, usa Wilcoxon)   │
-  └──────────────────────────────────────────────────────────────────────────┘
+    ┌─ Comparar Una Muestra vs Valor de Referencia ───────────────────────────┐
+    │ Datos normales        → t_test_1sample (con popmean)                    │
+    │ Datos no normales     → t_test_1sample (con popmedian, usa Wilcoxon)    │
+    └─────────────────────────────────────────────────────────────────────────┘
 
-  ┌─ Comparar Dos Grupos Independientes ────────────────────────────────────┐
-  │ Datos normales        → t_test_2sample                                  │
-  │ Datos no normales     → mann_whitney_test                               │
-  └──────────────────────────────────────────────────────────────────────────┘
+    ┌─ Comparar Dos Grupos Independientes ────────────────────────────────────┐
+    │ Datos normales        → t_test_2sample                                  │
+    │ Datos no normales     → mann_whitney_test                               │
+    └─────────────────────────────────────────────────────────────────────────┘
 
-  ┌─ Comparar Dos Grupos Pareados ──────────────────────────────────────────┐
-  │ Datos normales        → t_test_paired                                   │
-  │ Datos no normales     → (use scipy.stats.wilcoxon directamente)        │
-  └──────────────────────────────────────────────────────────────────────────┘
+    ┌─ Comparar Dos Grupos Pareados ──────────────────────────────────────────┐
+    │ Datos normales        → t_test_paired                                   │
+    │ Datos no normales     → (use scipy.stats.wilcoxon directamente)         │
+    └─────────────────────────────────────────────────────────────────────────┘
 
-  ┌─ Comparar Múltiples Grupos ─────────────────────────────────────────────┐
-  │ Datos normales        → anova_oneway                                    │
-  │ Datos no normales     → kruskal_wallis_test                             │
-  └──────────────────────────────────────────────────────────────────────────┘
+    ┌─ Comparar Múltiples Grupos ─────────────────────────────────────────────┐
+    │ Datos normales        → anova_oneway                                    │
+    │ Datos no normales     → kruskal_wallis_test                             │
+    └─────────────────────────────────────────────────────────────────────────┘
 
-  ┌─ Probar Independencia entre Categóricas ────────────────────────────────┐
-  │ Variables categóricas → chi_square_test                                 │
-  └──────────────────────────────────────────────────────────────────────────┘
+    ┌─ Probar Independencia entre Categóricas ────────────────────────────────┐
+    │ Variables categóricas → chi_square_test                                 │
+    └─────────────────────────────────────────────────────────────────────────┘
 
 ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
 
 🎯 CARACTERÍSTICAS CLAVE:
 
-  ✓ Pruebas paramétricas y no paramétricas
-  ✓ Intervalos de confianza con múltiples métodos
-  ✓ Pruebas de normalidad completas
-  ✓ Interpretación automática de resultados
-  ✓ Manejo automático de valores faltantes
-  ✓ Salidas formateadas profesionales
-  ✓ Soporte para análisis bilateral y unilateral
+    ✓ Pruebas paramétricas y no paramétricas
+    ✓ Intervalos de confianza con múltiples métodos
+    ✓ Pruebas de normalidad completas
+    ✓ Interpretación automática de resultados
+    ✓ Manejo automático de valores faltantes
+    ✓ Salidas formateadas profesionales
+    ✓ Soporte para análisis bilateral y unilateral
 
 ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
 
 ⚠️  INTERPRETACIÓN DE RESULTADOS:
 
-  • Valor p < 0.05: Se rechaza H0 (evidencia significativa)
-  • Valor p ≥ 0.05: No se rechaza H0 (evidencia insuficiente)
-  • IC que no incluye el valor nulo: Evidencia contra H0
+    • Valor p < 0.05: Se rechaza H0 (evidencia significativa)
+    • Valor p ≥ 0.05: No se rechaza H0 (evidencia insuficiente)
+    • IC que no incluye el valor nulo: Evidencia contra H0
 
 ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
 
 📚 DOCUMENTACIÓN ADICIONAL:
-   Para más información sobre métodos específicos, use:
-   help(InferentialStats.nombre_metodo)
+    Para más información sobre métodos específicos, use:
+    help(InferentialStats.nombre_metodo)
 
 ╚════════════════════════════════════════════════════════════════════════════╝
     """
+            case "en-US":
+                help_text = """
+╔════════════════════════════════════════════════════════════════════════════╗
+║                   🔬 CLASS InferentialStats - COMPLETE HELP                ║
+╚════════════════════════════════════════════════════════════════════════════╝
+
+📝 DESCRIPTION:
+    Class for inferential statistics: hypothesis tests, intervals 
+    confidence and normality tests. Allows inferences to be made about 
+    populations from data samples.
+
+━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
+
+📋 MAIN METHODS:
+
+┌────────────────────────────────────────────────────────────────────────────┐
+│ 1. 📊 CONFIDENCE INTERVALS                                                 │
+└────────────────────────────────────────────────────────────────────────────┘
+
+    • .confidence_interval(column, confidence=0.95, statistic='mean')
+    
+        Calculate confidence intervals for different statistics
+        
+        Parameters:
+            column      : Column to analyze (str)
+            confidence  : Confidence level (float, default 0.95 = 95%)
+            statistic   : 'mean', 'median' or 'proportion'
+        
+        Return: (lower_bound, upper_bound, point_estimate)
+
+┌────────────────────────────────────────────────────────────────────────────┐
+│ 2. 🧪 HYPOTHESIS TESTING - A SAMPLE                                        │
+└────────────────────────────────────────────────────────────────────────────┘
+
+    • .t_test_1sample(column, popmean=None, popmedian=None, 
+                        alternative='two-sided')
+    
+        One sample t test (or Wilcoxon for median)
+    
+        Parameters:
+            column      : Column to analyze
+            popmean     : Hypothetical population mean (for t-test)
+            popmedian   : Hypothetical population median (for Wilcoxon)
+            alternative : 'two-sided', 'less', 'greater'
+
+┌────────────────────────────────────────────────────────────────────────────┐
+│ 3. 🧪 HYPOTHESIS TESTING - TWO SAMPLES                                     │
+└────────────────────────────────────────────────────────────────────────────┘
+
+    🔹 Parametric Tests:
+
+    • .t_test_2sample(column1, column2, equal_var=True, 
+                    alternative='two-sided')
+        Two independent samples t test
+
+    • .t_test_paired(column1, column2, alternative='two-sided')
+        Paired t test (dependent samples)
+
+    🔹 Non-Parametric Tests:
+
+    • .mann_whitney_test(column1, column2, alternative='two-sided')
+        Non-parametric alternative to the two-sample t-test
+
+    🔹 Extra Tests:
+    • .hypothesis_test(method='mean', column1=None, column2=None, 
+                    alpha=0.05, homoscedasticity='levene')
+    • .variance_test(column1, column2, method='levene', center='median')
+    
+
+┌────────────────────────────────────────────────────────────────────────────┐
+│ 4. 🧪 TESTING FOR MULTIPLE GROUPS                                          │
+└────────────────────────────────────────────────────────────────────────────┘
+
+    🔹 Parametric Tests:
+
+    • .anova_oneway(column, groups)
+        One-way ANOVA to compare multiple groups
+
+    🔹 Non-Parametric Tests:
+
+    • .kruskal_wallis_test(column, groups)
+        Non-parametric alternative to ANOVA
+
+┌────────────────────────────────────────────────────────────────────────────┐
+│ 5. 🧪 TESTS FOR CATEGORICAL VARIABLES                                      │
+└────────────────────────────────────────────────────────────────────────────┘
+
+    • .chi_square_test(column1, column2)
+        Chi-square test of independence between categorical variables
+
+┌────────────────────────────────────────────────────────────────────────────┐
+│ 6. 📈 NORMALITY TESTS                                                      │
+└────────────────────────────────────────────────────────────────────────────┘
+
+    • .normality_test(column, method='shapiro', test_statistic='mean')
+    
+        Tests whether the data follows a normal distribution
+        
+        Available methods:
+            'shapiro'      : Shapiro-Wilk (best for n ≤ 5000)
+            'ks'           : Kolmogorov-Smirnov
+            'anderson'     : Anderson-Darling
+            'jarque_bera'  : Jarque-Bera (based on skewness and kurtosis)
+            'all'          : Run all tests
+        
+        test_statistic: 'mean', 'median' o 'mode' to focus the distribution
+
+━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
+
+💡 EXAMPLES OF USE:
+
+    ┌─ Example 1: Confidence Intervals ───────────────────────────────────────┐
+    │ from inferential import InferentialStats                                │
+    │ import pandas as pd                                                     │
+    │                                                                         │
+    │ df = pd.read_csv('data.csv')                                            │
+    │ inf_stats = InferentialStats(df)                                        │
+    │                                                                         │
+    │ # CI for mean (95%)                                                     │
+    │ lower, upper, mean = inf_stats.confidence_interval(                     │
+    │     'salario',                                                          │
+    │     confidence=0.95,                                                    │
+    │     statistic='mean'                                                    │
+    │ )                                                                       │
+    │ print(f"IC 95%: [{lower:.2f}, {upper:.2f}]")                            │
+    │                                                                         │
+    │ # CI for the median (bootstrap)                                         │
+    │ lower, upper, median = inf_stats.confidence_interval(                   │
+    │     'edad',                                                             │
+    │     confidence=0.99,                                                    │
+    │     statistic='median'                                                  │
+    │ )                                                                       │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Example 2: One Sample t-test ──────────────────────────────────────────┐
+    │ # H0: μ = 50000 (the average salary is 50,000)                          │
+    │ # H1: μ ≠ 50000 (two-sided test)                                        │
+    │                                                                         │
+    │ result = inf_stats.t_test_1sample(                                      │
+    │     column='salary',                                                    │
+    │     popmean=50000,                                                      │
+    │     alternative='two-sided'                                             │
+    │ )                                                                       │
+    │                                                                         │
+    │ print(result)                                                           │
+    │ # Sample: t-statistic, p-value, interpretation                          │
+    │                                                                         │
+    │ # One-sided test                                                        │
+    │ result = inf_stats.t_test_1sample(                                      │
+    │     column='salary',                                                    │
+    │     popmean=50000,                                                      │
+    │     alternative='greater'  # H1: μ > 50000                              │
+    │ )                                                                       │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Example 3: Comparison of Two Groups ───────────────────────────────────┐
+    │ # Independent t test                                                    │ 
+    │ result = inf_stats.t_test_2sample(                                      │ 
+    │     column1='men_salary',                                               │ 
+    │     column2='women_salary',                                             │ 
+    │     equal_var=True,                                                     │ 
+    │     alternative='two-sided'                                             │ 
+    │ )                                                                       │ 
+    │ print(result)                                                           │ 
+    │                                                                         │ 
+    │ # Mann-Whitney test (non-parametric)                                    │ 
+    │     result = inf_stats.mann_whitney_test(                               │ 
+    │     column1='salary_group_a',                                           │ 
+    │     column2='salary_group_b',                                           │ 
+    │     alternative='two-sided'                                             │ 
+    │ )                                                                       │ 
+    │                                                                         │ 
+    │ # Paired t-test (before/after measurements)                             │ 
+    │ result = inf_stats.t_test_paired(                                       │ 
+    │     column1='weight_before',                                            │ 
+    │     column2='after_weight',                                             │ 
+    │     alternative='two-sided'                                             │ 
+    │)                                                                        │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Example 4: ANOVA and Kruskal-Wallis ───────────────────────────────────┐
+    │ # ANOVA to compare multiple groups                                      │ 
+    │ result = inf_stats.anova_oneway(                                        │ 
+    │     column='performance',                                               │ 
+    │     groups='department'                                                 │ 
+    │ )                                                                       │ 
+    │ print(result)                                                           │ 
+    │                                                                         │ 
+    │ # Kruskal-Wallis (non-parametric alternative)                           │ 
+    │ result = inf_stats.kruskal_wallis_test(                                 │ 
+    │     column='satisfaction',                                              │ 
+    │     groups='category'                                                   │ 
+    │)                                                                        │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Example 5: Chi-Square ─────────────────────────────────────────────────┐
+    │ # Test independence between categorical variables                       │ 
+    │ result = inf_stats.chi_square_test(                                     │ 
+    │     column1='gender',                                                   │ 
+    │     column2='product_preference'                                        │ 
+    │ )                                                                       │ 
+    │ print(result)                                                           │ 
+    │                                                                         │ 
+    │ # The result includes the contingency table                             │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Example 6: Normality Tests ────────────────────────────────────────────┐
+    │ # Shapiro-Wilk (recommended for n ≤ 5000)                               │ 
+    │ result = inf_stats.normality_test(                                      │ 
+    │     column='age',                                                       │ 
+    │     method='shapiro',                                                   │ 
+    │     test_statistic='mean'                                               │ 
+    │ )                                                                       │ 
+    │ print(result)                                                           │ 
+    │                                                                         │ 
+    │ # Kolmogorov-Smirnov                                                    │ 
+    │ result = inf_stats.normality_test(                                      │ 
+    │     column='salary',                                                    │ 
+    │     method='ks'                                                         │ 
+    │ )                                                                       │ 
+    │                                                                         │ 
+    │ # Run all tests                                                         │ 
+    │ results = inf_stats.normality_test(                                     │ 
+    │     column='income',                                                    │ 
+    │     method='all',                                                       │ 
+    │     test_statistic='median'                                             │ 
+    │ )                                                                       │ 
+    │                                                                         │ 
+    │ # Access each test                                                      │ 
+    │ print(results['shapiro'])                                               │ 
+    │ print(results['kolmogorov_smirnov'])                                    │ 
+    │ print(results['anderson_darling'])                                      │ 
+    │ print(results['jarque_bera'])                                           │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
+
+📊 GUÍA DE SELECCIÓN DE PRUEBAS:
+
+    ┌─ Compare A Sample vs Reference Value ───────────────────────────────────┐
+    │ Normal data           → t_test_1sample (with mean)                      │ 
+    │ Non-normal data       → t_test_1sample (with popmedian, uses Wilcoxon)  │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Compare Two Independent Groups ────────────────────────────────────────┐
+    │ Normal data           → t_test_2sample                                  │ 
+    │ Non-normal data       → mann_whitney_test                               │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Compare Two Paired Groups ─────────────────────────────────────────────┐
+    │ Normal data           → t_test_paired                                   │ 
+    │ Non-normal data       → (use scipy.stats.wilcoxon directly)             │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Compare Multiple Groups ───────────────────────────────────────────────┐
+    │ Normal data           → anova_oneway                                    │ 
+    │ Non-normal data       → kruskal_wallis_test                             │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+    ┌─ Testing Independence between Categories ───────────────────────────────┐
+    │ Categorical variables → chi_square_test                                 │
+    └─────────────────────────────────────────────────────────────────────────┘
+
+━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
+
+🎯 KEY FEATURES: 
+
+    ✓ Parametric and non-parametric tests 
+    ✓ Confidence intervals with multiple methods 
+    ✓ Complete normality tests 
+    ✓ Automatic interpretation of results 
+    ✓ Automatic handling of missing values 
+    ✓ Professional formatted outputs 
+    ✓ Support for bilateral and unilateral analysis
+
+━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
+
+⚠️  INTERPRETATION OF RESULTS: 
+
+    • P value < 0.05: H0 is rejected (significant evidence) 
+    • P value ≥ 0.05: H0 is not rejected (insufficient evidence) 
+    • CI that does not include the null value: Evidence against H0
+
+━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
+
+📚 ADDITIONAL DOCUMENTATION: 
+    For more information on specific methods, use: 
+    help(InferentialStats.method_name)
+
+╚════════════════════════════════════════════════════════════════════════════╝
+"""
         print(help_text)
 
 @dataclass
 class TestResult:
     """Clase para resultados de pruebas de hipótesis"""
     
-    def __init__(self, test_name: str, statistic: float, pvalue: float, 
-                 alternative: str, params: dict):
+    def __init__(self, test_name: str, statistic: float, alpha: float = 0.05,
+                    params: dict = None, pvalue: float = None, 
+                    alternative: str = None, critical_values=None, 
+                    significance_levels=None, homo_result=None):
         self.test_name = test_name
         self.statistic = statistic
         self.pvalue = pvalue
         self.alternative = alternative
         self.params = params
+        self.critical_values = critical_values
+        self.significance_levels = significance_levels
+        self.interpretation = "Aun no hay interpretacion"
+        self.homo_result = homo_result
+        self.alpha = alpha
+
+        if self.pvalue is not None:
+            if self.pvalue < self.alpha:
+                self.interpretation = "Se RECHAZA la hipótesis nula"
+            else:
+                self.interpretation = "Se RECHAZA la hipótesis alternativa"
         
     def __repr__(self):
         return self._format_output()
@@ -1014,28 +1398,83 @@ class TestResult:
         output.append(f"Fecha: {datetime.now().strftime('%Y-%m-%d %H:%M:%S')}")
         output.append(f"Hipótesis Alternativa: {self.alternative}")
         output.append("-" * 80)
-        
+            
         output.append("\nRESULTADOS:")
         output.append("-" * 80)
         output.append(f"{'Estadístico':<40} {self.statistic:>20.6f}")
-        output.append(f"{'Valor p':<40} {self.pvalue:>20.6e}")
-        
-        # Interpretación
-        alpha = 0.05
-        if self.pvalue < alpha:
-            interpretation = "❌ Se RECHAZA la hipótesis nula"
-        else:
-            interpretation = "✔️ No hay evidencia suficiente para rechazar la hipótesis nula"
-        
+
+        # Mostrar valores críticos o p-value
+        if self.critical_values is not None and self.significance_levels is not None:
+            output.append("Valores Críticos:")
+            for sl, cv in zip(self.significance_levels, self.critical_values):
+                output.append(f"  α = {sl:>6.3f} → {cv:.6f}")
+        elif self.pvalue is not None:
+            output.append(f"{'Valor p':<40} {self.pvalue:>20.6e}")
+
+        # -------------------------
+        # INTERPRETACIÓN
+        # -------------------------
         output.append("\nINTERPRETACIÓN:")
         output.append("-" * 80)
-        output.append(f"Alpha = {alpha}")
-        output.append(interpretation)
-        
-        output.append("\nPARÁMETROS:")
-        output.append("-" * 80)
-        for k, v in self.params.items():
-            output.append(f"{k:<40} {str(v):>20}")
-        
+
+        alpha = 0.05
+
+        # Caso tests con p-value
+        if self.pvalue is not None:
+            output.append(f"Alpha = {alpha}")
+
+            if self.pvalue < alpha:
+                output.append("❌ Se RECHAZA la hipótesis nula")
+            else:
+                output.append("✔️ No hay evidencia suficiente para rechazar la hipótesis nula")
+
+        # Caso tests con valores críticos (ej. Anderson-Darling)
+        else:
+            # Protección mínima
+            if self.significance_levels is None or self.critical_values is None:
+                output.append("Resultado no disponible")
+            else:
+                idx = min(
+                    range(len(self.significance_levels)),
+                    key=lambda i: abs(self.significance_levels[i] - alpha)
+                )
+
+                critical_value = self.critical_values[idx]
+
+                output.append(f"Nivel de significancia (α) = {alpha}")
+                output.append(f"Estadístico A² = {self.statistic:.4f}")
+                output.append(f"Valor crítico = {critical_value:.4f}")
+
+                if self.statistic > critical_value:
+                    output.append("❌ Se RECHAZA la hipótesis nula")
+                else:
+                    output.append("✔️ No hay evidencia suficiente para rechazar la hipótesis nula")
+
+        # -------------------------
+        # HOMOCEDASTICIDAD
+        # -------------------------
+        if isinstance(self.homo_result, dict):
+            homo = self.homo_result
+
+            if isinstance(homo, dict):
+                output.append("\nTEST DE HOMOCEDASTICIDAD:")
+                output.append(f"Método: {homo['method']}")
+                output.append(f"Estadístico: {homo['statistic']:.6f}")
+                output.append(f"Valor p: {homo['p_value']:.6e}")
+
+                if homo.get("equal_var") is True:
+                    output.append("✔️ Se asume igualdad de varianzas")
+                elif homo.get("equal_var") is False:
+                    output.append("❌ No se asume igualdad de varianzas")
+
+        # -------------------------
+        # PARÁMETROS
+        # -------------------------
+        if isinstance(self.params, dict):
+            output.append("\nPARÁMETROS:")
+            output.append("-" * 80)
+            for k, v in self.params.items():
+                output.append(f"{k:<40} {str(v):>20}")
+                
         output.append("=" * 80)
         return "\n".join(output)