用 Python 做问卷数据清洗:无效样本检测与处理实战
做问卷研究的人都知道,回收数据里一定有垃圾。直线作答、秒填、规律选择、前后矛盾......这些数据如果不清理掉,分析结果会被严重污染,信效度也过不了关。
这篇教程教你用 Python 从零写一套问卷数据清洗流程,覆盖最常见的无效样本类型。
适用场景: 毕业论文问卷、学术调研、市场研究、用户满意度调查等。
准备工作
你需要:
- Python 3.8+
pandas、numpy、matplotlib- 一份导出为 CSV 或 Excel 的问卷数据
bash
pip install pandas numpy matplotlib openpyxl
问卷星导出的数据一般是这个结构(每行一份问卷,每列一个题目):
| 序号 | 提交时间 | Q1 | Q2 | Q3 | Q4 | Q5 |
|---|---|---|---|---|---|---|
| 1 | 2025-03-15 14:30:02 | 满意 | 一般 | 满意 | 满意 | 满意 |
| 2 | 2025-03-15 14:30:05 | 非常不满意 | 非常不满意 | 非常不满意 | 非常不满意 | 非常不满意 |
不用多说,第 2 行明显是垃圾数据。但实际场景里无效样本通常没这么明显。
第一步:加载数据
python
import pandas as pd
import numpy as np
# 问卷星导出通常是 Excel 或 CSV
df = pd.read_excel('问卷数据.xlsx')
# 或
# df = pd.read_csv('问卷数据.csv')
print(f"原始数据量: {len(df)} 份")
print(df.head())
通常问卷星导出的列名比较乱,先整理一下:
python
# 找出题目列(通常包含 "Q" 或题号)
question_cols = [c for c in df.columns if c.startswith('Q') or '第' in str(c)]
time_col = [c for c in df.columns if '时间' in str(c) or 'time' in str(c).lower()]
print(f"题目列: {question_cols}")
print(f"时间列: {time_col}")
第二步:检测直线作答(Straight-lining)
直线作答是最常见的无效模式------所有题目选了同一个选项。
python
def detect_straight_lining(df, question_cols):
"""检测所有题目答案完全一致的样本"""
invalid_indices = []
for idx, row in df.iterrows():
answers = row[question_cols].dropna().tolist()
if len(answers) < 3:
continue
# 如果所有答案都相同
if len(set(answers)) == 1:
invalid_indices.append(idx)
return set(invalid_indices)
straight_lining = detect_straight_lining(df, question_cols)
print(f"直线作答: {len(straight_lining)} 份")
但直线作答也有误判的情况。比如 5 道题都是"满意",可能是真实感受。所以加一个更严格的判断------不要求完全一致,而是 80% 以上选项相同:
python
def detect_near_straight_lining(df, question_cols, threshold=0.8):
"""检测接近直线作答(超过 threshold 比例的答案相同)"""
invalid_indices = []
for idx, row in df.iterrows():
answers = row[question_cols].dropna().tolist()
if len(answers) < 3:
continue
# 找出现最多的答案
from collections import Counter
counts = Counter(answers)
most_common_ratio = counts.most_common(1)[0][1] / len(answers)
if most_common_ratio >= threshold:
invalid_indices.append(idx)
return set(invalid_indices)
near_straight = detect_near_straight_lining(df, question_cols)
print(f"接近直线作答: {len(near_straight)} 份")
第三步:检测秒填(Speeders)
如果有人在 10 秒内"填完"了一份 20 题的问卷,基本可以确定是无效数据。
python
def detect_speeders(df, time_col, min_seconds=10):
"""检测提交时间异常短的样本"""
if not time_col:
return set()
time_data = pd.to_datetime(df[time_col[0]])
# 按时间排序,计算相邻提交的时间差
time_sorted = time_data.sort_values()
time_diffs = time_sorted.diff().dt.total_seconds()
# 方法1:如果数据里有"答题时长"字段,直接用
duration_col = [c for c in df.columns if '时长' in str(c) or '用时' in str(c)]
if duration_col:
durations = pd.to_numeric(df[duration_col[0]], errors='coerce')
speeder_indices = df[durations < min_seconds].index
return set(speeder_indices)
# 方法2:没有时长字段,用提交间隔估算(粗略)
return set()
speeders = detect_speeders(df, time_col, min_seconds=15)
print(f"秒填样本: {len(speeders)} 份")
如果问卷星导出的数据里包含答题时长列,直接按阈值过滤就行。一般规则:
- 10 题以内的问卷:< 10 秒视为秒填
- 10-20 题:< 20 秒视为秒填
- 20 题以上:< 30 秒视为秒填
第四步:检测规律作答(Pattern Response)
有些人会按规律选:1-2-3-4-5-1-2-3-4-5,或者 5-4-3-2-1-5-4-3-2-1。
python
def detect_pattern_response(df, question_cols):
"""检测规律性作答(递增、递减、交替等模式)"""
invalid_indices = []
for idx, row in df.iterrows():
answers = row[question_cols].dropna().tolist()
if len(answers) < 5:
continue
# 把文字选项转为数字(如果是文字的话需要先编码)
try:
numeric = pd.to_numeric(answers)
except:
# 文字选项,检查是否完全交替(A-B-A-B)
if len(answers) >= 6:
half = len(answers) // 2
first_half = answers[:half]
second_half = answers[half:half*2]
if first_half == second_half:
invalid_indices.append(idx)
continue
# 检测严格递增/递减
diffs = np.diff(numeric)
if np.all(diffs == diffs[0]) and diffs[0] != 0:
invalid_indices.append(idx)
continue
# 检测交替模式(1-3-1-3-1-3)
if len(numeric) >= 6:
even_vals = numeric[::2]
odd_vals = numeric[1::2]
if len(set(even_vals)) == 1 and len(set(odd_vals)) == 1:
invalid_indices.append(idx)
return set(invalid_indices)
patterns = detect_pattern_response(df, question_cols)
print(f"规律作答: {len(patterns)} 份")
第五步:检测矛盾回答
如果 Q1 问"您是否使用过本产品"选了"否",但 Q5 问"您对产品的哪些功能最满意"却填了具体内容,这就是矛盾回答。
python
def detect_contradictions(df, filter_col, filter_value, check_col):
"""
检测逻辑矛盾
filter_col: 筛选题列(如"是否使用过")
filter_value: 筛选值(如"否")
check_col: 应该为空的列(如后续满意度题)
"""
mask = df[filter_col] == filter_value
should_be_empty = df.loc[mask, check_col]
# 筛选了"否"的人,后续题不应该有内容
contradictions = should_be_empty[should_be_empty.notna() & (should_be_empty != '')]
return set(contradictions.index)
# 示例:选了"没使用过"但后面题有回答
# contradictions = detect_contradictions(df, 'Q1', '否', 'Q5')
这一步需要根据具体问卷逻辑来写,没有通用模板。
第六步:检测异常值(统计方法)
对于量表题(1-5 分或 1-7 分),可以用统计方法检测异常样本。
python
def detect_outliers_zscore(df, question_cols, threshold=3.0):
"""用 Z-score 检测每个样本的总分异常"""
# 把题目转为数值
numeric_df = df[question_cols].apply(pd.to_numeric, errors='coerce')
row_means = numeric_df.mean(axis=1)
z_scores = (row_means - row_means.mean()) / row_means.std()
outlier_indices = df[abs(z_scores) > threshold].index
return set(outlier_indices)
outliers = detect_outliers_zscore(df, question_cols)
print(f"统计异常: {len(outliers)} 份")
第七步:汇总所有无效样本
python
# 合并所有检测结果
all_invalid = straight_lining | near_straight | speeders | patterns | outliers
print(f"\n{'='*50}")
print(f"原始数据: {len(df)} 份")
print(f"直线作答: {len(straight_lining)} 份")
print(f"接近直线: {len(near_straight)} 份")
print(f"秒填样本: {len(speeders)} 份")
print(f"规律作答: {len(patterns)} 份")
print(f"统计异常: {len(outliers)} 份")
print(f"总无效: {len(all_invalid)} 份")
print(f"有效率: {(1 - len(all_invalid)/len(df))*100:.1f}%")
print(f"{'='*50}")
# 标记无效样本
df['is_invalid'] = df.index.isin(all_invalid)
df['invalid_reason'] = ''
df.loc[df.index.isin(straight_lining), 'invalid_reason'] += '直线作答;'
df.loc[df.index.isin(near_straight), 'invalid_reason'] += '接近直线;'
df.loc[df.index.isin(speeders), 'invalid_reason'] += '秒填;'
df.loc[df.index.isin(patterns), 'invalid_reason'] += '规律作答;'
df.loc[df.index.isin(outliers), 'invalid_reason'] += '统计异常;'
# 导出清洗结果
clean_df = df[~df['is_invalid']].copy()
clean_df.to_excel('问卷数据_已清洗.xlsx', index=False)
# 导出无效样本明细(方便人工复核)
invalid_df = df[df['is_invalid']].copy()
invalid_df.to_excel('无效样本明细.xlsx', index=False)
print(f"\n清洗后有效数据: {len(clean_df)} 份")
第八步:可视化清洗结果
python
import matplotlib.pyplot as plt
plt.rcParams['font.sans-serif'] = ['SimHei', 'Microsoft YaHei']
plt.rcParams['axes.unicode_minus'] = False
# 无效类型分布
reasons = {
'直线作答': len(straight_lining),
'接近直线': len(near_straight - straight_lining),
'秒填': len(speeders),
'规律作答': len(patterns),
'统计异常': len(outliers),
}
fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(12, 5))
# 左图:无效类型柱状图
ax1.bar(reasons.keys(), reasons.values(), color='#ff6b6b')
ax1.set_title('无效样本类型分布')
ax1.set_ylabel('份数')
for i, v in enumerate(reasons.values()):
ax1.text(i, v + 0.5, str(v), ha='center')
# 右图:有效率饼图
valid = len(clean_df)
invalid = len(all_invalid) - valid + len(all_invalid) # 避免双重计算
ax2.pie([len(clean_df), len(all_invalid)],
labels=['有效样本', '无效样本'],
autopct='%1.1f%%',
colors=['#51cf66', '#ff6b6b'])
ax2.set_title(f'样本有效率 (总计 {len(df)} 份)')
plt.tight_layout()
plt.savefig('数据清洗报告.png', dpi=150)
plt.show()
完整代码(一键运行)
把上面的函数整合成一个类,方便复用:
python
import pandas as pd
import numpy as np
from collections import Counter
class SurveyDataCleaner:
"""问卷数据清洗工具"""
def __init__(self, filepath: str, question_cols: list = None):
if filepath.endswith('.xlsx'):
self.df = pd.read_excel(filepath)
else:
self.df = pd.read_csv(filepath)
self.question_cols = question_cols or self._auto_detect_questions()
self.invalid_sets = {}
def _auto_detect_questions(self):
cols = [c for c in self.df.columns
if any(k in str(c) for k in ['Q', '第', '题'])]
return cols
def run_all(self, min_seconds=15, z_threshold=3.0, straight_threshold=0.8):
self.invalid_sets['直线作答'] = self._detect_straight()
self.invalid_sets['接近直线'] = self._detect_near_straight(straight_threshold)
self.invalid_sets['秒填'] = self._detect_speeders(min_seconds)
self.invalid_sets['规律作答'] = self._detect_pattern()
self.invalid_sets['统计异常'] = self._detect_outliers(z_threshold)
return self
def _detect_straight(self):
invalid = set()
for idx, row in self.df.iterrows():
ans = row[self.question_cols].dropna().tolist()
if len(ans) >= 3 and len(set(ans)) == 1:
invalid.add(idx)
return invalid
def _detect_near_straight(self, threshold):
invalid = set()
for idx, row in self.df.iterrows():
ans = row[self.question_cols].dropna().tolist()
if len(ans) >= 3:
counts = Counter(ans)
if counts.most_common(1)[0][1] / len(ans) >= threshold:
invalid.add(idx)
return invalid
def _detect_speeders(self, min_seconds):
duration_col = [c for c in self.df.columns
if '时长' in str(c) or '用时' in str(c)]
if not duration_col:
return set()
durations = pd.to_numeric(self.df[duration_col[0]], errors='coerce')
return set(self.df[durations < min_seconds].index)
def _detect_pattern(self):
invalid = set()
for idx, row in self.df.iterrows():
ans = row[self.question_cols].dropna().tolist()
if len(ans) < 5:
continue
try:
numeric = pd.array(ans, dtype='float')
diffs = np.diff(numeric)
if len(diffs) > 0 and np.all(diffs == diffs[0]) and diffs[0] != 0:
invalid.add(idx)
except:
pass
return invalid
def _detect_outliers(self, threshold):
numeric_df = self.df[self.question_cols].apply(
pd.to_numeric, errors='coerce')
row_means = numeric_df.mean(axis=1)
z_scores = (row_means - row_means.mean()) / row_means.std()
return set(self.df[abs(z_scores) > threshold].index)
@property
def all_invalid(self):
result = set()
for s in self.invalid_sets.values():
result |= s
return result
@property
def clean_data(self):
return self.df[~self.df.index.isin(self.all_invalid)].copy()
def report(self):
total = len(self.df)
invalid_count = len(self.all_invalid)
print(f"\n{'='*50}")
print(f"问卷数据清洗报告")
print(f"{'='*50}")
print(f"原始样本量: {total} 份")
for reason, indices in self.invalid_sets.items():
print(f" {reason}: {len(indices)} 份")
print(f"{'─'*50}")
print(f"无效样本合计: {invalid_count} 份")
print(f"有效样本: {total - invalid_count} 份")
print(f"有效率: {(1 - invalid_count/total)*100:.1f}%")
print(f"{'='*50}")
def export(self, output_path='清洗结果'):
clean = self.clean_data
clean.to_excel(f'{output_path}_有效数据.xlsx', index=False)
invalid_df = self.df[self.df.index.isin(self.all_invalid)].copy()
invalid_df['无效原因'] = ''
for reason, indices in self.invalid_sets.items():
invalid_df.loc[invalid_df.index.isin(indices), '无效原因'] += f'{reason};'
invalid_df.to_excel(f'{output_path}_无效样本.xlsx', index=False)
# 使用方法
if __name__ == '__main__':
cleaner = SurveyDataCleaner('问卷数据.xlsx')
cleaner.run_all(min_seconds=15, z_threshold=3.0)
cleaner.report()
cleaner.export()
关于样本质量的一个延伸
清洗完数据之后,很多人会发现一个问题:有效样本量不够了。
比如原本收了 300 份,洗完只剩 200 份,但你的论文需要 250 份以上的有效样本。这时候要么继续补数据,要么回头看看无效样本率是不是太高了。
如果无效样本率超过 20%,通常说明数据来源本身有问题------要么是分发渠道不合适,要么是激励方式导致乱填。
这也是为什么在做学术研究时,样本来源的选择很重要。相比在互填群里收回大量低质量数据,问卷鸭 这种真人分发渠道的无效样本率通常会更低,因为它走的是真实用户逐份填写,而不是批量自动化。价格大概是 0.2 元一份,如果你算一下重新收集数据的时间成本,其实并不贵。
不过不管数据来源是什么,拿到手之后该做的清洗一步都不能少。上面的代码覆盖了最常见的清洗场景,如果需要做信效度分析(Cronbach's Alpha、KMO 检验等),可以结合 factor_analyzer 库继续做。
相关资料
- 问卷星数据导出格式说明
- pandas 官方文档
- Cronbach's Alpha 计算:
pip install pingouin,然后用pg.cronbach_alpha(df[question_cols]) - KMO 检验:
pip install factor_analyzer,然后用calculate_kmo(df[question_cols])
数据清洗代码仓库:survey-data-cleaner(整理中)