来源:丁点帮你
作者:丁点helper
之前的文章讲了如何用R绘制箱形图,以此来帮助我们了解数据的整体分布情况、是否存在异常值。除此之外,箱形图还可以进行数据的组间比较。
多重线性回归,一般是指有多个自变量X,只有一个因变量Y。前面我们主要是以简单线性回归为例在介绍,两者的差距主要在于自变量X的数量,在只有一个X时,就称简单线性回归。
本次我们用到的是学生的课堂调查数据,包括了性别、年级、专业、身高、最喜欢的动物(讲数据清理时用的是这个变量,还记得吗)等变量。数据名:survey.csv,数据链接:
首先导入数据,存入survey这个数据集中:
survey <- read.csv("//Users//Desktop//survey.csv",
header = TRUE) # 获取数据中包含的变量名 names(survey)
[1] "ClassProb" "Status" "Year" "Division" "Gender" "HtCm" "Hand" "Haircut" "Exercise" [10] "Coursework" "Web" "TV" "Social" "Econ" "Animal" "Friends" "Pulse"
接下来我们以Gender作为分组依据,先来看看这个变量的情况。
table(survey$Gender) Choose not to answer Female Gender non-conforming Male 1 1 117 1 118
我们发现,除了female和male,有的同学回答了Choose not to answer,Gender non-conforming,还有同学什么都没填,空缺。
今天我们暂时将这三种特殊情况从数据中删去。
# 查看针对Gender这个变量,同学们有几类回答 levels(survey$Gender)
[1] "" "Choose not to answer" "Female" "Gender non-conforming" "Male"
在这五类回答中,我们想保留的是第3、第5类。也就是说,仅保留Gender为"Female" 或 "Male"的记录。
# 把更新后的数据存储在survey2这个对象中 survey2 <- survey[survey$Gender %in% levels(survey$Gender)[c(3,5)],]
这里,a %in%b的作用是,用a中的元素去匹配b中的任意元素,如果匹配成功,则返回结果为TRUE,反之,则结果为FALSE。
此时,上面的语句就简化为如下所示,c()里面是TRUE和FALSE的集合,是a中每个元素与b匹配的结果。
survey2 <- survey[c(),] # 这是为了便于理解写的简化语句,不能够运行的
survey2中仅保留了匹配结果为TRUE的记录:
table(survey2$Gender) Choose not to answer Female Gender non-conforming Male 0 0 117 0 118
哎?虽然记录被删了,但Gender中之前包含的五个类都还在,用下面的droplevels()这个函数删掉那些没有记录的类。
survey2$Gender <- droplevels(survey2$Gender)
table(survey2$Gender)
Female Male
117 118
多重线性回归,一般是指有多个自变量X,只有一个因变量Y。前面我们主要是以简单线性回归为例在介绍,两者的差距主要在于自变量X的数量,在只有一个X时,就称简单线性回归。
数据清理好之后,我们以身高HtCm这个变量为例,先用之前讲过的方法绘制箱形图,了解改变量的整体分布,然后对比性别之间的身高差异。
boxplot(survey2$HtCm, main="Boxplot of Ht in cm", col='orange', lwd=2)
一目了然,我们调查的是大学学生,却出现了身高小于100厘米的情况,不符合常理。现在去检查一下原始数据。
sort(survey2$HtCm) # 将身高从小到大排序
部分结果截图
实际操作中,大家要尽量核实那些极端身高数据的真实情况,修正数据。这里我们为便于教学,直接把那些小于100厘米的身高值记录为缺失。
然后利用整理后的身高数据绘制箱形图。
survey2$HtCm[survey2$HtCm < 100 ] <- NA
boxplot(survey2$HtCm, main="Boxplot of Ht in cm",
col='orange', lwd=2)
最后绘制不同性别学生的身高箱形图。
boxplot(survey2$HtCm~survey2$Gender,
main="Boxplot of Ht in cm",
col=c(2,3), lwd=2)
由图可知,男生的身高基本都高于女生。将两个箱形图放在一起,可以清晰地看到两组变量的大致情况,便于给两组做粗略的比较。
但是这男女生身高到底有没有统计学上的差异,肉眼是很难得出结论的,统计学上怎么做呢?
数据分析咨询请扫描二维码
若不方便扫码,搜微信号:CDAshujufenxi
大咖简介: 刘凯,CDA大咖汇特邀讲师,DAMA中国分会理事,香港金管局特聘数据管理专家,拥有丰富的行业经验。本文将从数据要素 ...
2025-04-22CDA持证人简介 刘伟,美国 NAU 大学计算机信息技术硕士, CDA数据分析师三级持证人,现任职于江苏宝应农商银行数据治理岗。 学 ...
2025-04-21持证人简介:贺渲雯 ,CDA 数据分析师一级持证人,互联网行业数据分析师 今天我将为大家带来一个关于用户私域用户质量数据分析 ...
2025-04-18一、CDA持证人介绍 在数字化浪潮席卷商业领域的当下,数据分析已成为企业发展的关键驱动力。为助力大家深入了解数据分析在电商行 ...
2025-04-17CDA持证人简介:居瑜 ,CDA一级持证人,国企财务经理,13年财务管理运营经验,在数据分析实践方面积累了丰富的行业经验。 一、 ...
2025-04-16持证人简介: CDA持证人刘凌峰,CDA L1持证人,微软认证讲师(MCT)金山办公最有价值专家(KVP),工信部高级项目管理师,拥有 ...
2025-04-15持证人简介:CDA持证人黄葛英,ICF国际教练联盟认证教练,前字节跳动销售主管,拥有丰富的行业经验。在实际生活中,我们可能会 ...
2025-04-14在 Python 编程学习与实践中,Anaconda 是一款极为重要的工具。它作为一个开源的 Python 发行版本,集成了众多常用的科学计算库 ...
2025-04-14随着大数据时代的深入发展,数据运营成为企业不可或缺的岗位之一。这个职位的核心是通过收集、整理和分析数据,帮助企业做出科 ...
2025-04-11持证人简介:CDA持证人黄葛英,ICF国际教练联盟认证教练,前字节跳动销售主管,拥有丰富的行业经验。 本次分享我将以教培行业为 ...
2025-04-11近日《2025中国城市长租市场发展蓝皮书》(下称《蓝皮书》)正式发布。《蓝皮书》指出,当前我国城市住房正经历从“增量扩张”向 ...
2025-04-10在数字化时代的浪潮中,数据已经成为企业决策和运营的核心。每一位客户,每一次交易,都承载着丰富的信息和价值。 如何在海量客 ...
2025-04-09数据是数字化的基础。随着工业4.0的推进,企业生产运作过程中的在线数据变得更加丰富;而互联网、新零售等C端应用的丰富多彩,产 ...
2025-04-094月7日,美国关税政策对全球金融市场的冲击仍在肆虐,周一亚市早盘,美股股指、原油期货、加密货币、贵金属等资产齐齐重挫,市场 ...
2025-04-08背景 3月26日,科技圈迎来一则重磅消息,苹果公司宣布向浙江大学捐赠 3000 万元人民币,用于支持编程教育。 这一举措并非偶然, ...
2025-04-07在当今数据驱动的时代,数据分析能力备受青睐,数据分析能力频繁出现在岗位需求的描述中,不分岗位的任职要求中,会特意标出“熟 ...
2025-04-03在当今数字化时代,数据分析师的重要性与日俱增。但许多人在踏上这条职业道路时,往往充满疑惑: 如何成为一名数据分析师?成为 ...
2025-04-02最近我发现一个绝招,用DeepSeek AI处理Excel数据简直太爽了!处理速度嘎嘎快! 平常一整天的表格处理工作,现在只要三步就能搞 ...
2025-04-01你是否被统计学复杂的理论和晦涩的公式劝退过?别担心,“山有木兮:统计学极简入门(Python)” 将为你一一化解这些难题。课程 ...
2025-03-31在电商、零售、甚至内容付费业务中,你真的了解你的客户吗? 有些客户下了一两次单就消失了,有些人每个月都回购,有些人曾经是 ...
2025-03-31