或许我们应该经常去泡馆、经常去图书馆借书看

或许我们应该经常去泡馆、经常去图书馆借书看

昨天熬夜写了那个查询我们学校同学的图书馆数据,今天让我们来一起用一下。

首先准备一份成绩单

这个成绩单是之前从数字暨大爬来的,出于隐私保护的目的我就不附上这份数据集的链接了,需要的朋友可以在留言去留言所要。这份成绩单有这些变量:

进馆次数和借阅图书的本数是后来生成的。

根据每个人的姓名和学号查询其进馆次数和借阅图书的本数

下面的Stata代码就能实现这一目标:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
clear all
cpd //设定工作目录为桌面
use cjd, clear
gen jgcs = ""
gen jsbs = ""
save temp, replace
forval i = 1(2)`=_N'{
use temp, clear
local xm = xm[`i']
local xh = xh[`i']
cap libdata `xm' `xh' , s(000) u(0)
local a = r(总进馆次数)
local b = r(总借阅图书的数量)
use temp, clear
local j = `i' + 1
replace jgcs = "`a'" in `i'/`j'
replace jsbs = "`b'" in `i'/`j'
save temp, replace
}
use temp, clear
order xm xh jgcs jsbs jd
destring, replace
label var jgcs "进馆次数"
label var jsbs "借阅图书的本数"
save cjd+tsg, replace

进行一些简单的相关分析

下面的图表分析虽然不能严格的说明因果关系,其中存在很多内生的干扰,例如双向因果。进馆次数多的人自然绩点更可能高,但是也可能是因为其绩点高所以才喜欢泡馆。因此下面不会过多的追究因果关系,仅仅是为了描述一下事实。

其次是数据的限制的问题,绩点、综测、思想测评成绩有两年的数据,但是进馆次数和借阅次数都是入校以来的总数。

综合上面的问题,下面的图表仅供调侃。

绩点与进馆次数

1
2
3
4
5
6
7
8
9
use cjd+tsg, clear
tw ///
sc jd jgcs if txmc == "2017年综合测评", mc(red*0.6) || ///
lfit jd jgcs if txmc == "2017年综合测评", lc(red*0.6) lp(solid) || ///
sc jd jgcs if txmc == "2016年综合测评", mc(green*0.6) || ///
lfit jd jgcs if txmc == "2016年综合测评", lc(green*0.6) lp(solid) ||, ///
xti("进馆次数") yti("绩点") ///
yla(, format(%6.2f)) leg(order(4 "2016年" 2 "2017年"))
gre g1


可见,泡馆果然有助于学习!

绩点与借阅次数

1
2
3
4
5
6
7
8
9
tw ///
sc jd jsbs if txmc == "2017年综合测评", mc(red*0.6) || ///
lfit jd jsbs if txmc == "2017年综合测评", lc(red*0.6) lp(solid) || ///
sc jd jsbs if txmc == "2016年综合测评", mc(green*0.6) || ///
lfit jd jsbs if txmc == "2016年综合测评", lc(green*0.6) lp(solid) ||, ///
yti("绩点") xti("借阅图书的本数") ///
yla(, format(%6.2f)) leg(order(4 "2016年" 2 "2017年"))
gre g2
`


多借点书,看不看不重要,只要借就有助于绩点的提升!

综测得分与进馆次数

1
2
3
4
5
6
7
8
9
use cjd+tsg, clear
tw ///
sc zc jgcs if txmc == "2017年综合测评", mc(red*0.6) || ///
lfit zc jgcs if txmc == "2017年综合测评", lc(red*0.6) lp(solid) || ///
sc zc jgcs if txmc == "2016年综合测评", mc(green*0.6) || ///
lfit zc jgcs if txmc == "2016年综合测评", lc(green*0.6) lp(solid) ||, ///
xti("进馆次数") yti("综测") ///
yla(, format(%6.2f)) leg(order(4 "2016年" 2 "2017年"))
gre g3


综测一直是一个神奇的东西,从趋势来看,爱泡馆的同学综测会提升,但是不爱泡馆综测会下滑!

综测得分与借阅次数

1
2
3
4
5
6
7
8
9
use cjd+tsg, clear
tw ///
sc zc jsbs if txmc == "2017年综合测评", mc(red*0.6) || ///
lfit zc jsbs if txmc == "2017年综合测评", lc(red*0.6) lp(solid) || ///
sc zc jsbs if txmc == "2016年综合测评", mc(green*0.6) || ///
lfit zc jsbs if txmc == "2016年综合测评", lc(green*0.6) lp(solid) ||, ///
xti("借阅次数") yti("综测") ///
yla(, format(%6.2f)) leg(order(4 "2016年" 2 "2017年"))
gre g4


不爱借书的同学,综测会大幅下滑!

思想测评成绩与借阅次数

1
2
3
4
5
6
7
8
9
use cjd+tsg, clear
tw ///
sc sxpd jsbs if txmc == "2017年综合测评", mc(red*0.6) || ///
lfit sxpd jsbs if txmc == "2017年综合测评", lc(red*0.6) lp(solid) || ///
sc sxpd jsbs if txmc == "2016年综合测评", mc(green*0.6) || ///
lfit sxpd jsbs if txmc == "2016年综合测评", lc(green*0.6) lp(solid) ||, ///
xti("借阅次数") yti("思想测评成绩") ///
yla(, format(%6.2f)) leg(order(4 "2016年" 2 "2017年"))
gre g5


不爱借书的同学思想会滑坡!

思想测评成绩与进馆次数

1
2
3
4
5
6
7
8
9
use cjd+tsg, clear
tw ///
sc sxpd jgcs if txmc == "2017年综合测评", mc(red*0.6) || ///
lfit sxpd jgcs if txmc == "2017年综合测评", lc(red*0.6) lp(solid) || ///
sc sxpd jgcs if txmc == "2016年综合测评", mc(green*0.6) || ///
lfit sxpd jgcs if txmc == "2016年综合测评", lc(green*0.6) lp(solid) ||, ///
xti("进馆次数") yti("思想测评成绩") ///
yla(, format(%6.2f)) leg(order(4 "2016年" 2 "2017年"))
gre g6


不爱泡馆的同学,思想严重滑坡!

# Stata

评论

程振兴

程振兴 @czxa.top
截止今天,我已经在本博客上写了607.9k个字了!

Your browser is out-of-date!

Update your browser to view this website correctly. Update my browser now

×