我正在使用R中的golub数据集(由AML和ALL分隔),并尝试对两个基因进行假设检验。对于AML患者组,我想找出与基因1000相比,900基因表达更高的患者所占的比例,然后我想看看在那些900基因表达更高的患者中,这个数字是否少于一半。我有一个大致的想法来做另一半,我在第一部分有类似的东西,但看到它的T/F,我试图将它切换到数字,它给0和1,但我想要实际的数字,而不是逻辑形式。
`gol.fac <- factor(golub.cl,levels=0:1, labels= c("ALL","AML"))
x <- golub[900,gol.fac=="AML"]
y <- golub[1000,gol.fac=="AML"]
z <-golub[900,gol.fac=="AML"] > golub[1000,gol.fac=="AML"]
k <- as.numeric(z)`发布于 2019-10-07 11:31:32
使用max
max(golub[900,gol.fac=="AML"], golub[1000,gol.fac=="AML"])或者,如果您有多个值,则使用pmax
pmax(golub[900,gol.fac=="AML"], golub[1000,gol.fac=="AML"])发布于 2019-10-08 00:34:13
不需要进行多个行的切片,只需通过一次子集设置来获得max
max(golub[900:1000, "AML"])https://stackoverflow.com/questions/58262943
复制相似问题