R 换班和拿到号码
我正在研究R中的golub数据集（由AML和ALL分隔），我试图对两个基因进行假设检验。对于AML患者组，我想找出基因900表达高于基因1000的患者比例，然后我想看看在基因900表达值较高的患者中，这个数字是否少于一半。我有一个大致的想法来做另一半，我在第一部分有类似的想法，但考虑到它的T/F，我尝试将它转换为数字，给出0和1，但我想要实际的数字，而不是逻辑形式R 换班和拿到号码,r,dataset,logic,numeric,hypothesis-test,R,Dataset,Logic,Numeric,Hypothesis Test,我正在研究R中的golub数据集（由AML和ALL分隔），我试图对两个基因进行假设检验。对于AML患者组，我想找出基因900表达高于基因1000的患者比例，然后我想看看在基因900表达值较高的患者中，这个数字是否少于一半。我有一个大致的想法来做另一半，我在第一部分有类似的想法，但考虑到它的T/F，我尝试将它转换为数字，给出0和1，但我想要实际的数字，而不是逻辑形式 `gol.fac <- factor(golub.cl,levels=0:1, labels= c("ALL","AML"))
`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)`
`gol.fac使用max
max(golub[900,gol.fac=="AML"], golub[1000,gol.fac=="AML"])
pmax(golub[900,gol.fac=="AML"], golub[1000,gol.fac=="AML"])
或者，如果您有多个值，则使用pmax
max(golub[900,gol.fac=="AML"], golub[1000,gol.fac=="AML"])
pmax(golub[900,gol.fac=="AML"], golub[1000,gol.fac=="AML"])
使用max
max(golub[900,gol.fac=="AML"], golub[1000,gol.fac=="AML"])
pmax(golub[900,gol.fac=="AML"], golub[1000,gol.fac=="AML"])
或者，如果您有多个值，则使用pmax
max(golub[900,gol.fac=="AML"], golub[1000,gol.fac=="AML"])
pmax(golub[900,gol.fac=="AML"], golub[1000,gol.fac=="AML"])
不要对行进行多个切片，只需通过一次子集设置来获取max
max(golub[900:1000, "AML"])
不要对行进行多个切片，只需通过一次子集设置来获取max
max(golub[900:1000, "AML"])