R中T检验
2022-05-01 本文已影响0人
pumpkinC
目的:使用R进行T检验分析。
data <- read.table("Pai.for.tTest", header=T)
t.test(data$Oilseed, data$sarson, alternative = "greater", paired = T, conf.level=0.99)
t.test(data$Oilseed, data$sarson, alternative = "greater", paired = T, conf.level=0.99)$p.value
文件格式(Pai.for.tTest):
Chr Pos Turniprape Oilseed sarson
A01 100000 0.3228234201884714 0.3490830380702748 0.1518795204597245
A01 105000 0.3245746994391083 0.3518123240508535 0.1546037695844976
A01 110000 0.3236503558450186 0.3521847444819548 0.1549785821320995
A01 115000 0.3190036423504874 0.3510108365304958 0.1576653434252903
A01 120000 0.3183541842914054 0.3512304957653968 0.1580963467984657
A01 125000 0.3195173160104653 0.3536353206979938 0.1580477810156511