prob2

library(ggplot2)

set.seed(100)

dados_f <- rf(10,1,9)
data <- data.frame(Valores = dados_f)
pdf_plot10 <- ggplot(data, aes(x = Valores)) +
  geom_density(alpha = 0.5) +
  labs(title = paste("PDF da Distribuição F (n =10 k=9)"),
       x = "Valor",
       y = "Densidade de Probabilidade") +
  theme_minimal()

print(pdf_plot10)

dados_t <- rt(10, 9)
z1 <- dados_t^2

datat1 <- data.frame(Valorest1 = z1)
pdf_plott2_10 <- ggplot(datat1, aes(x = Valorest1)) +
  geom_density(alpha = 0.5) +
  labs(title = paste("PDF3 da Distribuição F aprox (n =10 k=9)"),
       x = "Valor",
       y = "Densidade de Probabilidade") +
  theme_minimal()

print(pdf_plott2_10)

#####################################################################################


dados_f2 <- rf(100,1,99)
data2 <- data.frame(Valores2 = dados_f2)
pdf_plot100 <- ggplot(data2, aes(x = Valores2)) +
  geom_density(alpha = 0.5) +
  labs(title = paste("PDF2 da Distribuição F (n =100 k=99)"),
       x = "Valor",
       y = "Densidade de Probabilidade") +
  theme_minimal()

print(pdf_plot100)

dados_t2 <- rt(100, 99)
z2 <- dados_t2^2

datat2 <- data.frame(Valorest2 = z2)
pdf_plott2_100 <- ggplot(datat2, aes(x = Valorest2)) +
  geom_density(alpha = 0.5) +
  labs(title = paste("PDF3 da Distribuição F aprox (n =100 k=99)"),
       x = "Valor",
       y = "Densidade de Probabilidade") +
  theme_minimal()

print(pdf_plott2_100)

####################################################################################




dados_f3 <- rf(10000,1,9999)
data3 <- data.frame(Valores3 = dados_f3)
pdf_plot10000 <- ggplot(data3, aes(x = Valores3)) +
  geom_density(alpha = 0.5) +
  labs(title = paste("PDF3 da Distribuição F (n =10000 k=9999)"),
       x = "Valor",
       y = "Densidade de Probabilidade") +
  theme_minimal()

print(pdf_plot10000)

dados_t3 <- rt(10000, 9999)
z3 <- dados_t3^2

datat3 <- data.frame(Valorest3 = z3)
pdf_plott2_10000 <- ggplot(datat3, aes(x = Valorest3)) +
  geom_density(alpha = 0.5) +
  labs(title = paste("PDF3 da Distribuição F aprox (n =10000 k=9999)"),
       x = "Valor",
       y = "Densidade de Probabilidade") +
  theme_minimal()

print(pdf_plott2_10000)

hist(dados_f)

hist(dados_t)

hist(z1)

hist(dados_f2)

hist(dados_t2)

hist(z2)

hist(dados_f3)

hist(dados_t3)

hist(z3)