Zignoruj ​​NA w funkcji sapply

Question

Używam R i szukałem odpowiedzi, ale chociaż widziałem podobne pytania, to nie zadziałało dla mojego konkretnego problemu.

W moim zestawie danych próbuję użyć znaków NA jako symboli zastępczych, ponieważ zamierzam do nich wrócić po tym, jak otrzymam część mojej analizy, dlatego też chciałbym móc wykonać wszystkie moje obliczenia jako jeśli NA naprawdę tam nie było.

Oto mój problem z tabeli przykładowe dane

ROCA = c(1,3,6,2,1,NA,2,NA,1,NA,4,NA) 
ROCA <- data.frame (ROCA=ROCA)  # converting it just because that is the format of my original data 

#Now my function 
exceedes <- function (L=NULL, R=NULL, na.rm = T) 
{ 
    if (is.null(L) | is.null(R)) { 
     print ("mycols: invalid L,R.") 
     return (NULL)    
    } 
    test <-(mean(L, na.rm=TRUE)-R*sd(L,na.rm=TRUE)) 
    test1 <- sapply(L,function(x) if((x)> test){1} else {0}) 
    return (test1) 
} 
L=ROCA[,1] 
R=.5 
ROCA$newcolumn <- exceedes(L,R) 
names(ROCA)[names(ROCA)=="newcolumn"]="Exceedes1" 

otrzymuję błąd:

Error in if ((x) > test) { : missing value where TRUE/FALSE needed 

Jak znacie, to jest coś nie tak z funkcji sapply. Wszelkie pomysły na ignorowanie tych NA? Chciałbym spróbować na.omit, gdybym mógł go wstawić, aby wstawić wszystkie NA w miejscu, w którym były wcześniej, ale nie jestem pewien, jak to zrobić.

Answer 1

Stwierdzenie to jest dziwne:

test1 <- sapply(L,function(x) if((x)> test){1} else {0}) 

Spróbuj:

test1 <- ifelse(is.na(L), NA, ifelse(L > test, 1, 0)) 

Answer 2

Chcesz NA: S w wyniku? To znaczy, czy wiersze mają się składać?

wydaje się po prostu powracać L > test będzie działać wtedy. Dodanie kolumny również może być uproszczone (podejrzewam, że "Exeedes1" jest gdzieś w zmiennej).

exceedes <- function (L=NULL, R=NULL, na.rm = T) 
{ 
    if (is.null(L) | is.null(R)) { 
     print ("mycols: invalid L,R.") 
     return (NULL)    
    } 
    test <-(mean(L, na.rm=TRUE)-R*sd(L,na.rm=TRUE)) 

    L > test 
} 
L=ROCA[,1] 
R=.5 
ROCA[["Exceedes1"]] <- exceedes(L,R) 

Answer 3

Nie ma potrzeby sapply a ponieważ > anonimowa funkcja jest już wektorowy.

Również wydaje się dziwne określenie domyślnych wartości argumentów, które są nieprawidłowe. Domyślam się, że używasz tego jako kludge zamiast używać funkcji missing. Dobrą praktyką jest także zgłaszanie błędu, zamiast zwracania NULL, ponieważ nadal trzeba będzie próbować złapać, gdy funkcja zwróci NULL.

exceedes <- function (L, R, na.rm=TRUE) 
{ 
    if(missing(L) || missing(R)) { 
    stop("L and R must be provided") 
    } 
    test <- mean(L,na.rm=TRUE)-R*sd(L,na.rm=TRUE) 
    as.numeric(L > test) 
} 

ROCA <- data.frame(ROCA=c(1,3,6,2,1,NA,2,NA,1,NA,4,NA)) 
ROCA$Exceeds1 <- exceedes(ROCA[,1],0.5)