我有一个五列的数据框,如下所示:
我想对 A2/A1 与 B2/B1 进行 t 检验。我能想到的一种方法是首先计算比率并创建一个新列 C=A2/A1
和 D=B2/B1
。然后对 C 与 D 进行 t 检验。但这是正确的方法吗?或者是否有其他方法可以使用原始数据进行比较?
谢谢!
答案 0 :(得分:0)
这将执行配对 t 检验。
text="Gene A1 A2 B1 B2
gene1 100 20 2 20
gene2 10 200 2 200
gene3 10 2 2 2"
df=read.table(text=text, header=TRUE)
df[,"C"]=df[,"A2"]/df[,"A1"]
df[,"D"]=df[,"B2"]/df[,"B1"]
t.test(df[,"C"], df[,"D"], paired=TRUE)