NOTAS CPE I - TARDE

Como avisado NÃO CORRIGI as reposições das pessoas que já estavam passadas! Parabéns aos aprovados, os 2 alunos que foram pra final podem me escrever para remarcar a prova para quinta, caso não escrevam a prova final será dia terça.

Nota na final é quanto o aluno precisa tirar na final pra passar.

NOTAS CPE I - MANHÃ

Desculpa a demora pessoal, mas aí estão as notas. As pessoas que tiveram média abaixo de 4 PODERÃO fazer a prova final. A data pode ser mudada pra quinta de manhã, mas pra isso preciso que todos me escrevam até amanhã no meu email. Alyssonx@gmail.com


Nota na final é quanto o aluno precisa tirar na final pra passar.

Notas Estatística Vital

Os quadros amarelos são trabalhos ou listas que não foram entregues ou que foram entregues atrazadas na secretaria do Departamento de Estatística e ainda não pude pegar. Segunda feira entrego essas notas. (Para ver as notas basta clicar na imagem)

Notas de CPE - I Tarde

Devo informar aos alunos aprovados que não poderão repor a nota mais baixa, apenas os que não obteram média acima de 7 e os que faltaram provas poderam fazer reposição. Parabéns aos aprovados e espero que usem o conhecimento abtido no curso.

Análise de regressão simples no R

Existe uma grande quantidade de análises que podem ser feitas encima de uma análise de regressão, como por exemplo uma análise de resíduos. Porém meu intuito aqui é só mostrar a parte introdutória.

Para iniciar precisamos de duas variáveis, uma dependente ou resposta e a outra independente ou explicativa, que chamaremos de x e y e que serão gerados de forma aleatória. São gerados 20 valores atravéz da distribuição uniforme entre 1 e 10 para x e entre 10 e 20 para y, após isso esses valores são ordenados em ordem crecente pelo comando sort().

set.seed(1)
x=sort(runif(20,1,10))
set.seed(2)
y=sort(runif(20,10,20))

Note que usei o comando set.seed antes de cada comando para que os valores aleatórios sejam os mesmo em qualquer lugar em que eles forem gerados, ou seja, os resultados daqui serão iguais pra qualquer pessoa que queira reproduzir o exemplo. Anyway, seguimos começando a análise onde a primeira coisa que olhamos é a correlação das variáveis.

cor(x,y)
[1] 0.982508

Veja que é uma ótima correlação. Em seguida fazemos a análise de regressão com o comando lm() de onde construímos nossa reta y=a+bx. Note que o y~x é a explicação de como deve ser a equação, ou seja, x explica y, ou ainda, y=x.

regressao=lm(y~x) summary(regressao)

Call: lm(formula = y ~ x) Residuals: Min 1Q Median 3Q Max -0.9339 -0.4361 0.1189 0.2904 0.9063 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 8.14171 0.33723 24.14 3.65e-15 *** x 1.16125 0.05188 22.38 1.37e-14 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 0.5823 on 18 degrees of freedom Multiple R-squared: 0.9653, Adjusted R-squared: 0.9634 F-statistic: 501.1 on 1 and 18 DF, p-value: 1.367e-14

Com isso nossa reta fica y= 8.14171 + 1.16125 x onde tanto o coeficiente angular quanto o linear são segnificantes à 1%. Maravilha de regressão! E pra completar o grafico de dispersão com a reta de regressão.

plot(x,y)
abline(regressao)

Guia cantado dos episódios de TNG

Sabem o que eu mais gosto deste vídeo? É saber que eu já assisti todos os episódios de TNG, HAHAHAHAAHAA, ai ai ah esquece... sim TNG = Star Trek: The Next Generation.

Distribuição de probabildiade de GURPS

Uma das primeiras coisas que quis fazer quando comecei a estudar probabilidade era construir a distribuição de probabilidade dos jogos de RPG que conheço. Essa semana consegui a de GURPS. Para vocês que não jogam RPG entenderem o sistema funciona da seguinte forma: O personagem possui um nível de habilidade (NH) em certa pericia, por exemplo salto. Joga-se 3 dados de 6 faces (normalmente conhecido como d6) e compare-se os resultados com o NH, se o resultado for menor ou igual ao NH o personagem é bem sucedido. Os testes de GURPS sofrem a critica que um redutor no NH causa uma redução da chance de acerto diferente para diferentes NHs, por exemplo ao retirar 1 de um NH 12 você causará um redutor de 11,6% na chance de passar em um teste, já um -1 em um NH 17 causará apenas 1,4% de penalidade. Mas apesar disso ainda acho GURPS meu sistema favorito.

Consegui esses resultados usando os seguintes comandos no R

d1=d2=d3=c(1:6)

resultados=expand.grid(d1,d2,d3)

soma.resultados=resultados$Var1+resultados$Var2+resultados$Var3

table(soma.resultados)

Isso ai me deu a tabela de frequência, o resto eu fiz no Excel mesmo, :)

O meu maior desafio será montar a FP de Storyteller...