Publications by Adolfo Manoel
Exercícios TVE e Cópulas
Os dados para esses exercícios foram obtidos em https://finance.yahoo.com/quote/PBR/history?p=PBR . (Petróleo Brasileiro S.A. - Petrobras (PBR) Feb 09, 2020 - Feb 09, 2021). 1 Exercício 1 Estimação da estável - Método dos quantis 1.1 Pacotes Utilizados ## Instalando e Carregando os pacotes (bibliotecas) lbs<-c('tidyverse','lubridate',...
1127 sym R (9159 sym/43 pcs) 13 img 2 tbl
Atividade
Desempenho dos estimadores \(\hat{\beta}_0\) e \(\hat{\beta}_1\) # Vetor de Probabilidade, Gradiente e Hessiana Matriz<-function(b0,b1,X,y,P,W){ beta<-c(b0,b1) P <- 1 / ( 1+exp(-X %*% beta) ) Grad<-t(X) %*% (y-P) Hess<- -t(X) %*% W %*% X return(list(prob=P, grad=Grad, hess=Hess)) } ...
93 sym R (2032 sym/6 pcs) 1 tbl
Questão 06
Questão 01 (1o Modo) theta <- seq(0.01,0.99,0.02) var.est <- rep(0,length(theta)) n <- 10 licr <- ((1-2*theta)^2)*theta*(1-theta)/n na <- 50000 for (i in 1:length(theta)) { vec <- apply(matrix(rbinom(na*n,1,theta[i]),na,n),1,mean) estim <- n*vec*(1-vec)/(n-1) var.est[i] <- var(estim) } plot(theta,licr,type="l",ylim=c(0,0.007)) lines(t...
6581 sym R (5032 sym/11 pcs) 10 img
Implementação Prova 3
Questões 01 e 02 ## Geracao da amostra pela inversa da FDA theta<-2 # parametro u<-runif(1e4,0,1) x<- -log( 1 - u^(1/theta) ) ## Quantis da distribuicao Gama(n,1) alpha<-0.05 ; n<-length(x) q1<-qgamma(alpha/2,n,1) q2<-qgamma(1-alpha/2,n,1) ## Limites Inferior e Superior do IC li<- q1/sum(-log(1-exp(-x))) ls<- q2/sum(-log(1-exp(-x))) ...
90 sym R (2720 sym/20 pcs)
Questão 03 c)
Solução Item 3c) Dada a distribuição conjunta \[ \begin{align} f_{X,Y}(x,y)=\begin{cases} 4xy, & 0<x,y<1 \\\\ 0, & \mbox{c.c.}, \end{cases} \end{align} \] verifica-se, facilmente, que \[ \begin{align} f_X(x)=2x, \quad x\in(0,1) \quad \wedge \quad f_Y(y)=2y, \quad y\in(0,1). \end{align} ...
3829 sym R (729 sym/4 pcs)
Função Decisão do Teste
Questão 10 ### Funcao Decisao Item a) decisao<-function(amostra){ n<-length(amostra) if (sum(amostra) > qnorm(0.95)*sqrt(theta0*n) + n*theta0 ){ return("a1: Há fortes evidencias para rejeitar H0") }else{ return("a0: Não há evidencias suficientes para rejeitar H0") } } theta0<-2 ## H0: theta = 2 X H1: theta >...
107 sym R (3754 sym/36 pcs)
Questao Extra
Cálculo da probabilidade Como \(X\sim U(0,a)\) e \(Y\sim U(0,b)\), com \(0<a\leq b\); seque que $$\[\begin{align} F_Y(x)=\begin{cases} 0, & x<0,\\ \frac{x}{b}, & 0\leq x <b,\\ 1, & x\geq b. \end{cases} \qquad \wedge \qquad f_X(x)=\begin{cases} \frac{1}{a}, & x\in(0,a),\\ 0, & x\not\in (0,a)...
1071 sym R (535 sym/2 pcs)
Controle de Qualidade
Exemplo Um procedimento de controle de qualidade foi planejado para garantir um máximo de \(10\%\) de peças defeituosas na produção. A cada 5 minutos, sorteia-se uma amostra de 20 peças e, havendo mais de \(15\%\) defeituosas, encerra-se a produção para verificações. Objetivo Desejamos saber se a parada foi desnecessária e o momento da...
360 sym R (947 sym/5 pcs)
LimitesGraficos
Questão 90 Item a) ggplot()+ xlim(-1,1)+ geom_function(fun = function(x) sin(3*x)/(2*x), col="red") + geom_hline(yintercept = 1.5,col="blue")+ geom_vline(xintercept = 0,col="green")+ annotate("text", x = 0.6, y = 1.4, parse = TRUE, size = 4, label = "'Limite: ' * lim[x%->%0]~frac(sen~3*x,2*x)==frac(3,2)") Item b) gg...
332 sym R (7953 sym/24 pcs) 24 img
Exercício 98 b)
Geraçao de amostra da Estatística \(Y = X_{(n)}\). # Transformacao Integral da Probabilidade set.seed(24/11/2021) n <- 1e3 theta <- 2 u <- runif(n,0,1) y <- theta*u^(1/n) Convergência 1 # a ser feito Convergência 2 # a ser feito ...
95 sym R (150 sym/3 pcs)