目的:使用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
网友评论