计划执行与检查节点
executing-plans
按既定实施计划完成工作,保持过程检查和反馈。
ds-confidence-intervals【效率与协作】"Provides Constructs confidence intervals using bootstrap, analytical methods, and uncertainty quantification for parameter estimation"
独立 GitHub 快照 · 自动分类 · 以原作者文档为准
name: ds-confidence-intervals compatibility: opencode completeness: 95 content-types:
Comprehensive guide to confidence intervals in machine learning and data science workflows.
Confidence Intervals is a critical component of the machine learning workflow. This skill covers:
import numpy as np
import pandas as pd
from scipy import stats
def calculate_analytical_ci(data: np.ndarray, confidence: float = 0.95) -> tuple[float, float]:
"""Calculate analytical confidence interval using t-distribution."""
if not isinstance(data, np.ndarray):
data = np.asarray(data)
n = len(data)
if n < 2:
raise ValueError("Sample size must be at least 2")
mean = np.mean(data)
std_err = stats.sem(data)
margin = stats.t.ppf((1 + confidence) / 2, df=n - 1) * std_err
return (mean - margin, mean + margin)
def calculate_bootstrap_ci(data: np.ndarray, confidence: float = 0.95, n_resamples: int = 1000) -> tuple[float, float]:
"""Calculate bootstrap confidence interval via resampling."""
rng = np.random.default_rng(42)
bootstrap_means = np.empty(n_resamples)
for i in range(n_resamples):
sample = rng.choice(data, size=len(data), replace=True)
bootstrap_means[i] = np.mean(sample)
alpha = 1 - confidence
lower = np.percentile(bootstrap_means, 100 * alpha / 2)
upper = np.percentile(bootstrap_means, 100 * (1 - alpha / 2))
return (lower, upper)
import logging
import numpy as np
import pandas as pd
from typing import Any, Dict, Literal
logger = logging.getLogger(__name__)
class ConfidenceIntervalCalculator:
"""Production-grade calculator for confidence intervals."""
def __init__(self, confidence: float = 0.95, method: Literal['analytical', 'bootstrap'] = 'analytical', n_resamples: int = 1000):
self.confidence = confidence
self.method = method
self.n_resamples = n_resamples
logger.info(f"Initialized CI calculator with method={method}, confidence={confidence}")
def execute(self, data: pd.DataFrame, column: str) -> Dict[str, Any]:
"""Execute confidence interval calculation on specified column."""
if column not in data.columns:
raise KeyError(f"Column '{column}' not found in DataFrame")
values = data[column].dropna().values
if len(values) < 2:
raise ValueError("Insufficient non-null data for CI calculation")
try:
if self.method == 'analytical':
n = len(values)
mean = np.mean(values)
std_err = stats.sem(values)
margin = stats.t.ppf((1 + self.confidence) / 2, df=n - 1) * std_err
lower, upper = mean - margin, mean + margin
else:
rng = np.random.default_rng(42)
boots = np.array([np.mean(rng.choice(values, size=len(values), replace=True))
for _ in range(self.n_resamples)])
alpha = 1 - self.confidence
lower, upper = np.percentile(boots, 100 * alpha / 2), np.percentile(boots, 100 * (1 - alpha / 2))
logger.info(f"Calculated CI: [{lower:.4f}, {upper:.4f}]")
return {
'status': 'success'
'lower_bound': float(lower)
'upper_bound': float(upper)
'point_estimate': float(np.mean(values))
'method': self.method
'sample_size': int(n if self.method == 'analytical' else len(values))
}
except Exception as e:
logger.error(f"CI calculation failed: {e}")
return {'status': 'error', 'message': str(e)}
| Pitfall | Problem | Solution | |
Authoritative documentation links for this skill's domain. The model follows markdown links at load time to resolve external references and inline content.