diff --git a/Book/00-author.Rmd b/Book/00-author.Rmd index 1b4d51ae..19060d35 100644 --- a/Book/00-author.Rmd +++ b/Book/00-author.Rmd @@ -6,7 +6,7 @@ 本电子书由传鹏给南京师范大学心理学院研究生所开设的《R语言在心理学研究中的应用》课程转录而来。转录过程得到如下同学的大力帮助:\ -## 2024春季学期助教: +## 2024春季学期助教及电子书转录志愿者: - 陈钟奇,南京师范大学心理学院博士生; - 冯雨萌,香港城市大学社会与行为科学系硕士生; diff --git a/Book/1016-lesson16_2023.Rmd b/Book/1016-lesson16_2023.Rmd index ae6d3d22..76af3625 100644 --- a/Book/1016-lesson16_2023.Rmd +++ b/Book/1016-lesson16_2023.Rmd @@ -127,7 +127,7 @@ output : papaja::apa6_pdf (2)包和文件的调用: 在定义YAML头文件之后,需要进行检查是否该文件已经调用了需要使用的包,以及是否包含了需要的文件(例如.bib格式文件)。 (3)参考文献的引用方式 -如例子中**[@R-papaja]**,这是Papaja的引用格式,包含中括号[]、@符号以及参考文献标识。参考文献标识为.bib文件中@后面的内容。这种标识转换后为APA标准引用,以括号的形式引用文献。 +如例子中**[@R-papaja]**,这是`papaja`的引用格式,包含中括号[]、\@符号以及参考文献标识。参考文献标识为`.bib`文件中\@后面的内容。这种标识转换后为APA标准引用,以括号的形式引用文献。  还有我们通常使用的另一种引用,即在句子中直接引用文章作者,此时只需要略去大括号即可,如@Hu_2020 diff --git a/Book/_book/_main_files/figure-html/unnamed-chunk-255-1.png b/Book/_book/_main_files/figure-html/unnamed-chunk-255-1.png index a592fd91..59843991 100644 Binary files a/Book/_book/_main_files/figure-html/unnamed-chunk-255-1.png and b/Book/_book/_main_files/figure-html/unnamed-chunk-255-1.png differ diff --git a/Book/_book/_main_files/figure-html/unnamed-chunk-275-1.png b/Book/_book/_main_files/figure-html/unnamed-chunk-275-1.png index c6d3f2cb..259ff0c7 100644 Binary files a/Book/_book/_main_files/figure-html/unnamed-chunk-275-1.png and b/Book/_book/_main_files/figure-html/unnamed-chunk-275-1.png differ diff --git a/Book/_book/_main_files/figure-html/unnamed-chunk-277-1.png b/Book/_book/_main_files/figure-html/unnamed-chunk-277-1.png index f1fb41ac..d3b22dda 100644 Binary files a/Book/_book/_main_files/figure-html/unnamed-chunk-277-1.png and b/Book/_book/_main_files/figure-html/unnamed-chunk-277-1.png differ diff --git a/Book/_book/_main_files/figure-html/unnamed-chunk-278-1.png b/Book/_book/_main_files/figure-html/unnamed-chunk-278-1.png index a563a47b..7bc8452c 100644 Binary files a/Book/_book/_main_files/figure-html/unnamed-chunk-278-1.png and b/Book/_book/_main_files/figure-html/unnamed-chunk-278-1.png differ diff --git a/Book/_book/_main_files/figure-html/unnamed-chunk-288-1.png b/Book/_book/_main_files/figure-html/unnamed-chunk-288-1.png index 0352fd7b..d4d977e5 100644 Binary files a/Book/_book/_main_files/figure-html/unnamed-chunk-288-1.png and b/Book/_book/_main_files/figure-html/unnamed-chunk-288-1.png differ diff --git a/Book/_book/_main_files/figure-html/unnamed-chunk-289-1.png b/Book/_book/_main_files/figure-html/unnamed-chunk-289-1.png index 39df407a..8ad22ad5 100644 Binary files a/Book/_book/_main_files/figure-html/unnamed-chunk-289-1.png and b/Book/_book/_main_files/figure-html/unnamed-chunk-289-1.png differ diff --git a/Book/_book/index.html b/Book/_book/index.html index b4c910f1..cd2c4bf3 100644 --- a/Book/_book/index.html +++ b/Book/_book/index.html @@ -136,6 +136,29 @@ div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;} + @@ -162,8 +185,8 @@
2024年09月11日
+2024年09月12日
本电子书由传鹏给南京师范大学心理学院研究生所开设的《R语言在心理学研究中的应用》课程转录而来。转录过程得到如下同学的大力帮助:
## [1] "date" "Prac" "Sub" "Age" "Sex" "Hand" "Block"
-## [8] "Bin" "Trial" "Shape" "Label" "Match" "CorrResp" "Resp"
-## [15] "ACC" "RT"
+## [1] "date" "Prac" "Sub" "Age" "Sex" "Hand" "Block" "Bin" "Trial" "Shape" "Label"
+## [12] "Match" "CorrResp" "Resp" "ACC" "RT"
## 将全部列名都变成小写
df3 %>% dplyr::rename_with(
## 将字符向量全部变成小写; ~ 声明这是一个函数,.代表前面的数据(df3)传到.所在的位置
@@ -606,9 +628,8 @@ 7.1.1 批量导入数据 ## 即使用 tolower()对所有列名进行批量处理
##
) %>% colnames()
## [1] "date" "prac" "sub" "age" "sex" "hand" "block"
-## [8] "bin" "trial" "shape" "label" "match" "corrresp" "resp"
-## [15] "acc" "rt"
+## [1] "date" "prac" "sub" "age" "sex" "hand" "block" "bin" "trial" "shape" "label"
+## [12] "match" "corrresp" "resp" "acc" "rt"
skimr::skim(mt_raw) %>%
capture.output() %>%
.[1:12]
## [1] "── Data Summary ────────────────────────" " Values"
-## [3] "Name mt_raw" "Number of rows 25920 "
-## [5] "Number of columns 16 " "_______________________ "
-## [7] "Column type frequency: " " character 9 "
-## [9] " numeric 7 " "________________________ "
-## [11] "Group variables None " ""
+## [1] "── Data Summary ────────────────────────" " Values" "Name mt_raw"
+## [4] "Number of rows 25920 " "Number of columns 16 " "_______________________ "
+## [7] "Column type frequency: " " character 9 " " numeric 7 "
+## [10] "________________________ " "Group variables None " ""
skimr::skim()–2
-## [1] "── Variable type: character ──────────────────────────────────────────────────────────────"
-## [2] " skim_variable n_missing complete_rate min max empty n_unique whitespace"
-## [3] "1 Date 0 1 20 20 0 24362 0"
-## [4] "2 Prac 0 1 3 3 0 1 0"
-## [5] "3 Sex 0 1 1 6 0 4 0"
-## [6] "4 Hand 0 1 1 1 0 2 0"
-## [7] "5 Shape 0 1 9 12 0 4 0"
-## [8] "6 Label 0 1 9 12 0 4 0"
-## [9] "7 Match 0 1 5 8 0 2 0"
-## [10] "8 CorrResp 0 1 1 1 0 2 0"
-## [11] "9 Resp 658 0.975 1 5 0 9 0"
+## [1] "── Variable type: character ─────────────────────────────────────────────────────────────────────────────────────────────────────────"
+## [2] " skim_variable n_missing complete_rate min max empty n_unique whitespace"
+## [3] "1 Date 0 1 20 20 0 24362 0"
+## [4] "2 Prac 0 1 3 3 0 1 0"
+## [5] "3 Sex 0 1 1 6 0 4 0"
+## [6] "4 Hand 0 1 1 1 0 2 0"
+## [7] "5 Shape 0 1 9 12 0 4 0"
+## [8] "6 Label 0 1 9 12 0 4 0"
+## [9] "7 Match 0 1 5 8 0 2 0"
+## [10] "8 CorrResp 0 1 1 1 0 2 0"
+## [11] "9 Resp 658 0.975 1 5 0 9 0"
## [12] ""
skimr::skim()–3
-## [1] "── Variable type: numeric ────────────────────────────────────────────────────────────────"
-## [2] " skim_variable n_missing complete_rate mean sd p0 p25 p50"
-## [3] "1 Sub 0 1 8853. 9932. 7302 7313 7324 "
-## [4] "2 Age 0 1 20.8 2.48 18 19 20 "
-## [5] "3 Block 0 1 1.61 0.803 1 1 1 "
-## [6] "4 Bin 0 1 2.42 1.36 1 1 2 "
-## [7] "5 Trial 0 1 12.5 6.92 1 6.75 12.5 "
-## [8] "6 ACC 0 1 0.796 0.464 -1 1 1 "
-## [9] "7 RT 0 1 0.715 0.151 0.106 0.610 0.702"
-## [10] " p75 p100 hist "
-## [11] "1 7336 73370 ▇▁▁▁▁"
-## [12] "2 22 28 ▇▂▁▁▁"
-## [13] "3 2 3 ▇▁▃▁▃"
-## [14] "4 3 5 ▇▇▃▃▃"
-## [15] "5 18.2 24 ▇▇▆▇▇"
-## [16] "6 1 2 ▁▂▁▇▁"
-## [17] "7 0.805 1.18 ▁▂▇▅▁"
+## [1] "── Variable type: numeric ───────────────────────────────────────────────────────────────────────────────────────────────────────────"
+## [2] " skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist "
+## [3] "1 Sub 0 1 8853. 9932. 7302 7313 7324 7336 73370 ▇▁▁▁▁"
+## [4] "2 Age 0 1 20.8 2.48 18 19 20 22 28 ▇▂▁▁▁"
+## [5] "3 Block 0 1 1.61 0.803 1 1 1 2 3 ▇▁▃▁▃"
+## [6] "4 Bin 0 1 2.42 1.36 1 1 2 3 5 ▇▇▃▃▃"
+## [7] "5 Trial 0 1 12.5 6.92 1 6.75 12.5 18.2 24 ▇▇▆▇▇"
+## [8] "6 ACC 0 1 0.796 0.464 -1 1 1 1 2 ▁▂▁▇▁"
+## [9] "7 RT 0 1 0.715 0.151 0.106 0.610 0.702 0.805 1.18 ▁▂▇▅▁"
+## [10] NA
+## [11] NA
+## [12] NA
+## [13] NA
+## [14] NA
+## [15] NA
+## [16] NA
+## [17] NA
bruceR::Describe()
bruceR::Describe(mt_raw) %>%
capture.output() %>%
diff --git "a/Book/_book/\347\254\254\344\270\211\350\256\262git-rstudio\345\267\245\344\275\234\346\265\201.html" "b/Book/_book/\347\254\254\344\270\211\350\256\262git-rstudio\345\267\245\344\275\234\346\265\201.html"
index 83cb7895..02cd381f 100644
--- "a/Book/_book/\347\254\254\344\270\211\350\256\262git-rstudio\345\267\245\344\275\234\346\265\201.html"
+++ "b/Book/_book/\347\254\254\344\270\211\350\256\262git-rstudio\345\267\245\344\275\234\346\265\201.html"
@@ -136,6 +136,29 @@
div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
+
@@ -162,8 +185,8 @@
- 关于作者
- 1 第一讲:为什么要学习R
@@ -215,7 +238,7 @@
- 3 第三讲:Git & RStudio工作流
-- 3.1 Files and Folders System
+- 3.1 文件与文件夹结构系统
- 3.2 Git and Git Hub
- 3.2.1 Git
@@ -541,21 +564,16 @@
Chapter 3 第三讲:Git & RStudio工作流
-
-3.1 Files and Folders System
-第三节课将继续入门R语言的学习。本课程的核心目标是帮助大家开始使用R,即便是已经接触过R的同学,也可以通过这节课对R进行更系统的学习。去年我们的R课程中,Git的学习是在课程后期进行的,那时同学们已经接触了一定量的代码,意识到代码管理的重要性。但今年,我们选择提前介绍Git,因为我们的所有PPT和课件更新都将通过Git
-Hub进行。如果不熟悉Git
-Hub操作,同学们在获取最新课件时可能会面临反复下载和删除的问题,这不仅耗时,也不利于大家的学习。因此,我们决定将Git和Git
-Hub的学习放在课程的前面,让大家尽早掌握这些必要的技能。
-本次课的内容主要分为四个方面:首先是数据分析中如何更好地管理数据文件;其次是Git
-Hub的基本介绍和主要特点;接着是Git的基本功能和版本控制原理;最后是Git的远程文件版本控制和如何在远程代码托管平台或云平台上进行交互。这些技能对于团队合作和项目管理非常重要,尤其是在大型项目中。
+
+3.1 文件与文件夹结构系统
+第三节课将继续入门R语言的学习。本课程的核心目标是帮助大家开始使用R,即便是已经接触过R的同学,也可以通过这节课对R进行更系统的学习。去年我们的R课程中,在课程后期才开始学习Git。但今年我们选择提前介绍Git,因为所有PPT和课件的更新都会即时放在GitHub上。如果不熟悉GitHub操作,获取最新课件可能会需要反复下载文件的问题。这不仅耗时,也不利于大家的学习。因此,今年提前将Git和GitHub作为基础进行学习,让大家尽早掌握这些必要的技能。
+本讲的内容分为四个方面:首先是数据分析中如何更好地管理数据文件;其次是GitHub的基本介绍和主要特点;接着是Git的基本功能和版本控制原理;最后是Git的远程文件版本控制和如何在远程代码托管平台或云平台上进行交互。这些技能对于团队合作和项目管理非常重要,尤其是在大型项目中。
我们的目标有两个:一是提高大家的项目管理和版本控制技能,以增强学习效率、研究能力和团队协作能力;二是即便短期内不需要频繁编写代码的同学,也应该学会如何下载更新课件和完成第一次小作业。这些技能对于保持项目最新性和协作顺畅性至关重要。
在日常工作中,我们经常看到不规范的数据管理做法,如所有类型的文档和数据被随意放置在一个根目录下,没有分类。这种做法虽然能找到所需的数据和图表,但查找过程繁琐,效率低下。

-为了解决这个问题,建议对文件夹进行分类管理,并根据数据和文件的不同目的和使用方式,将它们分成不同的类别。此外,添加一个README文件,描述每个文件夹的功能和内容,以便未来的自己或其他研究者理解。
+为解决这个问题,我们建议对文件夹进行分类管理,并根据数据和文件的不同目的和使用方式,将它们分成不同的类别。此外,添加一个README文件,描述每个文件夹的功能和内容,以便未来的自己或其他研究者理解。

-我们的课题组采用了两种文件管理模式:针对实证研究的文件管理和针对元分析或元研究的文件管理。这些模式在我们的Git
-Hub仓库中有详细的模板,其他研究者可以参考。通过特定的网站,我们发现借鉴规范的文件夹结构模式对研究者有益,也可以向导师推荐。网址是:https://psych-transparency-guide.uni-koeln.de/folder-structure.html
+我们的课题组采用两种文件管理模式:针对实证研究的文件管理和针对元分析或元研究的文件管理。这些模式在我们的GitHub仓库中有详细的模板,大家可以参考。规范的文件夹结构模式对研究者是大大有益的,欢迎大家向自己课题组推荐。网址是:https://psych-transparency-guide.uni-koeln.de/folder-structure.html
在文件管理方面,我们发现一个常见问题是文件夹结构看起来整洁,但打开后发现文件组织混乱,包括文件命名问题。例如,对于毕业论文,应该在文件名中加上作者的名字,以便导师清楚地知道哪个文件是需要修改的。在数据分析过程中,版本控制非常重要,它可以帮助我们管理文件的各个版本,避免丢失重要的修改。
diff --git "a/Book/_book/\347\254\254\344\271\235\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\344\272\214\345\210\206\345\261\202\347\272\277\346\200\247\346\250\241\345\236\213.html" "b/Book/_book/\347\254\254\344\271\235\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\344\272\214\345\210\206\345\261\202\347\272\277\346\200\247\346\250\241\345\236\213.html"
index 3bfbe033..0207ff92 100644
--- "a/Book/_book/\347\254\254\344\271\235\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\344\272\214\345\210\206\345\261\202\347\272\277\346\200\247\346\250\241\345\236\213.html"
+++ "b/Book/_book/\347\254\254\344\271\235\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\344\272\214\345\210\206\345\261\202\347\272\277\346\200\247\346\250\241\345\236\213.html"
@@ -136,6 +136,29 @@
div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
+
@@ -162,8 +185,8 @@
- 关于作者
- 1 第一讲:为什么要学习R
@@ -215,7 +238,7 @@
- 3 第三讲:Git & RStudio工作流
-- 3.1 Files and Folders System
+- 3.1 文件与文件夹结构系统
- 3.2 Git and Git Hub
- 3.2.1 Git
@@ -670,8 +693,8 @@ 9.2 重复测量方差分析 regex = "(moral|immoral)(Self|Other)",
remove = FALSE)
[数据展示]
-
-
+
+
以match_raw.csv为例,一个2x2的被试内实验设计(Identity:Self vs.Other) x Valence:Moral vs.Immoral)),我们希望知道这两种条件之下被试的反应时是否存在显著差异然后我们就发现,这个数据实际上是有很多的,每个被试还有很多个试次,那么这种情况的话,大家通常做法是怎么做呢?如果通常我们用常规的这个重复测量方差分析,大家怎么做呢?比方说有40多个被试,然后有四种条件,我们最后会算出每个被试在每种条件下,比方说我们关心反应时间的话,就是反应时间的一个均值,那么最后我们得到比方说44个被试他们的这个反应时间的均值,然后呢我们就把它输入到SPSS里面,然后进行一个重复测量方差分析,然后把它对应好,这是我们常规的做法,它数理结构基本上就是这样的。每一个被试有自我有他人,自我和他人下面又有moral和immoral两个条件,other的也有moral和immoral。
这样的话,如果说我们在本科阶段我们讲方差分析的时候,我们就会告诉大家,这里面会进行方差的分解,我们把它分解为,不同条件之间的变异,或者不同的自变量引起的变异,以及这个被试的个体差异,我们主要关注的,比方说就是这个自变量它引起的变异在总体的变异当中的比值,然后根据这个比值去计算f值等等等等,计算MSE就是谁的,这是本科的或大概考研的时候,可能会涉及到内容,那么在SPSS里面,大家可能也看到过,类似的这个实现的方式,那么在r里面也很方便的进行实现.

@@ -682,8 +705,7 @@ 9.2 重复测量方差分析 dplyr::group_by(Sub,Identity,Valence) %>%
dplyr::summarise(RT = mean(RT)) %>%
dplyr::ungroup()
-## `summarise()` has grouped output by 'Sub', 'Identity'. You can override using the
-## `.groups` argument.
+## `summarise()` has grouped output by 'Sub', 'Identity'. You can override using the `.groups` argument.
## 本例为长数据
## RUN IN CONSOLE!
## 球形检验输出:
@@ -783,8 +805,8 @@ 9.3 分层线性模型/多层线
那么可能在有一些领域,它可能根据多层线性模型的思路,它发展出了一些新的,特异于处理某一些特定类型的数据的方法,那么它还会给它一些叠加一些其他的名字,最重要的就是你要看到它本质是什么,它是不是比方说用以线性模型,或是广义线性模型,广义线性模型我们下节课会讲以线性模型作为最核心的一个模型,然后去考虑他的这个数据的层级结构,或层次结构,如果他有做这样做,基本上你就可以确定,他的原理上,可能就是跟我们这里讲的是差不多的,这里跟多元回归是完全不一样,多元回归是说有多个自变量,但他是没有考虑这个层级的结构的,那么在这个多层线型模型或者层级线型模型当中,有两个很重要的概念,我们先给大家简单的说一下,因为我们这里的层级模型,它都是以回归模型,就是以正态分布为核心的,那么在这种现行的层级模型当中,一般我们会考虑,从截距和斜率这两个效应上面去考虑我们这些效应,那么层级模型当中,它最关注两个效应就是固定效应和随机效应,一个叫fixed effect,一个叫random effect,这个所谓的固定和随机,这个名字本身非常不好理解,那我这里也不去跟大家把这个做非常细致的展开,大家如果感兴趣的话,可以去看知乎上的一个博客,包寒吴霜老师写的,那个博客对随机效应做了一个比较详尽的一个梳理,那么我们就通过一个例子来给大家展示一下,什么叫做固定效应和随机效应,那比方说我们看一个非常简单明了的一个层级数据,比方说工龄和薪水之间的一个关系,我们想要调查高校老师的工龄和工资之间有没有关系,那么从某一个学校里面,随机抽取出5个学院,然后获得他们的工资和工龄之间的关系,这是一个网上的数据,这里有数据来源,大家如果拿到rmarkdown,可以点击这个来源,然后这个数据结构大家可以看到,它基本上是这么一个嵌套的结构,首先你整个大学或者你整个学校,然后它有不同的departments,它有不同的学院或者系,每个系下面有不同的人,你调查的时候,实际上就是在不同的学院下面,去搜不同的人,然后再去从这个每个人身上,找到两个数据,一个是他工作多少年,另外一个是他的这个工资,然后你关心的是工资,和他的工作年限之间关系。
在回归模型中一般会在截距和斜率上分别讨论固定效应和随机效应。
例如,关于研究教师的工龄(Experience)与薪水(Salary)之间是否存在关系。在某校随机抽取了5个学院的教师信息,具体数据如下:]
-
-
+
+
数据来源见(https://github.com/mkfreeman/hierarchical-models/blob/master/generate-data.R)
问题:是否可用工龄预测某个学校员工的工资?
数据结构
@@ -952,12 +974,9 @@
9.4 多层线性模型的应用
[固定效应]
-## [1] "Fixed effects:"
-## [2] " Estimate Std. Error t value"
-## [3] "(Intercept) 0.72328 0.00997 72.51"
-## [4] "IdentitySelf 0.01329 0.00624 2.13"
-## [5] "Valencemoral -0.00913 0.00744 -1.23"
-## [6] "IdentitySelf:Valencemoral -0.04150 0.01128 -3.68"
+## [1] "Fixed effects:" " Estimate Std. Error t value"
+## [3] "(Intercept) 0.72328 0.00997 72.51" "IdentitySelf 0.01329 0.00624 2.13"
+## [5] "Valencemoral -0.00913 0.00744 -1.23" "IdentitySelf:Valencemoral -0.04150 0.01128 -3.68"
比如说在这个模式里面,我们应该是可以把每一个被试在identity上的效应提取出来,看它在被试间是怎样一个变化,以及比方说valence的效应,或者它们之间的交互作用,我们也可以对它进行一个可视化,这个可视化是个非常简陋的可视化,这里全部采用的是默认值,比方说顺序就是他人和自我,道德也是moral和immoral,跟我们之前呈现的顺序不一样,在讲完了这个统计模型之后,我们会专门讲怎么在gPlot里面,对这种最后要呈现的结果进行精细的打磨。那么实际上在现在有很多包,包括像interaction这样的包,它能够帮助我们迅速地把一些我们关注的效应进行可视化,让我们看到它们之间是什么样的情况,比方说这里我们可以看到是一个比较明显的交互作用,在other的条件之下的immoral和moral的差别,肯定是要比在自我条件下的差别是要小的,所以在这里我们可以非常明显的看到交互作用的存在。
[交互效应的可视化]
## 一种快捷的方法
diff --git "a/Book/_book/\347\254\254\344\272\214\350\256\262\345\246\202\344\275\225\345\274\200\345\247\213\344\275\277\347\224\250r.html" "b/Book/_book/\347\254\254\344\272\214\350\256\262\345\246\202\344\275\225\345\274\200\345\247\213\344\275\277\347\224\250r.html"
index 6ad2ef90..b2453561 100644
--- "a/Book/_book/\347\254\254\344\272\214\350\256\262\345\246\202\344\275\225\345\274\200\345\247\213\344\275\277\347\224\250r.html"
+++ "b/Book/_book/\347\254\254\344\272\214\350\256\262\345\246\202\344\275\225\345\274\200\345\247\213\344\275\277\347\224\250r.html"
@@ -136,6 +136,29 @@
div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
+
@@ -162,8 +185,8 @@
- 关于作者
- 1 第一讲:为什么要学习R
@@ -215,7 +238,7 @@
- 3 第三讲:Git & RStudio工作流
-- 3.1 Files and Folders System
+- 3.1 文件与文件夹结构系统
- 3.2 Git and Git Hub
- 3.2.1 Git
diff --git "a/Book/_book/\347\254\254\345\205\253\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\344\270\200.html" "b/Book/_book/\347\254\254\345\205\253\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\344\270\200.html"
index b9f17573..ad43e522 100644
--- "a/Book/_book/\347\254\254\345\205\253\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\344\270\200.html"
+++ "b/Book/_book/\347\254\254\345\205\253\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\344\270\200.html"
@@ -136,6 +136,29 @@
div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
+
@@ -162,8 +185,8 @@
- 关于作者
- 1 第一讲:为什么要学习R
@@ -215,7 +238,7 @@
- 3 第三讲:Git & RStudio工作流
-- 3.1 Files and Folders System
+- 3.1 文件与文件夹结构系统
- 3.2 Git and Git Hub
- 3.2.1 Git
@@ -577,8 +600,8 @@ 8.2.1 独立样本t检
ALEX14 = case_when(TRUE ~ 6 - ALEX14),
ALEX16 = case_when(TRUE ~ 6 - ALEX16),
ALEX = rowSums(select(., starts_with("ALEX")))) # 反向计分后计算总分
-
-
+
+
要在R语言中进行t检验,我们可以使用自带的stats包中的t.test()
函数。这是一个非常常用的t检验函数。在使用t.test()
函数时,我们需要输入一些参数。第一个参数是经过筛选的数据框,第二个参数是我们感兴趣的自变量(如temp),第三个参数是分组变量(如romantic)。此外,我们还可以假定两组的方差是相等的。运行t.test()
函数后,我们可以得到结果,包括t值、自由度(df)和p值。在这个例子中,t值为0.34664,自由度为1425,而p值较大,表示在恋爱组和非恋爱组之间的体温差异不显著。从结果来看,恋爱状态对体温的影响似乎并不大。
stats::t.test(data = df.penguin, # 数据框
Temperature ~ romantic, # 因变量~自变量
@@ -768,8 +791,8 @@ 8.3.3 ANOVA代码实操|数据预
[tidy data]
-
-
+
+
8.3.4 代码实操|正态性检验
@@ -951,14 +974,10 @@ 8.4 线性回归res1 <- bruceR::MANOVA(data = df,
dv = "Temperature",
between = c("climate", "romantic"))
-## [1] "Anova Table (Type III tests)"
-## [2] ""
-## [3] "Response: Temperature"
-## [4] " Effect df MSE F ges p.value"
-## [5] "1 climate 2, 1421 0.19 50.88 *** .067 <.001"
-## [6] "2 romantic 1, 1421 0.19 1.00 <.001 .317"
-## [7] "3 climate:romantic 2, 1421 0.19 5.01 ** .007 .007"
-## [8] "---"
+## [1] "Anova Table (Type III tests)" ""
+## [3] "Response: Temperature" " Effect df MSE F ges p.value"
+## [5] "1 climate 2, 1421 0.19 50.88 *** .067 <.001" "2 romantic 1, 1421 0.19 1.00 <.001 .317"
+## [7] "3 climate:romantic 2, 1421 0.19 5.01 ** .007 .007" "---"
## [9] "Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '+' 0.1 ' ' 1"
在这个例子中,我们可以看到,我们对climate和romantic进行了方差分析,试图查看在不同的romantic条件下(即单身或非单身),climate是否有影响。这就是简单效应分析。
在R中,我们可以看到MANOVA
函数的输出结果。首先,它会输出描述性统计结果,包括climate和romantic的各个水平的均值、标准差和样本量。然后,它会给出方差分析表,包括climate的主效应、romantic的主效应和它们的交互作用的效应。我们可以看到,climate的主效应是显著的,这与我们的预期一致。即在不同气候带的人们体温存在差异,其p值非常小,F值也比较大。此外,我们也发现climate和romantic的交互作用是显著的,这也与我们的预期相符。
diff --git "a/Book/_book/\347\254\254\345\205\255\350\256\262\346\225\260\346\215\256\351\242\204\345\244\204\347\220\206.html" "b/Book/_book/\347\254\254\345\205\255\350\256\262\346\225\260\346\215\256\351\242\204\345\244\204\347\220\206.html"
index db191df0..1f9dc8b3 100644
--- "a/Book/_book/\347\254\254\345\205\255\350\256\262\346\225\260\346\215\256\351\242\204\345\244\204\347\220\206.html"
+++ "b/Book/_book/\347\254\254\345\205\255\350\256\262\346\225\260\346\215\256\351\242\204\345\244\204\347\220\206.html"
@@ -136,6 +136,29 @@
div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
+
@@ -162,8 +185,8 @@
- 关于作者
- 1 第一讲:为什么要学习R
@@ -215,7 +238,7 @@
- 3 第三讲:Git & RStudio工作流
-- 3.1 Files and Folders System
+- 3.1 文件与文件夹结构系统
- 3.2 Git and Git Hub
- 3.2.1 Git
@@ -631,64 +654,40 @@ 6.2.1 研究问题&数据情
读取数据之后,我们可以使用View
函数查看数据结构,包括站点、年龄等各种变量。我们还可以用colnames
查看变量,因为可能需要选择一些我们认为比较重要的变量,比如有一些我们可能不需要的变量,就可以直接将其筛选掉。
-
## [1] "Site" "age" "sex" "monogamous"
-## [5] "romantic" "health" "exercise" "eatdrink"
-## [9] "gluctot" "artgluctot" "smoke" "cigs"
-## [13] "avgtemp" "Temperature_t1" "Temperature_t2" "DEQ"
-## [17] "AvgHumidity" "mintemp" "language" "langfamily"
-## [21] "SNI1" "SNI2" "SNI3" "SNI4"
-## [25] "SNI5" "SNI6" "SNI7" "SNI8"
-## [29] "SNI9" "SNI10" "SNI11" "SNI12"
-## [33] "SNI13" "SNI14" "SNI15" "SNI16"
-## [37] "SNI17" "SNI18" "SNI19" "SNI20"
-## [41] "SNI21" "SNI22" "SNI23" "SNI24"
-## [45] "SNI25" "SNI26" "SNI27" "SNI28"
-## [49] "SNI29" "SNI30" "SNI31" "SNI32"
-## [53] "scontrol1" "scontrol2" "scontrol3" "scontrol4"
-## [57] "scontrol5" "scontrol6" "scontrol7" "scontrol8"
-## [61] "scontrol9" "scontrol10" "scontrol11" "scontrol12"
-## [65] "scontrol13" "stress1" "stress2" "stress3"
-## [69] "stress4" "stress5" "stress6" "stress7"
-## [73] "stress8" "stress9" "stress10" "stress11"
-## [77] "stress12" "stress13" "stress14" "phone1"
-## [81] "phone2" "phone3" "phone4" "phone5"
-## [85] "phone6" "phone7" "phone8" "phone9"
-## [89] "onlineid1" "onlineid2" "onlineid3" "onlineid4"
-## [93] "onlineid5" "onlineid6" "onlineid7" "onlineid8"
-## [97] "onlineid9" "onlineid10" "onlineide11" "ECR1"
-## [101] "ECR2" "ECR3" "ECR4" "ECR5"
-## [105] "ECR6" "ECR7" "ECR8" "ECR9"
-## [109] "ECR10" "ECR11" "ECR12" "ECR13"
-## [113] "ECR14" "ECR15" "ECR16" "ECR17"
-## [117] "ECR18" "ECR19" "ECR20" "ECR21"
-## [121] "ECR22" "ECR23" "ECR24" "ECR25"
-## [125] "ECR26" "ECR27" "ECR28" "ECR29"
-## [129] "ECR30" "ECR31" "ECR32" "ECR33"
-## [133] "ECR34" "ECR35" "ECR36" "HOME1"
-## [137] "HOME2" "HOME3" "HOME4" "HOME5"
-## [141] "HOME6" "HOME7" "HOME8" "HOME9"
-## [145] "SNS1" "SNS2" "SNS3" "SNS4"
-## [149] "SNS5" "SNS6" "SNS7" "ALEX1"
-## [153] "ALEX2" "ALEX3" "ALEX4" "ALEX5"
-## [157] "ALEX6" "ALEX7" "ALEX8" "ALEX9"
-## [161] "ALEX10" "ALEX11" "ALEX12" "ALEX13"
-## [165] "ALEX14" "ALEX15" "ALEX16" "KAMF1"
-## [169] "KAMF2" "KAMF3" "KAMF4" "KAMF5"
-## [173] "KAMF6" "KAMF7" "STRAQ_1" "STRAQ_2"
-## [177] "STRAQ_3" "STRAQ_4" "STRAQ_5" "STRAQ_6"
-## [181] "STRAQ_7" "STRAQ_8" "STRAQ_9" "STRAQ_10"
-## [185] "STRAQ_11" "STRAQ_12" "STRAQ_13" "STRAQ_14"
-## [189] "STRAQ_15" "STRAQ_16" "STRAQ_17" "STRAQ_18"
-## [193] "STRAQ_19" "STRAQ_20" "STRAQ_21" "STRAQ_22"
-## [197] "STRAQ_23" "STRAQ_24" "STRAQ_25" "STRAQ_26"
-## [201] "STRAQ_27" "STRAQ_28" "STRAQ_29" "STRAQ_30"
-## [205] "STRAQ_31" "STRAQ_32" "STRAQ_33" "STRAQ_34"
-## [209] "STRAQ_35" "STRAQ_36" "STRAQ_37" "STRAQ_38"
-## [213] "STRAQ_39" "STRAQ_40" "STRAQ_41" "STRAQ_42"
-## [217] "STRAQ_43" "STRAQ_44" "STRAQ_45" "STRAQ_46"
-## [221] "STRAQ_47" "STRAQ_48" "STRAQ_49" "STRAQ_50"
-## [225] "STRAQ_51" "STRAQ_52" "STRAQ_53" "STRAQ_54"
-## [229] "STRAQ_55" "STRAQ_56" "STRAQ_57" "socialdiversity"
+## [1] "Site" "age" "sex" "monogamous" "romantic" "health" "exercise"
+## [8] "eatdrink" "gluctot" "artgluctot" "smoke" "cigs" "avgtemp" "Temperature_t1"
+## [15] "Temperature_t2" "DEQ" "AvgHumidity" "mintemp" "language" "langfamily" "SNI1"
+## [22] "SNI2" "SNI3" "SNI4" "SNI5" "SNI6" "SNI7" "SNI8"
+## [29] "SNI9" "SNI10" "SNI11" "SNI12" "SNI13" "SNI14" "SNI15"
+## [36] "SNI16" "SNI17" "SNI18" "SNI19" "SNI20" "SNI21" "SNI22"
+## [43] "SNI23" "SNI24" "SNI25" "SNI26" "SNI27" "SNI28" "SNI29"
+## [50] "SNI30" "SNI31" "SNI32" "scontrol1" "scontrol2" "scontrol3" "scontrol4"
+## [57] "scontrol5" "scontrol6" "scontrol7" "scontrol8" "scontrol9" "scontrol10" "scontrol11"
+## [64] "scontrol12" "scontrol13" "stress1" "stress2" "stress3" "stress4" "stress5"
+## [71] "stress6" "stress7" "stress8" "stress9" "stress10" "stress11" "stress12"
+## [78] "stress13" "stress14" "phone1" "phone2" "phone3" "phone4" "phone5"
+## [85] "phone6" "phone7" "phone8" "phone9" "onlineid1" "onlineid2" "onlineid3"
+## [92] "onlineid4" "onlineid5" "onlineid6" "onlineid7" "onlineid8" "onlineid9" "onlineid10"
+## [99] "onlineide11" "ECR1" "ECR2" "ECR3" "ECR4" "ECR5" "ECR6"
+## [106] "ECR7" "ECR8" "ECR9" "ECR10" "ECR11" "ECR12" "ECR13"
+## [113] "ECR14" "ECR15" "ECR16" "ECR17" "ECR18" "ECR19" "ECR20"
+## [120] "ECR21" "ECR22" "ECR23" "ECR24" "ECR25" "ECR26" "ECR27"
+## [127] "ECR28" "ECR29" "ECR30" "ECR31" "ECR32" "ECR33" "ECR34"
+## [134] "ECR35" "ECR36" "HOME1" "HOME2" "HOME3" "HOME4" "HOME5"
+## [141] "HOME6" "HOME7" "HOME8" "HOME9" "SNS1" "SNS2" "SNS3"
+## [148] "SNS4" "SNS5" "SNS6" "SNS7" "ALEX1" "ALEX2" "ALEX3"
+## [155] "ALEX4" "ALEX5" "ALEX6" "ALEX7" "ALEX8" "ALEX9" "ALEX10"
+## [162] "ALEX11" "ALEX12" "ALEX13" "ALEX14" "ALEX15" "ALEX16" "KAMF1"
+## [169] "KAMF2" "KAMF3" "KAMF4" "KAMF5" "KAMF6" "KAMF7" "STRAQ_1"
+## [176] "STRAQ_2" "STRAQ_3" "STRAQ_4" "STRAQ_5" "STRAQ_6" "STRAQ_7" "STRAQ_8"
+## [183] "STRAQ_9" "STRAQ_10" "STRAQ_11" "STRAQ_12" "STRAQ_13" "STRAQ_14" "STRAQ_15"
+## [190] "STRAQ_16" "STRAQ_17" "STRAQ_18" "STRAQ_19" "STRAQ_20" "STRAQ_21" "STRAQ_22"
+## [197] "STRAQ_23" "STRAQ_24" "STRAQ_25" "STRAQ_26" "STRAQ_27" "STRAQ_28" "STRAQ_29"
+## [204] "STRAQ_30" "STRAQ_31" "STRAQ_32" "STRAQ_33" "STRAQ_34" "STRAQ_35" "STRAQ_36"
+## [211] "STRAQ_37" "STRAQ_38" "STRAQ_39" "STRAQ_40" "STRAQ_41" "STRAQ_42" "STRAQ_43"
+## [218] "STRAQ_44" "STRAQ_45" "STRAQ_46" "STRAQ_47" "STRAQ_48" "STRAQ_49" "STRAQ_50"
+## [225] "STRAQ_51" "STRAQ_52" "STRAQ_53" "STRAQ_54" "STRAQ_55" "STRAQ_56" "STRAQ_57"
+## [232] "socialdiversity"
另外,当我们需要进行问卷得分的计算时,我们可能需要选择与特定问卷相关的条目,然后对其进行预处理,以得到问卷得分。因此,我们必须要了解数据的当前结构,例如每个列的名称。我们得到的数据是比较干净的,每个名字基本上都可以反映出其大致信息,例如scontrol1
、scontrol2
等代表自我控制问卷的各个条目。这里的命名本身就是比较规范的。如果数据的列名不是这么清晰,可能还需要进行额外的处理步骤,比如重新命名问卷等。我们就不再一一赘述了。
那么,我们在这里的主要目标是什么,假设我们选择这几个变量,作为数据预处理的练习,我们把身体温度选出来,然后求一个平均的温度。
●研究核心变量:
@@ -786,12 +785,10 @@ 6.2.1 研究问题&数据情
我们前面说在studio里面有可能显示不全,在课件里面也可能显示不全。但在这里,它会显示得非常全,我们可以拖着看,它有两页。我们也可以去查看它的列名,在这里我们可以非常全地看到它的列名。(前面的…需要再确认,每个代码都需要加上)
-## [1] "Temperature_t1" "Temperature_t2" "socialdiversity" "Site"
-## [5] "DEQ" "romantic" "ALEX1" "ALEX2"
-## [9] "ALEX3" "ALEX4" "ALEX5" "ALEX6"
-## [13] "ALEX7" "ALEX8" "ALEX9" "ALEX10"
-## [17] "ALEX11" "ALEX12" "ALEX13" "ALEX14"
-## [21] "ALEX15" "ALEX16" "Temperature_t1_new"
+## [1] "Temperature_t1" "Temperature_t2" "socialdiversity" "Site" "DEQ" "romantic"
+## [7] "ALEX1" "ALEX2" "ALEX3" "ALEX4" "ALEX5" "ALEX6"
+## [13] "ALEX7" "ALEX8" "ALEX9" "ALEX10" "ALEX11" "ALEX12"
+## [19] "ALEX13" "ALEX14" "ALEX15" "ALEX16" "Temperature_t1_new"
然后加入Tidyverse
的包。
@@ -1013,8 +1010,8 @@ 6.3.1 研究问题 & 数据
ACC: 被试的判断是否正确,1 = “正确”, 0 = “错误”, -1,
2表示未按键或按了两个键的情况,属于无效作答
RT: 被试做出判断的反应时,[200,1500]的反应时纳入分析
-
-
+
+
这些数据看起来是这样的,对吧?如果我们打开这些文件,我们现在需要计算自我优势效应,即在匹配效应下,好的我和好的他人之间的差异。在这种情况下,好的我和好的他人都是好的,但是我的反应速度比他人快,这表明我们对自身的信息有一个优势效应。因此,我们需要批量读取数据,但由于有44个文件,我们不可能手动输入每个文件。我们必须使用代码来提高效率。
接下来,我们需要拆分变量,了解如何对字符进行拆分,以及如何将长数据转换为宽数据。我们将稍后讨论这些操作。
@@ -1260,8 +1257,7 @@ 6.3.2 操作步骤# 分实验条件计算
df4 <- dplyr::group_by(df4, Sub, Shape, Label, Match)
df4 <- dplyr::summarise(df4, mean_ACC = mean(ACC), mean_RT = mean(RT))
-## `summarise()` has grouped output by 'Sub', 'Shape', 'Label'. You can override using the
-## `.groups` argument.
+## `summarise()` has grouped output by 'Sub', 'Shape', 'Label'. You can override using the `.groups` argument.
- 拆分变量[
extract
, filter
]
diff --git "a/Book/_book/\347\254\254\345\215\201\344\270\200\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\345\233\233\344\270\255\344\273\213\345\210\206\346\236\220.html" "b/Book/_book/\347\254\254\345\215\201\344\270\200\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\345\233\233\344\270\255\344\273\213\345\210\206\346\236\220.html"
index 6cb55162..03055236 100644
--- "a/Book/_book/\347\254\254\345\215\201\344\270\200\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\345\233\233\344\270\255\344\273\213\345\210\206\346\236\220.html"
+++ "b/Book/_book/\347\254\254\345\215\201\344\270\200\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\345\233\233\344\270\255\344\273\213\345\210\206\346\236\220.html"
@@ -136,6 +136,29 @@
div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
+
@@ -162,8 +185,8 @@
- 关于作者
- 1 第一讲:为什么要学习R
@@ -215,7 +238,7 @@
- 3 第三讲:Git & RStudio工作流
-- 3.1 Files and Folders System
+- 3.1 文件与文件夹结构系统
- 3.2 Git and Git Hub
- 3.2.1 Git
@@ -651,8 +674,8 @@ 11.3 中介分析
11.3.1 2.1 对于“机制”的表示——“图”
我们在做路径分析的时候,我们一般使用“图”来表征理论因果关系
-
-
+
+
图包括两部分:节点和边。节点表示具体变量,而箭头表示变量之间的关系;
对节点来说,在SEM中,观测变量用矩形表示,潜变量用椭圆表示。
@@ -660,12 +683,12 @@ 11.3.1 2.1 对于“机制”的
我们在做假设模型时,一般都不会假设出现下图这种情况,可以看到,
从x1到x2,x2到x3,x3又回来到x1了,形成了一个闭环。
-
-
+
+
我们一般来说更希望的是不存在循环的关系,如下图,在统计和机械学习领域,称其为有向无环图(Directed
Acyclic Graph, DAGs)
-
-
+
+
11.3.2 中介分析
@@ -675,8 +698,8 @@ 11.3.2 中介分析\(X \rightarrow M \rightarrow Y\) 。
-
-
+
+
这里会有一个问题,我们为什么会关心中介机制呢,或者是中介变量呢?
中介分析,它对我们现实生活到底有什么真实意义?
有一种可能的例子,比方说原生家庭是自变量X,心理健康因变量是Y,
@@ -689,20 +712,20 @@
11.3.2 中介分析\(e_1\)。这里就是一个简单的线性回归模型。
总方程:
\[Y = i_1 + cX + e_1\]
-
-
+
+
分解:
\[M = i_2 + aX + e_2\]
\[Y = i_3 + c'X + bM + e_3\]
当我们加入中介时,我们就假定有这么一个中介变量m,m是可以通过x来预测的:\(M = i_2 + aX + e_2\),且又能够再去预测y:\(Y = i_3 + c'X + bM + e_3\),
所以它最后形成这样的形式:
-
-
+
+
11.3.3 中介效应
-
-
+
+
\[ Y = i_1 + cX + e_1\]
\[ M = i_2 + aX + e_2\]
\[Y = i_3 + c'X + bM + e_3\]
@@ -775,8 +798,8 @@ 11.3.5 问题提出
-
+
+
这里我们不再去做探索分析,进行数据的预处理之后,直接计算各个变量的值,再去做中介分析。如果要从原始数据里面去计算社交网络的复杂度CSI,可以去仔细的去看一看这里面的代码,会涉及到比较多的反复的重新编码。
@@ -1007,9 +1030,9 @@ 11.3.7 PROCESS in bruceR()11.3.7 PROCESS in bruceR() capture.output() %>% .[27:43]
##
## Warning: nsim=1000 (or larger) is suggested!
-## [1] "Model Summary"
-## [2] ""
-## [3] "──────────────────────────────────────────────────"
-## [4] " (1) CBT (2) CSI (3) CBT "
-## [5] "──────────────────────────────────────────────────"
-## [6] "(Intercept) 36.426 *** 7.159 *** 36.426 ***"
-## [7] " (0.015) (0.050) (0.015) "
-## [8] "DEQ -0.001 0.029 *** -0.002 "
-## [9] " (0.001) (0.004) (0.001) "
-## [10] "CSI 0.046 ***"
-## [11] " (0.011) "
-## [12] "──────────────────────────────────────────────────"
-## [13] "R^2 0.001 0.082 0.024 "
-## [14] "Adj. R^2 -0.001 0.081 0.021 "
-## [15] "Num. obs. 763 763 763 "
-## [16] "──────────────────────────────────────────────────"
+## [1] "Model Summary" ""
+## [3] "──────────────────────────────────────────────────" " (1) CBT (2) CSI (3) CBT "
+## [5] "──────────────────────────────────────────────────" "(Intercept) 36.426 *** 7.159 *** 36.426 ***"
+## [7] " (0.015) (0.050) (0.015) " "DEQ -0.001 0.029 *** -0.002 "
+## [9] " (0.001) (0.004) (0.001) " "CSI 0.046 ***"
+## [11] " (0.011) " "──────────────────────────────────────────────────"
+## [13] "R^2 0.001 0.082 0.024 " "Adj. R^2 -0.001 0.081 0.021 "
+## [15] "Num. obs. 763 763 763 " "──────────────────────────────────────────────────"
## [17] "Note. * p < .05, ** p < .01, *** p < .001."
##
## Warning: nsim=1000 (or larger) is suggested!
@@ -1058,9 +1073,9 @@ 11.3.7 PROCESS in bruceR()
可以看到,bruceR::PROCESS
的结果输出很符合我们在论文中看到的结果呈现方式,
对我们研究生或者研究者来说,非常友好,而且跑出来的结果还可以输出到Word文档中去。
中介模型的数值本身并不能告诉我们因果关系,而在心理学中,利用中介变量去发现因果关系又是很重要的,尤其是当自变量是一些不可干预的变量时,通过干预中介来达到科学研究的目的是很好的方式。这就会涉及到我们如何去验证带有真正的因果的中介,当然,就如同RCT研究中想要得出因果都需要很复杂的假设,要去验证真正有因果的中介,是一件难上加难的事情。
这里,我们用一个国内的学者葛枭语老师,发表在JESP上面一篇文章为例子,来探讨如何验证中介中的因果。
- - + +假设:教材难度(X)通过焦虑(M)来影响努力程度(Y),可以穷举出在哪些情况下我们不能验证中介中的因果:
教材难度(X)不能影响焦虑(M)
• 测量 Y
对X进行操纵(如使用不同难度的教材),可以验证X对M的因果关系,但M与Y之间的因果关系并没有得到验证
- - + +但如果我们理论假设错误,测量的是焦虑(A),但实际上实验操纵引发的中介应为恐惧(M,即实际路径应为X - M - Y,而我们测量路径为X - A - Y),那么刚才的实验设计可能无法证伪,因此需要对A进行操纵:
@@ -1146,8 +1161,8 @@对X(如使用不同难度的教材)和A(控制组 vs 提供相关辅导以减轻焦虑)进行操纵,如果对A的操纵不能影响Y,则可以证明中介路径不合理
- - + +所以这里意味着我们要去进行带有因果性质的中介分析时,我们要做很多科学的思考, 在实验设计上,尽量去控制混淆变量,这样的话才能够去真正的把因果关系搞清楚。
如果去看其他的一些领域的文章, diff --git "a/Book/_book/\347\254\254\345\215\201\344\270\211\350\256\262\345\237\272\344\272\216\347\275\221\347\273\234\346\250\241\345\236\213\347\232\204\345\277\203\347\220\206\345\255\246\347\240\224\347\251\266.html" "b/Book/_book/\347\254\254\345\215\201\344\270\211\350\256\262\345\237\272\344\272\216\347\275\221\347\273\234\346\250\241\345\236\213\347\232\204\345\277\203\347\220\206\345\255\246\347\240\224\347\251\266.html" index 248b9b30..8a1b3abb 100644 --- "a/Book/_book/\347\254\254\345\215\201\344\270\211\350\256\262\345\237\272\344\272\216\347\275\221\347\273\234\346\250\241\345\236\213\347\232\204\345\277\203\347\220\206\345\255\246\347\240\224\347\251\266.html" +++ "b/Book/_book/\347\254\254\345\215\201\344\270\211\350\256\262\345\237\272\344\272\216\347\275\221\347\273\234\346\250\241\345\236\213\347\232\204\345\277\203\347\220\206\345\255\246\347\240\224\347\251\266.html" @@ -136,6 +136,29 @@ div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;} + @@ -162,8 +185,8 @@
接下来我们整合试次,得到每个被试在4中条件下的均值。
-## `summarise()` has grouped output by 'Sub', 'Identity'. You can override using the
-## `.groups` argument.
+## `summarise()` has grouped output by 'Sub', 'Identity'. You can override using the `.groups` argument.
-
-
+
+
因为我们在作图时想看到的是总体的交互作用,因此我们将所有被试的数据平均得到一个总体在4种条件下的RT和ACC均值。在下面的代码中,大家可以发现SD和SE,这是因为我们在作图时还需要在均值的基础上加上误差棒(error bar)。这里我们需要将误差棒和字母”T”进行区分。因为之前有个研究者拿字母”T”冒充误差棒,当然后来被人发现了,还在社交媒体上火了一把。这也提示我们作图时,拿最基础的柱图举例,不仅要画出集中趋势(平均值),还要画出离散趋势(误差项)。
df.match.sum <- df.match.subj %>%
@@ -605,11 +627,10 @@ 12.1.2 作图数据准备 SE_ACC = SD_ACC/sqrt(n()-1),
n = n()) %>%
dplyr::ungroup()
## `summarise()` has grouped output by 'Identity'. You can override using the `.groups`
-## argument.
+## `summarise()` has grouped output by 'Identity'. You can override using the `.groups` argument.
-
-
+
+
整体的均值加上误差,这样的图在传统意义上已经足够了。但是最近几年,大家发现如果只呈现总体趋势,对读者是有误导性的。如果我们呈现每个被试数据的话,读者会发现,整体的效应量实际上没有那么大。因此最近整体的趋势是鼓励大家在作图时不仅要可视化整体数据,还要将每个个体的数据也呈现在图上。这样能够保证读者在读图时不会过高地估计实验的效应量。
@@ -680,8 +701,8 @@此时,我们的数据框就变成了一个长形的数据,并使用DV这一列对数据因变量进行分类。合并与分类结束之后,我们就可以使用facet
来进行呈现了。在下面的代码中,papaja::theme_apa()
这一行之前其实和前面提到的代码是一致的,只是将y定义为合并之后的因变量grand_mean。这时我们需要新增的就是放入DV这一分类变量,scales = "free_y"
表示在不同面板上y轴的单位不固定,根据数据情况来自动调整。最后根据实际情况增加一些title和轴上的label。
p4 <- df.match.sum.long %>%
ggplot2::ggplot(.,
diff --git "a/Book/_book/\347\254\254\345\215\201\345\233\233\350\256\262\345\277\203\347\220\206\345\255\246\345\205\203\345\210\206\346\236\220\345\205\245\351\227\250.html" "b/Book/_book/\347\254\254\345\215\201\345\233\233\350\256\262\345\277\203\347\220\206\345\255\246\345\205\203\345\210\206\346\236\220\345\205\245\351\227\250.html"
index dc8287d3..68cd4686 100644
--- "a/Book/_book/\347\254\254\345\215\201\345\233\233\350\256\262\345\277\203\347\220\206\345\255\246\345\205\203\345\210\206\346\236\220\345\205\245\351\227\250.html"
+++ "b/Book/_book/\347\254\254\345\215\201\345\233\233\350\256\262\345\277\203\347\220\206\345\255\246\345\205\203\345\210\206\346\236\220\345\205\245\351\227\250.html"
@@ -136,6 +136,29 @@
div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
+
@@ -162,8 +185,8 @@
- 关于作者
- 1 第一讲:为什么要学习R
@@ -215,7 +238,7 @@
- 3 第三讲:Git & RStudio工作流
-- 3.1 Files and Folders System
+- 3.1 文件与文件夹结构系统
- 3.2 Git and Git Hub
- 3.2.1 Git
@@ -828,20 +851,13 @@ 14.2 元分析的实施df <- df[order(df$yi), ]
head(df)
##
-## id Information N ES VI Cohort Context Materials Valence
-## 12 12 Kühl and Münzer (2021), G2 60 -0.5278 0.0691 3 1 1 2
-## 10 10 Kühl and Zander (2017), S2 71 -0.4687 0.0583 3 1 1 2
-## 9 9 Kühl and Zander (2017), S1 77 -0.0640 0.0530 3 1 1 2
-## 18 18 Sinatra et al. (2016) 111 -0.0560 0.0541 3 1 1 1
-## 14 14 Li et al. (2021), S2 29 0.0477 0.0276 3 1 1 1
-## 4 4 d'Ailly et al. (1997) 100 0.2108 0.0096 1 1 1 1
-## yi vi
-## 12 -0.5278 0.0691
-## 10 -0.4687 0.0583
-## 9 -0.0640 0.0530
-## 18 -0.0560 0.0541
-## 14 0.0477 0.0276
-## 4 0.2108 0.0096
+## id Information N ES VI Cohort Context Materials Valence yi vi
+## 12 12 Kühl and Münzer (2021), G2 60 -0.5278 0.0691 3 1 1 2 -0.5278 0.0691
+## 10 10 Kühl and Zander (2017), S2 71 -0.4687 0.0583 3 1 1 2 -0.4687 0.0583
+## 9 9 Kühl and Zander (2017), S1 77 -0.0640 0.0530 3 1 1 2 -0.0640 0.0530
+## 18 18 Sinatra et al. (2016) 111 -0.0560 0.0541 3 1 1 1 -0.0560 0.0541
+## 14 14 Li et al. (2021), S2 29 0.0477 0.0276 3 1 1 1 0.0477 0.0276
+## 4 4 d'Ailly et al. (1997) 100 0.2108 0.0096 1 1 1 1 0.2108 0.0096
在元分析中,计算综合效应量是一个关键步骤,这一步骤实际上可以通过一行代码来完成。我们可以使用metafor软件包中的rma函数来计算。在这个函数中,我们需要输入的是每个研究的效应量及其变异性(variance)。此外,我们还需要指定分析使用的方法(method)。通常情况下,我们可以使用函数的默认方法,除非有特定的理由需要更改。
res <- metafor::rma(
yi,#每个研究的效应量
diff --git "a/Book/_book/\347\254\254\345\215\201\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\344\270\211\345\271\277\344\271\211\347\272\277\346\200\247\346\250\241\345\236\213.html" "b/Book/_book/\347\254\254\345\215\201\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\344\270\211\345\271\277\344\271\211\347\272\277\346\200\247\346\250\241\345\236\213.html"
index 5013174a..f63d875a 100644
--- "a/Book/_book/\347\254\254\345\215\201\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\344\270\211\345\271\277\344\271\211\347\272\277\346\200\247\346\250\241\345\236\213.html"
+++ "b/Book/_book/\347\254\254\345\215\201\350\256\262\345\233\236\345\275\222\346\250\241\345\236\213\344\270\211\345\271\277\344\271\211\347\272\277\346\200\247\346\250\241\345\236\213.html"
@@ -136,6 +136,29 @@
div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
+
@@ -162,8 +185,8 @@
- 关于作者
- 1 第一讲:为什么要学习R
@@ -215,7 +238,7 @@
- 3 第三讲:Git & RStudio工作流
-- 3.1 Files and Folders System
+- 3.1 文件与文件夹结构系统
- 3.2 Git and Git Hub
- 3.2.1 Git
@@ -548,20 +571,19 @@ 10.1 前章回顾和本章数据
level),层级模型可以提高统计检验力并且为研究者提供十分丰富的信息,因此也是学界越来越推荐的一种统计方法。
让我们先对本章数据进行预处理,这里用到的还是之前的认知实验的数据,我们先对正确率进行预处理。
这里需要注意一下我们对于正确率ACC的处理,实际上我们的认知实验数据当中包含了其他两种反应,这里我们直接删除了另外两种情况,只保留了正确和错误的反应,即0和1,当然也有一些实验会将1之外的所有反应归到0中做处理。另外我们筛选去除了反应在1500ms以上和200ms以下的反应时,因为这在经验上是不符合人类的反应速度的。
-
## `summarise()` has grouped output by 'Sub', 'Valence'. You can override using the
-## `.groups` argument.
+## `summarise()` has grouped output by 'Sub', 'Valence'. You can override using the `.groups` argument.
-
-
+
+
上面是我们通常所做的操作:将不同条件下的反应正确率做一个平均,然后进行方差分析。
(知识补充:easystats系统包是过去五六年快速发展起来的一个包系列,适用于统计分析,特别是心理学相关背景的统计分析,具体使用可以参考我们在B站上传的视频(链接如下)。https://www.bilibili.com/video/BV1rz421D7iJ/?spm_id_from=333.337.search-card.all.click)
让我们再来看一下正确率的原始数据。
-
-
+
+
可以发现其只存在0和1两种取值,这种分布显然不服从正态分布,因此我们不能简单地用前一章提到的一般线性模型进行处理。在传统的方差分析中,我们对正确率数据的处理是求出每个条件下的平均正确率再进行统计分析。这个平均正确率的取值作为一个连续数据,可以被放在坐标轴上形成一个分布,其与以0为原点,向两端无限延伸的标准正态分布也存在差异。因此我们需要对一般线性模型进行拓展,这也就是我们这一章所要讲的广义线性模型(Generalized
Linear Model, GLM)。
## [1] " AIC BIC logLik deviance df.resid "
-## [2] " 9378.8 9460.2 -4678.4 9356.8 11999 "
-## [3] "Random effects:"
-## [4] " Groups Name Variance Std.Dev. Corr "
-## [5] " Sub (Intercept) 1.53 1.24 "
-## [6] " IdentityOther 2.52 1.59 -0.86 "
-## [7] " Valenceimmoral 2.40 1.55 -0.85 0.83 "
-## [8] " IdentityOther:Valenceimmoral 3.33 1.83 0.69 -0.87 -0.82"
-## [9] "Number of obs: 12010, groups: Sub, 41"
-## [10] ""
-## [11] "Fixed effects:"
-## [12] " Estimate Std. Error z value Pr(>|z|) "
-## [13] "(Intercept) 2.014 0.114 17.7 <0.0000000000000002 ***"
+## [1] " 9378.8 9460.2 -4678.4 9356.8 11999 "
+## [2] ""
+## [3] " Groups Name Variance Std.Dev. Corr "
+## [4] " Sub (Intercept) 1.53 1.24 "
+## [5] " IdentityOther 2.52 1.59 -0.86 "
+## [6] " Valenceimmoral 2.40 1.55 -0.85 0.83 "
+## [7] " IdentityOther:Valenceimmoral 3.33 1.83 0.69 -0.87 -0.82"
+## [8] "Number of obs: 12010, groups: Sub, 41"
+## [9] ""
+## [10] "Fixed effects:"
+## [11] " Estimate Std. Error z value Pr(>|z|) "
+## [12] "(Intercept) 2.014 0.114 17.7 <0.0000000000000002 ***"
+## [13] "---"
接着我们建立一个随机效应只包含截距的模型。
#随机截距,固定斜率
mod <- lme4::glmer(data = df.match, #数据
@@ -698,24 +720,24 @@ 10.3.2 GLM代码实操#performance::model_performance(mod)
summary(mod) %>%
capture.output() %>% .[c(7:8,14:24,28:32)]
## [1] " AIC BIC logLik deviance df.resid "
-## [2] " 9639.0 9675.9 -4814.5 9629.0 12005 "
-## [3] "Random effects:"
-## [4] " Groups Name Variance Std.Dev."
-## [5] " Sub (Intercept) 0.237 0.487 "
-## [6] "Number of obs: 12010, groups: Sub, 41"
-## [7] ""
-## [8] "Fixed effects:"
-## [9] " Estimate Std. Error z value Pr(>|z|) "
-## [10] "(Intercept) 2.4964 0.1015 24.60 < 0.0000000000000002 ***"
-## [11] "IdentityOther -0.7160 0.0839 -8.53 < 0.0000000000000002 ***"
-## [12] "Valenceimmoral -0.9474 0.0818 -11.58 < 0.0000000000000002 ***"
-## [13] "IdentityOther:Valenceimmoral 0.8230 0.1086 7.58 0.000000000000034 ***"
-## [14] "Correlation of Fixed Effects:"
-## [15] " (Intr) IdnttO Vlncmm"
-## [16] "IdenttyOthr -0.519 "
-## [17] "Valencemmrl -0.533 0.641 "
-## [18] "IdnttyOth:V 0.401 -0.773 -0.754"
+## [1] " 9639.0 9675.9 -4814.5 9629.0 12005 "
+## [2] ""
+## [3] " Groups Name Variance Std.Dev."
+## [4] " Sub (Intercept) 0.237 0.487 "
+## [5] "Number of obs: 12010, groups: Sub, 41"
+## [6] ""
+## [7] "Fixed effects:"
+## [8] " Estimate Std. Error z value Pr(>|z|) "
+## [9] "(Intercept) 2.4964 0.1015 24.60 < 0.0000000000000002 ***"
+## [10] "IdentityOther -0.7160 0.0839 -8.53 < 0.0000000000000002 ***"
+## [11] "Valenceimmoral -0.9474 0.0818 -11.58 < 0.0000000000000002 ***"
+## [12] "IdentityOther:Valenceimmoral 0.8230 0.1086 7.58 0.000000000000034 ***"
+## [13] "---"
+## [14] " (Intr) IdnttO Vlncmm"
+## [15] "IdenttyOthr -0.519 "
+## [16] "Valencemmrl -0.533 0.641 "
+## [17] "IdnttyOth:V 0.401 -0.773 -0.754"
+## [18] NA
最后我们建立一个包含了所有固定效应和随机效应的全模型。
#随机截距,随机斜率
mod_full <- lme4::glmer(data = df.match, #数据
@@ -724,26 +746,26 @@ 10.3.2 GLM代码实操##performance::model_performance(mod_full)
summary(mod_full) %>%
capture.output() %>% .[c(6:8,13:18,21:26,30:34)]
## [1] ""
-## [2] " AIC BIC logLik deviance df.resid "
-## [3] " 9355.9 9459.4 -4664.0 9327.9 11996 "
-## [4] ""
-## [5] "Random effects:"
-## [6] " Groups Name Variance Std.Dev. Corr "
-## [7] " Sub (Intercept) 0.972 0.986 "
-## [8] " IdentityOther 1.771 1.331 -0.79 "
-## [9] " Valenceimmoral 1.028 1.014 -0.75 0.75 "
-## [10] ""
-## [11] "Fixed effects:"
-## [12] " Estimate Std. Error z value Pr(>|z|) "
-## [13] "(Intercept) 2.772 0.178 15.53 < 0.0000000000000002 ***"
-## [14] "IdentityOther -0.870 0.235 -3.71 0.00021 ***"
-## [15] "Valenceimmoral -1.150 0.190 -6.06 0.0000000013 ***"
-## [16] ""
-## [17] "Correlation of Fixed Effects:"
-## [18] " (Intr) IdnttO Vlncmm"
-## [19] "IdenttyOthr -0.801 "
-## [20] "Valencemmrl -0.783 0.741 "
+## [1] " AIC BIC logLik deviance df.resid "
+## [2] " 9355.9 9459.4 -4664.0 9327.9 11996 "
+## [3] ""
+## [4] "Random effects:"
+## [5] " Groups Name Variance Std.Dev. Corr "
+## [6] " Sub (Intercept) 0.972 0.986 "
+## [7] " IdentityOther 1.771 1.331 -0.79 "
+## [8] " Valenceimmoral 1.028 1.014 -0.75 0.75 "
+## [9] " IdentityOther:Valenceimmoral 2.306 1.518 0.54 -0.82 -0.74"
+## [10] "Fixed effects:"
+## [11] " Estimate Std. Error z value Pr(>|z|) "
+## [12] "(Intercept) 2.772 0.178 15.53 < 0.0000000000000002 ***"
+## [13] "IdentityOther -0.870 0.235 -3.71 0.00021 ***"
+## [14] "Valenceimmoral -1.150 0.190 -6.06 0.0000000013 ***"
+## [15] "IdentityOther:Valenceimmoral 0.988 0.272 3.63 0.00028 ***"
+## [16] "Correlation of Fixed Effects:"
+## [17] " (Intr) IdnttO Vlncmm"
+## [18] "IdenttyOthr -0.801 "
+## [19] "Valencemmrl -0.783 0.741 "
+## [20] "IdnttyOth:V 0.573 -0.821 -0.747"
我们在运行全模型的时候可以明显感受到运行时间相较于之前变长了。当我们的模型内参数越多,模型越复杂的时候,计算机就需要花更多时间去拟合模型,也会有些时候因为找不到合适的参数而导致模型无法拟合。习惯了SPSS的读者可能会难以忍受,但实际上我们在后面处理一些大数据或者跑机器学习的时候,等待会是一件很常见的事情。这就提示我们合理分配时间,把要运行的代码提前运行起来,然后去做别的工作。
在这里,我们也可以根据结果来判断R语言对我们自变量的编码方式,可以看到结果中除了截距外的第一项为”Identityother”,由此我们可以判断R将”Identityself”编码为了基线,并据此来计算相应的回归系数和估计值,下面的其他结果也类似。
接下来我们对上述三个模型进行比较。
@@ -761,18 +783,17 @@performance::compare_performance(mod_null, mod, mod_full, rank = TRUE, verbose = FALSE)
-## Some of the nested models seem to be identical and probably only vary in their
-## random effects.
+## Some of the nested models seem to be identical and probably only vary in their random effects.
让我们输出全模型的结果,尝试进行解读。
-## [1] ""
-## [2] "Fixed effects:"
-## [3] " Estimate Std. Error z value Pr(>|z|) "
-## [4] "(Intercept) 2.772 0.178 15.53 < 0.0000000000000002 ***"
-## [5] "IdentityOther -0.870 0.235 -3.71 0.00021 ***"
-## [6] "Valenceimmoral -1.150 0.190 -6.06 0.0000000013 ***"
-## [7] "IdentityOther:Valenceimmoral 0.988 0.272 3.63 0.00028 ***"
+## [1] "Fixed effects:"
+## [2] " Estimate Std. Error z value Pr(>|z|) "
+## [3] "(Intercept) 2.772 0.178 15.53 < 0.0000000000000002 ***"
+## [4] "IdentityOther -0.870 0.235 -3.71 0.00021 ***"
+## [5] "Valenceimmoral -1.150 0.190 -6.06 0.0000000013 ***"
+## [6] "IdentityOther:Valenceimmoral 0.988 0.272 3.63 0.00028 ***"
+## [7] "---"
在结果里面显示的估计值,并不直接等于p值,我们要根据前面转换函数的逆运算来讲其转换为p值,转换公式如下。
代入之后我们便可以求出不同实验条件下各自的p值。
@@ -805,12 +826,9 @@## [1] "Response: mean_ACC"
-## [2] " Effect df MSE F ges p.value"
-## [3] "1 Identity 1, 40 0.01 3.08 + .017 .087"
-## [4] "2 Valence 1, 40 0.01 16.26 *** .068 <.001"
-## [5] "3 Identity:Valence 1, 40 0.01 8.52 ** .038 .006"
-## [6] "---"
+## [1] "Response: mean_ACC" " Effect df MSE F ges p.value"
+## [3] "1 Identity 1, 40 0.01 3.08 + .017 .087" "2 Valence 1, 40 0.01 16.26 *** .068 <.001"
+## [5] "3 Identity:Valence 1, 40 0.01 8.52 ** .038 .006" "---"
我们可以得到一个f值,并且因为研究的被试量比较大,我们可以发现这里呈现的主效应以及交互作用,和后面用glm或者层级模型做出来的结果有一样的趋势。我们可以用EMMAMNS()函数来查看模型的一些具体值。
res %>%
bruceR::EMMEANS(effect = 'Valence', by = 'Identity') %>%
diff --git "a/Book/_book/\347\254\254\345\233\233\350\256\262\345\246\202\344\275\225\345\257\274\345\205\245\346\225\260\346\215\256.html" "b/Book/_book/\347\254\254\345\233\233\350\256\262\345\246\202\344\275\225\345\257\274\345\205\245\346\225\260\346\215\256.html"
index a8711447..e2c5c70b 100644
--- "a/Book/_book/\347\254\254\345\233\233\350\256\262\345\246\202\344\275\225\345\257\274\345\205\245\346\225\260\346\215\256.html"
+++ "b/Book/_book/\347\254\254\345\233\233\350\256\262\345\246\202\344\275\225\345\257\274\345\205\245\346\225\260\346\215\256.html"
@@ -136,6 +136,29 @@
div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
+
@@ -162,8 +185,8 @@
- 关于作者
- 1 第一讲:为什么要学习R
@@ -215,7 +238,7 @@
- 3 第三讲:Git & RStudio工作流
-- 3.1 Files and Folders System
+- 3.1 文件与文件夹结构系统
- 3.2 Git and Git Hub
- 3.2.1 Git
diff --git "a/Book/_book/\350\241\245\345\205\2051\345\246\202\344\275\225\350\277\233\350\241\214\345\237\272\346\234\254\347\232\204\346\225\260\346\215\256\345\210\206\346\236\220-\347\233\270\345\205\263\344\270\216\345\233\236\345\275\222.html" "b/Book/_book/\350\241\245\345\205\2051\345\246\202\344\275\225\350\277\233\350\241\214\345\237\272\346\234\254\347\232\204\346\225\260\346\215\256\345\210\206\346\236\220-\347\233\270\345\205\263\344\270\216\345\233\236\345\275\222.html"
index 6672c4e3..ba9ce37a 100644
--- "a/Book/_book/\350\241\245\345\205\2051\345\246\202\344\275\225\350\277\233\350\241\214\345\237\272\346\234\254\347\232\204\346\225\260\346\215\256\345\210\206\346\236\220-\347\233\270\345\205\263\344\270\216\345\233\236\345\275\222.html"
+++ "b/Book/_book/\350\241\245\345\205\2051\345\246\202\344\275\225\350\277\233\350\241\214\345\237\272\346\234\254\347\232\204\346\225\260\346\215\256\345\210\206\346\236\220-\347\233\270\345\205\263\344\270\216\345\233\236\345\275\222.html"
@@ -136,6 +136,29 @@
div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
+
@@ -162,8 +185,8 @@
- 关于作者
- 1 第一讲:为什么要学习R
@@ -215,7 +238,7 @@
- 3 第三讲:Git & RStudio工作流
-- 3.1 Files and Folders System
+- 3.1 文件与文件夹结构系统
- 3.2 Git and Git Hub
- 3.2.1 Git
@@ -609,8 +632,8 @@ 15.2 相关-代码实现 ) %>% # 根据子项目求综合平均
dplyr::select(sex, stress_mean, scontrol_mean)
使用head查看一下前五行
- - + +bruceR::Corr()
results.Corr <- capture.output({
bruceR::Corr(data = df.pg.corr[,c(2,3)],
diff --git "a/Book/_book/\350\241\245\345\205\2052\344\273\216\345\210\206\346\236\220\345\210\260\346\211\213\347\250\277.html" "b/Book/_book/\350\241\245\345\205\2052\344\273\216\345\210\206\346\236\220\345\210\260\346\211\213\347\250\277.html"
index b43d9e80..fc07b027 100644
--- "a/Book/_book/\350\241\245\345\205\2052\344\273\216\345\210\206\346\236\220\345\210\260\346\211\213\347\250\277.html"
+++ "b/Book/_book/\350\241\245\345\205\2052\344\273\216\345\210\206\346\236\220\345\210\260\346\211\213\347\250\277.html"
@@ -136,6 +136,29 @@
div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
+
@@ -162,8 +185,8 @@
- 关于作者
- 1 第一讲:为什么要学习R
@@ -215,7 +238,7 @@
- 3 第三讲:Git & RStudio工作流
-- 3.1 Files and Folders System
+- 3.1 文件与文件夹结构系统
- 3.2 Git and Git Hub
- 3.2.1 Git
@@ -634,7 +657,7 @@ 16.1.2 Part2: Papaja语
(2)包和文件的调用:
在定义YAML头文件之后,需要进行检查是否该文件已经调用了需要使用的包,以及是否包含了需要的文件(例如.bib格式文件)。
(3)参考文献的引用方式
-如例子中(R-papaja?),这是Papaja的引用格式,包含中括号[]、(符号以及参考文献标识?)。参考文献标识为.bib文件中@后面的内容。这种标识转换后为APA标准引用,以括号的形式引用文献。
+如例子中(Aust and Barth, n.d.),这是papaja
的引用格式,包含中括号[]、@符号以及参考文献标识。参考文献标识为.bib
文件中@后面的内容。这种标识转换后为APA标准引用,以括号的形式引用文献。
还有我们通常使用的另一种引用,即在句子中直接引用文章作者,此时只需要略去大括号即可,如@Hu_2020

@@ -679,6 +702,12 @@ 16.1.4.1 小作业
+
+
+Aust, Frederik, and Marius Barth. n.d. “papaja: Create APA Manuscripts with R Markdown.”
+
diff --git "a/Book/_book/\350\241\245\345\205\2053\346\225\210\345\272\224\351\207\217\345\222\214\345\205\203\345\210\206\346\236\220.html" "b/Book/_book/\350\241\245\345\205\2053\346\225\210\345\272\224\351\207\217\345\222\214\345\205\203\345\210\206\346\236\220.html"
index b7ea8c82..6cf4e8c9 100644
--- "a/Book/_book/\350\241\245\345\205\2053\346\225\210\345\272\224\351\207\217\345\222\214\345\205\203\345\210\206\346\236\220.html"
+++ "b/Book/_book/\350\241\245\345\205\2053\346\225\210\345\272\224\351\207\217\345\222\214\345\205\203\345\210\206\346\236\220.html"
@@ -136,6 +136,29 @@
div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
+
@@ -162,8 +185,8 @@
- 关于作者
- 1 第一讲:为什么要学习R
@@ -215,7 +238,7 @@
- 3 第三讲:Git & RStudio工作流
-- 3.1 Files and Folders System
+- 3.1 文件与文件夹结构系统
- 3.2 Git and Git Hub
- 3.2.1 Git
diff --git a/Book/_main_files/figure-html/unnamed-chunk-255-1.png b/Book/_main_files/figure-html/unnamed-chunk-255-1.png
index a592fd91..59843991 100644
Binary files a/Book/_main_files/figure-html/unnamed-chunk-255-1.png and b/Book/_main_files/figure-html/unnamed-chunk-255-1.png differ
diff --git a/Book/_main_files/figure-html/unnamed-chunk-275-1.png b/Book/_main_files/figure-html/unnamed-chunk-275-1.png
index c6d3f2cb..259ff0c7 100644
Binary files a/Book/_main_files/figure-html/unnamed-chunk-275-1.png and b/Book/_main_files/figure-html/unnamed-chunk-275-1.png differ
diff --git a/Book/_main_files/figure-html/unnamed-chunk-277-1.png b/Book/_main_files/figure-html/unnamed-chunk-277-1.png
index f1fb41ac..d3b22dda 100644
Binary files a/Book/_main_files/figure-html/unnamed-chunk-277-1.png and b/Book/_main_files/figure-html/unnamed-chunk-277-1.png differ
diff --git a/Book/_main_files/figure-html/unnamed-chunk-278-1.png b/Book/_main_files/figure-html/unnamed-chunk-278-1.png
index a563a47b..7bc8452c 100644
Binary files a/Book/_main_files/figure-html/unnamed-chunk-278-1.png and b/Book/_main_files/figure-html/unnamed-chunk-278-1.png differ
diff --git a/Book/_main_files/figure-html/unnamed-chunk-288-1.png b/Book/_main_files/figure-html/unnamed-chunk-288-1.png
index 0352fd7b..d4d977e5 100644
Binary files a/Book/_main_files/figure-html/unnamed-chunk-288-1.png and b/Book/_main_files/figure-html/unnamed-chunk-288-1.png differ
diff --git a/Book/_main_files/figure-html/unnamed-chunk-289-1.png b/Book/_main_files/figure-html/unnamed-chunk-289-1.png
index 39df407a..8ad22ad5 100644
Binary files a/Book/_main_files/figure-html/unnamed-chunk-289-1.png and b/Book/_main_files/figure-html/unnamed-chunk-289-1.png differ
diff --git a/Book/_main_files/figure-latex/Outlier-1.pdf b/Book/_main_files/figure-latex/Outlier-1.pdf
index e581d300..5341a0dc 100644
Binary files a/Book/_main_files/figure-latex/Outlier-1.pdf and b/Book/_main_files/figure-latex/Outlier-1.pdf differ
diff --git a/Book/_main_files/figure-latex/corr2-1.pdf b/Book/_main_files/figure-latex/corr2-1.pdf
index 878ff12c..45671ba3 100644
Binary files a/Book/_main_files/figure-latex/corr2-1.pdf and b/Book/_main_files/figure-latex/corr2-1.pdf differ
diff --git a/Book/_main_files/figure-latex/forest plot-1.pdf b/Book/_main_files/figure-latex/forest plot-1.pdf
index d994acbd..6b29feff 100644
Binary files a/Book/_main_files/figure-latex/forest plot-1.pdf and b/Book/_main_files/figure-latex/forest plot-1.pdf differ
diff --git a/Book/_main_files/figure-latex/funnel plot-1.pdf b/Book/_main_files/figure-latex/funnel plot-1.pdf
index 519dd9ec..ec72b191 100644
Binary files a/Book/_main_files/figure-latex/funnel plot-1.pdf and b/Book/_main_files/figure-latex/funnel plot-1.pdf differ
diff --git a/Book/_main_files/figure-latex/leave 1 out method-1.pdf b/Book/_main_files/figure-latex/leave 1 out method-1.pdf
index 611b56fd..46b0cf44 100644
Binary files a/Book/_main_files/figure-latex/leave 1 out method-1.pdf and b/Book/_main_files/figure-latex/leave 1 out method-1.pdf differ
diff --git a/Book/_main_files/figure-latex/plot.lm-1.pdf b/Book/_main_files/figure-latex/plot.lm-1.pdf
index c6d5bedb..325054e7 100644
Binary files a/Book/_main_files/figure-latex/plot.lm-1.pdf and b/Book/_main_files/figure-latex/plot.lm-1.pdf differ
diff --git a/Book/_main_files/figure-latex/result.Corr-1.pdf b/Book/_main_files/figure-latex/result.Corr-1.pdf
index 9cbeb75b..49a95f37 100644
Binary files a/Book/_main_files/figure-latex/result.Corr-1.pdf and b/Book/_main_files/figure-latex/result.Corr-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-113-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-113-1.pdf
index db708aae..5f8ad2ce 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-113-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-113-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-114-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-114-1.pdf
index 45097e98..7fde2838 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-114-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-114-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-115-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-115-1.pdf
index dd6de1e7..41c1ef43 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-115-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-115-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-116-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-116-1.pdf
index ccbcc791..3d90e5d0 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-116-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-116-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-118-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-118-1.pdf
index b73f1599..fee37626 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-118-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-118-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-120-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-120-1.pdf
index c332f812..20b5821d 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-120-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-120-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-122-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-122-1.pdf
index 7afdfa0f..2068c623 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-122-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-122-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-124-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-124-1.pdf
index 9cad35c3..bcab1099 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-124-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-124-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-126-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-126-1.pdf
index 5e8d78b1..2c4f2a8b 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-126-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-126-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-129-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-129-1.pdf
index 4d3f0e1e..3d8b1f2b 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-129-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-129-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-130-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-130-1.pdf
index adb278ed..9eaebd38 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-130-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-130-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-131-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-131-1.pdf
index 382b7039..93f62842 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-131-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-131-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-132-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-132-1.pdf
index 4daef0df..4f6c5ec3 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-132-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-132-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-133-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-133-1.pdf
index 759b847d..a8cd902a 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-133-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-133-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-134-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-134-1.pdf
index d5c8b0fc..a0b199e3 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-134-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-134-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-135-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-135-1.pdf
index c4f1c7f5..80d2c626 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-135-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-135-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-136-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-136-1.pdf
index 94e9134f..1e6c37f5 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-136-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-136-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-156-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-156-1.pdf
index 00c4907e..e2b1e315 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-156-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-156-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-157-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-157-1.pdf
new file mode 100644
index 00000000..b2af4851
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-157-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-186-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-186-1.pdf
new file mode 100644
index 00000000..4a50e5d5
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-186-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-187-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-187-1.pdf
new file mode 100644
index 00000000..58142892
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-187-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-188-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-188-1.pdf
new file mode 100644
index 00000000..4ec9c305
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-188-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-189-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-189-1.pdf
new file mode 100644
index 00000000..d0d7a2eb
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-189-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-190-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-190-1.pdf
index 92c6ab2a..bcb75a8e 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-190-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-190-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-193-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-193-1.pdf
new file mode 100644
index 00000000..1f20885a
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-193-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-203-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-203-1.pdf
new file mode 100644
index 00000000..be5372d6
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-203-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-214-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-214-1.pdf
new file mode 100644
index 00000000..3bf133f7
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-214-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-215-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-215-1.pdf
new file mode 100644
index 00000000..e9a005ed
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-215-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-216-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-216-1.pdf
new file mode 100644
index 00000000..5994f2fd
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-216-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-224-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-224-1.pdf
new file mode 100644
index 00000000..ac0adb7d
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-224-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-238-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-238-1.pdf
index e5bae779..5be347b0 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-238-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-238-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-239-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-239-1.pdf
index 6736d14e..94cfefc1 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-239-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-239-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-241-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-241-1.pdf
index 0f19a6eb..097f6a2c 100644
Binary files a/Book/_main_files/figure-latex/unnamed-chunk-241-1.pdf and b/Book/_main_files/figure-latex/unnamed-chunk-241-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-255-1.png b/Book/_main_files/figure-latex/unnamed-chunk-255-1.png
new file mode 100644
index 00000000..f484308d
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-255-1.png differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-256-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-256-1.pdf
new file mode 100644
index 00000000..71732a35
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-256-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-267-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-267-1.pdf
new file mode 100644
index 00000000..2ba99ffb
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-267-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-268-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-268-1.pdf
new file mode 100644
index 00000000..86598105
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-268-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-269-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-269-1.pdf
new file mode 100644
index 00000000..8fbed969
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-269-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-271-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-271-1.pdf
new file mode 100644
index 00000000..0057fbe6
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-271-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-272-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-272-1.pdf
new file mode 100644
index 00000000..779e4948
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-272-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-273-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-273-1.pdf
new file mode 100644
index 00000000..65648341
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-273-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-274-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-274-1.pdf
new file mode 100644
index 00000000..4a3d7992
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-274-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-275-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-275-1.pdf
new file mode 100644
index 00000000..fda07d15
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-275-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-277-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-277-1.pdf
new file mode 100644
index 00000000..0ca85b0a
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-277-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-278-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-278-1.pdf
new file mode 100644
index 00000000..0726620d
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-278-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-279-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-279-1.pdf
new file mode 100644
index 00000000..2ece3a2a
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-279-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-280-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-280-1.pdf
new file mode 100644
index 00000000..1895d16e
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-280-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-281-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-281-1.pdf
new file mode 100644
index 00000000..b14dba76
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-281-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-285-1.pdf b/Book/_main_files/figure-latex/unnamed-chunk-285-1.pdf
new file mode 100644
index 00000000..cf10ab38
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-285-1.pdf differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-288-1.png b/Book/_main_files/figure-latex/unnamed-chunk-288-1.png
new file mode 100644
index 00000000..ff5376b9
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-288-1.png differ
diff --git a/Book/_main_files/figure-latex/unnamed-chunk-289-1.png b/Book/_main_files/figure-latex/unnamed-chunk-289-1.png
new file mode 100644
index 00000000..bbb2a214
Binary files /dev/null and b/Book/_main_files/figure-latex/unnamed-chunk-289-1.png differ
diff --git a/Book/mybib.bib b/Book/mybib.bib
index 3d171582..fc50d5f3 100644
--- a/Book/mybib.bib
+++ b/Book/mybib.bib
@@ -32,3 +32,8 @@ @Book{QGCZ2011:StochProc
year = {2011},
publisher = {高等教育出版社},
}
+
+@article{R-papaja,
+ author = {Frederik Aust and Marius Barth},
+ title = {{papaja}: {Create} {APA} manuscripts with {R Markdown}}
+}
\ No newline at end of file
diff --git a/Book/pic/chp12/p1.pdf b/Book/pic/chp12/p1.pdf
index 87eb8dd6..53ce121e 100644
Binary files a/Book/pic/chp12/p1.pdf and b/Book/pic/chp12/p1.pdf differ
diff --git a/Book/rsconnect/documents/index.Rmd/bookdown.org/hcp4715/R4PsyBook.dcf b/Book/rsconnect/documents/index.Rmd/bookdown.org/hcp4715/R4PsyBook.dcf
index 31539e7a..dbc714aa 100644
--- a/Book/rsconnect/documents/index.Rmd/bookdown.org/hcp4715/R4PsyBook.dcf
+++ b/Book/rsconnect/documents/index.Rmd/bookdown.org/hcp4715/R4PsyBook.dcf
@@ -5,6 +5,6 @@ account: hcp4715
server: bookdown.org
hostUrl: https://bookdown.org/__api__
appId: 13774
-bundleId: 100533
+bundleId: 100556
url: https://bookdown.org/hcp4715/R4PsyBook/
version: 1
diff --git a/README.md b/README.md
index 0d33014f..af8e5821 100644
--- a/README.md
+++ b/README.md
@@ -49,3 +49,18 @@ Shield: [.
[](http://creativecommons.org/licenses/by/4.0/)
+
+