Dados faltantes e outras dúvidas

Olá Marcos,

Estou me aventurando com meus dados no R e tem algumas dúvidas que estão surgindo.

Primeiro, estou tendo problemas para calcular parâmetros quando tenho algum dado NA. Teria uma outra forma de marcar dados faltantes? Porque o programa não está reconhecendo o NA como algo a ser ignorado para aquele cálculo, ele tenta calcular o NA junto.

Outra dúvida, preciso fazer um gráfico de barras que tenha quantos lagartos ativos em cada hora, assim como a variação da temperatura naqueles horários. Os dados são os seguintes: Estacao Data Hora Numero Sexo Idade CRC Tsb Tar Atividade

Como posso fazer pro R calcular quantos lagartos ativos tem no intervalo de hora de 0900 a 0959, por exemplo e depois fazer um gráfico de barras por horários com esses valores? É possível? E se eu quiser filtrar essa quantidade de lagartos ativos por sexo ou idade, eu posso também?

Obrigada,

Rê

Oi, Renata!

Vários comandos do R reconhecem o argumento na.rm, de forma que você pode dizer na.rm=T para indicar que os NAs devem ser ignorados. Então fica assim:

mean(variável, na.rm=T)

Não vai funcionar em todas as funções, mas deve dar certo na maioria.

Daqui a pouco monto um exemplo para responder a segunda pergunta e posto aqui, ok?

Abraços!

Ah, blz, vou tentar!

Thanks

Bom, vamos lá continuar. Nesta primeira parte, vou mostrar como podemos contar o número de lagartos em atividade dentro de uma faixa de horário. Depois continuo, avançando para as outras questões.

Os dados eu inventei, imaginando mais ou menos como devem ser os seus:

Código:


dados<-structure(list(Hora = c(900L, 904L, 910L, 935L, 936L, 938L, 942L,
945L, 947L, 949L, 951L, 957L, 1001L, 1010L, 1030L, 1045L), Sexo = structure(c(2L,
1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L), .Label = c("feminino",
"masculino"), class = "factor"), Idade = c(2L, 6L, 4L, 6L, 3L,
3L, 2L, 6L, 2L, 3L, 4L, 3L, 2L, 3L, 2L, 2L), Atividade = structure(c(2L,
1L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 2L), .Label = c("não",
"sim"), class = "factor")), .Names = c("Hora", "Sexo", "Idade",
"Atividade"), class = "data.frame", row.names = c(NA, -16L))

attach(dados)

dados

#Quantos lagartos ativos entre 09:00 e 09:59:

tabela.9<-table(Hora[Hora<1000], Atividade[Hora<1000])
tabela.9
colSums(tabela.9)

#Se precisar do valor direto, para usar em outras funções, por exemplo:

somas.9<-colSums(tabela.9)
somas.9["sim"]

#E se quiser um gráfico de barras simples:

barplot(somas.9, beside=T)
#Vamos ver mais sobre este tipo de gráfico no curso, com detalhes de como editar, colocar legenda e etc.

Dando continuidade, com os mesmos dados:

Código:


#Filtrando por outras variáveis:

#Vamos começar filtrando pelos ativos no mesmo intervalo e que são machos:
tabela.9.masc<-table(Hora[Hora<1000 & Sexo=="masculino"], Atividade[Hora<1000& Sexo=="masculino"])

#Agora os ativos das 09:30 em diante e com idade maior do que 3:
tabela.930.idmais3<-table(Hora[Hora>=930 & Idade>3], Atividade[Hora>=930 & Idade>3])

#E assim por diante.

Oi Marcos, consegui rodar teu exemplo, mas ainda está um pouco confuso sobre como fazer isso nos meus dados. Vou aguardar as novas aulas pra ver se fica mais claro. Obrigada

Renata Cardoso Vieira escreveu:Oi Marcos, consegui rodar teu exemplo, mas ainda está um pouco confuso sobre como fazer isso nos meus dados. Vou aguardar as novas aulas pra ver se fica mais claro. Obrigada

Então, Renata, o que eu fiz aqui são coisas que ao longo do curso vão aparecendo, pois em exemplos ao longo do curso nós vamos trabalhar com seleção de parte dos dados. Siga mais um pouco, mas qualquer coisa volte aqui para pedir ajuda, certo?

Abraços!