上次給業(yè)務(wù)方交付用戶分層報告時我基于默認(rèn)參數(shù)畫了一版小提琴圖12 個用戶群約 40 萬條活躍時長樣本圖倒是“一鍵生成”毫無波瀾。結(jié)果業(yè)務(wù)方拎著圖來找我連問三句這個鼓包是系統(tǒng) bug 嗎中間那組為什么有兩個腰最右邊只有三十個樣本怎么能畫出這么光滑的曲線我當(dāng)場意識到默認(rèn)參數(shù)的小提琴圖在真實數(shù)據(jù)面前遠(yuǎn)遠(yuǎn)不夠用。這篇是系列第三篇我默認(rèn)你會畫基礎(chǔ)小提琴圖——hue、split、inner、寬度比例這些基礎(chǔ)操作自己翻前兩篇即可。第三篇想聊的是比“畫出來”更重要的階段它畫的是不是事實多組比較怎么排布才不誤導(dǎo)人統(tǒng)計標(biāo)注怎么加才經(jīng)得起追問大數(shù)據(jù)量怎么畫才不會跑十分鐘這些都是我從交付現(xiàn)場踩出來的問題直接按坑順序講。1. 面對真實數(shù)據(jù)的第一道坎小提琴圖不是“萬能分布圖”1.1 樣本量不足時你的圖正在“無中生有”先給一個反直覺的結(jié)論默認(rèn)參數(shù)下一組只有 5~10 個樣本點(diǎn)的時候小提琴圖體現(xiàn)出來的分布是誤導(dǎo)性最強(qiáng)的。小提琴圖的形狀來自核密度估計KDE原理是在每個觀測點(diǎn)上放一個平滑核函數(shù)疊加后歸一化。核函數(shù)的寬度帶寬由數(shù)據(jù)自動估算但樣本量少時沒有數(shù)據(jù)的區(qū)域幾乎完全由這個“核寬度”決定而不是由真實樣本事實決定。5 個點(diǎn)可以走出一個“雙峰”10 個點(diǎn)可以走出一條“厚尾”但這些峰和尾換個隨機(jī)種子就可能消失。也就是說你在小樣本下看到的小提琴輪廓大概率是噪聲的重建結(jié)果。我在交付項目里定的經(jīng)驗規(guī)則如下樣本量區(qū)間建議做法原因n 10直接用抖動散點(diǎn)圖 箱線圖小提琴輪廓基本與真實分布無關(guān)10 ≤ n 50可以用小提琴但必須疊加原始數(shù)據(jù)點(diǎn)提供原始證據(jù)防過度解讀輪廓50 ≤ n 500小提琴圖的最佳使用區(qū)間KDE 參數(shù)趨于穩(wěn)定形態(tài)可信n ≥ 500放心使用但注意帶寬選擇和渲染性能誤導(dǎo)風(fēng)險低性能負(fù)擔(dān)上升重點(diǎn)在小樣本疊加原始點(diǎn)這個操作。為什么因為小提琴圖的“平滑”本身就是一種視覺修辭它在小數(shù)量的支撐下一旦被放大看很容易讓讀者相信“形狀背后有大量連續(xù)的數(shù)據(jù)”。疊加點(diǎn)之后讀者一眼就能看出“原來只是 7 個點(diǎn)撐出來的形狀”誤讀概率大幅下降。實操細(xì)節(jié)疊加散點(diǎn)時固定隨機(jī)種子點(diǎn)的大小控制在 3~5 之間alpha 設(shè)為 0.3 左右抖動幅度不要超過小提琴寬度的一半。這樣既能看到原始數(shù)據(jù)又不會把小提琴輪廓蓋掉。1.2 零膨脹、雙峰和截斷分布的三類典型失真小樣本只是第一道坎更隱蔽的是數(shù)據(jù)本身的形狀問題。零膨脹數(shù)據(jù)。典型場景是用戶行為數(shù)據(jù)活躍時長、點(diǎn)擊次數(shù)、購買金額大量樣本集中在 0 或極小值區(qū)域。默認(rèn) KDE 會在 0 這個“尖峰”周圍產(chǎn)生一個細(xì)長的鼓包同時因為平滑作用0 的信號會被拖到旁邊的小正值上看起來像一條光滑曲線從 0 緩慢“生長”出來。但事實上0 和 0.5 可能根本不來自同一個分布——前者是“沒發(fā)生”后者是“發(fā)生了但很少”?；煸谝黄甬媯乌厔菥统鰜砹?。處理思路先做描述性拆分——報告“零值占比”非零部分單獨(dú)畫小提琴圖。如果需要看整體分布用零膨脹兩成分模型擬合后再畫預(yù)測分布而不是拿核密度直接硬畫。雙峰數(shù)據(jù)。默認(rèn)帶寬在雙峰數(shù)據(jù)下極易把兩個峰“糊”成一個寬峰。原因很直接R 中bw.nrd0Silverman 規(guī)則和 Python 默認(rèn)帶寬都是為近似正態(tài)數(shù)據(jù)設(shè)計的遇到雙峰數(shù)據(jù)帶寬估計往往偏大平滑過頭兩個峰中間的谷被填平。診斷方法很簡單同一組數(shù)據(jù)把帶寬依次設(shè)為默認(rèn)值的 0.5 倍、1 倍、2 倍畫三張圖并排對比。如果三張圖的輪廓結(jié)構(gòu)發(fā)生了質(zhì)變雙峰變單峰說明當(dāng)前帶寬選擇不可靠需要手動下調(diào)。截斷數(shù)據(jù)。百分比、評分、概率這類有天然上下限的數(shù)據(jù)KDE 會把核函數(shù)“伸”出邊界。比如 0~100 的評分100 分處的核函數(shù)會向 110 方向延伸圖上就出現(xiàn)一條從 100 延出的虛尾。業(yè)務(wù)方看到會愣難道有人打了 110 分解決辦法是邊界修正或者至少設(shè)置cut0/trimTRUE切斷尾部別讓讀者誤讀物理范圍。1.3 什么場景該換成直方圖、箱線圖或 ECDF踩過這些坑之后我把自己對“小提琴圖生態(tài)位”的理解整理成了一個判斷框架如果目標(biāo)是判斷“是否存在雙峰 / 多峰”直方圖 適當(dāng) binwidth 比小提琴圖誠實直方圖的箱寬對應(yīng)真實數(shù)據(jù)區(qū)間沒有平滑修飾。如果目標(biāo)是“保留每個點(diǎn)的完整信息”ECDF累積分布或抖動散點(diǎn)圖更好沒有任何信息被平滑掉。如果目標(biāo)是“組間中位數(shù)和四分位距的直接比較”箱線圖更緊湊讀取成本更低。如果目標(biāo)是“分布形狀 組間統(tǒng)計對比 多分組復(fù)合展示”小提琴圖配合統(tǒng)計標(biāo)注是目前唯一順手的選擇。總結(jié)成一句話小提琴圖是分布比較的“集成展示器”但不是分布估計的“顯微鏡”。使用前先分清你要它演哪個角色。2. 把核密度估計調(diào)準(zhǔn)小提琴圖才有統(tǒng)計底氣2.1 帶寬是靈魂自動規(guī)則的兩處失效R 的bw.nrd0、Python 的scott/silverman這些默認(rèn)帶寬公式都有一個隱含前提數(shù)據(jù)近似正態(tài)分布。Silverman 規(guī)則的經(jīng)驗公式大致是bw 0.9 * min(sd, IQR / 1.349) * n^(-1/5)直觀翻譯一下它把數(shù)據(jù)的離散程度濃縮成一個數(shù)標(biāo)準(zhǔn)差和四分位距取較小者再根據(jù)樣本量打折扣。問題在于數(shù)據(jù)偏態(tài)嚴(yán)重時標(biāo)準(zhǔn)差會被長尾放大四分位距會被密集中心主導(dǎo)兩者取最小往往不能反映局部分布的變化。比如右偏分布中IQR可能很小而sd很大兩者取 min 后帶寬被壓得過窄小提琴圖就會變成鋸齒狀。我在項目里的做法是三層判斷全局帶寬審計同一組數(shù)據(jù)帶寬乘 0.5、1、2 各畫一版如果形狀結(jié)構(gòu)峰、谷、肩變化不大說明這個帶寬大致可靠如果結(jié)構(gòu)變了說明默認(rèn)帶寬不能直接信任。分組帶寬是否統(tǒng)一組間樣本量差異超過 10 倍或組間離散程度差異超過 5 倍時建議放棄統(tǒng)一帶寬改為按組計算帶寬否則統(tǒng)一帶寬才有助于公平比較。在圖上標(biāo)注帶寬值圖注里寫清楚bw 0.23讓讀者知道曲線的“平滑尺度”是什么。這個習(xí)慣能減少大量誤讀尤其是交付給數(shù)據(jù)團(tuán)隊以外的人時。2.2 邊界截斷別讓分布尾巴延伸出數(shù)據(jù)范圍默認(rèn) KDE 的 x 軸范圍會向數(shù)據(jù)范圍之外擴(kuò)展幾個帶寬這是數(shù)學(xué)上正常的但對業(yè)務(wù)圖是致命的。ggplot2 的geom_violin有trim參數(shù)可以裁剪延伸部分Python 的seaborn.violinplot也有cut參數(shù)設(shè)置為cut0即可切斷超出數(shù)據(jù)范圍的尾部。實操代碼sns.violinplot(datadf, xgroup, yvalue, cut0, bw_method0.2, innerquart)ggplot(df, aes(x group, y value)) geom_violin(trim TRUE, bw 0.2) stat_summary(fun.data median_iqr, geom pointrange)注意cut0后圖頂部會顯得“平頭”這是正常的因為數(shù)據(jù)范圍之外本來就沒有信息。如果有人質(zhì)疑“怎么砍掉了尾巴”正好借機(jī)解釋 KDE 的邊界效應(yīng)。2.3 偏態(tài)數(shù)據(jù)的三條路線log 變換、分位數(shù)變換與零膨脹兩成分建模偏態(tài)數(shù)據(jù)是業(yè)務(wù)數(shù)據(jù)的主力形態(tài)我總結(jié)出三條路線log 變換對數(shù)據(jù)先取 log 再做 KDE最后坐標(biāo)軸按對數(shù)刻度呈現(xiàn)。這樣做小提琴輪廓描述的是“對數(shù)空間下的分布”比較偏態(tài)數(shù)據(jù)的相對差異更合理。有個細(xì)節(jié)必須注意變換后軸標(biāo)簽要標(biāo)注清楚比如log10(金額)否則業(yè)務(wù)方容易拿原尺度的預(yù)期來讀圖。分位數(shù)變換把原始值映射到標(biāo)準(zhǔn)正態(tài)的分位數(shù)上rank-based inverse normal transformation畫出來一定是接近對稱的。這種圖適合“只看形狀”的內(nèi)部探索基本不適合對外交付因為坐標(biāo)軸失去了原始尺度語義。零膨脹兩成分建模先畫出“非零部分”的小提琴圖再把零值占比單獨(dú)標(biāo)注出來。具體做法是小提琴圖只承載非零數(shù)據(jù)的分布零值占比用圖注或附加一個獨(dú)立的窄條符號表示。效果比把所有零值都塞進(jìn) KDE 干凈得多也更能解釋業(yè)務(wù)規(guī)律。我在實際交付中遇到“70% 的零值 30% 的長尾”這種數(shù)據(jù)首選第三條路線因為它不糊弄數(shù)據(jù)也不糊弄讀者。3. 多組比較的高級形態(tài)半小提琴、雨云圖與分組小提琴3.1 配對數(shù)據(jù)用半小提琴和雨云圖信息密度翻倍當(dāng)數(shù)據(jù)是配對的前測 / 后測、實驗組 / 對照組、同一用戶兩個時間點(diǎn)完整的小提琴圖反而有問題兩組密度對稱繪制浪費(fèi)了大量水平空間而且配對關(guān)系完全沒有體現(xiàn)。這時候半小提琴圖half-violin和雨云圖raincloud plot是更好的形態(tài)。雨云圖的經(jīng)典結(jié)構(gòu)下方是箱線圖或點(diǎn)圖上方半側(cè)小提琴展示分布中間再用連線體現(xiàn)配對的個體軌跡。R 里用gghalves包可以一行實現(xiàn)geom_half_violin()Python 可以用ptitprince庫。library(ggplot2) library(gghalves) ggplot(df, aes(x time, y score, fill time)) geom_half_violin(side l, position position_nudge(x -0.1)) geom_point(aes(color subject), size 1, alpha 0.5, position position_jitter(width 0.05)) geom_line(aes(group subject), color grey70, alpha 0.4)import ptitprince as pt ax pt.RainCloud(xtime, yscore, datadf, paletteSet2, bw0.2, width_violin0.6, orienth)這類圖中最關(guān)鍵的不是小提琴本身而是中間的連接線——它讓“每位用戶都有一條軌跡”這個信息浮出水面比單純比較兩塊分布有說服力得多。匯報時我會額外把配對差異的中位數(shù)差標(biāo)注出來比如“Δ中位數(shù) 3.2”讀圖效率明顯高于完整小提琴圖。3.2 組內(nèi)嵌套比較分組小提琴圖的布局與寬度邏輯第三種高頻需求是一個 X 軸上有兩大類每類下面再分小類。例如對照組 / 實驗組 × 三個地區(qū)共 6 個分布。直接畫 6 條獨(dú)立小提琴圖會丟失嵌套結(jié)構(gòu)正確做法是用hue分組。Python 里新版 seaborn 有個非常有用的gap參數(shù)sns.violinplot(datadf, xregion, yvalue, huearm, splitFalse, innerquart, gap0.05, density_normwidth)gap控制同一 X 刻度下不同 hue 組之間的間隔。默認(rèn) 0 時兩側(cè)緊貼組數(shù)一多容易被擠成“面條”我一般設(shè) 0.05~0.1。density_normwidth會讓所有組縮放到相同寬度組間比較形狀更公平如果改用count寬度會反映樣本量差異適合強(qiáng)調(diào)樣本量失衡的場景。R 側(cè)的對應(yīng)寫法ggplot(df, aes(x region, y value, fill arm)) geom_violin(position position_dodge(width 0.8), scale width)這里有個容易犯的錯一旦 hue 有 3 個以上水平顏色區(qū)分度會急劇下降。我習(xí)慣在顏色之外增加線型或 alpha 分層比如某一組用虛線邊框、另一組用高透明度避免單純依賴顏色通道。寬度語義需要反復(fù)權(quán)衡默認(rèn)scalewidth會把所有組拉平形狀可比性優(yōu)先scalecount保留樣本量差異但會產(chǎn)生額外的視覺噪聲。我的默認(rèn)選擇是width樣本量信息交給圖注和表格請記住這個決策邏輯。3.3 拆分小提琴圖兩組對比時最緊湊但注意讀數(shù)陷阱拆分小提琴圖split violin把兩組密度各畫一半共享同一個 X 軸位置。Seaborn 一行搞定sns.violinplot(datadf, xregion, yvalue, huearm, splitTrue, innerbox)但我必須提醒一個非常隱蔽的陷阱拆分時兩組的中位數(shù)線、四分位線在交界處交織如果用innerquart畫四分位數(shù)短線很容易讓讀者把“A 組的 Q1 線”錯看成“B 組的中位線”。所以我用 split 時幾乎只用innerbox讓箱線圖作為內(nèi)部參照二分位和四分位信息靠箱體結(jié)構(gòu)表達(dá)而不是靠容易混淆的短線。R 沒有原生的geom_split_violin需要引入see包或手工鏡像路徑。如果你的項目沒有強(qiáng)烈約束用 R我建議這種情況直接用 Python seaborn內(nèi)置實現(xiàn)穩(wěn)定不折騰。3.4 實戰(zhàn)案例A/B 實驗 × 三個區(qū)域的八萬樣本交付圖把上面的技術(shù)串起來說一個實際案例。上個月做活動彈窗實驗評估三組策略對照組、策略 A、策略 B三個區(qū)域共 8 萬用戶點(diǎn)擊時長數(shù)據(jù)。交付圖用的是分組小提琴圖并疊加了點(diǎn)層和統(tǒng)計標(biāo)注?？磮D拆解主軸是區(qū)域3 個刻度× hue 是策略3 個水平gap0.08分組形狀層每組的箱線圖直接嵌在小提琴內(nèi)部innerbox讀者一眼看到中位數(shù)、IQR 和分布形態(tài)散點(diǎn)層每組抽樣 10%最多 500 點(diǎn)繪制抖動散點(diǎn)透明度 0.2作用有三個——交代樣本量量級、暴露極端值、提升信息顆粒度標(biāo)注層只在華東區(qū)域的“策略 A vs 對照組”顯著對比處畫括號和星號校正后 p 0.01。業(yè)務(wù)方拿到這張圖說的原話是“一次能讀出三層信息”。這正是小提琴圖進(jìn)階的核心價值把描述統(tǒng)計、分布推斷和樣本證據(jù)放進(jìn)同一個“容器”。4. 統(tǒng)計標(biāo)注的正確疊加方式中心指標(biāo)、顯著性檢驗與置信區(qū)間4.1 中心趨勢標(biāo)注均值、中位數(shù)、眾數(shù)不能混著畫小提琴圖不像箱線圖自帶中位數(shù)線默認(rèn)畫完讀者對“中心在哪”的感覺來自密度峰。問題是偏態(tài)數(shù)據(jù)下密度峰、中位數(shù)、均值是三個不同的位置不標(biāo)注就默認(rèn)傳遞了“峰中心”的信息但那個位置的業(yè)務(wù)含義往往最弱。我的標(biāo)注規(guī)則均值 / 中位數(shù)點(diǎn)圖上一律用圓點(diǎn)或菱形點(diǎn)疊加均值配誤差線中位數(shù)通常不配誤差線除非做 bootstrap。四分位數(shù)線innerquart會畫出 Q1 / 中位數(shù) / Q3 三條短線注意這些是原始數(shù)據(jù)的分位數(shù)不是 KDE 的分位數(shù)。密度輪廓的視覺重心偶爾會和它們不重合這是正常的不要強(qiáng)行對齊。峰值位置如果要強(qiáng)調(diào)“用戶活躍時長集中在 18 分鐘”用虛線單獨(dú)標(biāo)記峰值坐標(biāo)并配上數(shù)值文本。原則一張圖最多同時標(biāo)注兩種位置指標(biāo)例如“中位數(shù)點(diǎn) IQR 線”超過兩種讀者就會完全混亂。4.2 顯著性星號先整體檢驗再校正兩兩比較給多組小提琴圖疊加顯著性標(biāo)注時最常見的錯誤是做了 n 組兩兩檢驗把原始 p 值直接標(biāo)星完全不控制多重比較。組數(shù) k5 時兩兩組合有 10 次檢驗純隨機(jī)數(shù)據(jù)也幾乎必然出現(xiàn)至少一個 p0.05。這種圖一旦被數(shù)據(jù)團(tuán)隊發(fā)現(xiàn)整張圖的可信度直接崩盤。我的代碼流程from scipy.stats import kruskal from scikit_posthocs import posthoc_dunn # 1. 先做整體檢驗 p_overall kruskal(*[group[value].values for _, group in df.groupby(group)]).pvalue # 2. 整體顯著后再做校正的兩兩比較 if p_overall 0.05: dunn_df posthoc_dunn(df, val_colvalue, group_colgroup, p_adjustbonferroni) sig_pairs [(i, j) for i in dunn_df.columns for j in dunn_df.index if dunn_df.loc[i, j] 0.05]括號高度要按該組小提琴的頂部位置動態(tài)計算不要用全局固定值。各組密度峰值高度不同固定括號高度會要么戳進(jìn)小提琴內(nèi)部要么離得太遠(yuǎn)。用密度峰值乘以 1.05 作為括號高度起點(diǎn)逐級往上堆疊。R 側(cè)有g(shù)gpubr::stat_compare_means但有個坑它默認(rèn)做兩兩wilcox.test的原始 p 值沒有自動多重比較校正。用這個包時必須顯式傳入p.adjust.method bonferroni之類的參數(shù)否則就是給自己埋雷。4.3 置信區(qū)間疊加bootstrap 優(yōu)先于正態(tài)近似正式報告經(jīng)常要求疊加置信區(qū)間但小提琴圖本身已包含形狀信息再疊加基于正態(tài)近似的誤差線等于把“分布假設(shè)”塞進(jìn)“分布展示”里兩邊不搭。更穩(wěn)妥的做法是 bootstrap對均值重采樣 1000 次計算均值置信區(qū)間用誤差線標(biāo)注圖注里寫明“誤差線為 95% bootstrap CI”。對中位數(shù)可用分位數(shù) bootstrap或者 R 的DescTools::MedianCI不要直接套 t 分布公式——中位數(shù)的抽樣分布根本不是 t 分布硬套會得出偏窄的區(qū)間。我的一般做法是把“統(tǒng)計推斷”和“分布展示”分開。分布展示用小提琴圖推斷結(jié)論用森林圖或簡單誤差線圖兩張圖并列放在報告里。如果必須合并比如一張圖就要講完故事標(biāo)注務(wù)必克制最多兩組比較括號三組以上根本畫不干凈。5. 大數(shù)據(jù)量下的性能優(yōu)化與正確姿勢5.1 四十萬行數(shù)據(jù)卡死的計算根源KDE 的默認(rèn)實現(xiàn)要為每個樣本點(diǎn)計算對每個網(wǎng)格點(diǎn)的核密度貢獻(xiàn)復(fù)雜度大約是 O(n × m)其中 n 是樣本數(shù)m 是網(wǎng)格點(diǎn)數(shù)默認(rèn)往往 100~500 個。40 萬行 × 200 網(wǎng)格點(diǎn) 8000 萬次核函數(shù)評估再加上排序、帶寬自動估計、多組重復(fù)計算卡頓幾乎是必然的。另一個消耗容易被忽略帶寬自動估算本身就是計算量組數(shù)越多重復(fù)估算越多。而渲染階段矢量路徑節(jié)點(diǎn)數(shù)量會爆炸性增長導(dǎo)出的 PDF / SVG 動輒幾十 MB排版系統(tǒng)直接報警。5.2 三個優(yōu)化方案等量抽樣、直方圖聚合、降低網(wǎng)格精度方案 A等比例抽樣。最常見也最直觀。如果 40 萬行來自 12 個組每組保留 300~500 個樣本抽樣后畫出的小提琴形狀與全量幾乎無差別。關(guān)鍵細(xì)節(jié)是保持組間抽樣率平衡別讓某組保留 10000、另一組只剩 100。sampled df.groupby(group, group_keysFalse).apply( lambda g: g.sample(min(500, len(g)), random_state42) )方案 B基于直方圖求密度。把數(shù)據(jù)先分箱得到頻數(shù)100~200 bins再用頻數(shù)做平滑能大幅降低參與計算的觀測數(shù)量。用numpy.histogramscipy.ndimage.gaussian_filter可以手寫輪廓與直接 KDE 幾乎一致但計算量小一個數(shù)量級。方案 C降低網(wǎng)格精度。R 側(cè)geom_violin內(nèi)部對應(yīng)stat_density可以把 KDE 網(wǎng)格點(diǎn)數(shù)降到 50~100n參數(shù)。Python 側(cè)bw_method設(shè)為固定值避免分組自動估算帶寬的重復(fù)計算。輸出矢量圖前ggsave用ragg::agg_png或柵格化參數(shù)把圖形按 300dpi 輸出PDF 體積能縮小 90%。5.3 Python / R / Plotly 三端性能對照環(huán)境40 萬行 12 組的實測感受優(yōu)化建議Python seaborn (matplotlib)計算約 2~5 秒PDF 導(dǎo)出體積暴增抽樣 cut0 必要時柵格化R ggplot2計算約 3~8 秒SVG 巨大geom_violin(n50) 輸出 PNG 而非矢量Plotly交互流暢但 hover 點(diǎn)太多會卡render_modewebgl只保留匯總統(tǒng)計的 hover實際項目里我 90% 用 R / Python 靜態(tài)圖做對外交付5% 用 Plotly 交互版本用于在線報表只 hover 匯總統(tǒng)計不做全量散點(diǎn)互動剩下 5% 做駕駛艙看板走預(yù)聚合 受控交互的思路。6. 進(jìn)階自查清單六個最常見也最隱蔽的坑6.1 因子順序與圖例順序悄悄不一致pandas 讀入數(shù)據(jù)后分組列如果是文本型排序按字母序R 里factor的默認(rèn) levels 同樣不一定符合業(yè)務(wù)口徑。畫圖前不顯式設(shè)置pd.Categorical(..., categories[...])或factor(..., levels...)圖序和圖例順序就會和業(yè)務(wù)方的口徑打架。排查方式很簡單畫圖命令前打印一下分組列的 categories 或 levels確認(rèn)順序。6.2 統(tǒng)一帶寬抹平了真實差異當(dāng)組間離散程度差異很大時統(tǒng)一帶寬會把窄分布組畫得過于平滑把寬分布組畫得過于毛糙。處理思路不是“選哪個更好”而是先明確目的目標(biāo)是識別差異選“按組帶寬”目標(biāo)是展示差異大小選“統(tǒng)一帶寬”并在圖注寫明。最怕的是不假思索地讓繪圖庫用默認(rèn)值然后拿著結(jié)果硬解釋。6.3 離散數(shù)據(jù)畫出偽連續(xù)曲線5 星評分、整數(shù)年齡這類離散數(shù)據(jù)畫小提琴圖KDE 會把一個又一個離散柱子平滑成看似連續(xù)的曲線但圖中完全看不出“中間沒有 2.5 分的人”這個事實。處理方式有兩種一是在圖上疊加原始點(diǎn)陣或直方條讓讀者意識到底層是離散的二是干脆用直方圖或柱狀圖替代。小提琴圖的平滑性非常有欺騙性這一點(diǎn)在業(yè)務(wù)交付時要格外警惕。6.4 物理量被畫出負(fù)值尾巴價格、時長、溫度這類有物理下界的數(shù)據(jù)經(jīng)常被默認(rèn) KDE 畫出負(fù)值尾巴。修復(fù)順序先判斷是否需要 log 變換偏態(tài)時再設(shè)置cut0/trimTRUE如果還是有邊界溢出就要用反射法邊界核密度R 的logspline或 Python 的statsmodels都有邊界校正實現(xiàn)?？?X 軸最小刻度是快速發(fā)現(xiàn)這個問題的第一道關(guān)卡。6.5 分組太擠變成“面條圖”分組緊湊時默認(rèn)寬度會按組數(shù)量壓縮小提琴圖變成細(xì)面條。處理density_normwidth保證各組寬度可比同時把gap設(shè)大一點(diǎn)0.05~0.1。如果你想讓寬度真的反映樣本量差異用density_normcount但這是明確的視覺決策不是默認(rèn)行為的意外結(jié)果。6.6 小樣本上的顯著性星號是統(tǒng)計濫用重災(zāi)區(qū)n 15 時無論小提琴圖的形狀看起來多“顯著”第一匯報對象永遠(yuǎn)是原始數(shù)據(jù)點(diǎn)加效應(yīng)量Cliffs delta 或 Cohens d。我把一條原則寫進(jìn)了團(tuán)隊自檢卡有 3 個以上分組或 100 個以下樣本時任何 p 值都必須配合效應(yīng)量一起出現(xiàn)否則寧可不標(biāo)。這條規(guī)則擋掉過至少三次會被審稿人打回的初稿。最后分享一個我自己的畫圖流程習(xí)慣現(xiàn)在做任何小提琴圖我都會先在本地跑一遍“三視圖測試”——同一數(shù)據(jù)分別用默認(rèn)帶寬、0.5 倍帶寬、1.5 倍帶寬畫三張預(yù)覽確認(rèn)形狀結(jié)構(gòu)穩(wěn)定之后才談配色、圖例和最終交付。這個習(xí)慣幫我擋住了多次“看起來很美、數(shù)據(jù)上站不住”的返工。小提琴圖第三篇比起參數(shù)和 API我更希望你記住一件事情分布圖的第一要務(wù)是忠實地呈現(xiàn)數(shù)據(jù)結(jié)構(gòu)而不是提前給結(jié)論找一個漂亮的形狀。