изменение класса и получение номеров

Question

Я работаю с набором данных Голуб в R (разделенных AML и ALL), и я пытаюсь выполнить проверку гипотезы в отношении двух генов. Что касается группы пациентов с ОМЛ, я хочу выяснить, какая доля пациентов с более высокой экспрессией гена 900 по сравнению с геном 1000, а затем я хочу увидеть, если из тех, кто имеет более высокую величину экспрессии для гена 900,число меньше половины. У меня есть общая идея сделать вторую половину, и у меня было что-то подобное для первой части, но, учитывая, что его T / FI пытался переключить его на числовое значение, которое дало 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)`

Answer 1

Вместо того, чтобы делать несколько срезов строк, просто получите max путем подмножества один раз

max(golub[900:1000, "AML"])

Answer 2

Используйте max

max(golub[900,gol.fac=="AML"], golub[1000,gol.fac=="AML"])

Или, если у вас есть несколько значений, используйте pmax

pmax(golub[900,gol.fac=="AML"], golub[1000,gol.fac=="AML"])