当前位置: 代码网 > it编程>前端脚本>Python > 使用Python做问卷数据清洗的完整流程与完整代码

使用Python做问卷数据清洗的完整流程与完整代码

2026年09月01日 Python 我要评论
做问卷研究的人都知道,回收数据里一定有垃圾。直线作答、秒填、规律选择、前后矛盾……这些数据如果不清理掉,分析结果会被严重污染,信效度也过不了关。这篇教程教你用 python

做问卷研究的人都知道,回收数据里一定有垃圾。直线作答、秒填、规律选择、前后矛盾……这些数据如果不清理掉,分析结果会被严重污染,信效度也过不了关。

这篇教程教你用 python 从零写一套问卷数据清洗流程,覆盖最常见的无效样本类型。代码可以直接拿去用。

适用场景: 毕业论文问卷、学术调研、市场研究、用户满意度调查等。

准备工作

你需要:

  • python 3.8+
  • pandasnumpymatplotlib
  • 一份导出为 csv 或 excel 的问卷数据
pip install pandas numpy matplotlib openpyxl

问卷星导出的数据一般是这个结构(每行一份问卷,每列一个题目):

序号提交时间q1q2q3q4q5
12025-03-15 14:30:02满意一般满意满意满意
22025-03-15 14:30:05非常不满意非常不满意非常不满意非常不满意非常不满意

不用多说,第 2 行明显是垃圾数据。但实际场景里无效样本通常没这么明显。

第一步:加载数据

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())

通常问卷星导出的列名比较乱,先整理一下:

# 找出题目列(通常包含 "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)

直线作答是最常见的无效模式——所有题目选了同一个选项。

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% 以上选项相同

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)} 份")

把每一份问卷画成一条折线,无效模式其实肉眼就能认出来——直线作答是一条平的,
规律作答是锯齿:

上面那张图是拿 300 份合成数据跑本文的检测函数画出来的,蓝线(直线作答)和
橙线(1-2-3-4-5 循环)就是要剔的两类,绿红两条是正常样本的样子。

第三步:检测秒填(speeders)

如果有人在 10 秒内"填完"了一份 20 题的问卷,基本可以确定是无效数据。

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)} 份")

把答题时长画成直方图,秒填的那一簇会自己跳出来:

正常作答会落在一个比较宽的区间里(这份数据是 60~400 秒),秒填集中在最左端,
中间基本是空的——这个「断层」就是阈值该放的位置,比拍脑袋定 10 秒靠谱。

如果问卷星导出的数据里包含答题时长列,直接按阈值过滤就行。一般规则:

  • 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。

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 问"您对产品的哪些功能最满意"却填了具体内容,这就是矛盾回答。

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 分),可以用统计方法检测异常样本。

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)} 份")

第七步:汇总所有无效样本

# 合并所有检测结果
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)} 份")

第八步:可视化清洗结果

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()

跑出来长这样:

这份 300 样本的合成数据里,直线作答 14 份、规律作答 23 份、秒填 31 份,
去重后无效样本共 58 份,有效率 80.7%。注意三类是有重叠的——
一个人可能既秒填又直线作答,所以柱状图三根加起来会大于饼图里的无效数,
汇总时必须用集合取并集,不能直接相加。

完整代码(一键运行)

把上面的函数整合成一个类,方便复用:

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 库继续做。

以上就是使用python做问卷数据清洗的完整流程与完整代码的详细内容,更多关于python做问卷数据清洗的资料请关注代码网其它相关文章!

(0)

相关文章:

版权声明:本文内容由互联网用户贡献,该文观点仅代表作者本人。本站仅提供信息存储服务,不拥有所有权,不承担相关法律责任。 如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件至 2386932994@qq.com 举报,一经查实将立刻删除。

发表评论

验证码:
Copyright © 2017-2026  代码网 保留所有权利. 粤ICP备2024248653号
站长QQ:2386932994 | 联系邮箱:2386932994@qq.com