Skip to main content
100%

Gist 0383cb7786315e4c8a11

✓ Published0🌍 Public
Aakeshavan
Last edited Nov 25, 2014
Created on Nov 25, 2014

This example demonstrates fitting a Gaussian mixture model to multidimensional data using the expectation-maximization (EM) algorithm. It shows how the model estimates means, standard deviations, and mixture proportions for each dimension of the dataset. The code uses pandas to read the CSV file, NumPy for numerical operations, and matplotlib for inline plotting. It applies the EM algorithm separately to each column of the data, then prints the fitted parameters for each dimension.

AI-generated description

Similar vizzes

Loading thumbnail…

Gist d4b2a37f3ec9a90469c8

This visualization simulates longitudinal brain-imaging measurements across multiple subjects and time points, showing how percentile values (e.g., MTR_25) vary and are modeled with ordinary least squares regression per subject. It first computes subject means, average standard deviations, and true beta coefficients from the input CSV, then generates synthetic data for three time points using normal distributions. The code relies on pandas for data manipulation, statsmodels formula API for OLS fitting, numpy for random sampling, and seaborn for plotting.

Aakeshavan
61% match
Loading thumbnail…

Gist 5c928c271ee8dd212aef9236c69ed525

The example shows a PCA-based analysis of neuroimaging and clinical data from a multiple sclerosis dataset, displaying how imaging-derived features relate to demographic and clinical variables. It computes principal components from diffusion and magnetization transfer metrics, then explores correlations with clinical measures like EDSS and timed walk. The notebook uses pandas for data handling, scikit-learn's PCA and RandomForestRegressor for modeling, and statsmodels for statistical testing, with seaborn and matplotlib providing visualization support.

Aakeshavan
61% match
Loading thumbnail…

Gist c432d79c438c772e0332

This visualization demonstrates a regression analysis workflow from a Julia notebook, showing how a random feature map (phi) transforms 2D samples into a 250-dimensional space for least-squares fitting. The code computes the optimal weight vector w using normal equations with ridge regularization, evaluating the objective function across different lambda values. It uses Julia’s linear algebra operations (`transpose`, `\` for matrix division, `eye`) and `randn` for Gaussian sampling, with SVG rendering and animation to display the fitted model’s error surface and weight trajectories.

Aakeshavan
58% match
Loading thumbnail…

Gist a253f619b0b60f855f45

The notebook explores resting-state EPIC data by loading subcortical time-series from a CSV file into a pandas DataFrame, displaying the first rows to inspect voxel-level signal values alongside spatial coordinates. It uses `pd.read_csv` to parse the tabular file, which contains columns for voxel labels and time points, and `%pylab inline` to enable matplotlib plotting within the notebook. The description highlights raw time-series data before any averaging or smoothing, showing the initial structure of the dataset for a single subject.

Aakeshavan
57% match
Loading thumbnail…

Gist b6be7f1d29a9c01f70ec

This visualization shows the results of a permutation test for a Partial Least Squares (PLS) regression model predicting EDSS scores. It displays two side-by-side kernel density estimate plots of the distribution of $R^2$ values and root mean PRESS values across 1,000 random permutations, with vertical dashed lines marking the observed model performance (R² = 0.1282, rmPRESS = 1.6954). The empirical p-values computed from the permutation distributions are annotated directly on each plot. The chart uses matplotlib and seaborn, removes chartjunk, and displays the density plots with shaded colors, axis labels, and titles for each panel.

Aakeshavan
57% match
Loading thumbnail…

Gist 6857a0243ad24b95924f

This example uses pandas and Matplotlib within an IPython notebook to aggregate and compare FreeSurfer brain segmentation statistics across multiple CSV files. It processes volume, area, thickness, and aseg measurements, grouping them by subject ID and averaging values to produce normalized and unnormalized output tables. The code merges separate CSV files, collapses left/right hemisphere labels into combined metrics, and then saves the resulting aggregated data for each measurement type.

Aakeshavan
57% match
Loading thumbnail…

Gist e5affea97473980182b59a9d4efb945e

This example visualizes the relationship between brain connectivity measures and disability in multiple sclerosis patients. It fits ordinary least squares regression models predicting the EDSS score from diffusion imaging metrics and demographic covariates across each brain connection, then displays the resulting p-values and R-squared values in a table. The code uses pandas for data handling, statsmodels' `ols` API for the regression models, and seaborn with pylab inline for plotting.

Aakeshavan
56% match
Loading thumbnail…

Gist 9674847

This example simulates multi-site MRI data to illustrate how brain volume relates to age and disease severity across patient and control groups. It compares three statistical models for combining site data, including a linear calibration transformation. The notebook uses Python with NumPy, Matplotlib, seaborn, and statsmodels to generate synthetic data and fit ordinary least squares and robust linear models.

Aakeshavan
56% match