abelhas/Morfometria.R at main · ogilmar/abelhas · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
#Selecionando o diretorio
setwd("~/Documentos_pessoais/Estudos_abelha_solitaria")


#Importando os dados
dados1=read.table("dados_1.txt",h=T)

#Indentificando as vari?veis
Repeticao=dados1[,1]
Abelhas=dados1[,2]
Sexo=dados1[,3]
Nascimento=dados1[,4]
Pluviosidade=dados1[,5]
Peso=dados1[,6]
Comprimento=dados1[,7]
MedidaItertecular=dados1[,8]
TemperaturaExterna=dados1[,9]
UmidadeExterna=dados1[,10]
TemperaturaInterna=dados1[,11]
UmidadeInterna=dados1[,12]

#Analisando as variaveis morfometricas

## Peso corporal

# Medidas de tendencia central e dispersao
resumoPeso <- c( media = mean(Peso),
                mediana = median(Peso),
                desvio_padrao = sd(Peso),
                minimo = min(Peso),
                maximo = max(Peso))
print(resumoPeso)

# Histograma dos pesos
ggplot(dados1, aes(x = Peso)) +
  geom_histogram(binwidth = 0.05, fill = "blue", color = "black", alpha = 0.7) +
  labs(title = "Histograma do Peso das Abelhas",
       x = "Peso (gramas)",
       y = "Frequ?ncia") +
  theme_minimal()

modeloPeso<-(glm(Peso~Pluviosidade+Sexo*Pluviosidade+TemperaturaExterna+TemperaturaInterna+UmidadeExterna+UmidadeInterna))
summary(modeloPeso)

# Grafico de densidade (distribui??o suavizada)
ggplot(dados1, aes(x = Peso, fill = NULL)) +
  geom_density(alpha = 0.5) +
  labs(title = "Densidade do Peso",
       x = "Peso (gramas)",
       y = "Densidade") +
  theme_minimal()

#ANOVA do peso

modeloPeso=glm(Peso~Pluviosidade+Sexo*Pluviosidade+TemperaturaExterna+TemperaturaInterna+UmidadeExterna+UmidadeInterna)
anova(modeloPeso,test = "F")

modeloComprimento=glm(Comprimento~Pluviosidade+Sexo*Pluviosidade+TemperaturaExterna+TemperaturaInterna+UmidadeExterna+UmidadeInterna)
anova(modeloComprimento,test = "F")

modeloMedIntertecular=glm(MedidaItertecular~Pluviosidade+Sexo*Pluviosidade+TemperaturaExterna+TemperaturaInterna+UmidadeExterna+UmidadeInterna)
anova(modeloMedIntertecular,test="F")