R语言 配对t检验,对子变量到底填在哪?t.test paired=TRUE???

使用本地数据 鸢尾花(yuān wěi huā)做配对t检验 t.test示例

提取两个物种的处理介绍看之前提到的空因子处理

#加载数据和包
data(iris)
library(tidyverse)
#提取这两个物种做t检验
iris2=filter(iris,Species=='setosa'|Species=='versicolor')
#Error in t.test.formula(iris$Species ~ iris$Sepal.Length) : 
#分组因子只能有两层
iris2$Species=droplevels(iris2$Species)
table(iris2$Species)
iris$block <-rep(1:50,3)
str(iris$block)
table(iris$block)
iris$block=as.factor(iris$block)

#秩和检验的区组是有设置 区组的变量 block

friedman.test(iris$Sepal.Length,iris$Species,iris$block)

而配对t检验中没有相应的变量填入

那么它到底是怎么查找的呢?只是设置了paired=True

t.test( iris2$Sepal.Length ~ iris2$Species,var.equal = T)
t.test(iris2$Sepal.Length ~ iris2$Species,equal = T,paired=TRUE)

可以看到p_value是不同的数值

在这里插入图片描述

下面是一个小实验,证明结果

#配对t检验的原理=对子相减,然后做差值的单样本t检验,与总体0相比

data2=filter(iris,Species=='setosa') 
data3=filter(iris,Species=='versicolor') 
data4=left_join(data2,data3,by=c('block'))
data5=mutate(data4,difference=data4$Sepal.Length.x-data4$Sepal.Length.y)
t.test(data5$difference,mu=0)

结果一致

在这里插入图片描述

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值