TURMA 4.1 - Como indexar 2 condições para o mesmo fator de exploração de dados

Boa noite!

Tenho uma dúvida de como indexar 2 condições para gerar resultados de média, variância, sd, se....
Ex;
Tenho 2 variáveis explicativas qualitativas (Variedade e Região) e tenho 1 variável resposta quantitativa (Oxigenio 18/16)

UA VARIEDADE O18 REGIAO
AMOSTRA_1 Cabernet_Sauvignon 1,56 reg1
AMOSTRA_2 Merlot 1,32 reg1
AMOSTRA_3 Merlot 1,57 reg2
AMOSTRA_4 Cabernet_Sauvignon 2,16 reg2
AMOSTRA_5 Cabernet_Sauvignon 2,27 reg3
AMOSTRA_6 Merlot 2,00 reg3
....
....

Como se faz para indexar no sentido de obter valor de erro padrão para O18, na região1 e para cultivar Merlot?
Vi no curso somente a indexação para uma variável qualitativa.
se(O18[REGIAO=="reg1"])
E para reg1 + VARIEDADE "Merlot"?

Alguém pode ajudar?

Obrigado!!
Leonardo Reffatti

Olá, Leonardo!

Oba, vinhos!

Vamos lá, vou mostrar por exemplo:

Código:

#Os dados, para o R:
vinhos<-structure(list(UA = structure(1:6, .Label = c("AMOSTRA_1", "AMOSTRA_2",
"AMOSTRA_3", "AMOSTRA_4", "AMOSTRA_5", "AMOSTRA_6"), class = "factor"),
    VARIEDADE = structure(c(1L, 2L, 2L, 1L, 1L, 2L), .Label = c("Cabernet_Sauvignon",
    "Merlot"), class = "factor"), O18 = c(1.56, 1.32, 1.57, 2.16,
    2.27, 2), REGIAO = structure(c(1L, 1L, 2L, 2L, 3L, 3L), .Label = c("reg1",
    "reg2", "reg3"), class = "factor")), class = "data.frame", row.names = c(NA,
-6L))

#Conferindo:
summary(vinhos)

library(sciplot)  #Para a função se

#O que você já fez:
se(vinhos$O18[vinhos$REGIAO=="reg1"])

#O que você queria:
se(vinhos$O18[vinhos$REGIAO=="reg1" & vinhos$VARIEDADE=="Merlot"])

#Mas o R não conseguiu calcular, pois tem apenas um Merlot na região 1!
#Caso a planilha fosse maior, daria certo tranquilamente